direct assignment of source terms to nodes

23 views
Skip to first unread message

cb

unread,
Dec 8, 2016, 10:09:49 AM12/8/16
to ogs-...@googlegroups.com
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: christo...@gpi.uni-kiel.de
home: www.uni-kiel.de/ghm
_______________________________________

Marc Walther

unread,
Dec 8, 2016, 10:23:45 AM12/8/16
to ogs-...@googlegroups.com, Christoph Beyer
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
--


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

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

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

Wenqing Wang

unread,
Dec 9, 2016, 4:16:12 AM12/9/16
to ogs-...@googlegroups.com
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

Thomas Fischer

unread,
Dec 12, 2016, 12:46:46 AM12/12/16
to ogs-...@googlegroups.com
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. Thomas Fischer

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

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

cb

unread,
Dec 12, 2016, 2:57:07 AM12/12/16
to ogs-...@googlegroups.com
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
>
>

--
Reply all
Reply to author
Forward
0 new messages