density changes (for density driven flow processes)

Georg_Kosakowski · June 21, 2012, 1:29pm

Hi all,
who is working on density driven flow processes, especially Richards_flow or liquid, groundwater flow and knows/understands where Density is called for these processes?
I would like to pass fluid densities from my chemical solver GEMS to the corresponding routines and wonder how to do this in an efficient way.
You may contact me directly via email (georg.kosakowski@psi.ch).
Thanks and greetings, Georg

jens-olaf_Delfs · June 21, 2012, 2:17pm

Hi,
If think you need density model 3.

If you take the Boussinesq approximation, you have to pass the concentration in the function
void CFiniteElementStd::Assemble_Gravity()
You have to specify dens_arg[2] = concentration and call Density(dens_arg)

(maybe it calls Density() in your version, then you have to include the densarg).

The *.mfp input looks like:

#FLUID_PROPERTIES
$FLUID_TYPE
LIQUID
$DENSITY
3 1e+3 0 0.03
$VISCOSITY

1 1e-3
#STOP

Better ask Ashok to confirm of this, because he is more familiar with this stuff than I am

···

2012/6/21 Georg Kosakowski georg.kosakowski@psi.ch

Hi all,
who is working on density driven flow processes, especially Richards_flow or liquid, groundwater flow and knows/understands where Density is called for these processes?

I would like to pass fluid densities from my chemical solver GEMS to the corresponding routines and wonder how to do this in an efficient way.
You may contact me directly via email (georg.kosakowski@psi.ch).

Thanks and greetings, Georg

wenqing.wang · June 21, 2012, 4:05pm

Hi Georg,

If the density is calculated at points and it is not needed to be

calculated again in the local assembly, I suggest to do as follows
for each global assembly:

calculate the density for each element at the centroid of the

element

save all values of element in a vector.

make CFluidProperties::Density(double* variables) return the density

value for an element from the vector.

Best regards,

Wenqing

···

On 21.06.2012 16:17, jens-olaf Delfs wrote:

Hi,

  If think you need density model 3.



  If you take the Boussinesq approximation, you have to pass the

concentration in the function

  void CFiniteElementStd::Assemble_Gravity()

  You have to specify dens_arg[2] = concentration and call

Density(dens_arg)

  (maybe it calls Density() in your version, then you have to

include the densarg).

  The *.mfp input looks like:



  #FLUID_PROPERTIES

   $FLUID_TYPE

    LIQUID

   $DENSITY

    3 1e+3 0 0.03

   $VISCOSITY

    1 1e-3

  #STOP



  Better ask Ashok to confirm of this, because he is more familiar

with this stuff than I am

2012/6/21 Georg Kosakowski georg.kosakowski@psi.ch

      Hi
all,
      who is working on density driven flow processes, especially
Richards_flow or liquid, groundwater flow and
knows/understands where Density is called for these processes?
      I would like to pass fluid densities from my chemical solver
GEMS to the corresponding routines and wonder how to do this
in an efficient way.
      You may contact me directly via email (georg.kosakowski@psi.ch).

      Thanks and greetings, Georg

Georg_Kosakowski · June 22, 2012, 11:27am

Gruetzi Wenqing,
thanks for the hint. I guess that is the way to go!
Greetings, Georg

···

Am Donnerstag, 21. Juni 2012 18:05:15 UTC+2 schrieb Wenqing Wang:

2012/6/21 Georg Kosakowski georg.kosakowski@psi.ch
      Hi
all,
      who is working on density driven flow processes, especially
Richards_flow or liquid, groundwater flow and
knows/understands where Density is called for these processes?
      I would like to pass fluid densities from my chemical solver
GEMS to the corresponding routines and wonder how to do this
in an efficient way.
      You may contact me directly via email (georg.kosakowski@psi.ch).

      Thanks and greetings, Georg

Hi,

  If think you need density model 3.



  If you take the Boussinesq approximation, you have to pass the