Commit ebacce5b authored by Lambert Theisen's avatar Lambert Theisen 🚀

Add idea to have modified inflow BC

- Replace un-unw = eps chi ((p-pw)+sigmann)
- By un-unw = eps chi ((p-pw)-sigmann)
- Should be allowed because kappa1 is arbitrary following Torrilhon2017 (see sec 5.2)
- Resulting system is then anti-symm
- But can be made fully sym by partial row inversion
- Add esol for change inflowBC but not used for now
parent 0e34b354
......@@ -579,7 +579,7 @@ class Solver:
(tt(psi_) + (1/2) * nn(psi_))
)
+ (1/xi_tilde) * nt(sigma_) * nt(psi_)
) * df.ds + sum([
) * df.ds + sum([ # TODO: Fix inflow with minus
bcs[bc]["epsilon_w"] * nn(sigma_) * nn(psi_) * df.ds(bc)
for bc in bcs.keys()
])
......@@ -628,6 +628,7 @@ class Solver:
lhs = [None] * 5
rhs = [None] * 5
# 1) Left-hand sides
# TODO: Fix inflow with minus in lhs4 before f
lhs[0] = +1*a(s, r) -b(theta, r)-c(r, sigma) +0 +0
lhs[1] = +1*b(kappa, s)+0 +0 +0 +0
lhs[2] = +1*c(s, psi) +0 +d(sigma, psi)-e(u, psi)+f(p, psi)
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment