Skip to content
Snippets Groups Projects
Commit 4901ea14 authored by Johannes Boog's avatar Johannes Boog
Browse files

[test] add bm files for test OGS6_msh

parent c45c43bc
No related branches found
No related tags found
1 merge request!50Add support for mesh files in MSH format
Hide whitespace changes
Inline Side-by-side
inst/extdata/benchmarks/t1_1Dsource/t1_1Dsource.gml 0 → 100755
+ 40
− 0
View file @ 4901ea14
<?xml version="1.0" encoding="utf-8"?>
<OpenGeoSysGLI>
<name>geometry</name>
<points>
<point id="0" x="100" y="0" z="0" name="POINT0"/>
<point id="1" x="100" y="10" z="0" name="POINT1"/>
<point id="2" x="100" y="10" z="10" name="POINT2"/>
<point id="3" x="100" y="0" z="10" name="POINT3"/>
<point id="4" x="0" y="0" z="0" name="POINT4"/>
<point id="5" x="0" y="10" z="0" name="POINT5"/>
<point id="6" x="0" y="10" z="10" name="POINT6"/>
<point id="7" x="0" y="0" z="10" name="POINT7"/>
</points>
<polylines>
<polyline id="0" name="POLYLINE1">
<pnt>0</pnt>
<pnt>1</pnt>
<pnt>2</pnt>
<pnt>3</pnt>
<pnt>0</pnt>
</polyline>
<polyline id="1" name="POLYLINE2">
<pnt>4</pnt>
<pnt>5</pnt>
<pnt>6</pnt>
<pnt>7</pnt>
<pnt>4</pnt>
</polyline>
</polylines>
<surfaces>
<surface id="0" name="SURFACE1">
<element p1="0" p2="2" p3="1"/>
<element p1="2" p2="0" p3="3"/>
</surface>
<surface id="1" name="SURFACE2">
<element p1="4" p2="6" p3="5"/>
<element p1="6" p2="4" p3="7"/>
</surface>
</surfaces>
</OpenGeoSysGLI>
inst/extdata/benchmarks/t1_1Dsource/t1_1Dsource.msh 0 → 100755
+ 62
− 0
View file @ 4901ea14
#FEM_MSH
$PCS_TYPE
NO_PCS
$NODES
44
0 10 10 10
1 0 10 10
2 0 0 10
3 10 0 10
4 10 10 0
5 0 10 0
6 0 0 0
7 10 0 0
8 20 10 10
9 20 0 10
10 20 10 0
11 20 0 0
12 30 10 10
13 30 0 10
14 30 10 0
15 30 0 0
16 40 10 10
17 40 0 10
18 40 10 0
19 40 0 0
20 50 10 10
21 50 0 10
22 50 10 0
23 50 0 0
24 60 10 10
25 60 0 10
26 60 10 0
27 60 0 0
28 70 10 10
29 70 0 10
30 70 10 0
31 70 0 0
32 80 10 10
33 80 0 10
34 80 10 0
35 80 0 0
36 90 10 10
37 90 0 10
38 90 10 0
39 90 0 0
40 100 10 10
41 100 0 10
42 100 10 0
43 100 0 0
$ELEMENTS
10
0 0 hex 4 5 6 7 0 1 2 3
1 0 hex 10 4 7 11 8 0 3 9
2 0 hex 14 10 11 15 12 8 9 13
3 0 hex 18 14 15 19 16 12 13 17
4 1 hex 22 18 19 23 20 16 17 21
5 1 hex 26 22 23 27 24 20 21 25
6 1 hex 30 26 27 31 28 24 25 29
7 1 hex 34 30 31 35 32 28 29 33
8 1 hex 38 34 35 39 36 32 33 37
9 1 hex 42 38 39 43 40 36 37 41
#STOP
inst/extdata/benchmarks/t1_1Dsource/t1_1Dsource.prj 0 → 100644
+ 213
− 0
View file @ 4901ea14
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh>t1_1Dsource.msh</mesh>
<geometry>t1_1Dsource.gml</geometry>
<processes>
<process>
<name>HeatConduction</name>
<type>HEAT_CONDUCTION</type>
<integration_order>2</integration_order>
<process_variables>
<process_variable>temperature</process_variable>
</process_variables>
<secondary_variables>
</secondary_variables>
</process>
</processes>
<media>
<medium id="0">
<phases/>
<properties>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>rho</parameter_name>
</property>
<property>
<name>heat_capacity</name>
<type>Parameter</type>
<parameter_name>c</parameter_name>
</property>
<property>
<name>thermal_conductivity</name>
<type>Parameter</type>
<parameter_name>lambda</parameter_name>
</property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
</properties>
</medium>
<medium id="1">
<phases/>
<properties>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>rho</parameter_name>
</property>
<property>
<name>heat_capacity</name>
<type>Parameter</type>
<parameter_name>c</parameter_name>
</property>
<property>
<name>thermal_conductivity</name>
<type>Parameter</type>
<parameter_name>lambda</parameter_name>
</property>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="HeatConduction">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1.0e-9</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial> 0.0 </t_initial>
<t_end> 1 </t_end>
<timesteps>
<pair>
<repeat> 1 </repeat>
<delta_t> 1 </delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>t1_1Dsource</prefix>
<data_mode>Ascii</data_mode>
<compress_output>false</compress_output>
<timesteps>
<pair>
<repeat> 1 </repeat>
<each_steps> 1 </each_steps>
</pair>
</timesteps>
<variables>
<variable> temperature </variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<parameter>
<name>lambda</name>
<type>Group</type>
<group_id_property>MaterialIDs</group_id_property>
<index_values>
<index>0</index>
<value>100</value>
</index_values>
<index_values>
<index>1</index>
<value>300</value>
</index_values>
</parameter>
<parameter>
<name>c</name>
<type>Group</type>
<group_id_property>MaterialIDs</group_id_property>
<index_values>
<index>0</index>
<value>0</value>
</index_values>
<index_values>
<index>1</index>
<value>0</value>
</index_values>
</parameter>
<parameter>
<name>rho</name>
<type>Group</type>
<group_id_property>MaterialIDs</group_id_property>
<index_values>
<index>0</index>
<value>0</value>
</index_values>
<index_values>
<index>1</index>
<value>0</value>
</index_values>
</parameter>
<parameter>
<name>initial_temperature</name>
<type>Constant</type>
<value>1.5</value>
</parameter>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>1.0</value>
</parameter>
<parameter>
<name>h0</name>
<type>Constant</type>
<value>1.5</value>
</parameter>
</parameters>
<process_variables>
<process_variable>
<name>temperature</name>
<components>1</components>
<order>1</order>
<initial_condition>initial_temperature</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>geometry</geometrical_set>
<geometry>SURFACE1</geometry>
<type>Neumann</type>
<parameter>h0</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>geometry</geometrical_set>
<geometry>SURFACE2</geometry>
<type>Dirichlet</type>
<parameter>T0</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>10</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i cg -p jacobi -tol 1e-10 -maxiter 10000</lis>
<eigen>
<solver_type>CG</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-10</error_tolerance>
</eigen>
<petsc>
<prefix>T</prefix>
<parameters>-T_ksp_type cg -T_pc_type bjacobi -T_ksp_rtol 1e-10 -T_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment