Temperature-dependent fluid properties

Mailing List Archives ogs-users
1

Dear all,
I’m using coupled Richards Flow and Heat transport for modeling heat propagation in a partially saturated porous medium.
I need OGS to change fluid properties (both water and air) in function of temperature during the simulation.
Does anyone know how to do that?
May I write a curve f(T, property value) in the .rfd file? or it works only for time-dependent curves?

thanks a lot
Nicolò

2

Thanks for quick reply Marc,
I will look for and take a look to Thomas Nagel option.
But actually I'd like to change thermal conductivity only because I'm working on T range 20-80 °C.
In unsaturated media conductivity changes because water vapour diffusion contributes to heat propagation,
this way air conductivity can change by 1 order of magnitude.
I'd like to check whether or not my experimental data are confirmed by numerics without considering vapour flow.
The output of interest is only temperature at given control points, no flow output is needed.

Cheers,
Nicolò

···

Il 06/03/2014 10.02, Marc Walther ha scritto:

Hi Nicolò,

Richards-Flow only considers the fluid, not the air. If you want air and
fluid to be assessed independently, you will need to take a look at some
multiphase/gas flow application.
Also, we generally assume local thermal equilibrium (see this eg.:
Thermodynamic equilibrium - Wikipedia ) in heat transport. I know, that Thomas Nagel added an option for non-equilibrium, if you need this.

I assume, you first want to change density or viscosity. For this, you
can take a look at models 4 (both fluid properties) in the
rf_mfp_new.cpp:

https://svn.ufz.de/ogs/browser/trunk/sources/FEM/rf_mfp_new.cpp#L918

https://svn.ufz.de/ogs/browser/trunk/sources/FEM/rf_mfp_new.cpp#L1459

and add these to your .mfp file.

Concerning the .rfd file, I know that it should work, but I am not sure
how - sorry.

Best,
Marc

On Mi, 2014-03-05 at 08:10 -0800, Nicolò Giordano wrote:

Return-Path: <ogs-users
[email protected]>
Received: from server-30 (server-30-tg [10.30.8.30]) by cyrus-sol
(Cyrus
  v2.3.16) with LMTPSA (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA
  bits=256/256 verify=YES); Wed, 05 Mar 2014 17:10:17 +0100
X-Sieve: CMU Sieve 2.3
Received: from mailin1.zih.tu-dresden.de ([141.30.67.64]) by server-30
with
  esmtp (Exim 4.76) (envelope-from
  <[email protected]>) id
  1WLEOf-0005cu-KZ for [email protected]; Wed, 05 Mar
2014
  17:10:17 +0100
Received: from mail-ve0-f189.google.com ([209.85.128.189]) by
  mailin1.zih.tu-dresden.de with esmtp (Exim 4.63) (envelope-from
  <[email protected]>) id
  1WLEOT-0007S2-Uu for [email protected]; Wed, 05 Mar 2014
17:10:08
  +0100
Received: by mail-ve0-f189.google.com with SMTP id jx11sf302776veb.6
for
  <[email protected]>; Wed, 05 Mar 2014 08:10:04 -0800 (PST)
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
d=googlegroups.com;
  s=20120806;
  h=date:from:to:message-id:subject:mime-version:x-original-sender
  :reply-to:precedence:mailing-list:list-id:list-post:list-help
  :list-archive:sender:list-subscribe:list-unsubscribe:content-type;
  bh=uA74G+d+yrdYvydbad5NQvfTvP4tmglTCqmOHXogOQI=;
  b=MkWn9RIMZyEy+tQd3CLQD9/SqBGz6xc4tyXb0a+agHPKSbN4+OwEMinFKmpAbAIjxG
  CbpRr4IFlDZKaZgQXgBXnEr/LzoFaD0eBF7XsD7jrBMygZba0x51ZiQBDE+KMgsHJlgt
  THeJywVVbtMCd6IKArvEOTWFtZ78EotUOvAha+Y6+RVnd27oyUpGUILNb8+x1KjbYx5J
  dnhqCHOSB3kx8HO2lwEYsvCSuXccyfw1u437da8lsTefskrDtp+3N5K4QcAzo1nnHbo7
  xEs0ouoYtNvWHRNjC/gpkXbYsLW7aY5crkSlF8NUXkgWra4zdTEfBOyAnfbTcVWxTHJy
ZXoQ==
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com;
  s=20120113;
  h=date:from:to:message-id:subject:mime-version:x-original-sender
  :reply-to:precedence:mailing-list:list-id:list-post:list-help
  :list-archive:sender:list-subscribe:list-unsubscribe:content-type;
  bh=uA74G+d+yrdYvydbad5NQvfTvP4tmglTCqmOHXogOQI=;
  b=VV2vOi3pUlbB145FRtnATlPuoHEf8SauwHfksbly65ZGPkIxPUL//tVPpTJzYUnfHq
  vkKBx6ZT0wOSmlh0cYeRGyQB4HKo1UOQrohszA+ZM9dwTpze1k4o2mJEII7K8z4iRITz
  G4E3AaRX87awGV7FnaBZKRofCqmzyu6I30rxSkZKoGvs8EKTHNltzqCPYoch/K7JZMV/
  V042l2iY+lFbK2WYKQ6K4l35W2YSj1nSBOByg5ECHRR9uxkNNz3IOSQceoF5q7NwjclK
  216JHE2IE9pK3ruE1HbIQtWH7ak7P0jq+F6v1cr5/w//S12XhcpViRrJf5L3LvbqgoUk
v+Pg==
X-Received: by 10.140.108.6 with SMTP id
i6mr143009qgf.0.1394035804772;
  Wed, 05 Mar 2014 08:10:04 -0800 (PST)
X-BeenThere: [email protected]
Received: by 10.140.91.120 with SMTP id y111ls129658qgd.0.gmail; Wed,
05
  Mar 2014 08:10:04 -0800 (PST)
X-Received: by 10.140.27.43 with SMTP id
40mr20349qgw.27.1394035804425;
  Wed, 05 Mar 2014 08:10:04 -0800 (PST)
Date: Wed, 5 Mar 2014 08:10:03 -0800 (PST)
From: =?UTF-8?Q?Nicol=C3=B2_Giordano?= <[email protected]>
To: [email protected]
Message-Id: <[email protected]>
Subject: [ogs-users] Temperature-dependent fluid properties
MIME-Version: 1.0
X-Original-Sender: [email protected]
Reply-To: [email protected]
Precedence: list
Mailing-list: list [email protected]; contact
  [email protected]
List-ID: <ogs-users.googlegroups.com>
X-Google-Group-Id: 896866772606
List-Post: <http://groups.google.com/group/ogs-users/post&gt;,
  <mailto:[email protected]>
List-Help: <Google Groups Help,
  <mailto:[email protected]>
List-Archive: <http://groups.google.com/group/ogs-users&gt;
Sender: [email protected]
List-Subscribe: <http://groups.google.com/group/ogs-users/subscribe&gt;,
  <mailto:[email protected]>
List-Unsubscribe:
<http://groups.google.com/group/ogs-users/subscribe&gt;,
  <mailto:googlegroups-manage+896866772606
[email protected]>
Content-Type: multipart/alternative;
X-TUD-Virus-Scanned: mailin1.zih.tu-dresden.de
X-TUD-Spam-Status: No, hits=0.0 required=5 tests=[HTML_MESSAGE=0.001]
X-TUD-Spam-Level: /
X-Evolution-Source: 1389691061.2636.2@envinf133l

------=_Part_830_5760096.1394035803732
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: quoted-printable

Dear all,
I'm using coupled Richards Flow and Heat transport for modeling
heat=20
propagation in a partially saturated porous medium. =20
I need OGS to change fluid properties (both water and air) in function
of=
=20
temperature during the simulation.
Does anyone know how to do that?
May I write a curve f(T, property value) in the .rfd file? or it works
only=
=20
for time-dependent curves?

thanks a lot
Nicol=C3=B2

--=20
You received this message because you are subscribed to the Google
Groups "=
ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an e=
mail to [email protected].
For more options, visit https://groups.google.com/groups/opt_out\.

------=_Part_830_5760096.1394035803732
Content-Type: text/html; charset=UTF-8
Content-Transfer-Encoding: quoted-printable

<div dir=3D"ltr">Dear all,<br>I'm using coupled Richards Flow and Heat
tran=
sport for modeling heat propagation in a partially saturated porous
medium.=
&nbsp; <br>I need OGS to change fluid properties (both water and air)
in fu=
nction of temperature during the simulation.<br>Does anyone know how
to do =
that?<br>May I write a curve f(T, property value) in the .rfd file? or
it w=
orks only for time-dependent curves?<br><br>thanks a
lot<br>Nicol=C3=B2<br>=
</div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google
Groups &=
quot;ogs-users&quot; group.<br />
To unsubscribe from this group and stop receiving emails from it, send
an e=
mail to <a href=3D"mailto:ogs-users
[email protected]">ogs-users=
[email protected]</a>.<br />
For more options, visit <a
href=3D"https://groups.google.com/groups/opt_out=
">https://groups.google.com/groups/opt_out&lt;/a&gt;\.&lt;br />

------=_Part_830_5760096.1394035803732--

3

Dear Nicolò,

in the benchmark
\trunk\benchmarks\THM\deco_2015_b2s1\

you will find an example, how temperature dependent viscosity, density and compressibility (not activitated in this BM but explained principally) can be used within Richards flow.

Best regards,

Jobst

···

2014-03-06 11:08 GMT+01:00 Nicolò Giordano [email protected]:

Thanks for quick reply Marc,

I will look for and take a look to Thomas Nagel option.

But actually I’d like to change thermal conductivity only because I’m working on T range 20-80 °C.

In unsaturated media conductivity changes because water vapour diffusion contributes to heat propagation,

this way air conductivity can change by 1 order of magnitude.

I’d like to check whether or not my experimental data are confirmed by numerics without considering vapour flow.

The output of interest is only temperature at given control points, no flow output is needed.

Cheers,

Nicolò

Il 06/03/2014 10.02, Marc Walther ha scritto:

Hi Nicolò,

Richards-Flow only considers the fluid, not the air. If you want air and

fluid to be assessed independently, you will need to take a look at some

multiphase/gas flow application.

Also, we generally assume local thermal equilibrium (see this eg.:

http://en.wikipedia.org/wiki/Thermodynamic_equilibrium#Local_and_global_equilibrium ) in heat transport. I know, that Thomas Nagel added an option for non-equilibrium, if you need this.

I assume, you first want to change density or viscosity. For this, you

can take a look at models 4 (both fluid properties) in the

rf_mfp_new.cpp:

https://svn.ufz.de/ogs/browser/trunk/sources/FEM/rf_mfp_new.cpp#L918

https://svn.ufz.de/ogs/browser/trunk/sources/FEM/rf_mfp_new.cpp#L1459

and add these to your .mfp file.

Concerning the .rfd file, I know that it should work, but I am not sure

how - sorry.

Best,

Marc

On Mi, 2014-03-05 at 08:10 -0800, Nicolò Giordano wrote:

Return-Path: <ogs-users

[email protected]>

Received: from server-30 (server-30-tg [10.30.8.30]) by cyrus-sol

(Cyrus

v2.3.16) with LMTPSA (version=TLSv1/SSLv3 cipher=DHE-RSA-AES256-SHA

bits=256/256 verify=YES); Wed, 05 Mar 2014 17:10:17 +0100

X-Sieve: CMU Sieve 2.3

Received: from mailin1.zih.tu-dresden.de ([141.30.67.64]) by server-30

with

esmtp (Exim 4.76) (envelope-from

[email protected]) id

1WLEOf-0005cu-KZ for [email protected]; Wed, 05 Mar

2014

17:10:17 +0100

Received: from mail-ve0-f189.google.com ([209.85.128.189]) by

mailin1.zih.tu-dresden.de with esmtp (Exim 4.63) (envelope-from

[email protected]) id

1WLEOT-0007S2-Uu for [email protected]; Wed, 05 Mar 2014

17:10:08

+0100

Received: by mail-ve0-f189.google.com with SMTP id jx11sf302776veb.6

for

[email protected]; Wed, 05 Mar 2014 08:10:04 -0800 (PST)

DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;

d=googlegroups.com;

s=20120806;

h=date:from:to:message-id:subject:mime-version:x-original-sender

:reply-to:precedence:mailing-list:list-id:list-post:list-help

:list-archive:sender:list-subscribe:list-unsubscribe:content-type;

bh=uA74G+d+yrdYvydbad5NQvfTvP4tmglTCqmOHXogOQI=;

b=MkWn9RIMZyEy+tQd3CLQD9/SqBGz6xc4tyXb0a+agHPKSbN4+OwEMinFKmpAbAIjxG

CbpRr4IFlDZKaZgQXgBXnEr/LzoFaD0eBF7XsD7jrBMygZba0x51ZiQBDE+KMgsHJlgt

THeJywVVbtMCd6IKArvEOTWFtZ78EotUOvAha+Y6+RVnd27oyUpGUILNb8+x1KjbYx5J

dnhqCHOSB3kx8HO2lwEYsvCSuXccyfw1u437da8lsTefskrDtp+3N5K4QcAzo1nnHbo7

xEs0ouoYtNvWHRNjC/gpkXbYsLW7aY5crkSlF8NUXkgWra4zdTEfBOyAnfbTcVWxTHJy

ZXoQ==

DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com;

s=20120113;

h=date:from:to:message-id:subject:mime-version:x-original-sender

:reply-to:precedence:mailing-list:list-id:list-post:list-help

:list-archive:sender:list-subscribe:list-unsubscribe:content-type;

bh=uA74G+d+yrdYvydbad5NQvfTvP4tmglTCqmOHXogOQI=;

b=VV2vOi3pUlbB145FRtnATlPuoHEf8SauwHfksbly65ZGPkIxPUL//tVPpTJzYUnfHq

vkKBx6ZT0wOSmlh0cYeRGyQB4HKo1UOQrohszA+ZM9dwTpze1k4o2mJEII7K8z4iRITz

G4E3AaRX87awGV7FnaBZKRofCqmzyu6I30rxSkZKoGvs8EKTHNltzqCPYoch/K7JZMV/

V042l2iY+lFbK2WYKQ6K4l35W2YSj1nSBOByg5ECHRR9uxkNNz3IOSQceoF5q7NwjclK

216JHE2IE9pK3ruE1HbIQtWH7ak7P0jq+F6v1cr5/w//S12XhcpViRrJf5L3LvbqgoUk

v+Pg==

X-Received: by 10.140.108.6 with SMTP id

i6mr143009qgf.0.1394035804772;

Wed, 05 Mar 2014 08:10:04 -0800 (PST)

X-BeenThere: [email protected]

Received: by 10.140.91.120 with SMTP id y111ls129658qgd.0.gmail; Wed,

05

Mar 2014 08:10:04 -0800 (PST)

X-Received: by 10.140.27.43 with SMTP id

40mr20349qgw.27.1394035804425;

Wed, 05 Mar 2014 08:10:04 -0800 (PST)

Date: Wed, 5 Mar 2014 08:10:03 -0800 (PST)

From: =?UTF-8?Q?Nicol=C3=B2_Giordano?= [email protected]

To: [email protected]

Message-Id: [email protected]

Subject: [ogs-users] Temperature-dependent fluid properties

MIME-Version: 1.0

X-Original-Sender: [email protected]

Reply-To: [email protected]

Precedence: list

Mailing-list: list [email protected]; contact

[email protected]

List-ID: <ogs-users.googlegroups.com>

X-Google-Group-Id: 896866772606

List-Post: <http://groups.google.com/group/ogs-users/post>,

mailto:[email protected]

List-Help: <http://groups.google.com/support/>,

mailto:[email protected]

List-Archive: <http://groups.google.com/group/ogs-users>

Sender: [email protected]

List-Subscribe: <http://groups.google.com/group/ogs-users/subscribe>,

mailto:[email protected]

List-Unsubscribe:

<http://groups.google.com/group/ogs-users/subscribe>,

<mailto:googlegroups-manage+896866772606

[email protected]>

Content-Type: multipart/alternative;

X-TUD-Virus-Scanned: mailin1.zih.tu-dresden.de

X-TUD-Spam-Status: No, hits=0.0 required=5 tests=[HTML_MESSAGE=0.001]

X-TUD-Spam-Level: /

X-Evolution-Source: 1389691061.2636.2@envinf133l

------=_Part_830_5760096.1394035803732

Content-Type: text/plain; charset=UTF-8

Content-Transfer-Encoding: quoted-printable

Dear all,

I’m using coupled Richards Flow and Heat transport for modeling

heat=20

propagation in a partially saturated porous medium. =20

I need OGS to change fluid properties (both water and air) in function

of=

=20

temperature during the simulation.

Does anyone know how to do that?

May I write a curve f(T, property value) in the .rfd file? or it works

only=

=20

for time-dependent curves?

thanks a lot

Nicol=C3=B2

–=20

You received this message because you are subscribed to the Google

Groups "=

ogs-users" group.

To unsubscribe from this group and stop receiving emails from it, send

an e=

mail to [email protected].

For more options, visit https://groups.google.com/groups/opt_out.

------=_Part_830_5760096.1394035803732

Content-Type: text/html; charset=UTF-8

Content-Transfer-Encoding: quoted-printable

Dear all,
I'm using coupled Richards Flow and Heat

tran=

sport for modeling heat propagation in a partially saturated porous

medium.=

 
I need OGS to change fluid properties (both water and air)

in fu=

nction of temperature during the simulation.
Does anyone know how

to do =

that?
May I write a curve f(T, property value) in the .rfd file? or

it w=

orks only for time-dependent curves?

thanks a

lot
Nicol=C3=B2
=

–

You received this message because you are subscribed to the Google

Groups &=

quot;ogs-users" group.

To unsubscribe from this group and stop receiving emails from it, send

an e=

mail to <a href=3D"mailto:ogs-users

[email protected]">ogs-users=

[email protected].

For more options, visit <a

href=3D"https://groups.google.com/groups/opt_out=

">https://groups.google.com/groups/opt_out.

------=_Part_830_5760096.1394035803732–

–

You received this message because you are subscribed to the Google Groups “ogs-users” group.

To unsubscribe from this group and stop receiving emails from it, send an email to [email protected].

For more options, visit https://groups.google.com/groups/opt_out.