Poroelasticity Calculations

NairnMPM can do poroelasticity calculations coupled with stresses and strains through tracking of pore pressure of saturated media. This feature is only available in OSParticulas.

Introduction

Poroelasticity is modeled by extension of method first described by . In brief, each particle tracks a pore pressure that models pressure within a liquid for a material saturated in the material. As the material volume changes, the pore pressure increases or decreases. The pore pressure can transport between particles based on Darcy tensor for permittivity of the material for the given liquid. Some new poroelasticity material properties control coupling between stresses and strains and the pore pressure.

The extension from Biot are to model poroelasticity in anisotropic materials and to allow the pore pressure to drop below zero.

Activating Poroelasticity

In scripted files, diffusion is activated with the command

Diffusion (YesOrNo),<(refConc)>

In XML input files, diffusion is activated with the <Diffusion> command, which must be within the <<MPMHeader> element:

<Diffusion reference = '(refConc)'/>

where

• (YesOrNo) must be "Yes" or "No" to activate or not activate diffusion calculations. In XML input files, the presence of a <Diffusion> command activates diffusion. The default is "No".
• (refConc) is used to set a reference concentration potential (between 0 and 1) that corresponds to zero strain. All diffusion calculations are done in terms of a concentration potential from 0 to 1 where 1 is the saturation concentration of a material type. The default (refConc) is 0.

Poroelasticity Material Properties

Concencentration changes are coupled to stress and strains through concentration expansion coefficients defined for the materials. By default, all moisture expansion coefficients are zero which decouples diffusion and strains. By entering non-zero values, the coupling will occur. Isotropic materials have a single solvent expansion coefficient while anisotropic will have two or three solvent expansion coefficients.

The rate of diffusion is controlled by the solvent diffusion constant defined for each material. Isotropic materials have a single solvent diffusion constant while anisotropic will have two or three solvent diffusion constants.

To be able to model diffusion in composite materials where different phases may absorb different amounts of solvent, all diffusion calculations are done in terms of a chemical potential for the solvent in the material, where chemical or concentration potential, μ, is approximate by

      [math]\displaystyle{ \mu= {c\over c_{sat}} }[/math]

where c is the weight fraction of solvent absorbed in the material and c_sat is the saturation solvent weight fraction for that material (which is specified in the material definition). This concentration potential varies from 0 to 1 and equilibrium conditions corresponds to all particles being at the same concentration potential.

Archived Concentrations

When calculated concentrations and concentration gradients are archived, they are converted to actual concentration in weight fraction using the material's saturation concentration setting:

      [math]\displaystyle{ c = c_{sat}\mu }[/math]

This conversion applies both to particle archives and to global archiving.

Poroelasticity Boundary Conditions

When diffusion is activated, you can set, the possible concentration boundary conditions are:

• You can set concentration on the grid.
• Rigid particles can provide moving concentration boundary conditions.
• You can set a concentration flux on particle surfaces.

Note that setting initial particle concentrations different than the reference concentration will cause strains to immediately evolve toward the changed state. The net effect will be an instantaneous "impact" of concentration change that might cause undesirable dynamic effects.