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 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

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`

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`

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

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:`