Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA full Eulerian finite difference approach for solving fluid-structure coupling problems
501
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid using a finite difference scheme. A volume-of-fluid formulation (Hirt and Nichols (1981, J. Comput. Phys., 39, 201)), which has been widely used for multiphase flow simulations, is applied to describing the multi-component geometry. The temporal change in the solid deformation is described in the Eulerian frame by updating a left Cauchy-Green deformation tensor, which is used to express constitutive equations for nonlinear Mooney-Rivlin materials. In this paper, various verifications and validations of the present full Eulerian method, which solves the fluid and solid motions on a fixed grid, are demonstrated, and the numerical accuracy involved in the fluid-structure coupling problems is examined.
rate research
Read More
Particle-in-Cell (PIC) methods are widely used computational tools for fluid and kinetic plasma modeling. While both the fluid and kinetic PIC approaches have been successfully used to target either kinetic or fluid simulations, little was done to combine fluid and kinetic particles under the same PIC framework. This work addresses this issue by proposing a new PIC method, PolyPIC, that uses polymorphic computational particles. In this numerical scheme, particles can be either kinetic or fluid, and fluid particles can become kinetic when necessary, e.g. particles undergoing a strong acceleration. We design and implement the PolyPIC method, and test it against the Landau damping of Langmuir and ion acoustic waves, two stream instability and sheath formation. We unify the fluid and kinetic PIC methods under one common framework comprising both fluid and kinetic particles, providing a tool for adaptive fluid-kinetic coupling in plasma simulations.
Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newtons equations of motion are often sufficient to solve for the motion of cylindrical particles with low density ratios, for more complex particles - such as a body with a protrusion - they become unstable. We present an implicit formulation of the coupling between rigid body dynamics and fluid dynamics within the framework of the immersed boundary projection method. Similarly to previous work on this method, the resulting matrix equation in the present approach is solved using a block-LU decomposition. Each step of the block-LU decomposition is modified to incorporate the rigid body dynamics. We show that our method achieves second-order accuracy in space and first-order in time (third-order for practical settings), only with a small additional computational cost to the original method. Our implicit coupling yields stable solution for density ratios as low as $10^{-4}$. We also consider the influence of fictitious fluid located inside the rigid bodies on the accuracy and stability of our method.
Soft solids in fluids find wide range of applications in science and engineering, especially in the study of biological tissues and membranes. In this study, an Eulerian finite volume approach has been developed to simulate fully resolved incompressible hyperelastic solids immersed in a fluid. We have adopted the recently developed reference map technique (RMT) by Valkov et. al (J. Appl. Mech., 82, 2015) and assessed multiple improvements for this approach.These modifications maintain the numerical robustness of the solver and allow the simulations without any artificial viscosity in the solid regions (to stabilize the solver). This has also resulted in eliminating the striations (wrinkles) of the fluid-solid interface that was seen before and hence obviates the need for any additional routines to achieve a smooth interface. An approximate projection method has been used to project the velocity field onto a divergence free field. Cost and accuracy improvements of the modifications on the method have also been discussed.
We present a new nonlinear mode decomposition method to visualize the decomposed flow fields, named the mode decomposing convolutional neural network autoencoder (MD-CNN-AE). The proposed method is applied to a flow around a circular cylinder at $Re_D=100$ as a test case. The flow attributes are mapped into two modes in the latent space and then these two modes are visualized in the physical space. Because the MD-CNN-AEs with nonlinear activation functions show lower reconstruction errors than the proper orthogonal decomposition (POD), the nonlinearity contained in the activation function is considered the key to improve the capability of the model. It is found by applying POD to each field decomposed using the MD-CNN-AE with hyperbolic tangent activation that a single nonlinear MD-CNN-AE mode contains multiple orthogonal bases, in contrast to the linear methods, i.e., POD and the MD-CNN-AE with linear activation. We further assess the proposed MD-CNN-AE by applying it to a transient process of a circular cylinder wake in order to examine its capability for flows containing high-order spatial modes. The present results suggest a great potential for the nonlinear MD-CNN-AE to be used for feature extraction of flow fields in lower dimension than POD, while retaining interpretable relationships with the conventional POD modes.
Numerical simulations of cardiovascular mass transport pose significant challenges due to the wide range of Peclet numbers and backflow at Neumann boundaries. In this paper we present and discuss several numerical tools to address these challenges in the context of a stabilized finite element computational framework. To overcome numerical instabilities when backflow occurs at Neumann boundaries, we propose an approach based on the prescription of the total flux. In addition, we introduce a consistent flux outflow boundary condition and demonstrate its superior performance over the traditional zero diffusive flux boundary condition. Lastly, we discuss discontinuity capturing (DC) stabilization techniques to address the well-known oscillatory behavior of the solution near the concentration front in advection-dominated flows.We present numerical examples in both idealized and patient-specific geometries to demonstrate the efficacy of the proposed procedures. The three contributions dis-cussed in this paper enable to successfully address commonly found challenges when simulating mass transport processes in cardiovascular flows.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
