all terms for additional CTM figures and variable-sized angled braces for averaged terms

This commit is contained in:
Yeongdo Park 2021-07-13 03:26:04 +09:00
parent 940d1d17b7
commit 2e20f720f2
2 changed files with 92 additions and 1 deletions

View file

@ -610,7 +610,7 @@ class AveragedField (FieldBase):
tfield = fdict[tgt]
self.fset = tfield.checkFluctuation()
self.latex = r"\langle {} \rangle".format(tfield.latex)
self.latex = r"\left\langle {} \right\rangle".format(tfield.latex)
if not self.fset:
self.tgt = tgt

View file

@ -0,0 +1,91 @@
# fsd weighted (vn+sd) correlations
[u, v, w, y]
c(latex="c") = 1.0 - y
wrate (latex="\dot{\omega}") = rxn_rate(c)
f (latex="\Sigma_f'")= sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))
uc = u' * c'
lapc(latex="\Delta c") = d2dx(c) + d2dy(c) + d2dz(c)
sd(latex="S_d") = ($rod * lapc + wrate) / f
nx(latex="n_x") = - ddx(c) / f
ny(latex="n_y") = - ddy(c) / f
nz(latex="n_z") = - ddz(c) / f
nx_lvl = sqr(nx')
ny_lvl = sqr(ny')
nz_lvl = sqr(nz')
nx2 = sqr(nx)
ny2 = sqr(ny)
nz2 = sqr(nz)
nn = nx2 + ny2 + nz2
unx = u * nx
vny = v * ny
wnz = w * nz
unx_fluc = u' * nx'
vny_fluc = v' * ny'
wnz_fluc = w' * nz'
vn(latex="v_n") = unx + vny + wnz
divn(latex="(\nabla \cdot \vec{n})") = ddx(nx) + ddy(ny) + ddz(nz)
divn2 = sqr(divn)
# grad(n) : grad(n)
gradn_gradn(latex="\nabla \vec{n} : \nabla \vec{n}") = sqr(ddx(nx)) + sqr(ddy(ny)) + sqr(ddz(nz)) + 2*(ddx(ny)*ddy(nx)) + 2*(ddx(nz)*ddz(nx)) + 2*(ddy(nz)*ddz(ny))
gradn_gradn_lvl = gradn_gradn' * gradn_gradn'
absk = abs(divn)
absk_lvl = sqr(absk')
vnsd = vn+sd
vnsd_nx = vnsd * nx
# u * nx_lvl
# nx * unx_fluc
sd_nx = sd' * nx'
u_nx_nx = u' * sqr(nx')
v_ny_nx = v' * ny' * nx'
w_nz_nx = w' * nz' * nx'
vnsd_divn = vnsd * divn
nx_divn = nx' * divn'
u_divn = u' * divn'
vn_divn (latex="\vec{v}'' \cdot \vec{n}'' (\nabla\cdot \vec{n})''") = (u' * nx' + v' * ny' + w' * nz') * divn'
sd_divn = sd' * divn'
avg { u, v, w, c, f, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c, uc }
avg f {
u, v, w,
nx, ny, nz,
nx2, ny2, nz2,
nn,
sd,
absk, divn,
absk_lvl,
nx_lvl, ny_lvl, nz_lvl,
unx, vny, wnz,
unx_fluc, vny_fluc, wnz_fluc,
vnsd_nx,
sd_nx,
u_nx_nx ,
v_ny_nx ,
w_nz_nx ,
vnsd_divn,
nx_divn,
u_divn,
vn_divn,
sd_divn
}