MAMMA  1.0
Conduitsolver
Modules | Functions/Subroutines | Variables
constitutive.f90 File Reference

Go to the source code of this file.

Modules

module  constitutive
 Constitutive equations.
 

Functions/Subroutines

subroutine constitutive::initialize_models
 
subroutine constitutive::allocate_phases_parameters
 Initialization of relaxation flags. More...
 
subroutine constitutive::eos
 Equation of state. More...
 
subroutine constitutive::sound_speeds (C_mix, mach)
 Local sound speeds. More...
 
subroutine constitutive::eval_densities
 Phases densities. More...
 
subroutine constitutive::solve_cubic (a, b, c, y1, y2, y3)
 Solution of a cubic. More...
 
subroutine constitutive::f_xdis_eq
 Equilibrium Dissolved gas. More...
 
subroutine constitutive::f_beta_eq
 Equilibrium Crystal content. More...
 
subroutine constitutive::f_growth_rate
 This subroutine compute the growth rates for the different crystal phases. More...
 
subroutine constitutive::f_nucleation_rate
 Lithostatic pressure. More...
 
subroutine constitutive::lithostatic_pressure
 Lithostatic pressure. More...
 
subroutine constitutive::hydrostatic_pressure
 Hydrostatic pressure. More...
 
subroutine constitutive::press_relax_term (pressure_relaxation)
 Pressure relaxation term. More...
 
subroutine constitutive::vel_relax_term (velocity_relaxation)
 Velocity relaxation term. More...
 
subroutine constitutive::f_permkc
 Magma permeability. More...
 
subroutine constitutive::mixture_viscosity
 Mixture viscosity. More...
 
subroutine constitutive::f_viscliq
 Magma viscosity. More...
 
subroutine constitutive::f_viscmelt
 Melt viscosity. More...
 
subroutine constitutive::f_theta
 Crystal relative viscosity. More...
 
subroutine constitutive::f_bubbles
 Exsolved gas relative viscosity. More...
 
subroutine constitutive::f_alfa (xtot, xmax, r_beta, r_rho_md, r_rho_2, r_alfa_2)
 Bottom exsolved gas. More...
 
subroutine constitutive::f_alfa3 (r_p_2, xtot, r_beta, r_rho_md, r_rho_g, r_alfa_g)
 Bottom exsolved gas. More...
 

Variables

real *8 constitutive::cv_2
 exsolved gas specific heat capacity at constant volume More...
 
real *8 constitutive::cv_m
 melt specific heat capacity at constant volume More...
 
real *8, dimension(:), allocatable constitutive::cv_c
 crystals specific heat capacity at constant volume More...
 
real *8, dimension(:), allocatable constitutive::cv_d
 dissolved gas heat capacity at constant volume More...
 
real *8, dimension(:), allocatable constitutive::cv_g
 exsolved gas heat capacity at constant volume More...
 
real *8 constitutive::c0_m
 melt sound speed at atmospheric conditions More...
 
real *8, dimension(:), allocatable constitutive::c0_c
 crystals sound speed at atm conditions More...
 
real *8, dimension(:), allocatable constitutive::c0_d
 dissolved gas sound speed at atm conditions More...
 
real *8, dimension(:), allocatable constitutive::c0_g
 exsolved gas sound speed at atm conditions More...
 
real *8 constitutive::gamma_2
 exsolved gas adiabatic exponent More...
 
real *8 constitutive::gamma_m
 melt adiabatic exponent More...
 
real *8, dimension(:), allocatable constitutive::gamma_c
 crystals adiabatic exponent More...
 
real *8, dimension(:), allocatable constitutive::gamma_d
 dissolved gas adiabatic exponent More...
 
real *8, dimension(:), allocatable constitutive::gamma_g
 exsolved gas adiabatic exponent More...
 
real *8 constitutive::rho0_2
 exsolved gas reference density More...
 
real *8 constitutive::rho0_m
 melt reference density More...
 
real *8, dimension(:), allocatable constitutive::rho0_c
 crystals reference density More...
 
real *8, dimension(:), allocatable constitutive::rho0_d
 dissolved gas reference density More...
 
real *8, dimension(:), allocatable constitutive::rho0_g
 exsolved gas reference density More...
 
real *8 constitutive::t0_2
 exsolved gas reference temperature More...
 
real *8 constitutive::t0_m
 melt reference temperature More...
 
real *8, dimension(:), allocatable constitutive::t0_c
 crystals gas reference temperature More...
 
real *8, dimension(:), allocatable constitutive::t0_d
 dissolved gas reference temperature More...
 
real *8, dimension(:), allocatable constitutive::t0_g
 exsolved gas reference temperature More...
 
real *8 constitutive::p0_m
 melt reference pressure More...
 
real *8, dimension(:), allocatable constitutive::p0_c
 crystals reference pressure More...
 
real *8, dimension(:), allocatable constitutive::p0_d
 dissolved gas reference pressure More...
 
real *8, dimension(:), allocatable constitutive::p0_g
 exsolved gas reference pressure More...
 
real *8 constitutive::bar_e_2
 exsolved gas formation energy More...
 
real *8 constitutive::bar_e_m
 melt formation energy More...
 
real *8, dimension(:), allocatable constitutive::bar_e_c
 crystals formation energy More...
 
real *8, dimension(:), allocatable constitutive::bar_e_d
 dissolved gas formation energy More...
 
real *8, dimension(:), allocatable constitutive::bar_e_g
 exsolved gas formation energy More...
 
real *8 constitutive::bar_p_m
 melt cohesion pressure More...
 
real *8, dimension(:), allocatable constitutive::bar_p_c
 crystals cohesion pressure More...
 
real *8, dimension(:), allocatable constitutive::bar_p_d
 dissolved gas cohesion pressure More...
 
real *8, dimension(:), allocatable constitutive::bar_p_g
 exsolved gas cohesion pressure More...
 
real *8 constitutive::s0_m
 melt reference entropy More...
 
real *8, dimension(:), allocatable constitutive::s0_c
 crystals reference entropy More...
 
real *8, dimension(:), allocatable constitutive::s0_d
 dissolved gas reference entropy More...
 
real *8, dimension(:), allocatable constitutive::s0_g
 exsolved gas reference entropy More...
 
real *8, dimension(:), allocatable constitutive::pc_g
 critical gas pressure More...
 
real *8, dimension(:), allocatable constitutive::tc_g
 critical gas temperature More...
 
real *8, dimension(:), allocatable constitutive::a_g
 parameter for the VDW EOS More...
 
real *8, dimension(:), allocatable constitutive::b_g
 parameter for the VDW EOS More...
 
character *20 constitutive::gas_law
 equation of state for gas
More...
 
complex *16 constitutive::cv_1
 dis.gas+melt+crystals specific heat capacity at constant volume More...
 
complex *16 constitutive::e_mix
 total internal energy More...
 
complex *16 constitutive::e_1
 local specific internal energy of the melt-crystals phase More...
 
complex *16 constitutive::e_2
 local specific internal energy of the exsolved gas More...
 
complex *16 constitutive::e_m
 local specific internal energy of the melt More...
 
complex *16, dimension(:), allocatable constitutive::e_c
 local specific internal energy of the crystals More...
 
complex *16, dimension(:), allocatable constitutive::e_d
 local specific internal energy of the dissolved gas More...
 
complex *16, dimension(:), allocatable constitutive::e_g
 local specific internal energy of the exsolved gas More...
 
complex *16 constitutive::p_1
 local pressure of the melt-crystals phase More...
 
complex *16 constitutive::p_2
 local pressure of the exsolved gas More...
 
complex *16 constitutive::s_1
 local specific entropy of the melt-crystals phase More...
 
complex *16 constitutive::s_2
 local specific entropy of the exsolved gas More...
 
complex *16 constitutive::s_m
 local specific entropy of the melt More...
 
complex *16, dimension(:), allocatable constitutive::s_c
 local specific entropy of the crystals More...
 
complex *16, dimension(:), allocatable constitutive::s_d
 local specific entropy of the dissolved gas More...
 
complex *16, dimension(:), allocatable constitutive::s_g
 local specific entropy of the exsolved gas More...
 
complex *16 constitutive::mu_1
 free Gibbs energy of the melt-crystals phase More...
 
complex *16 constitutive::mu_2
 free Gibbs energy of the exsolved gas More...
 
complex *16 constitutive::mu_m
 free Gibbs energy of the melt More...
 
complex *16, dimension(:), allocatable constitutive::mu_c
 free Gibbs energy of the crystals More...
 
complex *16, dimension(:), allocatable constitutive::mu_d
 free Gibbs energy of the dissolved gas More...
 
complex *16, dimension(:), allocatable constitutive::mu_g
 free Gibbs energy of the exsolved gas More...
 
complex *16 constitutive::rho_1
 dis_gas+melt+crystals phase local density More...
 
complex *16 constitutive::rho_2
 exsolved gas local density More...
 
complex *16 constitutive::rho_m
 melt local density More...
 
complex *16, dimension(:), allocatable constitutive::rho_c
 crystals local density More...
 
complex *16, dimension(:), allocatable constitutive::rho_d
 dissolved gas local density More...
 
complex *16, dimension(:), allocatable constitutive::rho_g
 exsolved gas local density More...
 
complex *16 constitutive::rho_md
 dis_gas+melt local density More...
 
complex *16, dimension(:), allocatable constitutive::rhob_c
 crystals bulk density More...
 
complex *16 constitutive::rhob_m
 melt bulk density More...
 
complex *16 constitutive::alfa_1
 dis_gas+melt+crystals phase local volume fraction More...
 
complex *16 constitutive::alfa_2
 total exsolved gas local volume fraction More...
 
complex *16, dimension(:), allocatable constitutive::alfa_g
 exsolved gas phases local volume fraction More...
 
complex *16 constitutive::alfarho_2
 bulk density of the exsolved gas More...
 
complex *16 constitutive::alfa_m_1
 melt volume fraction in phase 1 More...
 
complex *16, dimension(:), allocatable constitutive::alfa_d_1
 dissolved gas volume fractions in phase 1 More...
 
complex *16, dimension(:), allocatable constitutive::alfa_g_2
 exsolved gas volume fractions in phase 2 More...
 
complex *16, dimension(:), allocatable constitutive::beta
 crystal volume fraction in the melt-crystals phase More...
 
complex *16, dimension(:), allocatable constitutive::beta_eq
 equil. cry. volume fraction in the melt-crystals phase More...
 
complex *16, dimension(:,:), allocatable constitutive::mom_cry
 moments of the crystal referred to the melt-crystals phase More...
 
complex *16, dimension(:), allocatable constitutive::growth_rate
 growth rate for the crystals More...
 
complex *16, dimension(:), allocatable constitutive::nucleation_rate
 nulceation rate for the crystals More...
 
complex *16 constitutive::u_1
 melt-crystals phase local velocity More...
 
complex *16 constitutive::u_2
 exsolved gas local velocity More...
 
complex *16 constitutive::x_1
 melt-crystals phase local mass fraction More...
 
complex *16 constitutive::x_2
 exsolved gas local mass fraction More...
 
complex *16, dimension(:), allocatable constitutive::x_c
 crystals mass fraction (with respect to the mixture) More...
 
complex *16 constitutive::x_m
 melt mass fraction (with respect to the mixture) More...
 
complex *16, dimension(:), allocatable constitutive::x_d
 dissolved gas mass fraction (with respect to the mixture) More...
 
complex *16, dimension(:), allocatable constitutive::x_g
 exsolved gas mass fraction (with respect to the mixture) More...
 
complex *16, dimension(:), allocatable constitutive::x_d_md
 dissolved gas mass fraction in the melt+dis.gas phase More...
 
complex *16, dimension(:), allocatable constitutive::x_d_md_eq
 equil. dis. gas mass fraction in the melt+dis.gas phase More...
 
complex *16, dimension(:), allocatable constitutive::x_c_1
 cristal mass fractions in phase 1 More...
 
complex *16 constitutive::x_m_1
 melt mass fraction in phase 1 More...
 
complex *16 constitutive::x_md_1
 melt+dis.gas mass fraction in phase 1 More...
 
complex *16, dimension(:), allocatable constitutive::x_d_1
 dissolved gas mass fractions in phase 1 More...
 
complex *16, dimension(:), allocatable constitutive::x_g_2
 exsolved gas mass fractions in phase 2 More...
 
complex *16 constitutive::t
 mixture local temperature More...
 
complex *16 constitutive::rho_mix
 mixture local density More...
 
complex *16 constitutive::u_mix
 mixture velocity More...
 
complex *16 constitutive::cv_mix
 mixture specific heat capacity at constant volume More...
 
complex *16 constitutive::visc_mix
 mixture viscosity More...
 
complex *16 constitutive::s
 mixture entropy More...
 
complex *16 constitutive::c_1
 first phase local sound speed More...
 
complex *16 constitutive::c_2
 second phase local sound speed More...
 
integer constitutive::n_drag_models
 
character(len=30), dimension(20) constitutive::available_drag_models
 
character *30 constitutive::drag_funct_model
 drag function model
More...
 
complex *16 constitutive::drag_funct
 drag function for the relative velocity:
More...
 
real *8 constitutive::drag_funct_coeff
 coefficient for the drag function for the relative velocity:
More...
 
real *8 constitutive::bubble_number_density
 bubble number density More...
 
real *8 constitutive::log10_bubble_number_density
 
real *8 constitutive::tortuosity_factor
 tortuosity factor More...
 
real *8 constitutive::throat_bubble_ratio
 throat bubble ratio More...
 
real *8 constitutive::friction_coefficient
 
real *8 constitutive::c_d
 
real *8 constitutive::r_a
 
character *20 constitutive::p_relax_model
 pressure relaxation model
More...
 
complex *16 constitutive::tau_p
 pressure relaxation rate More...
 
real *8 constitutive::tau_p_coeff
 pressure relaxation coefficient:
More...
 
real *8, dimension(:), allocatable constitutive::tau_c
 crystallization parameter:
More...
 
real *8, dimension(:), allocatable constitutive::beta0
 chamber (equilibrium) crystal volume fraction More...
 
real *8, dimension(:), allocatable constitutive::beta_max
 maximum crystal volume fraction More...
 
character *20 constitutive::crystallization_model
 Model for the equilibrium crystal volume fraction:
. More...
 
real *8 constitutive::frag_eff
 index of fragmentation in the interval [0;1] More...
 
integer constitutive::fragmentation_model
 Parameter to choose the fragmentation model:
. More...
 
real *8 constitutive::frag_thr
 Threshold for the fragmentation. More...
 
real *8 constitutive::tau_frag_exp
 Fragmentation exponent. More...
 
complex *16 constitutive::tau_frag
 fragmentation rate More...
 
real *8 constitutive::tau_frag_coeff
 fragmentation coefficient:
More...
 
real *8, dimension(:), allocatable constitutive::tau_d
 exsolution parameter:
More...
 
character *20 constitutive::exsol_model
 String for exsolution model:
. More...
 
real *8, dimension(:), allocatable constitutive::solub
 Solubility parameter for the Henry's law. More...
 
real *8, dimension(:), allocatable constitutive::solub_exp
 
real *8, dimension(:), allocatable constitutive::x_ex_dis_in
 
real *8 constitutive::grav
 gravitational acceleration More...
 
real *8 constitutive::k_cr
 country rock permeability More...
 
real *8 constitutive::rho_cr
 contry rock density More...
 
real *8 constitutive::visc_2
 gas viscosity More...
 
complex *16 constitutive::visc_melt
 melt viscosity More...
 
complex *16 constitutive::visc_1
 melt+crystal viscosity More...
 
complex *16 constitutive::visc_rel_bubbles
 relative viscosity due to bubbles More...
 
integer constitutive::n_theta_models
 
character(len=30), dimension(20) constitutive::available_theta_models
 
character *30 constitutive::theta_model
 Parameter to choose the model for the influence of crystal on the mixture: 'Lejeune_and_Richet1995' 'Dingwell1993' 'Melnik_and_Sparks1999' 'Costa2005' 'Melnik_and_Sparks2005' 'Vona_et_al2011' 'Vona_et_al2011_mod' 'Vona_et_al2013_eq19' 'Vona_et_al2013_eq20' 'Vona_et_al2013_eq21' 'Fixed_value' More...
 
complex *16 constitutive::theta
 Relative viscosity of the crystals. More...
 
real *8 constitutive::theta_fixed
 Fixed value for the relative viscosity of the crystals. More...
 
real *8 constitutive::c1
 Coefficients for the relative viscosity models. More...
 
real *8 constitutive::c2
 
real *8 constitutive::c3
 
real *8 constitutive::p_lith
 Lithostatic pressure. More...
 
real *8 constitutive::p_hydro
 Hydrostatic pressure. More...
 
real *8 constitutive::zeta_lith
 Elevation above the bottom for the evaluation of the lithostatic pressure. More...
 
integer constitutive::n_bubble_models
 
character(len=30), dimension(20) constitutive::available_bubble_models
 
character *20 constitutive::bubbles_model
 Parameter to choose the model for the influence of the bubbles on the mixture:
. More...
 
logical constitutive::explosive
 Flag to choose the eruptive style:
. More...
 
complex *16 constitutive::permkc
 Magma permeability. More...
 
complex *16 constitutive::inv_permkc
 
real *8 constitutive::alfa_switch
 
real *8 constitutive::a_2nd
 
real *8 constitutive::b_2nd
 
real *8 constitutive::perm0
 
integer constitutive::n_visc_melt_models
 
character(len=30), dimension(20) constitutive::available_visc_melt_models
 
character *30 constitutive::visc_melt_model
 Parameter to select the melt viscosity (bubbles and crystal-free) model:
. More...
 
real *8, dimension(12) constitutive::wt_init
 
real *8 constitutive::rho_w
 Water density. More...
 
real *8 constitutive::xa
 
real *8 constitutive::xb
 
real *8 constitutive::xc