Hello all,
I am working on a coastal aquifer 3D model, for wich before getting to the model itself, I am doing some previous simpler models in order to gradually understand the process before getting to a more complex model.
I am having problems to model a very simple variable density 3D model. It consists of a rectangular mesh with 6 Layers (1000x2400x45; length perpendicular to the coastline x length along the coastline x depth). Everything is fine if I run the simulation with a density ratio of 0 (i.e. without considering density variability), but once I apply a value of 0.025, the head and concentration distribution calculated are senseless (i.e. orders of magnitude without any sense like e+12 for the head ).
The BCs I apply are:
- Constant Saltwater Head at the seaside boundary
- Constant Flow (Neumann) at the inland boundary.
- Constant concentration (35g/L) at the seaside boundary. I apply a min-mass flow constraint of 0, in order to desactivate the BC where water flows out of the model (i.e. the shallow part of the boundary)
- Constant concentration (0 g/L) at the inland boundary.
Thank you