A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems

This is the peer reviewed version of the following article: [de-Pouplana, I., and Oñate, E. (2017) A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems. Int. J. Numer. Anal. Meth. Geomech., 41: 110–134. doi: 10.1002/nag.2550], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nag.2550/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving."A new mixed displacement-pressure element for solving solid–pore fluid interaction problems is presented. In the resulting coupled system of equations, the balance of momentum equation remains unaltered, while the mass balance equation for the pore fluid is stabilized with the inclusion of higher-order terms multiplied by arbitrary dimensions in space, following the finite calculus (FIC) procedure. The stabilized FIC-FEM formulation can be applied to any kind of interpolation for the displacements and the pressure, but in this work, we have used linear elements of equal order interpolation for both set of unknowns. Examples in 2D and 3D are presented to illustrate the accuracy of the stabilized formulation for solid–pore fluid interaction problems.Peer ReviewedPostprint (author's final draft

Finite element modelling of fracture propagation in saturated media using quasi-zero-thickness interface elements

A new computational technique for the simulation of 2D and 3D fracture propagation processes in saturated porous media is presented. A non-local damage model is conveniently used in conjunction with interface elements to predict the degradation pattern of the domain and insert new fractures followed by remeshing. FIC-stabilized elements of equal order interpolation in the displacement and the pore pressure have been successfully used under complex conditions near the undrained-incompressible limit. A bilinear cohesive fracture model describes the mechanical behaviour of the joints. A formulation derived from the cubic law models the fluid flow through the crack. Examples in 2D and 3D, using 3-noded triangles and 4-noded tetrahedra respectively, are presented to illustrate the accuracy and robustness of the proposed methodology.Peer ReviewedPostprint (author's final draft

Low-order stabilized finite element for the full Biot formulation in soil mechanics at finite strain

This article presents a novel finite element formulation for the Biot equation using low-order elements. Additionally, an extra degree of freedom is introduced to treat the volumetric locking steaming from the effective response of the medium; its balance equation is also stabilized. The accuracy of the proposed formulation is demonstrated by means of numerical analyses.Peer ReviewedPostprint (author's final draft

Desarrollo de nuevos métodos para el análisis fluido-estructura mediante PFEM

El objetivo de la tesina será avanzar en la implementación de métodos numéricos para la modelización de la interacción entre fluido-estructura y fluido-partícula utilizando PFEM, elementos sólidos en tres dimensiones y partículas discretas (DEM). Con este tipo de metodología se pretende estudiar los posibles efectos producidos por desastres naturales como avenidas, inundaciones, desprendimientos, etc. en las que un fluido puede arrastrar objetos de distinto tamañ

PFEM–DEM for particle-laden flows with free surface

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-019-00244-1This work proposes a fully Lagrangian formulation for the numerical modeling of free-surface particle-laden flows. The fluid phase is solved using the particle finite element method (PFEM), while the solid particles embedded in the fluid are modeled with the discrete element method (DEM). The coupling between the implicit PFEM and the explicit DEM is performed through a sub-stepping staggered scheme. This work only considers suspended spherical particles that are assumed not to affect the fluid motion. Several tests are presented to validate the formulation. The PFEM–DEM results show very good agreement with analytical solutions, laboratory tests and numerical results from alternative numerical methods.Peer ReviewedPostprint (author's final draft

A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems

This is the peer reviewed version of the following article: [de-Pouplana, I., and Oñate, E. (2017) A FIC-based stabilized mixed finite element method with equal order interpolation for solid–pore fluid interaction problems. Int. J. Numer. Anal. Meth. Geomech., 41: 110–134. doi: 10.1002/nag.2550], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nag.2550/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving."A new mixed displacement-pressure element for solving solid–pore fluid interaction problems is presented. In the resulting coupled system of equations, the balance of momentum equation remains unaltered, while the mass balance equation for the pore fluid is stabilized with the inclusion of higher-order terms multiplied by arbitrary dimensions in space, following the finite calculus (FIC) procedure. The stabilized FIC-FEM formulation can be applied to any kind of interpolation for the displacements and the pressure, but in this work, we have used linear elements of equal order interpolation for both set of unknowns. Examples in 2D and 3D are presented to illustrate the accuracy of the stabilized formulation for solid–pore fluid interaction problems.Peer Reviewe

Enhanced continuum damage modeling of mechanical failure in ice and rocks

Modeling fracture in geomaterials is essential to the understanding of many physical phenomenon which may posses natural hazards e.g. landslides, faults and iceberg calving or man-made processes e.g. hydraulic fracture and excavations. Continuum Damage Mechanics (CDM) models the crack as a solid region with a degraded stiffness. This continuum definition of cracks in CDM allows more feasible coupling with other forms of material non-linearity and eliminates the need to track complicated crack geometry. Using CDM to analyze fracture for the modeling of fracture in geomaterials encounters several challenges e.g.: 1) the need to model the multiple physical processes occurring in geomaterials, typically: coupled fluid flow and solid deformation, 2) the need to consider non-local damage and transport in order to capture the underlying long range interactions and achieve mesh-independent finite element solutions and 3) the elevated computational cost associated with non-linear mixed finite element formulations. The research presented in this thesis aims at improving the CDM formulations for modeling fracture geomaterials. This research can be divided into three main parts. The first is the introduction of a novel non-local damage transport formulation for modeling fracture in poroelastic media. The mathematical basis of the formulation are derived from thermodynamic equilibrium that considers non-local processes and homogenization principles. The non-local damage transport model leads to two additional regularization equations, one for non-local damage and the other for non-local transport which is reduced to non-local permeability. We consider two options for the implementation of the derived non-local transport damage model. The first option is the four-field formulation which extends the (u/P) formulation widely used in poroelasticity to include the non-local damage and transport phenomena. The second option is the three-field formulation, which is based on the coupling of the regularization equations under the assumptions of similar damage and permeability length scales and similar driving local stress/strain for the evolution of the damage and permeability. The three-field formulation is computationally cheaper but it degrades the physical modeling capabilities of the model. For each of these formulations, a non-linear mixed-finite element solution is developed and the Jacobian matrix is derived analytically. The developed formulations are used in the analysis of hydraulic fracture and consolidation examples. In the second part, a novel approach for CDM modeling of hydraulic fracture of glaciers is pretended. The presence of water-filled crevasses is known to increase the penetration depth of crevasses and this has been hypothesized to play an important role controlling iceberg calving rate. Here, we develop a continuum damage-based poro-mechanics formulation that enables the simulation of water-filled basal and/or surface crevasse propagation. The formulation incorporates a scalar isotropic damage variable into a Maxwell-type viscoelastic constitutive model for glacial ice and the effect of the water pressure on fracture propagation using the concept of effective solid stress. We illustrate the model by simulating quasi-static hydro-fracture in idealized rectangular slabs of ice in contact with the ocean. Our results indicate that water-filled basal crevasses only propagate when the water pressure is sufficiently large and that the interaction between simultaneously propagating water-filled surface and basal crevasses can have a mutually positive influence leading to deeper crevasse propagation which can critically affect glacial stability. In the third part, we propose a coupled Boundary Element Method (BEM) and Finite Element Method (FEM) for modeling localized damage growth in structures. BEM offers the flexibility of modeling large domains efficiently while the nonlinear damage growth is accurately accounted by a local FEM mesh. An integral-type nonlocal continuum damage mechanics with adapting FEM mesh is used to model multiple damage zones and follow their propagation in the structure. Strong form coupling, BEM hosted, is achieved using Lagrange multipliers. Since the non-linearity is isolated in the FEM part of the system of equations, the system size is reduced using Schur complement approach, then, the solution is obtained by a monolithic Newton method that is used to solve both domains simultaneously. The method is applied to multiple fractures growth benchmark problems and shows good agreement with the literature

Research and Technology, 1998

This report selectively summarizes the NASA Lewis Research Center's research and technology accomplishments for the fiscal year 1998. It comprises 134 short articles submitted by the staff scientists and engineers. The report is organized into five major sections: Aeronautics, Research and Technology, Space, Engineering and Technical Services, and Commercial Technology. A table of contents and an author index have been developed to assist readers in finding articles of special interest. This report is not intended to he a comprehensive summary of all the research and technology work done over the past fiscal year. Most of the work is reported in Lewis-published technical reports, journal articles, and presentations prepared by Lewis staff and contractors. In addition, university grants have enabled faculty members and graduate students to engage in sponsored research that is reported at technical meetings or in journal articles. For each article in this report, a Lewis contact person has been identified, and where possible, reference documents are listed so that additional information can be easily obtained. The diversity of topics attests to the breadth of research and technology being pursued and to the skill mix of the staff that makes it possible. At the time of publication, NASA Lewis was undergoing a name change to the NASA John H. Glenn Research Center at Lewis Field