Skip to content

Draft: Apply the known solutions to the initial condition for the Picard scheme

wenqing requested to merge wenqing/ogs:apply_the_known_solutions_to_init into master

This MR presents a change to apply the Dirichlet boundary conditions to the initial solutions for the Picard scheme. Without this change, the first Picard iteration of the first time step may fails for some problems due to the large difference between the solutions of the previous and the current Picard iteration. In other words, the change enhances the convergence of the Picard iterations. This MR is a part of MR !3095 (closed).

Due to the changes, the results of all benchmarks that use the Picard method are changed. I have checked the results of the affected benchmarks via ParaView against the old reference results, and found the current results look reasonable. For the affected benchmarks that have slight changes in results, the comparison tolerance were increased. For the affected benchmarks that have large changes in results (in the initial or first several time steps), the reference result files were replaced with the present result files.

The affected benchmarks are

	  9 - ogs-2D_ComponentTransport_ConcentrationDiffusionAndStorage-time-vtkdiff (Failed)
	 11 - ogs-2D_ComponentTransport_StaggeredScheme_ConcentrationDiffusionAndStorage-time-vtkdiff (Failed)
	 13 - ogs-2D_ComponentTransport_DiffusionAndStorageAndAdvection-time-vtkdiff (Failed)
	 15 - ogs-2D_ComponentTransport_StaggeredScheme_DiffusionAndStorageAndAdvection-time-vtkdiff (Failed)
	 25 - ogs-2D_ComponentTransport_DiffusionAndStorageAndAdvectionAndDispersion-time-vtkdiff (Failed)
	 27 - ogs-2D_ComponentTransport_OpenBoundaryBC-time-vtkdiff (Failed)
	 29 - ogs-2D_ComponentTransport_Advective_and_NonAdvective_comparison-time-vtkdiff (Failed)
	 31 - ogs-2D_ComponentTransport_StaggeredScheme_DiffusionAndStorageAndAdvectionAndDispersion-time-vtkdiff (Failed)
	 33 - ogs-2D_MultiComponentTransport_StaggeredScheme_DiffusionAndStorageAndAdvectionAndDispersion-time-vtkdiff (Failed)
	 35 - ogs-2D_ComponentTransport_DiffusionAndStorageAndAdvectionAndDecay-time-vtkdiff (Failed)
	 37 - ogs-2D_ComponentTransport_StaggeredScheme_DiffusionAndStorageAndAdvectionAndDecay-time-vtkdiff (Failed)
	 39 - ogs-2D_ComponentTransport_DiffusionAndStorageAndAdvectionAndDispersionHalf-time-vtkdiff (Failed)
	 41 - ogs-2D_ComponentTransport_StaggeredScheme_DiffusionAndStorageAndAdvectionAndDispersionHalf-time-vtkdiff (Failed)
	 43 - ogs-3D_ComponentTransport_surfaceflux-time-vtkdiff (Failed)
	 47 - ogs-3D_StaggeredScheme_ComponentTransport_surfaceflux-time-vtkdiff (Failed)
	 59 - ogs-ComponentTransport_ConTracer_1d-time-vtkdiff (Failed)
	 61 - ogs-ComponentTransport_ConTracer_2d-time-vtkdiff (Failed)
	 62 - ogs-Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcite (Failed)
	 63 - ogs-Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac (Failed)
	 64 - ogs-Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_flag_formula (Failed)
	 67 - ogs-1D_ReactiveMassTransport_Wetland-time-vtkdiff (Failed)
	 75 - ogs-HT_SimpleSynthetics_PressureParabolicTemperatureParabolic-time-vtkdiff (Failed)
	 77 - ogs-HT_SimpleSynthetics_CoupledPressureParabolicTemperatureParabolic-time-vtkdiff (Failed)
	 85 - ogs-HT_a_DECOVALEX_THMC_based_Example-time-vtkdiff (Failed)
	 93 - ogs-HT_SimpleSynthetics_PressureParabolicTemperatureParabolicStaggered-time-vtkdiff (Failed)
	 95 - ogs-HT_SimpleSynthetics_CoupledPressureParabolicTemperatureParabolicStaggered-time-vtkdiff (Failed)
	 97 - ogs-HT_FaultedCube_rev0-time-vtkdiff (Failed)
	 98 - ogs-Parabolic/HT/SimpleSynthetics/deactivated_subdomain/HT_DeactivatedSubdomain (Failed)
	100 - ogs-1D_HeatConduction_dirichlet-vtkdiff (Failed)
	139 - ogs-HeatTransportBHE_1U_3D_beier_sandbox-time-vtkdiff (Failed)
	141 - ogs-HeatTransportBHE_coaxial_pipe_3D_deep_BHE_CXA-time-vtkdiff (Failed)
	143 - ogs-HeatTransportBHE_coaxial_pipe_3D_deep_BHE_CXC-time-vtkdiff (Failed)
	147 - ogs-HeatTransportBHE_3D_BHE_groundwater_advection-time-vtkdiff (Failed)
	226 - ogs-LiquidFlow_Dupuit_BC_BC_Storage-time-vtkdiff (Failed)
	228 - ogs-LiquidFlow_TimeDependentHeterogeneousDirichletBCs-time-vtkdiff (Failed)
	230 - ogs-LiquidFlow_TimeDependentHeterogeneousSourceTerm-time-vtkdiff (Failed)
	232 - ogs-LiquidFlow_Flux_3D_Hex-time-vtkdiff (Failed)
	234 - ogs-LiquidFlow_Flux_3D_Pyramid-time-vtkdiff (Failed)
	236 - ogs-LiquidFlow_Flux_2D_Quads-time-vtkdiff (Failed)
	238 - ogs-LiquidFlow_Flux_2D_Tris-time-vtkdiff (Failed)
	240 - ogs-Parabolic/LiquidFlow/SimpleSynthetics/PrimaryVariableConstraintDirichletBC/cuboid_1x1x1_hex_1000_Dirichlet_Dirichlet_2 (Failed)
	269 - ogs-LIE_M_cohesive_zone_modeI-time-vtkdiff (Failed)
	280 - ogs-Parabolic/RichardsComponentTransport/Padilla/Padilla_NaCl1/Padilla_NaCl1 (Failed)
	281 - ogs-Parabolic/RichardsComponentTransport/Padilla/Padilla_NaCl1/Padilla_NaCl1_quadratic (Failed)
	282 - ogs-Parabolic/RichardsComponentTransport/Padilla/Padilla_NaCl6/Padilla_NaCl6 (Failed)
	290 - ogs-2D_RichardsFlow_h_us_quad_small_PID_adaptive_dt-vtkdiff (Failed)

Merge request reports