Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStabilized material point methods for coupled large deformation and fluid flow in porous materials
61
0
0.0
(
0
)
Added by
Jinhyun Choo
Publication date
2019
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The material point method (MPM) has been increasingly used for the simulation of large deformation processes in fluid-infiltrated porous materials. For undrained poromechanical problems, however, standard MPMs are numerically unstable because they use low-order interpolation functions that violate the inf-sup stability condition. In this work, we develop stabilized MPM formulations for dynamic and quasi-static poromechanics that permit the use of standard low-order interpolation functions notwithstanding the drainage condition. For the stabilization of both dynamic and quasi-static formulations, we utilize the polynomial pressure projection method whereby a stabilization term is augmented to the balance of mass. The stabilization term can be implemented with both the original and generalized interpolation material point (GIMP) methods, and it is compatible with existing time-integration methods. Here we use fully-implicit methods for both dynamic and quasi-static poromechanical problems, aided by a block-preconditioned Newton-Krylov solver. The stabilized MPMs are verified and investigated through several numerical examples under dynamic and quasi-static conditions. Results show that the proposed MPM formulations allow standard low-order interpolation functions to be used for both the solid displacement and pore pressure fields of poromechanical formulations, from undrained to drained conditions, and from dynamic to quasi-static conditions.
rate research
Read More
In this article, we present new random walk methods to solve flow and transport problems in unsaturated/saturated porous media, including coupled flow and transport processes in soils, heterogeneous systems modeled through random hydraulic conductivity and recharge fields, processes at the field and regional scales. The numerical schemes are based on global random walk algorithms (GRW) which approximate the solution by moving large numbers of computational particles on regular lattices according to specific random walk rules. To cope with the nonlinearity and the degeneracy of the Richards equation and of the coupled system, we implemented the GRW algorithms by employing linearization techniques similar to the $L$-scheme developed in finite element/volume approaches. The resulting GRW $L$-schemes converge with the number of iterations and provide numerical solutions that are first-order accurate in time and second-order in space. A remarkable property of the flow and transport GRW solutions is that they are practically free of numerical diffusion. The GRW solutions are validated by comparisons with mixed finite element and finite volume solutions in one- and two-dimensional benchmark problems. They include Richards equation fully coupled with the advection-diffusion-reaction equation and capture the transition from unsaturated to saturated flow regimes. For completeness, we also consider decoupled flow and transport model problems for saturated aquifers.
In this paper, we study a model for the transport of an external component, e.g., a surfactant, in variably saturated porous media. We discretize the model in time and space by combining a backward Euler method with the linear Galerkin finite elements. The Newton method and the L-Scheme are employed for the linearization and the performance of these schemes is studied numerically. A special focus is set on the effects of dynamic capillarity on the transport equation.
We study several iterative methods for fully coupled flow and reactive transport in porous media. The resulting mathematical model is a coupled, nonlinear evolution system. The flow model component builds on the Richards equation, modified to incorporate nonstandard effects like dynamic capillarity and hysteresis, and a reactive transport equation for the solute. The two model components are strongly coupled. On one hand, the flow affects the concentration of the solute; on the other hand, the surface tension is a function of the solute, which impacts the capillary pressure and, consequently, the flow. After applying an Euler implicit scheme, we consider a set of iterative linearization schemes to solve the resulting nonlinear equations, including both monolithic and two splitting strategies. The latter include a canonical nonlinear splitting and an alternate linearized splitting, which appears to be overall faster in terms of numbers of iterations, based on our numerical studies. The (time discrete) system being nonlinear, we investigate different linearization methods. We consider the linearly convergent L-scheme, which converges unconditionally, and the Newton method, converging quadratically but subject to restrictions on the initial guess. Whenever hysteresis effects are included, the Newton method fails to converge. The L-scheme converges; nevertheless, it may require many iterations. This aspect is improved by using the Anderson acceleration. A thorough comparison of the different solving strategies is presented in five numerical examples, implemented in MRST, a toolbox based on MATLAB.
We demonstrate the effective use of randomized methods for linear algebra to perform network-based analysis of complex vortical flows. Network theoretic approaches can reveal the connectivity structures among a set of vortical elements and analyze their collective dynamics. These approaches have recently been generalized to analyze high-dimensional turbulent flows, for which network computations can become prohibitively expensive. In this work, we propose efficient methods to approximate network quantities, such as the leading eigendecomposition of the adjacency matrix, using randomized methods. Specifically, we use the Nystrom method to approximate the leading eigenvalues and eigenvectors, achieving significant computational savings and reduced memory requirements. The effectiveness of the proposed technique is demonstrated on two high-dimensional flow fields: two-dimensional flow past an airfoil and two-dimensional turbulence. We find that quasi-uniform column sampling outperforms uniform column sampling, while both feature the same computational complexity.
This work proposes a new stabilized $P_1times P_0$ finite element method for solving the incompressible Navier--Stokes equations. The numerical scheme is based on a reduced Bernardi--Raugel element with statically condensed face bubbles and is pressure-robust in the small viscosity regime. For the Stokes problem, an error estimate uniform with respect to the kinematic viscosity is shown. For the Navier--Stokes equation, the nonlinear convection term is discretized using an edge-averaged finite element method. In comparison with classical schemes, the proposed method does not require tuning of parameters and is validated for competitiveness on several benchmark problems in 2 and 3 dimensional space.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
