StokesContinuity <: CollectionContainer

Data structure containing parameter and array objects for Stokes-continuity solver.

Fields

• parameters::Parameters: Parameter groups for Stokes solver configuration
• arrays::Arrays: Array groups for velocity, pressure, and stress fields

Constructor

StokesContinuity(ynum::Int, xnum::Int)::StokesContinuity

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Parameters <: AbstractParameterCollection

Parameter collection for Stokes-continuity solver.

Fields

• itype_density::ParameterInt: Density model ID
• stype_density::ParameterStr: Density model name
• picard::Picard: Picard iteration parameters
• solution_norms::SolutionNormsStokes: Solution norm tracking
• residual_norms::ResidualNormsStokes: Residual norm tracking
• velocity_calc_option::VelocityCalcOptions: Velocity calculation options
• output::OutputStokes: Output parameters
• build::BuildSysStokes: Sparse matrix build parameters

Constructor

Parameters(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Picard <: AbstractParameterGroup

Parameter group for Picard iteration parameters.

Fields

• nglobal::ParameterInt: Maximum number of global plasticity iterations (a.k.a. Picard iterations)
• tolerance_picard::ParameterFloat: Convergence criterion for the global plasticity iterations (a.k.a. Picard iterations)
• iconverge::ParameterInt: Convergence flag
• iglobal::ParameterInt: Iteration counter for global stokes loop
• itype_global::ParameterInt: Type of global loop: 0 = marker-base; 1 = node-based
• stype_global::ParameterStr: Global loop option name
• obj_list::Vector{Union{ParameterFloat, ParameterInt}}: List of parameter objects

Nested Dot Access

• nglobal = model.stokes_continuity.parameters.picard.nglobal.value
• tolerance_picard = model.stokes_continuity.parameters.picard.tolerance_picard.value
• iconverge = model.stokes_continuity.parameters.picard.iconverge.value
• iglobal = model.stokes_continuity.parameters.picard.iglobal.value
• itype_global = model.stokes_continuity.parameters.picard.itype_global.value
• stype_global = model.stokes_continuity.parameters.picard.stype_global.value

Constructor

Picard()

Create a new Picard parameter group with default values.

Returns

• Picard: New Picard parameter group with initialized values

SolutionNormsStokes <: AbstractParameterGroup

Parameter group for Stokes solution norm tracking.

Fields

• dsoluv1_abs_inf::ParameterFloat: Infinity norm of solution change
• dsoluv1_rel_inf::ParameterFloat: Relative infinity norm of solution change
• dsoluv1_abs_L2::ParameterFloat: L2 norm of solution change
• dsoluv1_rel_L2::ParameterFloat: Relative L2 norm of solution change
• dvx1_abs_L2::ParameterFloat: L2 norm of vx change
• dvx1_rel_L2::ParameterFloat: Relative L2 norm of vx change
• dvy1_abs_L2::ParameterFloat: L2 norm of vy change
• dvy1_rel_L2::ParameterFloat: Relative L2 norm of vy change
• dpr1_abs_L2::ParameterFloat: L2 norm of pressure change
• dpr1_rel_L2::ParameterFloat: Relative L2 norm of pressure change
• dvxy_abs_inf::ParameterFloat: Infinity norm of velocity solution
• dvxy_rel_inf::ParameterFloat: Relative infinity norm of velocity solution
• dvxy_abs_L2::ParameterFloat: L2 norm of velocity solution
• dvxy_rel_L2::ParameterFloat: Relative L2 norm of velocity solution
• global_yield_error::ParameterFloat: Global yield stress error
• obj_list::Vector{ParameterFloat}: List of parameter objects

Nested Dot Access

• dsoluv1_abs_inf = model.stokes_continuity.parameters.solution_norms.dsoluv1_abs_inf.value
• dsoluv1_abs_L2 = model.stokes_continuity.parameters.solution_norms.dsoluv1_abs_L2.value
• dvx1_abs_L2 = model.stokes_continuity.parameters.solution_norms.dvx1_abs_L2.value

Constructor

SolutionNormsStokes()

Returns

• SolutionNormsStokes: New SolutionNormsStokes parameter group with initialized values

ResidualNormsStokes <: AbstractParameterGroup

Parameter group for Stokes residual norm tracking.

Fields

• resnl_L2_ini::ParameterFloat: Initial L2 norm of Stokes system
• resnl_L2::ParameterFloat: Current L2 norm of Stokes system
• resnl_rel_L2::ParameterFloat: Relative L2 norm of Stokes system
• resx_L2::ParameterFloat: L2 norm of vx residual
• resy_L2::ParameterFloat: L2 norm of vy residual
• resc_L2::ParameterFloat: L2 norm of pressure residual
• resx_L2_ini::ParameterFloat: Initial L2 norm of vx residual
• resy_L2_ini::ParameterFloat: Initial L2 norm of vy residual
• resc_L2_ini::ParameterFloat: Initial L2 norm of pressure residual
• resnlx_L2::ParameterFloat: Non-linear L2 norm of vx residual
• resnly_L2::ParameterFloat: Non-linear L2 norm of vy residual
• resnlc_L2::ParameterFloat: Non-linear L2 norm of pressure residual
• resnlx_rel_L2::ParameterFloat: Relative non-linear L2 norm of vx residual
• resnly_rel_L2::ParameterFloat: Relative non-linear L2 norm of vy residual
• resnlc_rel_L2::ParameterFloat: Relative non-linear L2 norm of pressure residual
• resnlx_L2_ini::ParameterFloat: Initial non-linear L2 norm of vx residual
• resnly_L2_ini::ParameterFloat: Initial non-linear L2 norm of vy residual
• resnlc_L2_ini::ParameterFloat: Initial non-linear L2 norm of pressure residual
• obj_list::Vector{ParameterFloat}: List of parameter objects

Nested Dot Access

• resnl_L2_ini = model.stokes_continuity.parameters.residual_norms.resnl_L2_ini.value
• resnl_L2 = model.stokes_continuity.parameters.residual_norms.resnl_L2.value
• resnl_rel_L2 = model.stokes_continuity.parameters.residual_norms.resnl_rel_L2.value
• resx_L2 = model.stokes_continuity.parameters.residual_norms.resx_L2.value
• resy_L2 = model.stokes_continuity.parameters.residual_norms.resy_L2.value
• resc_L2 = model.stokes_continuity.parameters.residual_norms.resc_L2.value

Constructor

ResidualNormsStokes()

Returns

• ResidualNormsStokes: New ResidualNormsStokes parameter group with initialized values

VelocityCalcOptions <: AbstractParameterGroup

Parameter group for velocity calculation options.

Fields

• itype_velocity::ParameterInt: Velocity calculation options: -1 all zeros; 0 solve equations; 1 solid body rotation
• stype_velocity::ParameterStr: Velocity calculation option name
• obj_list::Vector{ParameterInt}: List of parameter objects

Nested Dot Access

• itype_velocity = model.stokes_continuity.parameters.velocity_calc_option.itype_velocity.value
• stype_velocity = model.stokes_continuity.parameters.velocity_calc_option.stype_velocity.value

Constructor

VelocityCalcOptions()

Returns

• VelocityCalcOptions: New VelocityCalcOptions parameter group with initialized values

OutputStokes <: AbstractParameterGroup

Parameter group for Stokes solver output parameters.

Fields

• outtest::ParameterInt: 0 off; 1 output grid, vis., RHS, and bc as text files
• outtest2::ParameterInt: 1 output sparse matrix, rhs and solution vector as text files
• outtest3::ParameterInt: 1 output sparse matrix, rhs and solution vector as npy files
• obj_list::Vector{ParameterInt}: List of parameter objects

Nested Dot Access

• outtest = model.stokes_continuity.parameters.output.outtest.value
• outtest2 = model.stokes_continuity.parameters.output.outtest2.value
• outtest3 = model.stokes_continuity.parameters.output.outtest3.value

Constructor

OutputStokes()

Returns

• OutputStokes: New OutputStokes parameter group with initialized values

BuildSysStokes <: AbstractParameterGroup

Parameter group for sparse matrix build parameters for Stokes system.

Fields

• ibuild::ParameterInt: 0 define full system, 1 only define non-zero elements
• hshift_to_vxR::ParameterInt: Horizontal shift index for system building equal to (ynum-1)*3
• N::ParameterInt: Number of rows/columns in NxN matrix equal to (xnum-1)(ynum-1)3
• nonzero_max::ParameterInt: Maximum number of non-zero values allowed set equal to xnumynum31
• pscale::ParameterFloat: Coefficient for scaling pressure
• iuse_interface_stabilization::ParameterInt: Flag for interface stabilization: 0 off; 1 on
• obj_list::Vector{Union{ParameterFloat, ParameterInt}}: List of parameter objects

Nested Dot Access

• ibuild = model.stokes_continuity.parameters.build.ibuild.value
• hshift_to_vxR = model.stokes_continuity.parameters.build.hshift_to_vxR.value
• N = model.stokes_continuity.parameters.build.N.value
• nonzero_max = model.stokes_continuity.parameters.build.nonzero_max.value
• pscale = model.stokes_continuity.parameters.build.pscale.value
• iuse_interface_stabilization = model.stokes_continuity.parameters.build.iuse_interface_stabilization.value

Constructor

BuildSysStokes(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• BuildSysStokes: New BuildSysStokes parameter group with initialized values

Arrays <: AbstractArrayCollection

Data structure containing array groups for Stokes-continuity solver.

Fields

• rhs::RhsStokes: Right-hand side arrays for Stokes equations
• staggered_grid_velocity::StaggeredGridVelocity: Velocity on staggered grid
• basic_grid_velocity::BasicGridVelocity: Velocity on basic grid
• velocity_solution::VelocitySolution: Velocity solution arrays
• stokes_solution::StokesSolution: Combined Stokes solution arrays
• strain_rate_and_spin::StrainRateAndSpin: Strain rate and spin tensor
• stress::Stress: Stress tensor components
• stress_change::StressChange: Stress change arrays
• pressure::Pressure: Pressure field arrays
• density::Density: Density field arrays
• shear_modulus::ShearModulus: Shear modulus arrays
• viscosity::Viscosity: Viscosity field arrays
• plastic_def::PlasticDef: Plastic deformation arrays
• residuals::ResidualStokes: Residual arrays for Stokes equations

Constructor

Arrays(ynum::Int, xnum::Int)::Arrays

Create a new Arrays collection with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

RhsStokes <: AbstractArrayGroup

Array group containing right-hand side arrays for Stokes equations.

Fields

• RX1::ScalarArray2DState: (ynum + 1, xnum) : Right-hand side values on vx grid for x-Stokes equation.
• RY1::ScalarArray2DState: (ynum, xnum + 1) : Right-hand side values on vy grid for y-Stokes equation.
• RC1::ScalarArray2DState: (ynum - 1, xnum - 1) : Right-hand side values on pressure grid for continuity equation.
• RHS::RhsStokesArray1DState: Discretized right-hand side array for Stokes-continuity equation

Nested Dot Access

• RX1 = model.stokes_continuity.arrays.rhs.RX1.array
• RY1 = model.stokes_continuity.arrays.rhs.RY1.array
• RC1 = model.stokes_continuity.arrays.rhs.RC1.array
• RHS = model.stokes_continuity.arrays.rhs.RHS.array

Constructor

RhsStokes(ynum::Int, xnum::Int)

Create a new RhsStokes array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• RhsStokes: New RhsStokes array group with initialized arrays

StaggeredGridVelocity <: AbstractArrayGroup

Array group containing velocity arrays on staggered grid.

Fields

• vy1_old::ScalarArray2DState: (ynum, xnum + 1) : Old y-component of velocity on staggered vy grid in m/s.
• vy1::ScalarArray2DState: (ynum, xnum + 1) : Y-component of velocity on staggered vy grid in m/s.
• vx1_old::ScalarArray2DState: (ynum + 1, xnum) : Old x-component of velocity on staggered vx grid in m/s.
• vx1::ScalarArray2DState: (ynum + 1, xnum) : X-component of velocity on staggered vx grid in m/s.

Nested Dot Access

• vy1_old = model.stokes_continuity.arrays.staggered_grid_velocity.vy1_old.array
• vy1 = model.stokes_continuity.arrays.staggered_grid_velocity.vy1.array
• vx1_old = model.stokes_continuity.arrays.staggered_grid_velocity.vx1_old.array
• vx1 = model.stokes_continuity.arrays.staggered_grid_velocity.vx1.array

Constructor

StaggeredGridVelocity(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• StaggeredGridVelocity: New StaggeredGridVelocity array group with initialized arrays

BasicGridVelocity <: AbstractArrayGroup

Array group containing velocity arrays on basic grid.

Fields

• vxb::ScalarArray2DState: (ynum, xnum) : x-component of velocity interpolated to basic grid in m/s.
• vyb::ScalarArray2DState: (ynum, xnum) : y-component of velocity interpolated to basic grid in m/s.

Nested Dot Access

• vxb = model.stokes_continuity.arrays.basic_grid_velocity.vxb.array
• vyb = model.stokes_continuity.arrays.basic_grid_velocity.vyb.array

Constructor

BasicGridVelocity(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• BasicGridVelocity: New BasicGridVelocity array group with initialized arrays

VelocitySolution <: AbstractArrayGroup

Array group containing velocity solution arrays.

Fields

• vxy_old::SolutionArray1DState: Old solution array with vx and vy solutions in m/s.
• vxy::SolutionArray1DState: Solution array with vx and vy solutions in m/s.

Nested Dot Access

• vxy_old = model.stokes_continuity.arrays.velocity_solution.vxy_old.array
• vxy = model.stokes_continuity.arrays.velocity_solution.vxy.array

Constructor

VelocitySolution(ynum::Int, xnum::Int)

Create a new VelocitySolution array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• VelocitySolution: New VelocitySolution array group with initialized arrays

StokesSolution <: AbstractArrayGroup

Array group containing combined Stokes solution arrays.

Fields

• soluv1_old::SolutionArray1DState: Old Stokes-continuity solution array with velocity and pressure in m/s.
• soluv1::SolutionArray1DState: Stokes-continuity solution array with velocity and pressure in m/s or Pa.

Nested Dot Access

• soluv1_old = model.stokes_continuity.arrays.stokes_solution.soluv1_old.array
• soluv1 = model.stokes_continuity.arrays.stokes_solution.soluv1.array

Constructor

StokesSolution(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• StokesSolution: New StokesSolution array group with initialized arrays

StrainRateAndSpin <: AbstractArrayGroup

Array group containing strain rate and spin tensor arrays.

Fields

• exy::ScalarArray2DState: (ynum, xnum) : Shear strain rate in 1/s.
• exx::ScalarArray2DState: (ynum - 1, xnum - 1) : Normal strain rate in 1/s.
• eii::ScalarArray2DState: (ynum - 1, xnum - 1) : Second invariant of deviatoric strain rate in 1/s.
• esp::ScalarArray2DState: (ynum, xnum) : Spin in 1/s.
• eii_plastic_basic::ScalarArray2DState: (ynum, xnum) : Second invariant of plastic strain rate on basic grid in 1/s.
• eii_plastic_pressure::ScalarArray2DState: (ynum - 1, xnum - 1) : Second invariant of plastic strain rate on pressure grid in 1/s.

Nested Dot Access

• exy = model.stokes_continuity.arrays.strain_rate_and_spin.exy.array
• exx = model.stokes_continuity.arrays.strain_rate_and_spin.exx.array
• eii = model.stokes_continuity.arrays.strain_rate_and_spin.eii.array
• esp = model.stokes_continuity.arrays.strain_rate_and_spin.esp.array

Constructor

StrainRateAndSpin(ynum::Int, xnum::Int)

Create a new StrainRateAndSpin array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• StrainRateAndSpin: New StrainRateAndSpin array group with initialized arrays

Stress <: AbstractArrayGroup

Array group containing stress tensor component arrays.

Fields

• sxy0::ScalarArray2DState: (ynum, xnum) : Old deviatoric shear stress array on basic grid in Pa used for: (1) storing copy of visco-elasto-plastic shear stress transport array (sxy1) and (2) storing product of old visco-elasto-plastic shear stress transport array (sxy1) and viscoelastic factor used in RHS terms of visco-elasto-plastic Stokes equations.
• sxy1::ScalarArray2DState: (ynum, xnum) : Visco-elasto-plastic deviatoric shear stress in Pa interpolated to basic grid from advected markers.
• sxy2::ScalarArray2DState: (ynum, xnum) : Visco-elasto-plastic deviatoric shear stress in Pa on basic grid calculated using updated deviatoric strain rates from Stokes velocity solution.
• sxx0::ScalarArray2DState: (ynum - 1, xnum - 1) : Old deviatoric normal stress array on basic grid in Pa used for: (1) storing copy of visco-elasto-plastic normal stress transport array (sxx1) and (2) storing product of old visco-elasto-plastic normal stress transport array (sxx1) and viscoelastic factor used in RHS terms of visco-elasto-plastic Stokes equations.
• sxx1::ScalarArray2DState: (ynum - 1, xnum - 1) : Visco-elasto-plastic deviatoric normal stress interpolated to basic grid from advected markers in Pa.
• sxx2::ScalarArray2DState: (ynum - 1, xnum - 1) : Visco-elasto-plastic deviatoric normal stress on pressure grid calculated using updated deviatoric strain rates from Stokes velocity solution in Pa.

Nested Dot Access

• sxy0 = model.stokes_continuity.arrays.stress.sxy0.array
• sxy1 = model.stokes_continuity.arrays.stress.sxy1.array
• sxy2 = model.stokes_continuity.arrays.stress.sxy2.array
• sxx0 = model.stokes_continuity.arrays.stress.sxx0.array

Constructor

Stress(ynum::Int, xnum::Int)

Create a new Stress array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• Stress: New Stress array group with initialized arrays

StressChange <: AbstractArrayGroup

Array group containing stress change arrays.

Fields

• dsxy::ScalarArray2DState: (ynum, xnum) : Dual purpose array: First stores post-Stokes viscoelastic stress change (sxy2-sxy1) and then updated to final (remaining) shear stress change by subtracting relaxed grid nodal-marker stress differences (dsxyn) in Pa.
• dsxx::ScalarArray2DState: (ynum - 1, xnum - 1) : Dual purpose array: First stores post-Stokes normal viscoelastic stress change (sxx2-sxx1) and then updated to final (remaining) normal stress change by subtracting relaxed grid nodal-marker stress differences (dsxxn) in Pa.
• dsxyn::ScalarArray2DState: (ynum, xnum) : Subgrid nodal-marker shear stress change on basic grid interpolated from markers in Pa.
• dsxxn::ScalarArray2DState: (ynum - 1, xnum - 1) : Subgrid nodal-marker normal stress change on basic grid interpolated from markers in Pa.

Nested Dot Access

• dsxy = model.stokes_continuity.arrays.stress_change.dsxy.array
• dsxx = model.stokes_continuity.arrays.stress_change.dsxx.array
• dsxyn = model.stokes_continuity.arrays.stress_change.dsxyn.array
• dsxxn = model.stokes_continuity.arrays.stress_change.dsxxn.array

Constructor

StressChange(ynum::Int, xnum::Int)

Create a new StressChange array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• StressChange: New StressChange array group with initialized arrays

Pressure <: AbstractArrayGroup

Array group containing pressure field arrays.

Fields

• pr1_old::ScalarArray2DState: (ynum - 1, xnum - 1) : Old pressure solution on pressure grid in Pa.
• pr1::ScalarArray2DState: (ynum - 1, xnum - 1) : Pressure solution on pressure grid in Pa.

Nested Dot Access

• pr1_old = model.stokes_continuity.arrays.pressure.pr1_old.array
• pr1 = model.stokes_continuity.arrays.pressure.pr1.array

Constructor

Pressure(ynum::Int, xnum::Int)

Create a new Pressure array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• Pressure: New Pressure array group with initialized arrays

Density <: AbstractArrayGroup

Array group containing density field arrays.

Fields

• rho0::ScalarArray2DState: (ynum, xnum) : Backup of transport array with density on basic grid in kg/m^3.
• rho1::ScalarArray2DState: (ynum, xnum) : Transport array with density on basic grid in kg/m^3.
• rho0_vy::ScalarArray2DState: (ynum, xnum + 1) : Backup of transport array with density on vy grid in kg/m^3.
• rho1_vy::ScalarArray2DState: (ynum, xnum + 1) : Transport array with density on vy grid in kg/m^3.

Nested Dot Access

• rho0 = model.stokes_continuity.arrays.density.rho0.array
• rho1 = model.stokes_continuity.arrays.density.rho1.array
• rho0_vy = model.stokes_continuity.arrays.density.rho0_vy.array
• rho1_vy = model.stokes_continuity.arrays.density.rho1_vy.array

Constructor

Density(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• Density: New Density array group with initialized arrays

ShearModulus <: AbstractArrayGroup

Array group containing shear modulus arrays.

Fields

• mus0::ScalarArray2DState: (ynum, xnum) : Backup of shear modulus for shear stress on basic grid in Pa.
• mus1::ScalarArray2DState: (ynum, xnum) : Shear elastic modulus interpolated to basic grid from markers in Pa.
• mun0::ScalarArray2DState: (ynum - 1, xnum - 1) : Backup of shear modulus for normal stress on pressure grid in Pa.
• mun1::ScalarArray2DState: (ynum - 1, xnum - 1) : Normal elastic modulus interpolated to pressure grid from markers in Pa.

Nested Dot Access

• mus0 = model.stokes_continuity.arrays.shear_modulus.mus0.array
• mus1 = model.stokes_continuity.arrays.shear_modulus.mus1.array
• mun0 = model.stokes_continuity.arrays.shear_modulus.mun0.array
• mun1 = model.stokes_continuity.arrays.shear_modulus.mun1.array

Constructor

ShearModulus(ynum::Int, xnum::Int)

Create a new ShearModulus array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• ShearModulus: New ShearModulus array group with initialized arrays

Viscosity <: AbstractArrayGroup

Array group containing viscosity field arrays.

Fields

• etas0::ScalarArray2DState: (ynum, xnum) : Old viscoplastic shear viscosity array in Pa.s used for: (1) backup of viscoplastic shear viscosity (etas1) interpolated from markers used for problematic interpolation and (2) visco-elasto-plastic viscosity coefficients used in visco-elastic Stokes equations.
• etas1::ScalarArray2DState: (ynum, xnum) : Viscoplastic shear viscosity in Pa.s interpolated to basic grid from markers.
• etan0::ScalarArray2DState: (ynum - 1, xnum - 1) : Old viscoplastic normal viscosity array in Pa.s used for: (1) backup of viscoplastic normal viscosity (etan1) interpolated from markers used for problematic interpolation and (2) visco-elasto-plastic viscosity coefficients used in visco-elastic Stokes equations.
• etan1::ScalarArray2DState: (ynum - 1, xnum - 1) : Viscoplastic normal viscosity in Pa.s interpolated to pressure grid from markers.
• eta_flow::ScalarArray2DState: (ynum, xnum) : Flow viscosity in Pa.s interpolated to basic grid from markers.
• eta_flow0::ScalarArray2DState: (ynum, xnum) : Backup of flow viscosity on basic grid in Pa.s used for problematic interpolation.

Nested Dot Access

• etas0 = model.stokes_continuity.arrays.viscosity.etas0.array
• etas1 = model.stokes_continuity.arrays.viscosity.etas1.array
• etan0 = model.stokes_continuity.arrays.viscosity.etan0.array
• etan1 = model.stokes_continuity.arrays.viscosity.etan1.array

Constructor

Viscosity(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• Viscosity: New Viscosity array group with initialized arrays

PlasticDef <: AbstractArrayGroup

Array group containing plastic deformation arrays.

Fields

• plastics::ScalarArray2DState: (ynum, xnum) : Plastic deformation indicator on basic grid.
• plastics0::ScalarArray2DState: (ynum, xnum) : Backup of plastic deformation indicator on basic grid.
• plasticn::ScalarArray2DState: (ynum - 1, xnum - 1) : Plastic deformation indicator on pressure grid.
• plasticn0::ScalarArray2DState: (ynum - 1, xnum - 1) : Backup of plastic deformation indicator on pressure grid.
• plastic_yield::ScalarArray2DState: (ynum, xnum) : Plastic yielding marker on basic grid.
• yield_error::ScalarArray2DState: (ynum, xnum) : Plastic yielding stress error on basic grid.
• cohesion_grid::ScalarArray2DState: (ynum, xnum) : Cohesion on basic grid interpolated from markers in Pa.
• cohesion_grid0::ScalarArray2DState: (ynum, xnum) : Backup of cohesion on basic grid interpolated from markers in Pa.
• fric_degrees_grid::ScalarArray2DState: (ynum, xnum) : Friction angle in degrees on basic grid interpolated from markers in Degrees.
• fric_degrees_grid0::ScalarArray2DState: (ynum, xnum) : Backup of friction angle in degrees on basic grid in Degrees.
• dilatation_grid::ScalarArray2DState: (ynum - 1, xnum - 1) : Dilatation angle on pressure grid interpolated from markers in Degrees.
• dilatation_grid0::ScalarArray2DState: (ynum - 1, xnum - 1) : Backup of dilatation angle on pressure grid interpolated from markers.

Nested Dot Access

• plastics = model.stokes_continuity.arrays.plastic_def.plastics.array
• cohesion_grid = model.stokes_continuity.arrays.plastic_def.cohesion_grid.array

Constructor

PlasticDef(ynum::Int, xnum::Int)

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• PlasticDef: New PlasticDef array group with initialized arrays

ResidualStokes <: AbstractArrayGroup

Array group containing residual arrays for Stokes equations.

Fields

• resx1::ScalarArray2DState: (ynum + 1, xnum) : x-Stokes equation residual on staggered vx grid in m/s.
• resnlx::ScalarArray2DState: (ynum + 1, xnum) : Non-linear x-Stokes equation residual on staggered vx grid in m/s.
• resy1::ScalarArray2DState: (ynum, xnum + 1) : y-Stokes equation residual on staggered vy grid in m/s.
• resnly::ScalarArray2DState: (ynum, xnum + 1) : Non-linear y-Stokes equation residual on staggered vy grid in m/s.
• resc1::ScalarArray2DState: (ynum - 1, xnum - 1) : Continuity equation residual on pressure grid in Pa.
• resnlc::ScalarArray2DState: (ynum - 1, xnum - 1) : Non-linear continuity equation residual on pressure grid in Pa.
• resnl::SolutionArray1DState: Non-linear stokes-continuity system residual.

Nested Dot Access

• resx1 = model.stokes_continuity.arrays.residuals.resx1.array
• resy1 = model.stokes_continuity.arrays.residuals.resy1.array
• resc1 = model.stokes_continuity.arrays.residuals.resc1.array
• resnl = model.stokes_continuity.arrays.residuals.resnl.array

Constructor

ResidualStokes(ynum::Int, xnum::Int)

Create a new ResidualStokes array group with the given grid dimensions.

Arguments

• ynum::Int: Number of grid points in y-direction
• xnum::Int: Number of grid points in x-direction

Returns

• ResidualStokes: New ResidualStokes array group with initialized arrays