96 lines
1.7 KiB
Text
96 lines
1.7 KiB
Text
# fsd weighted (vn+sd) correlations
|
|
|
|
[u, v, w, y]
|
|
|
|
c = 1.0 - y
|
|
wrate = rxn_rate(c)
|
|
fsd = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))
|
|
|
|
sd = ((d2dx(c) + d2dy(c) + d2dz(c)) * $rod + wrate) / fsd
|
|
|
|
nx = - ddx(c) / fsd
|
|
ny = - ddy(c) / fsd
|
|
nz = - ddz(c) / fsd
|
|
|
|
unx = u * nx
|
|
vny = v * ny
|
|
wnz = w * nz
|
|
vn = unx + vny + wnz
|
|
|
|
divn = ddx(nx) + ddy(ny) + ddz(nz)
|
|
|
|
divn2 = sqr(divn)
|
|
|
|
# grad(n) : grad(n)
|
|
gradn_gradn = sqr(ddx(nx)) + sqr(ddy(ny)) + sqr(ddz(nz)) + 2*(ddx(ny)*ddy(nx)) + 2*(ddx(nz)*ddz(nx)) + 2*(ddy(nz)*ddz(ny))
|
|
|
|
absk = abs(divn)
|
|
|
|
lapc = d2dx(c) + d2dy(c) + d2dz(c)
|
|
|
|
vnsd = vn+sd
|
|
|
|
divn_nx = divn*nx
|
|
|
|
gradn_vnsd = ddx(vnsd) * nx + ddy(vnsd) * ny + ddz(vnsd) * nz
|
|
|
|
gradt_vnsd1 = ddx(vnsd) - nx * gradn_vnsd
|
|
|
|
divn2_dndn = divn2 - gradn_gradn
|
|
|
|
sign_divn = $dsign(1.0, divn)
|
|
|
|
sign_divn2_dndn = sign_divn * divn2_dndn
|
|
|
|
vnsd_divn = vnsd * divn
|
|
|
|
vnsd_divn_cor = vnsd' * divn'
|
|
|
|
vnsd_divn_nx_cor = vnsd' * divn_nx'
|
|
|
|
vnsd_divn2_dndn_cor = vnsd' * divn2_dndn'
|
|
|
|
vnsd_sign_divn2_dndn_cor = vnsd' * sign_divn2_dndn'
|
|
|
|
vnsd_divn_nx = vnsd * divn_nx
|
|
|
|
vnsd_divn2_dndn = vnsd * divn2_dndn
|
|
|
|
vnsd_sign_divn2_dndn = vnsd * sign_divn2_dndn
|
|
|
|
vnsd_divn_nxfluc = vnsd * divn * nx'
|
|
|
|
vnsd_divn_divnfluc = vnsd * divn * divn'
|
|
|
|
vnsd_divn_abskfluc = vnsd * divn * absk'
|
|
|
|
vnsd_dndn = vnsd * gradn_gradn
|
|
|
|
sign_vnsd_dndn = sign_divn * vnsd * gradn_gradn
|
|
|
|
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,
|
|
unx, vny, wnz,
|
|
vn, sd,
|
|
gradt_vnsd1,
|
|
vnsd,
|
|
vnsd_divn,
|
|
vnsd_divn_cor,
|
|
divn_nx,
|
|
vnsd_divn_nx,
|
|
vnsd_divn_nx_cor,
|
|
divn2_dndn,
|
|
vnsd_divn2_dndn,
|
|
vnsd_divn2_dndn_cor,
|
|
sign_divn2_dndn,
|
|
vnsd_sign_divn2_dndn,
|
|
vnsd_sign_divn2_dndn_cor,
|
|
|
|
vnsd_divn_nxfluc,
|
|
vnsd_divn_divnfluc,
|
|
vnsd_divn_abskfluc,
|
|
vnsd_dndn,
|
|
sign_vnsd_dndn
|
|
|
|
}
|