diff --git a/code/code_gen/post.py b/code/code_gen/post.py index fdd3aa5..8914860 100644 --- a/code/code_gen/post.py +++ b/code/code_gen/post.py @@ -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 diff --git a/code/termspecs/termspec.20210708 b/code/termspecs/termspec.20210708 new file mode 100644 index 0000000..2322675 --- /dev/null +++ b/code/termspecs/termspec.20210708 @@ -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 +}