search

gerlero / Fronts.jl

Nonlinear diffusion problems in Julia
https://gerlero.github.io/Fronts.jl
MIT License
7 stars 0 forks source link

How to set gravitation potential #190

Open kongdd opened 2 months ago

kongdd commented 2 months ago

The Richards Equation is:

$$ \frac{\partial \theta}{\partial t}=\frac{\partial}{\partial z} [D \frac{\partial \theta}{\partial z}+K ] $$

It seems that Fronts.jl ignored the gravitation potential term $\frac{\partial K}{\partial z}$. Any suggestion how to solve this problem?

gerlero commented 2 months ago

Indeed Fronts implements the horizontal case of the Richards equation, where such a term is neglected (which is the case in pretty much all of our problems in microfluidics).

I'm not against adding support for such a term, but it'd need to dispatch to the FiniteDifference() algorithm because the default BoltzmannODE() algorithm is not compatible with that term.