OpenIFEM:A High Performance Modular Open-Source Software of the Immersed Finite Element Method for Fluid-Structure Interactions

We present a high performance modularly-built open-source software-OpenIFEM.OpenIFEM is a C++implementation of the modified immersed finite element method(mIFEM)to solve fluid-structure interaction(FSI)problems.This software is modularly built to perform multiple tasks including fluid dynamics(incompressible and slightly compressible fluid models),linear and nonlinear solid mechanics,and fully coupled fluid-structure interactions.Most of open-source software packages are restricted to certain discretization methods;some are under-tested,under-documented,and lack modularity as well as extensibility.OpenIFEM is designed and built to include a set of generic classes for users to adapt so that any fluid and solid solvers can be coupled through the FSI algorithm.In addition,the package utilizes well-developed and tested libraries.It also comes with standard test cases that serve as software and algorithm validation.The software can be built on cross-platform,i.e.,Linux,Windows,and Mac OS,using CMake.Efficient parallelization is also implemented for high-performance computing for large-sized problems.OpenIFEM is documented using Doxygen and publicly available to download on GitHub.It is expected to benefit the future development of FSI algorithms and be applied to a variety of FSI applications.

Parallelization strategies for resolved simulations of fluid-structure-particle interactions

Fluid-structure-particle interactions in three spatial dimensions happen in many environmental and engineering flows.This paper presents the parallel algorithms for the hybrid diffuse and sharp interface immersed boundary(IB)method developed in our previous work.For the moving structure modeled using the sharp interface IB method,a recursive box method is developed for efficiently classifying the background grid nodes.For the particles modeled using the diffuse interface IB method,a‘master-slave’approach is adopted.For the particle-particle interaction(PPI)and particle-structure interaction(PSI),a fast algorithm for classifying the active and inactive Lagrangian points,which discretize the particle surface,is developed for the‘dry’contact approach.The results show that the proposed recursive box method can reduce the classifying time from 52seconds to 0.3 seconds.Acceptable parallel efficiency is obtained for cases with different particle concentrations.Furthermore,the lubrication model is utilized when a particle approaches a wall,enabling an accurate simulation of the rebounding phenomena in the benchmark particle-wall collision problem.At last,the capability of the proposed computational framework is demonstrated by simulating particle-laden turbulent channel flows with rough walls.

Blade Wrap Angle Impact on Centrifugal Pump Performance:Entropy Generation and Fluid-Structure Interaction Analysis

The centrifugal pump is a prevalent power equipment widely used in different engineering patterns,and the impeller blade wrap angle significantly impacts its performance.A numerical investigation was conducted to analyze the influence of the blade wrap angle on flow characteristics and energy distribution of a centrifugal pump evaluated as a low specific speed with a value of 69.This study investigates six impellermodels that possess varying blade wrap angles(95°,105°,115°,125°,135°,and 145°)that were created while maintaining the same volute and other geometrical characteristics.The investigation of energy loss was conducted to evaluate the values of total and entropy generation rates(TEG,EGR).The fluid-structure interaction was considered numerically using the software tools ANSYS Fluent and ANSYSWorkbench.The elastic structural dynamic equation was used to estimate the structural response,while the shear stress transport k–ωturbulence model was utilized for the fluid domain modeling.The findings suggest that the blade wrap angle has a significant influence on the efficiency of the pump.The impeller featuring a blade wrap angle of 145°exhibits higher efficiency,with a notable increase of 3.76%relative to the original model.Variations in the blade wrap angle impact the energy loss,shaft power,and pump head.The model with a 145°angle exhibited a maximum equivalent stress of 14.8MPa and a total deformation of 0.084 mm.The results provide valuable insights into the intricate flow mechanism of the centrifugal pump,particularly when considering various blade wrap angles.

Tailoring Light–Matter Interactions in Overcoupled Resonator for Biomolecule Recognition and Detection

Plasmonic nanoantennas provide unique opportunities for precise control of light–matter coupling in surface-enhanced infrared absorption(SEIRA)spectroscopy,but most of the resonant systems realized so far suffer from the obstacles of low sensitivity,narrow bandwidth,and asymmetric Fano resonance perturbations.Here,we demonstrated an overcoupled resonator with a high plasmon-molecule coupling coefficient(μ)(OC-Hμresonator)by precisely controlling the radiation loss channel,the resonator-oscillator coupling channel,and the frequency detuning channel.We observed a strong dependence of the sensing performance on the coupling state,and demonstrated that OC-Hμresonator has excellent sensing properties of ultra-sensitive(7.25%nm^(−1)),ultra-broadband(3–10μm),and immune asymmetric Fano lineshapes.These characteristics represent a breakthrough in SEIRA technology and lay the foundation for specific recognition of biomolecules,trace detection,and protein secondary structure analysis using a single array(array size is 100×100μm^(2)).In addition,with the assistance of machine learning,mixture classification,concentration prediction and spectral reconstruction were achieved with the highest accuracy of 100%.Finally,we demonstrated the potential of OC-Hμresonator for SARS-CoV-2 detection.These findings will promote the wider application of SEIRA technology,while providing new ideas for other enhanced spectroscopy technologies,quantum photonics and studying light–matter interactions.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

T cell interactions with microglia in immune-inflammatory processes of ischemic stroke

The primary mechanism of secondary injury after cerebral ischemia may be the brain inflammation that emerges after an ischemic stroke,which promotes neuronal death and inhibits nerve tissue regeneration.As the first immune cells to be activated after an ischemic stroke,microglia play an important immunomodulatory role in the progression of the condition.After an ischemic stroke,peripheral blood immune cells(mainly T cells)are recruited to the central nervous system by chemokines secreted by immune cells in the brain,where they interact with central nervous system cells(mainly microglia)to trigger a secondary neuroimmune response.This review summarizes the interactions between T cells and microglia in the immune-inflammatory processes of ischemic stroke.We found that,during ischemic stroke,T cells and microglia demonstrate a more pronounced synergistic effect.Th1,Th17,and M1 microglia can co-secrete proinflammatory factors,such as interferon-γ,tumor necrosis factor-α,and interleukin-1β,to promote neuroinflammation and exacerbate brain injury.Th2,Treg,and M2 microglia jointly secrete anti-inflammatory factors,such as interleukin-4,interleukin-10,and transforming growth factor-β,to inhibit the progression of neuroinflammation,as well as growth factors such as brain-derived neurotrophic factor to promote nerve regeneration and repair brain injury.Immune interactions between microglia and T cells influence the direction of the subsequent neuroinflammation,which in turn determines the prognosis of ischemic stroke patients.Clinical trials have been conducted on the ways to modulate the interactions between T cells and microglia toward anti-inflammatory communication using the immunosuppressant fingolimod or overdosing with Treg cells to promote neural tissue repair and reduce the damage caused by ischemic stroke.However,such studies have been relatively infrequent,and clinical experience is still insufficient.In summary,in ischemic stroke,T cell subsets and activated microglia act synergistically to regulate inflammatory progression,mainly by secreting inflammatory factors.In the future,a key research direction for ischemic stroke treatment could be rooted in the enhancement of anti-inflammatory factor secretion by promoting the generation of Th2 and Treg cells,along with the activation of M2-type microglia.These approaches may alleviate neuroinflammation and facilitate the repair of neural tissues.

Smoothed particle hydrodynamics and its applications in fluid-structure interactions

In ocean engineering, the applications are usually related to a free surface which brings so many interesting physical phenomena （e.g. water waves, impacts, splashing jets, etc.）. To model these complex free surface flows is a tough and challenging task for most computational fluid dynamics （CFD） solvers which work in the Eulerian framework. As a Lagrangian and meshless method, smoothed particle hydrodynamics （SPH） offers a convenient tracking for different complex boundaries and a straightforward satisfaction for different boundary conditions. Therefore SPH is robust in modeling complex hydrodynamic problems characterized by free surface boundaries, multiphase interfaces or material discontinuities. Along with the rapid development of the SPH theory, related numerical techniques and high-performance computing technologies, SPH has not only attracted much attention in the academic community, but also gradually gained wide applications in industrial circles. This paper is dedicated to a review of the recent developments of SPH method and its typical applications in fluid-structure interactions in ocean engineering. Different numerical techniques for improving numerical accuracy, satisfying different boundary conditions, improving computational efficie- ncy, suppressing pressure fluctuations and preventing the tensile instability, etc., are introduced. In the numerical results, various typical fluid-structure interaction problems or multiphase problems in ocean engineering are described, modeled and validated. The prospective developments of SPH in ocean engineering are also discussed.

Smoothed particle hydrodynamics(SPH) for modeling fluid-structure interactions

Fluid-structure interaction(FSI) is a class of mechanics-related problems with mutual dependence between the fluid and structure parts and it is observable nearly everywhere, in natural phenomena to many engineering systems. The primary challenges in developing numerical models with conventional grid-based methods are the inherent nonlinearity and timedependent nature of FSI, together with possible large deformations and moving interfaces. Smoothed particle hydrodynamics(SPH) method is a truly Lagrangian and meshfree particle method that conveniently treats large deformations and naturally captures rapidly moving interfaces and free surfaces. Since its invention, the SPH method has been widely applied to study different problems in engineering and sciences, including FSI problems. This article presents a review of the recent developments in SPH based modeling techniques for solving FSI-related problems. The basic concepts of SPH along with conventional and higher order particle approximation schemes are first introduced. Then, the implementation of FSI in a pure SPH framework and the hybrid approaches of SPH with other grid-based or particle-based methods are discussed. The SPH models of FSI problems with rigid, elastic and flexible structures, with granular materials, and with extremely intensive loadings are demonstrated. Some discussions on several key techniques in SPH including the balance of accuracy, stability and efficiency, the treatment of material interface, the coupling of SPH with other methods, and the particle regularization and adaptive particle resolution are provided as concluding marks.

Numerical simulations of the discontinuous progression of cerebral aneurysms based on fluid-structure interactions study

Investigations into the characteristics of hemodynamics will provide a better understanding of the pathology of cerebral aneurysms for clinicians.In this work,a steady state discontinuous-growth model of the cerebral aneurysms was proposed.With the assumption of the fluid-structure interaction between the wall of blood vessel and blood,a fluid-structure coupling numerical simulation for this model was built using software ANSYS and CFX.The simulation results showed that as the aneurysm volume increased,a blood flow vortex came forth,the vortex region became asymptotically larger,and eddy density became gradually stronger in it.After the emergence of the vortex region,the blood flow in the vicinity of the downstream in the aneurysms volume turned into bifurcated flow,and the location of the flow bifurcated point was shifted with the aneurysm volume growing while directions of the shear stress applied to two sides of the bifurcated point were opposite.The Von Mises stress distribution along the wall of aneurysm volume decreased in the prior period and increased gradually in the later period.The maximum stress was in the neck of the volume and the minimum was on the distal end in the whole process of growth.It was shown that as the aneurysm increased the maximum deformation location of the aneurysm,vertical to the streamline,was transferred from the distal end of the aneurysm to its neck,then back to its distal end of the aneurysm again.

An improved algorithm for fluid-structure interaction of high-speed trains under crosswind

Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome （OSAHS）, three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.

Hybrid algorithm for modeling of fluid-structure interaction in incompressible, viscous flows

The objective of this paper is to present and to validate a new hybrid coupling （HC） algorithm for modeling of fluid-structure interaction （FSI） in incompressible, viscous flows. The HC algorithm is able to avoid numerical instability issues associated with artificial added mass effects, which are often encountered by standard loosely coupled （LC） and tightly coupled （TC） algorithms, when modeling the FSI response of flexible structures in incompressible flow. The artificial added mass effect is caused by the lag in exchange of interfacial displacements and forces between the fluid and solid solvers in partitioned algorithms. The artificial added mass effect is much more prominent for light/flexible struc- tures moving in water, because the fluid forces are in the same order of magnitude as the solid forces, and because the speed at which numerical errors propagate in an incom- pressible fluid. The new HC algorithm avoids numerical instability issues associated with artificial added mass effects by embedding Theodorsen＇s analytical approximation of the hydroelastic forces in the solution process to obtain better initial estimates of the displacements. Details of the new HC algorithm are presented. Numerical validation studies are shown for the forced pitching response of a steel and a plastic hydrofoil. The results show that the HC algorithm is able to converge faster, and is able to avoid numerical insta- bility issues, compared to standard LC and TC algorithms, when modeling the transient FSI response of a plastic hydrofoil. Although the HC algorithm is only demonstrated for a NACA0009 hydrofoil subject to pure pitching motion, the method can be easily extended to model general 3-D FSI response and stability of complex, flexible structures in turbulent, incompressible, multiphase flows.

NUMERICAL SOLUTION OF FLUID-STRUCTURE INTERACTION IN LIQUID-FILLED PIPES BY METHOD OF CHARACTERISTICS

Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.

Fluid-structure interaction simulation of three-dimensional flexible hydrofoil in water tunnel

The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.

3D numerical simulation on fluid-structure interaction of structure subjected to underwater explosion with cavitation

In the underwater-shock environment, cavitation occurs near the structural surface. The dynamic response of fluid-structure interactions is influenced seriously by the cavitation effects. It is also the difficulty in the field of underwater explosion. With the traditional boundary element method and the finite element method （FEM）, it is difficult to solve the nonlinear problem with cavitation effects subjected to the underwater explosion. To solve this problem, under the consideration of the cavitation effects and fluid compressibility, with fluid viscidity being neglected, a 3D numerical model of transient nonlinear fluid-structure interaction subjected to the underwater explosion is built. The fluid spectral element method （SEM） and the FEM are adopted to solve this model. After comparison with the FEM, it is shown that the SEM is more precise than the FEM, and the SEM results are in good coincidence with benchmark results and experiment results. Based on this, combined with ABAQUS, the transient fluid-structure interaction mechanism of the 3D submerged spherical shell and ship stiffened plates subjected to the underwater explosion is discussed, and the cavitation region and its influence on the structural dynamic responses are presented. The paper aims at providing references for relevant research on transient fluid-structure interaction of ship structures subjected to the underwater explosion.

Dynamic Adaptive Finite Element Analysis of Acoustic Wave Propagation Due to Underwater Explosion for Fluid-structure Interaction Problems

In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.

Effects of renal artery stenosis on realistic model of abdominalaorta and renal arteries incorporating fluid-structureinteraction and pulsatile non-Newtonian blood flow

The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.

Fluid-Structure Interaction in Problems of Patient Specific Transcatheter Aortic Valve Implantation with and Without Paravalvular Leakage Complication

Paravalvular Leakage(PVL)has been recognized as one of the most dangerous complications in relation to Transcathether Aortic Valve Implantation(TAVI)activities.However,data available in the literature about Fluid Structure Interaction(FSI)for this specific problem are relatively limited.In the present study,the fluid and structure responses of the hemodynamics along the patient aorta model and the aortic wall deformation are studied with the aid of numerical simulation taking into account PVL and 100%TAVI valve opening.In particular,the aorta without valve(AWoV)is assumed as the normal condition,whereas an aorta with TAVI 26 mm for 100%Geometrical Orifice Area(GOA)is considered as the patient aorta with PVL complication.A 3D patient-specific aorta model is elaborated using the MIMICS software.Implantation of the identical TAVI valve of Edward SAPIEN XT 26(Edwards Lifes ciences,Irvine,California)is considered.An undersized 26 mm TAVI valve with 100%valve opening is selected to mimic the presence of PVL at the aortic annulus.The present research indicates that the existence of PVL can increase the blood velocity,pressure drop and WSS in comparison to normal conditions,thereby paving the way to the development of recirculation flow,thrombus formation,aorta wall collapse,aortic rupture and damage of endothelium.

Improved frequency modeling and solution for parallel liquid-filled pipes considering both fluid-structure interaction and structural coupling

The dynamic characteristics of a single liquid-filled pipe have been broadly studied in the previous literature.The parallel liquid-filled pipe(PLFP)system is also widely used in engineering,and its structure is more complex than that of a single pipe.However,there are few reports about the dynamic characteristics of the PLFPs.Therefore,this paper proposes improved frequency modeling and solution for the PLFPs,involving the logical alignment principle and coupled matrix processing.The established model incorporates both the fluid-structure interaction(FSI)and the structural coupling of the PLFPs.The validity of the established model is verified by modal experiments.The effects of some unique parameters on the dynamic characteristics of the PLFPs are discussed.This work provides a feasible method for solving the FSI of multiple pipes in parallel and potential theoretical guidance for the dynamic analysis of the PLFPs in engineering.

DYNAMIC SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS INVOLVING LARGE-AMPLITUDE SLOSHING

An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure is modeled as a rigid container supportedby a system consisting of springs and dashpots. The motion of the fluid is decomposed into twoparts: the large-displacement motion with the container and the large-amplitude sloshing relative tothe container. The former is conveniently dealt with by defining a container-fixed noninertiallocal frame, while the latter is easily handled by adopting an ALE kinematical description. Thisleads to an easy and accurate treatment of both the fluid-structure interface and the fluid freesurface without producing excessive distortion of the computational mesh. The coupling between thefluid and the structure is accomplished through the coupling matrices that can be easilyestablished. Two numerical examples, including a TLD-structure system and a simplified liquid-loadedvehicle system, are presented to demonstrate the effectiveness and reliability of the proposedmethod. The present work can also be applied to simulate fluid-structure problems incorporatingmultibody systems and several fluid domains.