research/incomp-flame-post
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

bugfix: FieldExporter construction, latex annotations to terms.vnsd_fluctuation

Browse source
This commit is contained in:
Yeongdo Park 2019-08-14 20:44:51 +09:00
parent a901ebe932
commit e9c61bf48c
2 changed files with 21 additions and 16 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
code/code_gen/post.py
View file

@ -431,7 +431,13 @@ class Field (FieldBase):
else:
self.latex = self.latex_given
self.exporter = self.attr.get("export")
self.exporter = None
try:
if self.attr["export"]:
self.exporter = FieldExporter(self.name, self.attr, self)
except KeyError:
pass
'''
for a in self.dep:

29
code/termspecs/termspec.vnsd_fluctuation
View file

@ -2,43 +2,42 @@
[u, v, w, y]
c = 1.0 - y
wrate = rxn_rate(c)
fsd = sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))
c(latex="c") = 1.0 - y
wrate (latex="\dot{\omega}") = rxn_rate(c)
fsd (latex="\Sigma_f'")= sqrt (sqr(ddx(c)) + sqr(ddy(c)) + sqr(ddz(c)))
sd = ((d2dx(c) + d2dy(c) + d2dz(c)) * $rod + wrate) / fsd
lapc(latex="\Delta c") = d2dx(c) + d2dy(c) + d2dz(c)
sd(latex="S_d") = ($rod * lapc + wrate) / fsd
nx = - ddx(c) / fsd
ny = - ddy(c) / fsd
nz = - ddz(c) / fsd
nx(latex="n_x") = - ddx(c) / fsd
ny(latex="n_y") = - ddy(c) / fsd
nz(latex="n_z") = - ddz(c) / fsd
unx = u * nx
vny = v * ny
wnz = w * nz
vn = unx + vny + wnz
vn(latex="v_n") = unx + vny + wnz
divn = ddx(nx) + ddy(ny) + ddz(nz)
divn(latex="\nabla \cdot \vec{n}") = 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))
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))
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
gradn_vnsd (latex="\nabla_{\mathrm{n}}(v_n+S_d)") = ddx(vnsd) * nx + ddy(vnsd) * ny + ddz(vnsd) * nz
gradt_vnsd1 = ddx(vnsd) - nx * gradn_vnsd
gradt_vnsd1 (latex="\nabla_{\top,x}(v_n+S_d)")= ddx(vnsd) - nx * gradn_vnsd
divn2_dndn = divn2 - gradn_gradn
sign_divn = $dsign(1.0, divn)
sign_divn(latex="\mathrm{sign}(\nabla \cdot \vec{n})") = $dsign(1.0, divn)
sign_divn2_dndn = sign_divn * divn2_dndn