22 #ifndef _aspect_parameters_h 23 #define _aspect_parameters_h 25 #include <deal.II/base/parameter_handler.h> 35 struct CompositionalFieldDescription;
38 namespace GeometryModel
81 no_Advection_no_Stokes
89 return input == NonlinearSolver::iterated_Advection_and_Newton_Stokes ||
90 input == NonlinearSolver::single_Advection_iterated_Newton_Stokes ||
91 input == NonlinearSolver::no_Advection_iterated_defect_correction_Stokes ||
92 input == NonlinearSolver::single_Advection_iterated_defect_correction_Stokes ||
93 input == NonlinearSolver::iterated_Advection_and_defect_correction_Stokes
108 net_translation_x = 0x1,
109 net_translation_y = 0x2,
110 net_translation_z = 0x4,
111 net_translation = 0x1+0x2+0x4,
112 linear_momentum_x = 0x8,
113 linear_momentum_y = 0x10,
114 linear_momentum_z = 0x20,
115 linear_momentum = linear_momentum_x+linear_momentum_y+linear_momentum_z,
117 angular_momentum = 0x80,
118 net_surface_rotation = 0x100,
119 any_translation = net_translation+linear_momentum,
120 any_rotation = net_rotation+angular_momentum+net_surface_rotation,
145 prescribed_field_with_diffusion
179 if (input ==
"isentropic compression")
180 return Formulation::isentropic_compression;
181 else if (input ==
"anelastic liquid approximation")
182 return Formulation::anelastic_liquid_approximation;
183 else if (input ==
"Boussinesq approximation")
184 return Formulation::boussinesq_approximation;
185 else if (input ==
"custom")
186 return Formulation::custom;
188 AssertThrow(
false, ExcNotImplemented());
223 if (input ==
"isentropic compression")
224 return Formulation::MassConservation::isentropic_compression;
225 else if (input ==
"hydrostatic compression")
226 return Formulation::MassConservation::hydrostatic_compression;
227 else if (input ==
"reference density profile")
228 return Formulation::MassConservation::reference_density_profile;
229 else if (input ==
"implicit reference density profile")
230 return Formulation::MassConservation::implicit_reference_density_profile;
231 else if (input ==
"incompressible")
232 return Formulation::MassConservation::incompressible;
233 else if (input ==
"projected density field")
234 return Formulation::MassConservation::projected_density_field;
235 else if (input ==
"ask material model")
236 return Formulation::MassConservation::ask_material_model;
238 AssertThrow(
false, ExcNotImplemented());
258 reference_density_profile
269 if (input ==
"real density")
270 return Formulation::TemperatureEquation::real_density;
271 else if (input ==
"reference density profile")
272 return Formulation::TemperatureEquation::reference_density_profile;
274 AssertThrow(
false, ExcNotImplemented());
300 if (input ==
"entropy viscosity")
301 return entropy_viscosity;
302 else if (input ==
"SUPG")
305 AssertThrow(
false, ExcNotImplemented());
312 return "entropy viscosity|SUPG";
331 return "block AMG|direct solver|block GMG";
337 if (input ==
"block AMG")
339 else if (input ==
"direct solver")
340 return direct_solver;
341 else if (input ==
"block GMG")
344 AssertThrow(
false, ExcNotImplemented());
364 return "GMRES|IDR(s)";
370 if (input ==
"GMRES")
372 else if (input ==
"IDR(s)")
375 AssertThrow(
false, ExcNotImplemented());
398 const MPI_Comm mpi_communicator);
426 void parse_parameters (ParameterHandler &prm,
427 const MPI_Comm mpi_communicator);
445 void parse_geometry_dependent_parameters (ParameterHandler &prm,
The NullspaceRemoval struct.
bool use_full_A_block_preconditioner
unsigned int min_grid_level
unsigned int n_expensive_stokes_solver_steps
std::vector< double > global_composition_min_preset
unsigned int composition_degree
static Kind parse(const std::string &input)
double linear_stokes_solver_tolerance
unsigned int n_cheap_stokes_solver_steps
static Kind parse(const std::string &input)
unsigned int initial_adaptive_refinement
double coarsening_fraction
std::vector< unsigned int > normalized_fields
double maximum_first_time_step
double global_temperature_max_preset
std::map< types::boundary_id, std::pair< std::string, std::string > > prescribed_traction_boundary_indicators
unsigned int initial_global_refinement
static const std::string pattern()
double AMG_aggregation_threshold
unsigned int stabilization_alpha
double linear_solver_A_block_tolerance
std::vector< double > additional_refinement_times
std::map< unsigned int, std::pair< std::string, unsigned int > > mapped_particle_properties
bool use_limiter_for_discontinuous_temperature_solution
bool use_locally_conservative_discretization
double nonlinear_tolerance
unsigned int stokes_gmres_restart_length
static Kind parse(const std::string &input)
unsigned int temperature_degree
unsigned int idr_s_parameter
bool mesh_deformation_enabled
Formulation::MassConservation::Kind formulation_mass_conservation
bool skip_solvers_on_initial_refinement
bool run_postprocessors_on_initial_refinement
bool use_discontinuous_composition_discretization
bool adapt_by_fraction_of_cells
AdvectionFieldMethod::Kind temperature_method
unsigned int AMG_smoother_sweeps
double temperature_solver_tolerance
NonlinearSolver::Kind nonlinear_solver
std::string world_builder_file
AdvectionStabilizationMethod::Kind advection_stabilization_method
unsigned int stokes_velocity_degree
void declare_parameters(ParameterHandler &prm)
bool enable_additional_stokes_rhs
NullspaceRemoval::Kind nullspace_removal
static const std::string pattern()
bool use_artificial_viscosity_smoothing
double discontinuous_penalty
bool volume_of_fluid_tracking_enabled
double reaction_time_step
double stabilization_gamma
double linear_solver_S_block_tolerance
bool use_direct_stokes_solver
bool use_limiter_for_discontinuous_composition_solution
bool use_operator_splitting
double refinement_fraction
unsigned int n_compositional_fields
MaterialModel::MaterialAveraging::AveragingOperation material_averaging
unsigned int reaction_steps_per_advection_step
bool use_conduction_timestep
double maximum_relative_increase_time_step
double composition_solver_tolerance
std::vector< std::string > names_of_compositional_fields
unsigned int max_nonlinear_iterations
bool run_postprocessors_on_nonlinear_iterations
bool use_discontinuous_temperature_discretization
static bool is_defect_correction(const typename NonlinearSolver::Kind &input)
static std::string get_options_string()
std::set< types::boundary_id > fixed_heat_flux_boundary_indicators
std::string output_directory
std::vector< double > stabilization_beta
std::vector< aspect::CompositionalFieldDescription > composition_descriptions
std::vector< double > stabilization_c_R
std::vector< std::string > compositional_fields_with_disabled_boundary_entropy_viscosity
std::string pressure_normalization
bool use_equal_order_interpolation_for_stokes
StokesSolverType::Kind stokes_solver_type
unsigned int advection_gmres_restart_length
bool enable_prescribed_dilation
std::string AMG_smoother_type
StokesKrylovType::Kind stokes_krylov_type
double global_temperature_min_preset
Formulation::TemperatureEquation::Kind formulation_temperature_equation
bool include_melt_transport
double adiabatic_surface_temperature
std::vector< unsigned int > chemical_composition_indices
unsigned int max_nonlinear_iterations_in_prerefinement
std::vector< double > global_composition_max_preset
unsigned int timing_output_frequency
std::vector< typename AdvectionFieldMethod::Kind > compositional_field_methods
Formulation::Kind formulation
unsigned int n_chemical_compositions
bool skip_setup_initial_conditions_on_initial_refinement
unsigned int adaptive_refinement_interval
double diffusion_length_scale