Output Flux through Boundary Surface

22 views
Skip to first unread message

Wenjie Xu

unread,
Apr 7, 2016, 9:30:06 PM4/7/16
to DevOGS

Dear All,

 

does anybody know how to output the mass flux through a boundary surface in OGS5?

 

Cheers

 

Wenjie

Norihiro Watanabe

unread,
Apr 8, 2016, 12:34:50 AM4/8/16
to ogs-...@googlegroups.com, Wenjie Xu
Setting $DAT_TYPE TOTAL_FLUX in OUT file might work, but I'm not sure. maybe someone from Kiel can help you.

best,
nori
> --
> You received this message because you are subscribed to the Google
> Groups "ogs-devs" group.
> To unsubscribe from this group and stop receiving emails from it, send
> an email to ogs-devs+u...@googlegroups.com
> <mailto:ogs-devs+u...@googlegroups.com>.
> For more options, visit https://groups.google.com/d/optout.

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


Marc Walther

unread,
Apr 8, 2016, 2:08:00 AM4/8/16
to ogs-...@googlegroups.com

You can also use Paraview to slice a plane of your domain and then use a filter called "Surface Flow" (or similar, I don't remember exactly).

Best,
Marc

- - - - - - - - - - - - - - - - - -

Sent from my mobile device, please excuse spelling mistakes.

Marc Walther
Hydrologist
JProf Dr. rer. nat.

Helmholtz Centre for Environmental Research - UFZ
Department of Environmental Informatics
Permoserstraße 15, 04318 Leipzig, Germany

Technische Universität Dresden
Faculty of Environmental Sciences
Department of Hydrosciences
Institute for Groundwater Management
01062 Dresden

Visiting Address UFZ
Permoserstraße 15, 04318 Leipzig
Building 7.2, Room E110

Visiting Address TUD
Bergstr. 66, 01069 Dresden
Neubau Chemie 2. Bauabschnitt, CHE-BA2, Room E60

Phone +49 341 235 10 54
Fax +49 341 235 19 39
Mobile +49 178 334 18 19

Email marc.w...@ufz.de
Web http:/​/​www​.​ufz​.​de
Personal http:/​/​www​.​ufz​.​de/​index​.​php​?​en=​32239
ResGate http:/​/​www​.​researchgate​.​net/​profile/​Marc_Walther

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ

Registered Office: Leipzig
Registration Office: Amtsgericht Leipzig
Trade Register Nr. B 4703
Chairman of the Supervisory Board: MinR Wilfried Kraus
Scientific Director: Prof. Dr. Georg Teutsch
Administrative Director: Dr. Heike Graßmann

- - - - - - - - - - - - - - - - - -

--
You received this message because you are subscribed to the Google Groups "ogs-devs" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ogs-devs+u...@googlegroups.com.

Christof Beyer

unread,
Apr 8, 2016, 2:19:34 AM4/8/16
to ogs-...@googlegroups.com
Hi Wenjie,

as Nori said, you can use the TOTAL_FLUX output for fluxes across
surfaces (or polylines in 2D):

#OUTPUT
$PCS_TYPE
HEAT_TRANSPORT
$GEO_TYPE
SURFACE SURF1
$DAT_TYPE
TOTAL_FLUX
$TIM_TYPE
STEPS 1

This works for heat & mass transport, liquid & groundwater flow.
Just set the correct $PCS_TYPE. For some boundary element types (tri,
tet), however, accuracy may not be perfect.

You can also print total heat, or mass stored in your model, e.g. for a
particular material group (MatGroup 0 in this case):

#OUTPUT
$PCS_TYPE
HEAT_TRANSPORT
$DAT_TYPE
CONTENT 0
$TIM_TYPE
STEPS 1

These functions are implemented in the KielTesting branch:
https://svn.ufz.de/svn/ogs/branches/kiel_testing

Jens Olaf Delfs can give you more information, as he implemented these
routines.

Schönen Tag noch, Christof
--
_______________________________________
Dr. Christof Beyer
Institute of Geosciences
Geohydromodelling
Christian-Albrechts-University Kiel
Ludewig-Meyn-Str. 10
24118 Kiel
Germany

phone: +49(0)431-8803172
fax: +49(0)431-8807606
mobile: +49(0)176-24297908
email: christo...@gpi.uni-kiel.de
home: http://www.ifg.uni-kiel.de/
_______________________________________

jens-olaf Delfs

unread,
Apr 8, 2016, 5:05:05 AM4/8/16
to ogs-...@googlegroups.com

--

Xuan Yu

unread,
Jan 25, 2018, 4:09:00 PM1/25/18
to ogs-devs
I just tried today with triangle mesh, it did not work. And I found a webpage here http://tk-opengeosys.blogspot.com/2016/04/ says that "there seems to be a bug for trinagles and tetrahedra"

Marc Walther

unread,
Jan 25, 2018, 6:01:14 PM1/25/18
to ogs-...@googlegroups.com, xuan...@gmail.com

Dear Xuan,

I suggest to use Paraview for getting a budget through any arbitrary planar surface. You can use the filters like this:

1) 'slice' (define your area where you want to know the flux),
2) if you have a 2d setup add 'linear extrusion' with length of unit 1 into the missing dimension here,
3) 'generate normals' on the new surface,
4) use a 'calculator' filter to compute the flux per node (type 'velocity.normals' to computer a scalar product, both variables should be vectors),
5) and finally 'integrate variables' to get the total flux through the area

Best,
Marc

- - - - - - - - - - - - - - - - - -

Sent from my mobile device, please excuse spelling mistakes.

Marc Walther
Hydrologist
JProf. Dr. rer. nat.

Helmholtz Centre for Environmental Research - UFZ
Department of Environmental Informatics
Permoserstraße 15, 04318 Leipzig, Germany

Technische Universität Dresden
Faculty of Environmental Sciences
Department of Hydrosciences
Institute for Groundwater Management
01062 Dresden

Visiting Address UFZ
Permoserstraße 15, 04318 Leipzig
Building 7.2, Room E110

Visiting Address TUD
Bergstr. 66, 01069 Dresden

Neubau Chemie 2. Bauabschnitt, CHE-BA2, Room P204


Web http:/​/​www​.​ufz​.​de
Personal http:/​/​www​.​ufz​.​de/​index​.​php​?​en=​32239
ResGate http:/​/​www​.​researchgate​.​net/​profile/​Marc_Walther

Helmholtz-Zentrum für Umweltforschung GmbH - UFZ
Registered Office: Leipzig
Registration Office: Amtsgericht Leipzig
Trade Register Nr. B 4703

Chairman of the Supervisory Board: MinDirig Wilfried Kraus
Scientific Director: Prof. Dr. Dr. hc. Georg Teutsch
Administrative Director: Prof. Dr. Heike Graßmann

- - - - - - - - - - - - - - - - - -

Reply all
Reply to author
Forward
0 new messages