OGS
|
<OpenGeoSysProject>
<mesh>th_decovalex.vtu</mesh>
<geometry>boundary.gml</geometry>
<processes>
<process>
<name>th_decovalex</name>
<type>HT</type>
<coupling_scheme>staggered</coupling_scheme>
<integration_order>2</integration_order>
<process_variables>
<temperature>T</temperature>
<pressure>p</pressure>
</process_variables>
<specific_body_force>0 -9.81</specific_body_force>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity" />
</secondary_variables>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>density</name>
<type>Linear</type>
<reference_value>1000.0</reference_value>
<independent_variable>
<variable_name>temperature</variable_name>
<reference_condition>293.0</reference_condition>
<slope>-4.3e-4</slope>
</independent_variable>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4000</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>18921600.0</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>storage</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>rho_solid</parameter_name>
</property>
<property>
<name>thermal_conductivity</name>
<type>Parameter</type>
<parameter_name>lambda</parameter_name>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Parameter</type>
<parameter_name>c_p</parameter_name>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>thermal_longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>6.3072e-14 0. 0. 6.3072e-14</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.01</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>EffectiveThermalConductivityPorosityMixing</type>
</property>
</properties>
</medium>
<medium id="1">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>density</name>
<type>Linear</type>
<reference_value>1000.0</reference_value>
<independent_variable>
<variable_name>temperature</variable_name>
<reference_condition>293.0</reference_condition>
<slope>-4.3e-4</slope>
</independent_variable>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4000</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>18921600.0</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>rho_solid</parameter_name>
</property>
<property>
<name>storage</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Parameter</type>
<parameter_name>lambda</parameter_name>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Parameter</type>
<parameter_name>c_p</parameter_name>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>thermal_longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_transversal_dispersivity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4000</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>3.1536e-10 0. 0. 3.1536e-10</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.41</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>EffectiveThermalConductivityPorosityMixing</type>
</property>
</properties>
</medium>
</media>
<time_loop>
<global_process_coupling>
<max_iter>6</max_iter>
<convergence_criteria>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1.e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1.e-14</reltol>
</convergence_criterion>
</convergence_criteria>
</global_process_coupling>
<processes>
<process ref="th_decovalex">
<nonlinear_solver>basic_picard_T</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>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>1e3</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>0.001</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>0.006</delta_t>
</pair>
<pair>
<repeat>9</repeat>
<delta_t>0.01</delta_t>
</pair>
<pair>
<repeat>9</repeat>
<delta_t>0.1</delta_t>
</pair>
<pair>
<repeat>39</repeat>
<delta_t>1.</delta_t>
</pair>
<pair>
<repeat>2</repeat>
<delta_t>2.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>16.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>30.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>40.0</delta_t>
</pair>
<pair>
<repeat>10</repeat>
<delta_t>86.0</delta_t>
</pair>
<pair>
<repeat>90</repeat>
<delta_t>100.0</delta_t>
</pair>
<pair>
<repeat>100</repeat>
<delta_t>900</delta_t>
</pair>
<pair>
<repeat>90</repeat>
<delta_t>1.e4</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="th_decovalex">
<nonlinear_solver>basic_picard_H</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>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>1.e3</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>0.001</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>0.006</delta_t>
</pair>
<pair>
<repeat>9</repeat>
<delta_t>0.01</delta_t>
</pair>
<pair>
<repeat>9</repeat>
<delta_t>0.1</delta_t>
</pair>
<pair>
<repeat>39</repeat>
<delta_t>1.</delta_t>
</pair>
<pair>
<repeat>2</repeat>
<delta_t>2.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>16.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>30.0</delta_t>
</pair>
<pair>
<repeat>1</repeat>
<delta_t>40.0</delta_t>
</pair>
<pair>
<repeat>10</repeat>
<delta_t>86.0</delta_t>
</pair>
<pair>
<repeat>90</repeat>
<delta_t>100.0</delta_t>
</pair>
<pair>
<repeat>100</repeat>
<delta_t>900</delta_t>
</pair>
<pair>
<repeat>90</repeat>
<delta_t>1.e4</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>th_decovalex</prefix>
<suffix>ts{:timestep}_t_{:time}</suffix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>10</each_steps>
</pair>
</timesteps>
<variables>
<variable>T</variable>
<variable>p</variable>
<variable>darcy_velocity</variable>
</variables>
</output>
</time_loop>
<parameters>
<parameter>
<name>p0</name>
<type>Constant</type>
<value>4.7e6</value>
</parameter>
<parameter>
<name>p_bc_top</name>
<type>Constant</type>
<value>4.3e+6</value>
</parameter>
<parameter>
<name>p_bc_bottom</name>
<type>Constant</type>
<value>4.7e+6</value>
</parameter>
<parameter>
<name>lambda</name>
<type>MeshElement</type>
<field_name>scaled_solid_thermal_conductivity</field_name>
</parameter>
<parameter>
<name>lambda_fluid</name>
<type>Constant</type>
<value>18921600.0</value>
</parameter>
<parameter>
<name>c_p</name>
<type>MeshElement</type>
<field_name>solid_heat_capacity</field_name>
</parameter>
<parameter>
<name>rho_solid</name>
<type>MeshElement</type>
<field_name>solid_density</field_name>
</parameter>
<parameter>
<name>n0</name>
<type>Constant</type>
<value>0.01</value>
</parameter>
<parameter>
<name>n1</name>
<type>Constant</type>
<value>0.41</value>
</parameter>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>298</value>
</parameter>
<parameter>
<name>TB_top</name>
<type>Constant</type>
<value>294.0</value>
</parameter>
<parameter>
<name>TB_bottom</name>
<type>Constant</type>
<value>310</value>
</parameter>
<parameter>
<name>mu0</name>
<type>Constant</type>
<values>6.3072e-14 0. 0. 6.3072e-14</values>
</parameter>
<parameter>
<name>mu1</name>
<type>Constant</type>
<values>3.1536e-10 0. 0. 3.1536e-10</values>
</parameter>
<parameter>
<name>hpower_on_arc</name>
<type>Constant</type>
<value>1.0</value>
</parameter>
<parameter>
<name>heat_power</name>
<type>CurveScaled</type>
<curve>time_heat_power</curve>
<parameter>hpower_on_arc</parameter>
</parameter>
<parameter>
<name>heat_power_c</name>
<type>Constant</type>
<value>2477507279.46192</value>
</parameter>
</parameters>
<curves>
<curve>
<name>time_heat_power</name>
<coords>0 1.00E-03 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 20 30 40 50 60 70 80 90 100 200 300 400 500 600 700 800 900 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 20000 30000 40000 50000 60000 70000 80000 90000 100000 200000 300000 400000 500000 600000 700000 800000 900000 1000000 1.00E+10</coords>
<values>2477507279.46192 2477459529.54288 2477411689.43088 2477363767.32528 2477315670.51024 2477267748.40464 2477219911.13088 2477171729.484 2477123892.21024 2477075710.56336 2477027873.2896 2472672907.16208 2467757091.56544 2462762401.59456 2457693385.056 2452552352.27712 2447343250.9344 2442069507.41424 2436734637.66528 2431341468.88992 2425894034.2272 2369163847.82544 2310711582.0096 2253930470.64432 2202214690.836 2158958419.37568 2112431166.82368 2078620739.1 2045550770.03376 2013204028.0392 1726032720.05136 1495088443.47744 1308713237.89344 1157694104.67264 1034749222.15824 934122524.20848 851262379.89264 782568686.43936 725193373.44672 448276352.8464 341492110.56144 277947936.256176 236509762.562928 209106613.956816 190948269.79944 178912492.01208 154295944.435872 143838731.939184 87693573.405456 68334737.77152 58916070.8136 52907269.428288 48548275.775136 45238544.066592 42688989.27792 41088819.808944 36699991.77528 21249719.02512 14044626.9792 9640538.832816 6944041.074528 5293033.113744 4282092.927936 3663114.947568 3284183.858832 3052109.993328 1972070.615664 1660867.028928 1472515.998336 1358407.99656 1289308.804704 1247386.360032 1222095.6864 1206749.827296 1197405.142848 1197405.142848</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>T</name>
<components>1</components>
<order>1</order>
<initial_condition>T0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>th_decovalex</geometrical_set>
<geometry>top</geometry>
<type>Dirichlet</type>
<parameter>TB_top</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>th_decovalex</geometrical_set>
<geometry>bottom</geometry>
<type>Dirichlet</type>
<parameter>TB_bottom</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>th_decovalex</geometrical_set>
<geometry>arc</geometry>
<type>Neumann</type>
<parameter>heat_power</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>p</name>
<components>1</components>
<order>1</order>
<initial_condition>p0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>th_decovalex</geometrical_set>
<geometry>top</geometry>
<type>Dirichlet</type>
<parameter>p_bc_top</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>th_decovalex</geometrical_set>
<geometry>bottom</geometry>
<type>Dirichlet</type>
<parameter>p_bc_bottom</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard_T</name>
<type>Picard</type>
<max_iter>10</max_iter>
<linear_solver>linear_solver_T</linear_solver>
</nonlinear_solver>
<nonlinear_solver>
<name>basic_picard_H</name>
<type>Picard</type>
<max_iter>10</max_iter>
<linear_solver>linear_solver_H</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>linear_solver_T</name>
<lis>-i cg -p jacobi -tol 1e-20 -maxiter 10000</lis>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-20</error_tolerance>
</eigen>
<petsc>
<prefix>T</prefix>
<parameters>-T_ksp_type bcgs -T_pc_type bjacobi -T_ksp_rtol 1e-16 -T_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
<linear_solver>
<name>linear_solver_H</name>
<lis>-i cg -p jacobi -tol 1e-20 -maxiter 10000</lis>
<eigen>
<solver_type>CG</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-20</error_tolerance>
</eigen>
<petsc>
<prefix>H</prefix>
<parameters>-H_ksp_type cg -H_pc_type bjacobi -H_ksp_rtol 1e-16 -H_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>