Search

The technique of identifying Lagrangian Coherent Structures (LCSs) has emerged as a powerful tool for studying incompressible flows. Yet, the discrepancies arising from incompressibility assumptions often compromise the accuracy of LCSs constructed using Lagrangian particle methods. In this study, we introduce a modified framework to compute Finite-Time Lyapunov Exponents (FTLEs), addressing the misalignment between the fully incompressible assumption of LCS theory and the inherent incompressibility loss in simulations of particle methods. We begin by examining the correlation between the minimum and maximum FTLEs. By incorporating the deformation gradient and Cauchy-Green strain tensor which account for the time-advancing errors of incompressibility based on continuum theory, we enhance the computational accuracy of FTLEs. Moreover, we introduce the modified FTLE algorithm to the incompressible particle-in-cell (PIC) method for resolving free surface flows and fluid-structure interaction problems. Finally, numerical examples including the Tayler-Green vortices, water sloshing with baffles, an eccentric box sinking in water, and three-dimensional shear-driven cavity problems with high Reynolds numbers are tested to validate the effectiveness of the modified FTLE algorithm and the improved LCSs. These results demonstrate that the proposed modification scheme adeptly counteracts the errors caused by incompressibility loss, enabling accurate computation of FTLEs and detection of LCSs. [ABSTRACT FROM AUTHOR]

In this paper, we establish a fluid-structure interaction (FSI) model to investigate the dynamic interactions between the underwater explosion (UNDEX) shock wave and a movable sphere near the free surface. We utilize the local discontinuous Galerkin (LDG) method to capture the propagation of the shock wave in the fluid domain and employ the pressure cutoff model to calculate cavitation effects. The fluid elements at the fluid-structure interface are directly coupled to the structural dynamic model, and the structural transient dynamic responses are coupled with fluid acoustic pressure at the fluid-structure interface in the governing equation. The validity of the present FSI model is verified by comparing with the continuous Galerkin method. Due to the advantage of the LDG method in capturing the discontinuous wave, the present model shows better properties than the traditional coupled acoustic-structural model. With the present FSI model, we investigate the interaction between the UNDEX shock wave and a submerged and floated sphere. Under the combined effects of the free surface and structure, the UNDEX shock and cavitation loading characteristics are analyzed, and the influences of complicated cavitation effects on dynamic responses of the sphere are discussed. [ABSTRACT FROM AUTHOR]

In recent years, the investigations on Lagrangian Coherent Structures (LCS) in complex flows have attracted increasing attention due to their ability to accurately describe flow field details. In this study, we propose a novel and accurate numerical technique based on the Lagrangian–Eulerian Stabilized Collocation Method (LESCM) for computing the Finite Time Lyapunov Exponents (FTLEs), which is essential for extracting LCSs in viscous incompressible flows. LESCM, previously employed for simulating fluid–structure interaction problems, was developed based on the Material Point Method (MPM). The hybrid Lagrangian–Eulerian​ description in LESCM enables explicit tracking of fluid flow trajectories throughout the simulation and direct evaluation of the FTLE field. Furthermore. The errors in FTLEs caused by the Particle Shifting Technique (PST) are completely avoided due to the Eulerian characteristic of LESCM as the deformation gradient is calculated based on fixed Eulerian nodes rather than fluid particles with unphysical shifting. Consequently, the novel technique based on LESCM surpasses the accuracy of pure Lagrangian particle methods and provides an accurate way of detecting complex LCSs in flow fields. Additionally, By harnessing the remarkable efficiency of LESCM, MATLAB can now handle up to 16 million particles with ease, eliminating the need for parallel computation techniques. Several numerical examples, such as 2D Taylor-Green​ vortices, flow passing a circular cylinder and a 3-dimensional shear driven cavity problem, have been tested to demonstrate the effectiveness of this novel approach under various conditions. [ABSTRACT FROM AUTHOR]

• A numerical model is developed to analyze underwater explosion of slender explosive and structure responses. • Slender explosive will generate a directional pressure field in the surrounding water. • The shape of explosive strongly influences shock wave dynamics in the near-field. • The direction effect of pressure and energy flux is significant. • The response of a plate subjected to near-field underwater explosion of a slender explosive is studied. The shape of a high-explosive is an important factor to consider in the near-field and contact explosion, but this has been rarely studied previously. In this paper, we first conduct experiments to study the directional effects of the underwater explosion shock wave of slender explosives. A numerical model based on the multi-material arbitrary Lagrangian Eulerian (MM–ALE) technique is then developed to study the influence of the explosive slender ratio on the generated shock waves. The numerical model is validated by experiment data. Then, tests are conducted with constant-mass cylindrical explosives with slender ratios ranging from 2.0 to 9.2. By comparing with the center-detonated spherical explosive, the peak pressure, pulse duration, impulse and energy flux of the underwater explosion of slender explosives are analyzed in detail. Furthermore, the effectiveness of the fluid-structure interaction algorithm is verified through the experimental data of a plate subjected to underwater contact explosion. Finally, the response of plates subjected to the underwater explosion of a slender explosive at different azimuths is studied. The interaction between the plate and the directional loads is analyzed. We conclude that the loads and structure responses of the near-field underwater explosion of slender explosives are different from those of spherical explosives. The analyses and results provide a reference for the near-field underwater explosion loads and structural response studies. Image, graphical abstract [ABSTRACT FROM AUTHOR]

An improved coupling method of smoothed particle hydrodynamics (SPH) and element bending group (EBG) is developed for modeling the interaction of viscous flows with free surface and flexible structures with free and fixed ends. SPH and EBG are both particle methods which are appealing in modeling problems with free surfaces, moving interfaces and large deformations. SPH is used to model viscous fluid, while EBG is used to model flexible structure. Structure particles are also used as moving boundary for SPH, and the interaction of flexible structure with fluid is therefore modeled through the interaction of structure particles and fluid particles. A fixed-end treatment is introduced for flexible structures. A free surface treatment and a surface tension model are used for free surface flow. The improved SPH–EBG method is applied to simulate problems of dam break flow on flexible structures. The good agreement of presented numerical results with existing experimental and numerical results clearly demonstrates the effectiveness of the SPH and EBG coupling approach in modeling fluid–flexible structure interactions. [ABSTRACT FROM AUTHOR]

Patient-specific simulation of the fluid–structure interaction (FSI) problems in cardiovascular systems plays an increasingly important role in fundamental research and clinical applications. However, modeling such problems is challenging, as they often involve non-trivial structural deformation, morphing flow domains, and complex interfaces. In this paper, we develop an incompressible SPH-total Lagrangian SPH (ISPH-TLSPH) framework well-adapted to cardiovascular FSI simulations. The matrix-free iterative ISPH method is used to simulate the hemodynamics and the stabilized TLSPH is used to simulate the dynamics of blood vessels. We propose a novel FSI coupling strategy to couple the ISPH with TLSPH, conforming to strict interface matching conditions. Moreover, we accurately incorporate the lumped parameter (0D) models into the 3D SPH framework to simulate the physiological effects of downstream vascular beds. Lastly, graphics processing unit (GPU) parallelization techniques are implemented in our framework to improve efficiency. The developed framework is first validated by investigating the pulse wave propagation in straight vessels under different boundary conditions. Then the FSI processes in the blood vessel with stenosis and the patient-specific aorta are modeled and investigated. The simulation results show that our framework is effective and efficient for the simulations of patient-specific blood vessels. • A framework that couples the matrix-free iterative ISPH with TLSPH is developed. • A novel coupling strategy for fluid–structure interaction is proposed. • A 3D-0D coupled model is introduced to the SPH framework. • GPU parallelization is implemented for accelerating SPH computations. [ABSTRACT FROM AUTHOR]

As a typical fluid–structure interaction (FSI) problem, water entry involves violent fluid flows and changing free surfaces, which presents great challenges for numerical modeling. Smoothed Particle Hydrodynamics (SPH) is a Lagrangian particle method that has natural advantages in modeling free surfaces and moving interfaces. However, SPH is computationally expensive due to the search of particle–particle interactions, and it causes great difficulties for performing large-scale simulations of 3D FSI problems. In this work, we present an accelerated SPH framework based on the Graphics Processing Unit (GPU) techniques to study water entry problems. The multi-threading programmed by Compute Unified Device Architecture (CUDA) is applied to enhance the computational performance in terms of efficiency and scale. Compared to the Single-CPU-based strategy, the newly presented GPU-accelerated SPH method is computationally more efficient with a speedup over hundreds times and enables a larger memory available for large-scale simulations of around ten million particles for three-dimensional cases. With the GPU-accelerated SPH method, the 3D water entry of a circular cylinder is investigated with some kinematic and dynamic characteristics explained. The results demonstrate that the rotational characteristic of a 3D cylinder in water entry is related to the dimensionless number γ defined as the ratio of the initial inclination angle to the initial velocity angle. The rotation of a cylinder changes from anticlockwise to clockwise with the increase in γ. A transition value of γ exists between the anticlockwise to clockwise rotation, which focuses on the range from 1.0 to 6.0. Meanwhile, the water entry of a 3D circular cylinder leads to a violent impact on the bottom of the cylinder, which causes a peak value of pressure being a maximum value at the early stage of the water entry. It is also indicated that the selection of the initial inclination angle has a great effect on the maximum pressure. • A GPU-accelerated SPH method for fluid–structure interaction problems is developed. • The GPU-accelerated SPH method can greatly improve the computational ability compared to CPU-based SPH. • The water entry characteristics of a 3D cylinder relate to the ratio of the initial inclination angle to velocity angle. • A cross-region that the trajectory of the centroid always passes through is determined by the initial velocity angle. [ABSTRACT FROM AUTHOR]

Numerical simulation of multiphase flows with large density ratios presents considerable challenges. Traditional Smoothed Particle Hydrodynamics (SPH) methods, while efficient in tracking multi-phase interfaces, often suffer from non-physical gaps near interfaces, compromising accuracy. In this paper, we develop an advanced multiphase model with enhanced interface continuity and stability, offering a reliable approach to both multiphase flows and fluid-structure interaction. The model employs equivalent continuity equations across different phases and refines the treatment of pressure and its gradients at interfaces through a smoothing technique, which effectively eliminates the non-physical gaps prevalent in the conventional SPH multiphase models. Moreover, an efficient and accurate method for calculating interface normal vectors and identifying interfaces is proposed. Based on the method, we introduce a series of techniques, including an improved interface force model, a particle displacement strategy, and a new particle penetration detection scheme, all of which significantly improve the interface continuity and the overall simulation accuracy of multiphase flows. Finally, the proposed multiphase SPH model is validated through several numerical examples, including two-dimensional static water, standing waves, dam break, oscillating multi-fluid, and wedge entry, as well as three-dimensional sphere entry. Comparative results with the traditional model indicate that the multiphase SPH model in this work ensures interface continuity and stability in long-period simulations, even under intense flow conditions. These results suggest that the smoothed-interface SPH multiphase model can eliminate non-physical gaps at interfaces, greatly enhancing interface continuity and stability, and highlighting its potential for accurate multiphase flow simulations. [ABSTRACT FROM AUTHOR]

• A novel coupling approach of ES-FEM-SPH is developed for solving thermal-fluid-structure interaction (TFSI) problems. • The updated Lagrangian ES-FEM is developed for solving the coupled thermal elastic problems. • The SPH method, after integrating with particle shifting technique and kernel gradient correction, is robust and effective in modeling thermal fluid flows. • The ghost particle coupling algorithm is developed for treating fluid structure conjugate heat transfer. Thermal-fluid-structure interaction (TFSI) problems are significant in science and engineering, and usually pose great challenges for numerical simulations due to the coupled effects of thermal convection, fluid flow and structure deformation. In this paper, a novel coupling approach of smoothed finite element method (ES-FEM) with an improved smoothed particle hydrodynamic (SPH) method is developed for TFSI problems. In the coupling approach, the edge based ES-FEM is used to model solid domain and the Lagrangian SPH is used to model fluid flow. In ES-FEM, the temperature and velocity gradient smoothing technique are applied over the edge-based smoothing domain for thermal structure coupling problems. In SPH, some state-of-art algorithms including kernel gradient correction (KGC) and particle shift technique (PST) are integrated to ensure computational accuracy for simulating thermal fluid flows. A ghost particle coupling algorithm is developed to handle fluid-structure interaction and fluid-structure conjugate heat transfer, and the kinematic condition, dynamics conditions and conservation of energy are satisfied. Four numerical examples are tested to demonstrate the effectiveness of the present coupling approach of ES-FEM-SPH for TFSI problems. Image, graphical abstract [ABSTRACT FROM AUTHOR]

Existing studies of fluid-structure interaction (FSI) in ocean engineering mainly focus on the interaction between Newtonian fluids and structure. The FSI problems involving non-Newtonian fluids, especially viscoplastic fluids, have rarely been studied while the inherent dynamic behavior is not clear. In this paper, an immersed boundary-lattice Boltzmann method (IB-LBM) is developed for numerical investigations on FSI problems involving viscoplastic fluids. The present IB-LBM is integrated with a hybrid multiple relaxation times (MRT) scheme where different diagonal relaxation matrices are used for modeling Newtonian and non-Newtonian fluids, and are combined in a hybrid manner using a step function to achieve smooth transition for Newtonian to non-Newtonian fluid behavior at the FSI area. Four benchmark problems are used to validate the IB-LBM with hybrid MRT scheme. It is demonstrated that the numerical model can avoid numerical instability when modeling viscoplastic fluid flow and reduce the numerical boundary slip in the IB-LBM. The numerical model is further used to study the viscoplastic fluid flow around a fixed and moving cylinder (or particle). We show that the present IB-LBM with the hybrid MRT scheme is effective in modeling FSI involving viscoplastic fluids while the obtained phenomena are quite different from those with Newtonian fluids. [ABSTRACT FROM AUTHOR]