Abstract

Abstract

We define stress and strain splittings appropriate to linearly elastic anisotropic materials with volumetric constraints. The treatment includes rigidtropic materials, which develop no strains under a stress pattern that is a null eigenvector of the compliance matrix. This model includes as special case incompressible materials, for which the eigenvector is hydrostatic stress. The main finding is that pressure and volumetric strain must be redefined as effective quantities. Using this idea, an energy decomposition that exactly separates deviatoric and volumetric energy follows.

Abstract

The Henry problem has played a key role in our understanding of seawater intrusion into coastal aquifers and in benchmarking density dependent flow codes. This paper seeks to modify Henry’s problem to ensure sensitivity to density variations and vertical salinity profiles that resemble field observations. In the proposed problem, the “dispersive Henry problem”, mixing is represented by means of the traditional Scheidegger dispersion tensor (dispersivity times water flux). Anisotropy in the hydraulic conductivity is acknowledged and Henry’s seaside boundary condition of prescribed salt concentration is replaced by a flux dependent boundary condition, which represents more realistically salt transport across the seaside boundary. This problem turns out to be very sensitive to density variations and its solution gets closer to reality. However, an improvement in the traditional Henry problem (gain in sensitivity and realism) can be also achieved if the value of the Peclet number is significantly reduced.

Although the dispersive problem lacks an analytical solution, it can shed light on flow in coastal aquifers. It provides significant information about the factors controlling seawater penetration, width of the mixing zone and influx of seawater. The width of the mixing zone depends basically on dispersion with longitudinal and transverse dispersion controlling different parts of the mixing zone but displaying similar overall effects. Toe penetration is mainly controlled by the horizontal permeability and by the geometric mean of the dispersivities. Finally, transverse dispersivity and the geometric mean of the hydraulic conductivity are the leading parameters controlling the amount of saltwater that enters the aquifer

Abstract

We present a local formulation for 2D Discrete Exterior Calculus (DEC) similar to that of the Finite Element Method (FEM), which allows a natural treatment of material heterogeneity (assigning material properties element by element). It also allows us to deduce, in a principled manner, anisotropic fluxes and the DEC discretization of the pullback of 1-forms by the anisotropy tensor, i.e. we deduce the the discrete action of the anisotropy tensor on primal 1-forms. Due to the local formulation, the computational cost of DEC is similar to that of the Finite Element Method with Linear interpolating functions (FEML). The numerical DEC solutions to the anisotropic Poisson equation show numerical convergence, are very close to those of FEML on fine meshes and are slightly better than those of FEML on coarse meshes.

Abstract

A geometrically non‐linear formulation for composites and the resulting explicit dynamic finite element algorithm are presented. The proposed formulation assumes that small elastic and large plastic strains, being the anisotropy considered using tensors which map the model variables at each time step into an equivalent isotropic space, where the integration of the rate constitutive equations is performed. The evolution of the internal variables is calculated in the auxiliary spaces, taking into account the material non‐linear behaviour, and the results mapped back to the real stress space. The updating of the mapping tensors for each new spatial configuration allows the treatment of general anisotropic materials under large strain and can be extended to treat multiphase composite materials using the mixing theory. The behaviour of the composite is dictated by the mechanical response of each substance, and the resultant model allows a fully non‐linear analysis combining different material models, such as damage in one compounding substance, elastoplastic behaviour in the other, while a third substance behaves elastically. 

Abstract

We study three “incompressibility flavors” of linearly-elastic anisotropic solids that exhibit volumetric constraints: isochoric, hydroisochoric and rigidtropic. An isochoric material deforms without volume change under any stress system. An hydroisochoric material does so under hydrostatic stress. A rigidtropic material undergoes zero deformations under a certain stress pattern. Whereas the three models coalesce for isotropic materials, important differences appear for anisotropic behavior. We find that isochoric and hydroisochoric models under certain conditions may be hampered by unstable physical behavior. Rigidtropic models can represent semistable physical materials of arbitrary anisotropy while including isochoric and hydroisochoric behavior as special cases.

Abstract

We study three “incompressibility flavors” of linearly-elastic anisotropic solids that exhibit volumetric constraints: isochoric, hydroisochoric and rigidtropic. An isochoric material deforms without volume change under any stress system. An hydroisochoric material does so under hydrostatic stress. A rigidtropic material undergoes zero deformations under a certain stress pattern. Whereas the three models coalesce for isotropic materials, important differences appear for anisotropic behavior. We find that isochoric and hydroisochoric models under certain conditions may be hampered by unstable physical behavior. Rigidtropic models can represent semistable physical materials of arbitrary anisotropy while including isochoric and hydroisochoric behavior as special cases.

Abstract

We present a local formulation for 2D Discrete Exterior Calculus (DEC) similar to that of the Finite Element Method (FEM), which allows a natural treatment of material heterogeneity (assigning material properties element by element). It also allows us to deduce, in a principled manner, anisotropic fluxes and the DEC discretization of the pullback of 1-forms by the anisotropy tensor, i.e. we deduce the the discrete action of the anisotropy tensor on primal 1-forms. Due to the local formulation, the computational cost of DEC is similar to that of the Finite Element Method with Linear interpolating functions (FEML). The numerical DEC solutions to the anisotropic Poisson equation show numerical convergence, are very close to those of FEML on fine meshes and are slightly better than those of FEML on coarse meshes.