correct sign(divn)
This commit is contained in:
parent
7f87a663ec
commit
e5c2fe58e9
1 changed files with 9 additions and 3 deletions
|
|
@ -12,7 +12,10 @@ nx = - ddx(c) / fsd
|
||||||
ny = - ddy(c) / fsd
|
ny = - ddy(c) / fsd
|
||||||
nz = - ddz(c) / fsd
|
nz = - ddz(c) / fsd
|
||||||
|
|
||||||
vn = u * nx + v * ny + w * nz
|
unx = u * nx
|
||||||
|
vny = v * ny
|
||||||
|
wnz = w * nz
|
||||||
|
vn = unx + vny + wnz
|
||||||
|
|
||||||
divn = ddx(nx) + ddy(ny) + ddz(nz)
|
divn = ddx(nx) + ddy(ny) + ddz(nz)
|
||||||
|
|
||||||
|
|
@ -35,7 +38,7 @@ gradt_vnsd1 = ddx(vnsd) - nx * gradn_vnsd
|
||||||
|
|
||||||
divn2_dndn = divn2 - gradn_gradn
|
divn2_dndn = divn2 - gradn_gradn
|
||||||
|
|
||||||
sign_divn2_dndn = $dsign(divn2_dndn, divn)
|
sign_divn2_dndn = $dsign(1.0, divn) * divn2_dndn
|
||||||
|
|
||||||
vnsd_divn = vnsd * divn
|
vnsd_divn = vnsd * divn
|
||||||
|
|
||||||
|
|
@ -55,7 +58,10 @@ vnsd_sign_divn2_dndn = vnsd * sign_divn2_dndn
|
||||||
|
|
||||||
avg { c, u, v, w, fsd, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c }
|
avg { c, u, v, w, fsd, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c }
|
||||||
|
|
||||||
avg fsd { u, v, w, nx, ny, nz, absk, divn, sd, gradt_vnsd1,
|
avg fsd { u, v, w, nx, ny, nz, absk, divn,
|
||||||
|
unx, vny, wnz,
|
||||||
|
vn, sd,
|
||||||
|
gradt_vnsd1,
|
||||||
vnsd,
|
vnsd,
|
||||||
vnsd_divn,
|
vnsd_divn,
|
||||||
vnsd_divn_cor,
|
vnsd_divn_cor,
|
||||||
|
|
|
||||||
Loading…
Add table
Reference in a new issue