axisymmetric diffusion example

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion, or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

···

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Email: georg.kosakowski@psi.ch
homepage: http://www.psi.ch/les/georg-kosakowski

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough. X coordinate given in the mesh corresponds to a radial distance from origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

···

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion, or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch
Administrativer Geschäftsführer: Dr. Heike Graßmann

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a little
bit irritated that the books describes benchmarks which were obviously
removed or never added....That is not good.
Greetings, Georg

···

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion, or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Email: georg.kosakowski@psi.ch
homepage: http://www.psi.ch/les/georg-kosakowski

it seems the benchmark corresponds to benchmarks/C/diffusion/Diff_HTO_test. You can find a word "radial" in the RFD file. BUT there is no $AXISYMMETRY keyword written in its mesh file, which means OGS simulates it as a 1d problem. I guess someone forgot to add the keyword in the mesh file when current mesh file format was introduced.

n

···

On 05/13/2014 01:59 PM, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a little
bit irritated that the books describes benchmarks which were obviously
removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion, or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
Norihiro Watanabe, Dr.-Ing.
Department of Environmental Informatics (ENVINF)
Wissenschaftler

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Helmholtz Centre for Environmental Research GmbH - UFZ
Permoserstraße 15 / 04318 Leipzig / Germany

norihiro.watanabe@ufz.de / http://www.ufz.de
Telefon +49 341 235 1806

Sitz der Gesellschaft: Leipzig
Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703
Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus
Wissenschaftlicher Geschäftsführer: Prof. Dr. Georg Teutsch
Administrativer Geschäftsführer: Dr. Heike Graßmann

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One thing, elem->getVolume() for symmetrical element returns area of 2D cross section. Bastian will fix the volume calculation for symmetrical elements. The current volume calculation does affect the correction of results.

Best,

Wenqing

···

On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a little
bit irritated that the books describes benchmarks which were obviously
removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion, or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

Hi Wenqing,
your remark on the symmetric elements I do not completely understand. I
would like to use a 1D setup with linear elements. I guess that linear
elements are symmetrical, or?
I know I also use elem->getVolume() in the OGS-GEM coupling for
calculation of "node volumes"....I guess that would give now wrong
results, correct?

I guess it would be more save to create a 2d triangle mesh that
aproximates a segment of circle, or?
Greetings, georg

···

On 05/13/2014 02:11 PM, Wenqing Wang wrote:

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One
thing, elem->getVolume() for symmetrical element returns area of 2D
cross section. Bastian will fix the volume calculation for
symmetrical elements. The current volume calculation does affect the
correction of results.

Best,

Wenqing
On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a little
bit irritated that the books describes benchmarks which were obviously
removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find
this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion,
or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Email: georg.kosakowski@psi.ch
homepage: http://www.psi.ch/les/georg-kosakowski

Hi Georg,

For 1D element with asymmetrical assumption, the return value of elem->getVolume() should multiplied with 2*pi*r in the current version. You can try to correct 1D of that in elem-> computeVolume() while remove if(axisymmetry) conditions in calcLumpedMass() for your test.

For your case, 1D or 2D dose not matter except the correct volume calculation, if which is needed somewhere else.

Best,

Wenqing

···

On 13.05.2014 14:25, Georg Kosakowski wrote:

Hi Wenqing,
your remark on the symmetric elements I do not completely understand. I
would like to use a 1D setup with linear elements. I guess that linear
elements are symmetrical, or?
I know I also use elem->getVolume() in the OGS-GEM coupling for
calculation of "node volumes"....I guess that would give now wrong
results, correct?

I guess it would be more save to create a 2d triangle mesh that
aproximates a segment of circle, or?
Greetings, georg

  On 05/13/2014 02:11 PM, Wenqing Wang wrote:

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One
thing, elem->getVolume() for symmetrical element returns area of 2D
cross section. Bastian will fix the volume calculation for
symmetrical elements. The current volume calculation does affect the
correction of results.

Best,

Wenqing
On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a little
bit irritated that the books describes benchmarks which were obviously
removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a description
of a diffusion benchmark for an axisymmetric setup. Where do I find
this
benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion,
or do
I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

Hi Georg,

I have implemented the volume correction in my branche version. I am not sure yet how long it will take to get it back to the trunk version.

Do you need the volume correction?

Best regards,

Bastian

···

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 13. Mai 2014 14:11
An: ogs-users@googlegroups.com
Betreff: Re: [ogs-users] axisymmetric diffusion example

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One thing, elem->getVolume() for symmetrical element returns area of 2D cross section. Bastian will fix the volume calculation for
symmetrical elements. The current volume calculation does affect the
correction of results.

Best,

Wenqing
On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a
little bit irritated that the books describes benchmarks which were
obviously removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a
description of a diffusion benchmark for an axisymmetric setup.
Where do I find this benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion,
or do I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
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 ogs-users+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Hoi Bastian,
I guess I will need it. I am currently setting up a correct test
example, as for the HTO benchmark the setup and the description do not
agree...
Please just tell me in which branch (and which files) you implemented
the volume correction :wink:
Thanks, Georg

···

On 05/14/2014 10:27 AM, Graupner Bastian wrote:

Hi Georg,

I have implemented the volume correction in my branche version. I am not sure yet how long it will take to get it back to the trunk version.

Do you need the volume correction?

Best regards,

Bastian

-----Ursprüngliche Nachricht-----
Von: ogs-users@googlegroups.com [mailto:ogs-users@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 13. Mai 2014 14:11
An: ogs-users@googlegroups.com
Betreff: Re: [ogs-users] axisymmetric diffusion example

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One thing, elem->getVolume() for symmetrical element returns area of 2D cross section. Bastian will fix the volume calculation for
symmetrical elements. The current volume calculation does affect the
correction of results.

Best,

Wenqing
On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint...that is what I guessed. I am "only" still a
little bit irritated that the books describes benchmarks which were
obviously removed or never added....That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a
description of a diffusion benchmark for an axisymmetric setup.
Where do I find this benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion,
or do I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And...how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg

--
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 ogs-users+unsubscribe@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--
PD Dr. Georg Kosakowski * Laboratory for Waste Management
Paul Scherrer Institut * 5232 Villigen/PSI * Switzerland
Phone: +41-56-3104743 * Email: georg.kosakowski@psi.ch
homepage: http://www.psi.ch/les/georg-kosakowski

Hello Bastian,

I need this volume correction for 1d radial heat trasnport. Can you give me more info on this?

Regards,

Chan-Hee

···

On Wednesday, May 14, 2014 5:27:30 PM UTC+9, Bastian Graupner wrote:

Hi Georg,

I have implemented the volume correction in my branche version. I am not sure yet how long it will take to get it back to the trunk version.

Do you need the volume correction?

Best regards,

Bastian

-----Ursprüngliche Nachricht-----
Von: ogs-...@googlegroups.com [mailto:ogs-...@googlegroups.com] Im Auftrag von Wenqing Wang
Gesendet: Dienstag, 13. Mai 2014 14:11
An: ogs-...@googlegroups.com

Betreff: Re: [ogs-users] axisymmetric diffusion example

Hi Georg,

Nori is right: just put that keyword at beginning of mesh file. One thing, elem->getVolume() for symmetrical element returns area of 2D cross section. Bastian will fix the volume calculation for
symmetrical elements. The current volume calculation does affect the
correction of results.

Best,

Wenqing
On 13.05.2014 13:59, Georg Kosakowski wrote:

Hi Nori,
thanks for the hint…that is what I guessed. I am “only” still a
little bit irritated that the books describes benchmarks which were
obviously removed or never added…That is not good.
Greetings, Georg

On 05/13/2014 12:33 PM, Norihiro Watanabe wrote:

Hi Georg,

yes, setting $AXISYMMETRY keyword into the msh file should be enough.
X coordinate given in the mesh corresponds to a radial distance from
origin (r=0).

you might be able to get some hints from benchmarks/T/T_1D_axi which
is an axisymmetric model for heat transport in a hollow cylinder.

Best,
Nori

On 05/12/2014 11:30 AM, Georg Kosakowski wrote:

Hi all,
in addition to the forum I ask here:

In the (first) benchmark book there is in section 10.1.1 a
description of a diffusion benchmark for an axisymmetric setup.
Where do I find this benchmark in the benchmark folder?
/benchmarks/C/diffusion/Diff_HTO_test.xxx is for linear diffusion,
or do I miss something?
Would it be enough to put the $AXISYMMETRY keyword into the msh file?
And…how does this work then? Is 0,0,0 the center of symmetry?
Thanks, Georg


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 ogs-users+...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.