-
Notifications
You must be signed in to change notification settings - Fork 32
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
hypoellictic operator #45
Comments
Hi, what is the intended geometry ? Can you use a grid aligned with the coordinate directions ? |
Rectangle
Yes I can |
Ok, this is a good situation as otherwise the method as implemented now probably wouldn't converge.
You can use exact comparison here as the coordinates come from the arrays at the input. The you have convection-diffusion in the x direction which can be discretied by an upwind flux, and pure convection in y direction which can be discretized by a Godunov flux. I can give some more hints on these , just ask. Note that these fluxes introduce artificial diffusion which in particular would smear sharp fronts, so you should carefully judge the accuracy of the results. Central difference variants on finer grids may provide a more accurate alternative. If things don't work out properly with this package (which focuses on elliptic and parabolic problems), you may try Trixi.jl which is oriented at hyperbolic problems so it has more methods for the situation of vanishing diffusion. |
Hi,
I would like to use your package to solve an equation of the form
d/dt P = d2/dx2 P + d/dx( F(x,y) P) + d/dy (G(x,y) P)
with Dirichlet BC. However, I am struggling writing the flux function mainly because the "Laplacian" only applies to the x variable. Can you give me a hint please?
Thank you
The text was updated successfully, but these errors were encountered: