direct assignment of source terms to nodes

Dear all,

is there a source term method implemented in OGS to assign a specific flux of a flow pcs to a list of nodes?

As far as I understand, $DIS_TYPE DIRECT will assign a CONSTANT volume flux (i.e. m³/s) to all nodes listed in a text file, but I want to assing the flux in m³/s per m² of patch area for each node, i.e. a combination o DIRECT and CONSTANT_NEUMANN.

Cheers, 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: christof.beyer@gpi.uni-kiel.de
home: www.uni-kiel.de/ghm
_______________________________________

Hey Christoph,

I don't know, if I clearly understood you, but you can use the keyword
"$DISTRIBUTE_VOLUME_FLUX" to give a volume per time on a line or
surface, which will then distribute this flux over the element. If you
have your list of nodes as individual source terms, this might work.

Cheers,
Marc

···

On Thu, 2016-12-08 at 16:09 +0100, cb wrote:
Dear all,

is there a source term method implemented in OGS to assign a specific
flux of a flow pcs to a list of nodes?

As far as I understand, $DIS_TYPE DIRECT will assign a CONSTANT volume
flux (i.e. m³/s) to all nodes listed in a text file, but I want to
assing the flux in m³/s per m² of patch area for each node, i.e. a
combination o DIRECT and CONSTANT_NEUMANN.

Cheers, 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: christof.beyer@gpi.uni-kiel.de
home: www.uni-kiel.de/ghm
_______________________________________

--

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

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
Professorship for Contaminant Hydrology
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

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

Email marc.walther@ufz.de
Web http://www.ufz.de
Personal UFZ http://www.ufz.de/index.php?en=38484
Personal TUD https://tu-dresden.de/bu/umwelt/hydro/igw/das-institut/mitarbeiterinformationen/marc-walther
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. h.c. Georg Teutsch
Administrative Director: Prof. Dr. Heike Graßmann

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

Dear Christof,

There is no direct method for this kind of BC input. Maybe you can mesh the area first, then write a short script (or c++ code) to convert the 2D elements of mesh into surfaces of gli and to assign the Neumann BC to these surfaces.

Wenqing

···

On 12/08/2016 04:09 PM, cb wrote:

Dear all,

is there a source term method implemented in OGS to assign a specific flux of a flow pcs to a list of nodes?

As far as I understand, $DIS_TYPE DIRECT will assign a CONSTANT volume flux (i.e. m³/s) to all nodes listed in a text file, but I want to assing the flux in m³/s per m² of patch area for each node, i.e. a combination o DIRECT and CONSTANT_NEUMANN.

Cheers, Christof

I think the workflow [0] could be a good starting point. Maybe the last step of the workflow needs some modification/additional work.
For the converting of the surface mesh to an geometry it is possible to use the DataExplorer.

Best regards,
Tom

[0] https://docs.opengeosys.org/docs/tools/model-preparation/computesurfacenodeidsinpolygonalregion

···

On 2016-12-09 10:16, Wenqing Wang wrote:

There is no direct method for this kind of BC input. Maybe you can mesh
the area first, then write a short script (or c++ code) to convert the
2D elements of mesh into surfaces of gli and to assign the Neumann BC to
these surfaces.

--

Dr. Thomas Fischer

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

Phone +49 341 235 10 54
Fax +49 341 235 19 39

Email thomas.fischer@ufz.de
Web http://www.ufz.de

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. Georg Teutsch
Administrative Director: Dr. Heike Graßmann

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

Hi Thomas, Wenqing and Marc,

thanks for your suggestions. This looks like a pretty nice way to assign BCs or STs.

One more question: I guess the ExtractSurface tool will extract outer surfaces of a mesh. In case a ST has to be assigned to a surface *within* the 3D mesh (e.g. a fracture or similar) it would then be necessary to cut away the surrounding elements before surface extraction. In Gina, this can be achieved using the "delete by matgroup" tool. I just found the "Remove Mesh Elements" method description, this obviously does the job...

Best regards, Christof

···

Am 12.12.2016 um 06:46 schrieb Thomas Fischer:

On 2016-12-09 10:16, Wenqing Wang wrote:

There is no direct method for this kind of BC input. Maybe you can mesh
the area first, then write a short script (or c++ code) to convert the
2D elements of mesh into surfaces of gli and to assign the Neumann BC to
these surfaces.

I think the workflow [0] could be a good starting point. Maybe the last step of the workflow needs some modification/additional work.
For the converting of the surface mesh to an geometry it is possible to use the DataExplorer.

Best regards,
Tom

[0] https://docs.opengeosys.org/docs/tools/model-preparation/computesurfacenodeidsinpolygonalregion

--
_______________________________________
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: christof.beyer@gpi.uni-kiel.de
home: www.uni-kiel.de/ghm
_______________________________________