An iterative matrix-free method in implicit immersed boundary/continuum methods

2007 ◽  
Vol 85 (11-14) ◽  
pp. 739-748 ◽  
Author(s):  
X. Sheldon Wang
Keyword(s):  
Immersed Boundary ◽  
Matrix Free ◽  
Boundary Continuum

scholarly journals A Revisit of Implicit Monolithic Algorithms for Compressible Solids Immersed Inside a Compressible Liquid

Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 273
Author(s):  
Sheldon Wang
Keyword(s):  
Finite Element ◽  
Finite Element Methods ◽  
Immersed Boundary Method ◽  
Compressible Liquid ◽  
Working Fluid ◽  
Immersed Boundary ◽  
Constant Density ◽  
Liquid Density ◽  
Boundary Method ◽  
Matrix Free

With the development of mature Computational Fluid Dynamics (CFD) tools for fluids (air and liquid) and Finite Element Methods (FEM) for solids and structures, many approaches have been proposed to tackle the so-called Fluid–Structure Interaction or Fluid–Solid Interaction (FSI) problems. Traditional partitioned iterations are often used to link available FEM codes with CFD codes in the study of FSI systems. Although these procedures are convenient, fluid mesh adjustments according to the motion and finite deformation of immersed solids or structures can be challenging or even prohibitive. Moreover, complex dynamic behaviors of coupled FSI systems are often lost in these iterative processes. In this paper, the author would like to review the so-called monolithic approaches for the solution of coupled FSI systems as a whole in the context of the immersed boundary method. In particular, the focus is on the implicit monolithic algorithm for compressible solids immersed inside a compressible liquid. Notice here the main focus of this paper is on liquid or more precisely liquid phase of water as working fluid. Using the word liquid, the author would like to emphasize the consideration of the compressibility of the fluid and the assumption of constant density and temperature. It is a common practice to assume that the pressure variations are not strong enough to alter the liquid density in any significant fashion for acoustic fluid–solid interactions problems. Although the algorithm presented in this paper is not directly applicable to aerodynamics in which the density change is significant along with its relationship with the pressure and the temperature, the author did revisit his earlier work on merging immersed boundary method concepts with a fully-fledged compressible aerodynamic code based on high-order compact scheme and energy conservative form of governing equations. In the proposed algorithm, on top of a uniform background (ghost) mesh, a fully implicit immersed method is implemented with mixed finite element methods for compressible liquid as well as immersed compressible solids with a matrix-free Newton–Krylov iterative solution scheme. In this monolithic approach, with the simple modulo function, the immersed solid or structure points can be easily located and thus the displacement projections and force distributions stipulated in the immersed boundary method can be effectively and efficiently implemented. This feature coupled with the key concept of the immersed boundary method helps to avoid topologically challenging mesh adjustments and to incorporate parallel processing commands such as Message Passing Interface (MPI) and further vectorization of the numerical operation. Once these high-performance procedures are implemented coupled with the monolithic implicit matrix-free Newton–Krylov iterative scheme with immersed methods, effective and efficient reduced order modeling techniques can then be employed to explore phase and parametric spaces. The in-house developed programs are at the moment two-dimensional. Furthermore, based on the same approach implemented in one-dimensional test example with one continuum immersed in another continuum, such monolithic implicit matrix-free Newton–Krylov iterative approach can be extended for the study of composites with deformable aggregates and matrix.

Microstructure simulation of early paper forming using immersed boundary methods

TAPPI Journal ◽  
2011 ◽  
Vol 11 (11) ◽  
pp. 23-30 ◽  
Author(s):  
ANDREAS MARK ◽  
ERIK SVENNING ◽  
ROBERT RUNDQVIST ◽  
FREDRIK EDELVIK ◽  
ERIK GLATT ◽  
...  
Keyword(s):  
Early Paper ◽  
Contact Model ◽  
Fiber Suspension ◽  
Beam Equation ◽  
Immersed Boundary ◽  
Paper Machine ◽  
Immersed Boundary Methods ◽  
Element Discretization ◽  
Microstructure Simulation ◽  
Boundary Methods

Paper forming is the first step in the paper machine where a fiber suspension leaves the headbox and flows through a forming fabric. Complex physical phenomena occur as the paper forms, during which fibers, fillers, fines, and chemicals added to the suspension interact. Understanding this process is important for the development of improved paper products because the configuration of the fibers during this step greatly influences the final paper quality. Because the effective paper properties depend on the microstructure of the fiber web, a continuum model is inadequate to explain the process and the properties of each fiber need to be accounted for in simulations. This study describes a new framework for microstructure simulation of early paper forming. The simulation framework includes a Navier-Stokes solver and immersed boundary methods to resolve the flow around the fibers. The fibers were modeled with a finite element discretization of the Euler-Bernoulli beam equation in a co-rotational formulation. The contact model is based on a penalty method and includes friction and elastic and inelastic collisions. We validated the fiber model and the contact model against demanding test cases from the literature, with excellent results. The fluid-structure interaction in the model was examined by simulating an elastic beam oscillating in a cross flow. We also simulated early paper formation to demonstrate the potential of the proposed framework.

Sign in / Sign up

Export Citation Format

Share Document