ASPECT
Public Member Functions | Static Public Member Functions | Private Attributes | List of all members
aspect::MaterialModel::Rheology::DiffusionCreep< dim > Class Template Reference
Inheritance diagram for aspect::MaterialModel::Rheology::DiffusionCreep< dim >:
Inheritance graph
[legend]

Public Member Functions

 DiffusionCreep ()
 
void parse_parameters (ParameterHandler &prm, const std::unique_ptr< std::vector< unsigned int >> &expected_n_phases_per_composition=nullptr)
 
const DiffusionCreepParameters compute_creep_parameters (const unsigned int composition, const std::vector< double > &phase_function_values=std::vector< double >(), const std::vector< unsigned int > &n_phase_transitions_per_composition=std::vector< unsigned int >()) const
 
double compute_viscosity (const double pressure, const double temperature, const unsigned int composition, const std::vector< double > &phase_function_values=std::vector< double >(), const std::vector< unsigned int > &n_phase_transitions_per_composition=std::vector< unsigned int >()) const
 
double compute_viscosity (const double pressure, const double temperature, const double grain_size, const unsigned int composition, const std::vector< double > &phase_function_values=std::vector< double >(), const std::vector< unsigned int > &n_phase_transitions_per_composition=std::vector< unsigned int >()) const
 
std::pair< double, double > compute_strain_rate_and_derivative (const double stress, const double pressure, const double temperature, const DiffusionCreepParameters creep_parameters) const
 
std::pair< double, double > compute_strain_rate_and_derivative (const double stress, const double pressure, const double temperature, const double grain_size, const DiffusionCreepParameters creep_parameters) const
 
std::pair< double, double > compute_log_strain_rate_and_derivative (const double log_stress, const double pressure, const double temperature, const DiffusionCreepParameters creep_parameters) const
 
std::pair< double, double > compute_log_strain_rate_and_derivative (const double log_stress, const double pressure, const double temperature, const double grain_size, const DiffusionCreepParameters creep_parameters) const
 
- Public Member Functions inherited from aspect::SimulatorAccess< dim >
 SimulatorAccess ()
 
 SimulatorAccess (const Simulator< dim > &simulator_object)
 
virtual ~SimulatorAccess ()=default
 
virtual void initialize_simulator (const Simulator< dim > &simulator_object)
 
const Introspection< dim > & introspection () const
 
const Simulator< dim > & get_simulator () const
 
const Parameters< dim > & get_parameters () const
 
SimulatorSignals< dim > & get_signals () const
 
MPI_Comm get_mpi_communicator () const
 
TimerOutput & get_computing_timer () const
 
const ConditionalOStream & get_pcout () const
 
double get_time () const
 
double get_timestep () const
 
double get_old_timestep () const
 
unsigned int get_timestep_number () const
 
const TimeStepping::Manager< dim > & get_timestepping_manager () const
 
unsigned int get_nonlinear_iteration () const
 
const parallel::distributed::Triangulation< dim > & get_triangulation () const
 
double get_volume () const
 
const Mapping< dim > & get_mapping () const
 
std::string get_output_directory () const
 
bool include_adiabatic_heating () const
 
bool include_latent_heat () const
 
bool include_melt_transport () const
 
int get_stokes_velocity_degree () const
 
double get_adiabatic_surface_temperature () const
 
double get_surface_pressure () const
 
bool convert_output_to_years () const
 
unsigned int get_pre_refinement_step () const
 
unsigned int n_compositional_fields () const
 
double get_end_time () const
 
void get_refinement_criteria (Vector< float > &estimated_error_per_cell) const
 
void get_artificial_viscosity (Vector< float > &viscosity_per_cell, const bool skip_interior_cells=false) const
 
void get_artificial_viscosity_composition (Vector< float > &viscosity_per_cell, const unsigned int compositional_variable) const
 
const LinearAlgebra::BlockVector & get_current_linearization_point () const
 
const LinearAlgebra::BlockVector & get_solution () const
 
const LinearAlgebra::BlockVector & get_old_solution () const
 
const LinearAlgebra::BlockVector & get_old_old_solution () const
 
const LinearAlgebra::BlockVector & get_reaction_vector () const
 
const LinearAlgebra::BlockVector & get_mesh_velocity () const
 
const DoFHandler< dim > & get_dof_handler () const
 
const FiniteElement< dim > & get_fe () const
 
const LinearAlgebra::BlockSparseMatrix & get_system_matrix () const
 
const LinearAlgebra::BlockSparseMatrix & get_system_preconditioner_matrix () const
 
const MaterialModel::Interface< dim > & get_material_model () const
 
const GravityModel::Interface< dim > & get_gravity_model () const
 
const InitialTopographyModel::Interface< dim > & get_initial_topography_model () const
 
const GeometryModel::Interface< dim > & get_geometry_model () const
 
const AdiabaticConditions::Interface< dim > & get_adiabatic_conditions () const
 
bool has_boundary_temperature () const
 
const BoundaryTemperature::Manager< dim > & get_boundary_temperature_manager () const
 
const BoundaryHeatFlux::Interface< dim > & get_boundary_heat_flux () const
 
bool has_boundary_composition () const
 
const BoundaryComposition::Manager< dim > & get_boundary_composition_manager () const
 
const BoundaryTraction::Manager< dim > & get_boundary_traction_manager () const
 
std::shared_ptr< const InitialTemperature::Manager< dim > > get_initial_temperature_manager_pointer () const
 
const InitialTemperature::Manager< dim > & get_initial_temperature_manager () const
 
std::shared_ptr< const InitialComposition::Manager< dim > > get_initial_composition_manager_pointer () const
 
const InitialComposition::Manager< dim > & get_initial_composition_manager () const
 
const std::set< types::boundary_id > & get_fixed_temperature_boundary_indicators () const
 
const std::set< types::boundary_id > & get_fixed_heat_flux_boundary_indicators () const
 
const std::set< types::boundary_id > & get_fixed_composition_boundary_indicators () const
 
const std::set< types::boundary_id > & get_mesh_deformation_boundary_indicators () const
 
const BoundaryVelocity::Manager< dim > & get_boundary_velocity_manager () const
 
const HeatingModel::Manager< dim > & get_heating_model_manager () const
 
const MeshRefinement::Manager< dim > & get_mesh_refinement_manager () const
 
const MeltHandler< dim > & get_melt_handler () const
 
const VolumeOfFluidHandler< dim > & get_volume_of_fluid_handler () const
 
const NewtonHandler< dim > & get_newton_handler () const
 
const MeshDeformation::MeshDeformationHandler< dim > & get_mesh_deformation_handler () const
 
const LateralAveraging< dim > & get_lateral_averaging () const
 
const AffineConstraints< double > & get_current_constraints () const
 
bool simulator_is_past_initialization () const
 
double get_pressure_scaling () const
 
bool pressure_rhs_needs_compatibility_modification () const
 
bool model_has_prescribed_stokes_solution () const
 
TableHandler & get_statistics_object () const
 
const Postprocess::Manager< dim > & get_postprocess_manager () const
 
unsigned int n_particle_managers () const
 
const Particle::Manager< dim > & get_particle_manager (const unsigned int particle_manager_index) const
 
Particle::Manager< dim > & get_particle_manager (const unsigned int particle_manager_index)
 
bool is_stokes_matrix_free ()
 
const StokesMatrixFreeHandler< dim > & get_stokes_matrix_free () const
 
RotationProperties< dim > compute_net_angular_momentum (const bool use_constant_density, const LinearAlgebra::BlockVector &solution, const bool limit_to_top_faces=false) const
 

Static Public Member Functions

static void declare_parameters (ParameterHandler &prm)
 
- Static Public Member Functions inherited from aspect::SimulatorAccess< dim >
static void get_composition_values_at_q_point (const std::vector< std::vector< double >> &composition_values, const unsigned int q, std::vector< double > &composition_values_at_q_point)
 

Private Attributes

std::vector< double > prefactors_diffusion
 
std::vector< double > stress_exponents_diffusion
 
std::vector< double > grain_size_exponents_diffusion
 
std::vector< double > activation_energies_diffusion
 
std::vector< double > activation_volumes_diffusion
 
double fixed_grain_size
 

Detailed Description

template<int dim>
class aspect::MaterialModel::Rheology::DiffusionCreep< dim >

Definition at line 59 of file diffusion_creep.h.

Constructor & Destructor Documentation

§ DiffusionCreep()

template<int dim>
aspect::MaterialModel::Rheology::DiffusionCreep< dim >::DiffusionCreep ( )

Constructor.

Member Function Documentation

§ declare_parameters()

template<int dim>
static void aspect::MaterialModel::Rheology::DiffusionCreep< dim >::declare_parameters ( ParameterHandler &  prm)
static

Declare the parameters this function takes through input files.

§ parse_parameters()

template<int dim>
void aspect::MaterialModel::Rheology::DiffusionCreep< dim >::parse_parameters ( ParameterHandler &  prm,
const std::unique_ptr< std::vector< unsigned int >> &  expected_n_phases_per_composition = nullptr 
)

Read the parameters this class declares from the parameter file. If expected_n_phases_per_composition points to a vector of unsigned integers this is considered the number of phases for each compositional field and will be checked against the parsed parameters.

§ compute_creep_parameters()

template<int dim>
const DiffusionCreepParameters aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_creep_parameters ( const unsigned int  composition,
const std::vector< double > &  phase_function_values = std::vector< double >(),
const std::vector< unsigned int > &  n_phase_transitions_per_composition = std::vector< unsigned int >() 
) const

Compute the creep parameters for the diffusion creep law. If n_phase_transitions_per_composition points to a vector of unsigned integers this is considered the number of phase transitions for each compositional field and viscosity will be first computed on each phase and then averaged for each compositional field.

§ compute_viscosity() [1/2]

template<int dim>
double aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_viscosity ( const double  pressure,
const double  temperature,
const unsigned int  composition,
const std::vector< double > &  phase_function_values = std::vector< double >(),
const std::vector< unsigned int > &  n_phase_transitions_per_composition = std::vector< unsigned int >() 
) const

Compute the viscosity based on the diffusion creep law with the fixed grain size given in the input file. If n_phase_transitions_per_composition points to a vector of unsigned integers this is considered the number of phase transitions for each compositional field and viscosity will be first computed on each phase and then averaged for each compositional field.

§ compute_viscosity() [2/2]

template<int dim>
double aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_viscosity ( const double  pressure,
const double  temperature,
const double  grain_size,
const unsigned int  composition,
const std::vector< double > &  phase_function_values = std::vector< double >(),
const std::vector< unsigned int > &  n_phase_transitions_per_composition = std::vector< unsigned int >() 
) const

Compute the viscosity based on the diffusion creep law for the given grain_size. If n_phase_transitions_per_composition points to a vector of unsigned integers this is considered the number of phase transitions for each compositional field and viscosity will be first computed on each phase and then averaged for each compositional field.

§ compute_strain_rate_and_derivative() [1/2]

template<int dim>
std::pair<double, double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_strain_rate_and_derivative ( const double  stress,
const double  pressure,
const double  temperature,
const DiffusionCreepParameters  creep_parameters 
) const

Compute the strain rate and first stress derivative as a function of stress based on the diffusion creep law with the fixed grain size given in the input file.

§ compute_strain_rate_and_derivative() [2/2]

template<int dim>
std::pair<double, double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_strain_rate_and_derivative ( const double  stress,
const double  pressure,
const double  temperature,
const double  grain_size,
const DiffusionCreepParameters  creep_parameters 
) const

Compute the strain rate and first stress derivative as a function of stress based on the diffusion creep law for the given grain_size.

§ compute_log_strain_rate_and_derivative() [1/2]

template<int dim>
std::pair<double, double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_log_strain_rate_and_derivative ( const double  log_stress,
const double  pressure,
const double  temperature,
const DiffusionCreepParameters  creep_parameters 
) const

Compute the logarithm of strain rate and first derivative with respect to the logarithm of the stress based on the diffusion creep law with the fixed grain size given in the input file.

§ compute_log_strain_rate_and_derivative() [2/2]

template<int dim>
std::pair<double, double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::compute_log_strain_rate_and_derivative ( const double  log_stress,
const double  pressure,
const double  temperature,
const double  grain_size,
const DiffusionCreepParameters  creep_parameters 
) const

Compute the logarithm of strain rate and first derivative with respect to the logarithm of the stress based on the diffusion creep law for the given grain_size.

Member Data Documentation

§ prefactors_diffusion

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::prefactors_diffusion
private

List of diffusion creep prefactors A.

Definition at line 178 of file diffusion_creep.h.

§ stress_exponents_diffusion

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::stress_exponents_diffusion
private

List of diffusion creep stress exponents n (usually = 1).

Definition at line 183 of file diffusion_creep.h.

§ grain_size_exponents_diffusion

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::grain_size_exponents_diffusion
private

List of diffusion creep grain size exponents m.

Definition at line 188 of file diffusion_creep.h.

§ activation_energies_diffusion

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::activation_energies_diffusion
private

List of diffusion creep activation energies E.

Definition at line 193 of file diffusion_creep.h.

§ activation_volumes_diffusion

template<int dim>
std::vector<double> aspect::MaterialModel::Rheology::DiffusionCreep< dim >::activation_volumes_diffusion
private

List of diffusion creep activation volumes V.

Definition at line 198 of file diffusion_creep.h.

§ fixed_grain_size

template<int dim>
double aspect::MaterialModel::Rheology::DiffusionCreep< dim >::fixed_grain_size
private

Diffusion creep grain size d. This is read from the input file, and is only used by the functions that do not take the grain size as additional argument.

Definition at line 205 of file diffusion_creep.h.

The documentation for this class was generated from the following file: