MAMMA
1.0
Conduitsolver
|
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 |