Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
[T] GWF; Anisotropic test with local coord sys.
- Loading branch information
1 parent
6bf8029
commit baf3fe9
Showing
4 changed files
with
170 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
160 changes: 160 additions & 0 deletions
160
Tests/Data/Elliptic/cube_1x1x1_GroundWaterFlow/cube_1e4_anisotropic.prj
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,160 @@ | ||
<?xml version="1.0" encoding="ISO-8859-1"?> | ||
<OpenGeoSysProject> | ||
<mesh>cube_1x1x1_hex_1e4_material_groups.vtu</mesh> | ||
<geometry>cube_1x1x1.gml</geometry> | ||
<processes> | ||
<process> | ||
<name>GW23</name> | ||
<type>GROUNDWATER_FLOW</type> | ||
<integration_order>2</integration_order> | ||
<hydraulic_conductivity>K</hydraulic_conductivity> | ||
<process_variables> | ||
<process_variable>pressure</process_variable> | ||
</process_variables> | ||
<secondary_variables> | ||
<secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/> | ||
</secondary_variables> | ||
</process> | ||
</processes> | ||
<time_loop> | ||
<processes> | ||
<process ref="GW23"> | ||
<nonlinear_solver>basic_picard</nonlinear_solver> | ||
<convergence_criterion> | ||
<type>DeltaX</type> | ||
<norm_type>NORM2</norm_type> | ||
<abstol>1.e-6</abstol> | ||
</convergence_criterion> | ||
<time_discretization> | ||
<type>BackwardEuler</type> | ||
</time_discretization> | ||
<time_stepping> | ||
<type>SingleStep</type> | ||
</time_stepping> | ||
</process> | ||
</processes> | ||
<output> | ||
<type>VTK</type> | ||
<prefix>cube_1e4_anisotropic</prefix> | ||
<variables> | ||
<variable>pressure</variable> | ||
<variable>darcy_velocity</variable> | ||
</variables> | ||
</output> | ||
</time_loop> | ||
<local_coordinate_system> | ||
<basis_vector_0>e0</basis_vector_0> | ||
<basis_vector_1>e1</basis_vector_1> | ||
<basis_vector_2>e2</basis_vector_2> | ||
</local_coordinate_system> | ||
<parameters> | ||
<parameter> | ||
<name>e0</name> | ||
<type>Constant</type> | ||
<values>0.86602540378443864676 0 0.5</values> | ||
</parameter> | ||
<parameter> | ||
<name>e1</name> | ||
<type>Constant</type> | ||
<values>0 1 0</values> | ||
</parameter> | ||
<parameter> | ||
<name>e2</name> | ||
<type>Constant</type> | ||
<values>-0.5 0 0.86602540378443864676</values> | ||
</parameter> | ||
<parameter> | ||
<name>K</name> | ||
<type>Group</type> | ||
<group_id_property>MaterialIDs</group_id_property> | ||
<index_values> | ||
<index>0</index> | ||
<values>1 1 1</values> | ||
</index_values> | ||
<index_values> | ||
<index>1</index> | ||
<values>1 1 0.1</values> | ||
</index_values> | ||
<index_values> | ||
<index>2</index> | ||
<values>0.1 0.1 1</values> | ||
</index_values> | ||
<use_local_coordinate_system>true</use_local_coordinate_system> | ||
</parameter> | ||
<parameter> | ||
<name>p0</name> | ||
<type>Constant</type> | ||
<value>0</value> | ||
</parameter> | ||
<parameter> | ||
<name>p_Dirichlet_left</name> | ||
<type>Constant</type> | ||
<value>1</value> | ||
</parameter> | ||
<parameter> | ||
<name>p_Dirichlet_right</name> | ||
<type>Constant</type> | ||
<value>-1</value> | ||
</parameter> | ||
</parameters> | ||
<process_variables> | ||
<process_variable> | ||
<name>pressure</name> | ||
<components>1</components> | ||
<order>1</order> | ||
<initial_condition>p0</initial_condition> | ||
<boundary_conditions> | ||
<boundary_condition> | ||
<geometrical_set>cube_1x1x1_geometry</geometrical_set> | ||
<geometry>left</geometry> | ||
<type>Dirichlet</type> | ||
<parameter>p_Dirichlet_left</parameter> | ||
</boundary_condition> | ||
<boundary_condition> | ||
<geometrical_set>cube_1x1x1_geometry</geometrical_set> | ||
<geometry>right</geometry> | ||
<type>Dirichlet</type> | ||
<parameter>p_Dirichlet_right</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-16 -maxiter 10000</lis> | ||
<eigen> | ||
<solver_type>CG</solver_type> | ||
<precon_type>DIAGONAL</precon_type> | ||
<max_iteration_step>10000</max_iteration_step> | ||
<error_tolerance>1e-16</error_tolerance> | ||
</eigen> | ||
<petsc> | ||
<prefix>gw</prefix> | ||
<parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters> | ||
</petsc> | ||
</linear_solver> | ||
</linear_solvers> | ||
<test_definition> | ||
<vtkdiff> | ||
<file>cube_1e4_anisotropic_pcs_0_ts_1_t_1.000000.vtu</file> | ||
<field>pressure</field> | ||
<absolute_tolerance>1e-14</absolute_tolerance> | ||
<relative_tolerance>0</relative_tolerance> | ||
</vtkdiff> | ||
<vtkdiff> | ||
<file>cube_1e4_anisotropic_pcs_0_ts_1_t_1.000000.vtu</file> | ||
<field>darcy_velocity</field> | ||
<absolute_tolerance>1e-13</absolute_tolerance> | ||
<relative_tolerance>0</relative_tolerance> | ||
</vtkdiff> | ||
</test_definition> | ||
</OpenGeoSysProject> |
3 changes: 3 additions & 0 deletions
3
...s/Data/Elliptic/cube_1x1x1_GroundWaterFlow/cube_1e4_anisotropic_pcs_0_ts_1_t_1.000000.vtu
Git LFS file not shown
3 changes: 3 additions & 0 deletions
3
Tests/Data/Elliptic/cube_1x1x1_GroundWaterFlow/cube_1x1x1_hex_1e4_material_groups.vtu
Git LFS file not shown