Fluid-structure Interaction Approach Research Articles

Breakage risk analysis during the glass tempering process based on fluid–structure interaction approach

Breakage risk analysis during the glass tempering process based on fluid–structure interaction approach

Hydroelastic modelling of a deformable wave energy converter including power take-off

Given the advantages of flexible wave energy converters (FlexWECs), such as deformation-led energy harnessing and structural loading compliance, there has been a significant interest in FlexWECs in both academia and industries. To simulate the FlexWEC interaction with ocean surface waves, a 3D computational fluid-structure interaction approach is developed in this study. The fluid and solid governing equations are discretized using finite difference and finite element methods, respectively. An immersed boundary method is used to couple the two independent grid systems. A novel numerical technique is introduced to model the dielectric elastomer generator (DEG) as the power take-off (PTO). The wave energy capture performance is analysed for different PTO configurations and at various wave conditions. Based on the obtained results, the PTO damping coefficient and the relative wavelength range that maximizes the capture width ratio (CWR) are determined. The wavefield results also reveal the presence of wave-height enhancement and attenuation points around a single FlexWEC, providing potential site selection references when deploying multiple FlexWECs in an array.

Toward a physiological model of vascular wall vibrations in the arteriovenous fistula.

The mechanism behind hemodialysis arteriovenous fistula (AVF) failure remains poorly understood, despite previous efforts to correlate altered hemodynamics with vascular remodeling. We have recently demonstrated that transitional flow induces high-frequency vibrations in the AVF wall, albeit with a simplified model. This study addresses the key limitations of our original fluid-structure interaction (FSI) approach, aiming to evaluate the vibration response using a more realistic model. A 3D AVF geometry was generated from contrast-free MRI and high-fidelity FSI simulations were performed. Patient-specific inflow and pressure were incorporated, and a three-term Mooney-Rivlin model was fitted using experimental data. The viscoelastic effect of perivascular tissue was modeled with Robin boundary conditions. Prescribing pulsatile inflow and pressure resulted in a substantial increase in vein displacement ( %) and strain ( %), with a higher maximum spectral frequency becoming visible above -42 dB (from 200 to 500 Hz). Transitioning from Saint Venant-Kirchhoff to Mooney-Rivlin model led to displacement amplitudes exceeding 10 micrometers and had a substantial impact on strain ( %). Robin boundary conditions significantly damped high-frequency displacement ( %). Incorporating venous tissue properties increased vibrations by 91%, extending up to 700 Hz, with a maximum strain of 0.158. Notably, our results show localized, high levels of vibration at the inner curvature of the vein, a site known for experiencing pronounced remodeling. Our findings, consistent with experimental and clinical reports of bruits and thrills, underscore the significance of incorporating physiologically plausible modeling approaches to investigate the role of wall vibrations in AVF remodeling and failure.

Fluid-structure interaction simulation of mechanical aortic valves: a narrative review exploring its role in total product life cycle.

Over the last years computer modelling and simulation has emerged as an effective tool to support the total product life cycle of cardiovascular devices, particularly in the device preclinical evaluation and post-market assessment. Computational modelling is particularly relevant for heart valve prostheses, which require an extensive assessment of their hydrodynamic performance and of risks of hemolysis and thromboembolic complications associated with mechanically-induced blood damage. These biomechanical aspects are typically evaluated through a fluid-structure interaction (FSI) approach, which enables valve fluid dynamics evaluation accounting for leaflets movement. In this context, the present narrative review focuses on the computational modelling of bileaflet mechanical aortic valves through FSI approach, aiming to foster and guide the use of simulations in device total product life cycle. The state of the art of FSI simulation of heart valve prostheses is reviewed to highlight the variety of modelling strategies adopted in the literature. Furthermore, the integration of FSI simulations in the total product life cycle of bileaflet aortic valves is discussed, with particular emphasis on the role of simulations in complementing and potentially replacing the experimental tests suggested by international standards. Simulations credibility assessment is also discussed in the light of recently published guidelines, thus paving the way for a broader inclusion of in silico evidence in regulatory submissions. The present narrative review highlights that FSI simulations can be successfully framed within the total product life cycle of bileaflet mechanical aortic valves, emphasizing that credible in silico models evaluating the performance of implantable devices can (at least) partially replace preclinical in vitro experimentation and support post-market biomechanical evaluation, leading to a reduction in both time and cost required for device development.

Study on inertial migration behaviours under secondary flow patterns using a fluid-structure interaction approach

Study on inertial migration behaviours under secondary flow patterns using a fluid-structure interaction approach

Fluid–structure interaction of a marine rudder at incidence in the wake of a propeller

The structural response of a rudder in the wake of a marine propeller is investigated by one-way fluid–structure interaction approach. The unsteady pressure field gathered by detached eddy simulations is provided to a structural solver for the computation of deformations and stresses of the rudder. The study compares the structural response of the rudder at neutral and two equal and opposite rotations, which are representative of design conditions in straight motion and maneuvering conditions that are experienced under the action of the autopilot for course control or weak maneuvering. The analysis sheds light on the different structural behavior at the two opposite rotation angles, caused by asymmetrical variation along the span of the rudder of the angle of attack induced by the propeller slipstream, by considering the different role played by the tip and hub vortical systems. The test case consists of a rudder with a rectangular plane area and National Advisory Committee for Aeronautics 0015 sectional profile located past the E779A propeller. The propeller operates at low loading conditions, and the rudder is set at incidence δ=0∘,±4∘. The study shows that the response of the rudder is driven by flap and torsion and is asymmetric for the two and opposite rotations. The mean deformation and vibratory response are magnified for δ=−4∘ by at least 70% and 20% for the lateral and edgewise deflections, respectively, with respect to the opposite rudder incidence. In general, the excitation generated by the tip vortex is stronger than that of the hub vortex. In the most critical condition, at δ=−4∘, the excitation associated with the tip vortex is nearly double that of the hub vortex.

Hemodynamics and wall shear metrics in a pulmonary autograft: Comparing a fluid-structure interaction and computational fluid dynamics approach

ObjectiveIn young patients, aortic valve disease is often treated by placement of a pulmonary autograft (PA) which adapts to its new environment through growth and remodeling. To better understand the hemodynamic forces acting on the highly distensible PA in the acute phase after surgery, we developed a fluid-structure interaction (FSI) framework and comprehensively compared hemodynamics and wall shear-stress (WSS) metrics with a computational fluid dynamic (CFD) simulation. MethodsThe FSI framework couples a prestressed non-linear hyperelastic arterial tissue model with a fluid model using the in-house coupling code CoCoNuT. Geometry, material parameters and boundary conditions are based on in-vivo measurements. Hemodynamics, time-averaged WSS (TAWSS), oscillatory shear index (OSI) and topological shear variation index (TSVI) are evaluated qualitatively and quantitatively for 3 different sheeps. ResultsDespite systolic-to-diastolic volumetric changes of the PA in the order of 20 %, the point-by-point correlation of TAWSS and OSI obtained through CFD and FSI remains high (r > 0.9, p < 0.01) for TAWSS and (r > 0.8, p < 0.01) for OSI). Instantaneous WSS divergence patterns qualitatively preserve similarities, but large deformations of the PA leads to a decrease of the correlation between FSI and CFD resolved TSVI (r < 0.7, p < 0.01). Moderate co-localization between FSI and CFD is observed for low thresholds of TAWSS and high thresholds of OSI and TSVI. ConclusionFSI might be warranted if we were to use the TSVI as a mechano-biological driver for growth and remodeling of PA due to varying intra-vascular flow structures and near wall hemodynamics because of the large expansion of the PA.

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use.

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential. This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed. Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data. Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Adaptive Immersed Mesh Method (AIMM) for Fluid–Structure Interaction

Our paper proposes an innovative approach for modeling Fluid–Structure Interaction (FSI). Our method combines both traditional monolithic and partitioned approaches, creating a hybrid solution that facilitates FSI. At each time iteration, the solid mesh is immersed within a fluid–solid mesh, all while maintaining its independent Lagrangian hyperelastic solver. The Eulerian mesh encompasses both the fluid and solid components and accommodates various physical phenomena. We enhance the interaction between solid and fluid through anisotropic mesh adaptation and the Level-Set methods. This enables a more accurate representation of their interaction. Together, these components constitute the Adaptive Immersed Mesh Method (AIMM). For both solvers, we utilize the Variational Multi-Scale (VMS) method, mitigating potential spurious oscillations common with piecewise linear tetrahedral elements. The framework operates in 3D with parallel computing capabilities. Our method’s accuracy, robustness, and capabilities are assessed through a series of 2D numerical problems. Furthermore, we present various three-dimensional test cases and compare their results to experimental data.

Stability Evaluation of a Simplified Reclaimer Using Fluid–Structure Interaction

This paper describes the numerical modeling and simulation of the wind effects on an ore reclaimer structure using analytical and numerical methods. The physical model is a large machine with a height of 34 m and a width of 77 m, and due to its complexity, a simplified model was used. This study aims to investigate the influence of wind speed on the stability failure of the reclaimer and to provide a more efficient and precise stopping criterion. The simulations were performed using a two-way FSI (fluid–structure interaction) approach. An FSI analysis was performed to study the dynamic behavior of a numerical model consisting of two separate parts with contact constraints. This article also highlights the importance of FSI in improving the reliability of the stability failure. Finally, the numerical results showed differences compared to the analytical model, and the wind load limit to stability failure was observed at higher wind speeds. The structure was able to support wind velocities higher than suggested by the FEM (European Materials Handling Federation) standard without stability fails.

An Investigation on an SMA Valve in Micro-Channel Through FSI Approach

Micro-channels with built-in shape memory alloys (SMAs) can respond intelligently to channel temperature variations without applying an external force and control the flow rate. In this work, a 2D micro-channel with a shape memory alloy (SMA) valve is designed and studied using fluid–structure interaction (FSI) conditions. The behavior of the SMA valve, which includes bending under the effect of force loading and inlet fluid pressure, and the recovery of the shape of the valve, has been investigated in this study. The temperature at which the valve is completely recovered and the output fluid reaches the lowest value is called the closing temperature. This research investigates the impact of different geometrical and functional parameters on closing temperature, such as fluid inlet pressure and velocity, channel height and valve thickness. The outlet fluid flow rate of the SMA micro-channel is evaluated by employing a relevant constitutive equation to model the SMA valve behavior and implement it in Finite Element Method (FEM) multiphysics frameworks. The results show that the outlet flow rate relies on the shape recovery of the SMA valve and the fluid velocity field in the FSI simulations.

Investigating hermetic reciprocating compressor performance by using various machine learning methods

Due to their durability and efficiency, hermetic reciprocating compressors (HRCs) are used in refrigeration and air conditioning. Compressor performance and reliability concerns reduce system efficiency and raise maintenance costs. Machine learning (ML) is being used to improve hermetic reciprocating compressor performance, reliability, and energy economy. ML is used in hermetic reciprocating compressors for issue identification, performance improvement, predictive maintenance, and energy management. This research compared HRC performance factors such as mass flow rate, cooling capacity, compression power, coefficient of performance, exhaust line losses, and volumetric efficiency. Simple regression, probabilistic neural network, gradient boosted, polynomial regression, and random forest (RF) were used to examine and evaluate these parameters as outputs. Over three cycles, the Fluid-Structure Interaction (FSI) approach assessed compressor performance parameters. For compressor speeds of 1300, 2100, and 3000 rpm, mass flow rate, compression power, cooling efficiency coefficient, and exhaust line energy losses varied by 10%, 4%, 5%, and 6%. To gather ML algorithm inputs, the research used experimental, fluid-structure interaction, and ML methodologies. Experimental and FSI approaches produced 108 data points. These data points were randomly assigned, with 70% for learning and 30% for prediction. The mean convergence criterion for mass flow rate, cooling capacity, compression power, cooling efficiency coefficient, exhaust line energy losses, and volumetric efficiency parameters was 0.9966, 0.9969, 0.9572, 0.0561, 0.9925, and 0.4640 for all ML methods. Simple regression, probabilistic neural networks, gradient boosted, polynomial regression, and RF convergence criteria were 0.8978, 0.9999, 0.6016, 0.4439, and 0.7761.

A general FSI framework for an effective stress analysis on composite wind turbine blades

The fluid-structure interaction (FSI) technique has been extensively used and developed in the past decades. Commonly, the reduced-order models are used in FSI analyses to assure the numerical robustness and efficiency. However, due to the increasing demand for higher numerical resolutions in modern wind turbine composite blade applications, intrinsic limitations of reduced-order models, such as their inability to account for complex aerodynamic flow interactions, multi-motion couplings, and sophisticated composite properties, have become the weaknesses in existing reduced-order FSI approaches. In this study, we propose a general FSI framework, which combines the advantages of high-fidelity Computational Fluid Dynamics (CFD) and robust Multi-Body Dynamics (MBD) methods, and detailed Finite Element Analysis (FEA) for analysing the detailed stress distributions on the composite structures. The results of predicted dynamics and the Von Mises stress on the composite blade structures under given operation condition are compared and reasonably agreed with the literature results, with a significant computational cost reduction by nearly 25% is achieved. The proposed FSI framework can be a general approach to investigate the multi-physical interactions where the composite structure specifications are involved, coupling with complex dynamic motions in three-dimensional space.

Non-Newtonian pulsatile blood flow through the stenosed arteries: comparison between the viscoelastic and elastic arterial wall in response to the alterations

:In this study, we investigate the impact of aortic stenosis on the hemodynamics of pulsatile blood flow within a 3D aortic model. Employing a non-Newtonian Casson model with a hematocrit of 45%, our study introduces a preliminary hypothesis to simulate blood flow dynamics, incorporating both linear elastic and viscoelastic models to define the mechanical characteristics of the artery. Through simulations conducted with Ansys-Cfx (version 15), we utilize a 2-way fluid-structure interaction (FSI) approach, employing a Lagrangian-Eulerian formulation with second-order accuracy. We explore the influence of stenosis severity on variables including velocity profiles, pressure distribution, shear stress, wall displacement, and changes in the OSI parameter. Our investigation encompasses arteries with both elastic and viscoelastic walls. The key findings that arise from our results highlight the viscoelastic model’s demonstration of reduced radial wall displacement when compared to the linear elastic model. Additionally, we observe that elevated arterial stenosis percentages lead to the elongation of vortex length, heightened wall shear stress, and increased slope of velocity profiles downstream of the stenosed region. Furthermore, bulky obstruction of viscoelastic arteries as opposed to elastic, resulted in a maximum 5 percent increase in velocity profile and a 29.6% decrease in radial displacement. The zenith of shear stress occurs concomitantly with the velocity’s peak within the stenosed area. Viscoelastic arterial wall shear stress at the stenosis site escalates due to the rapid expansion of the stenosis. The viscoelastic wall, responding with a blend of viscous and elastic characteristics to applied stress, undergoes slight deformation in shape. Following stress reduction, the wall gradually reverts to its original form, thus alleviating some of the applied stress. In contrast, the elastic wall retains its altered shape due to stress preservation within the material. Additionally, we ascertain an augmentation in radial displacement corresponding with increased artery stenosis.

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ∼200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. If the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

Effect of stenotic shapes and arterial wall elasticity on the hemodynamics

The present study employs an arbitrary Lagrangian–Eulerian fluid–structure interaction approach to investigate pulsatile blood flow through a deformable stenosed channel. The flow is modeled by solving the incompressible continuity and momentum equations using finite element-based commercial solver COMSOL Multiphysics®. In this work, we explore the effects of different stenotic shapes—elliptical, round, and sinusoidal, degrees of stenosis (30%, 50%, and 70%), and arterial wall stiffnesses—0.5, 1.5, and 2.5 MPa on the velocity profile, pressure and wall shear stress distribution, and wall deformation. The oscillatory shear index (OSI) is analyzed to predict further plaque formation in the stenosed artery. We find that the flow velocity, wall shear stress, and pressure difference across the stenosed region increase with an increase in the stenotic severity and artery stiffness. The velocity profiles intersect at a radial location in the stenotic region termed critical radius, where relative magnitudes get reversed. With the increase in stenotic severity, the wall displacement decreases at the throat and increases at the upstream side. With the increase in wall stiffness, the wall deformation decreases, and shear stress increases, thereby increasing the pressure drop across the stenosed region. At a lower mass flow rate and a higher degree of stenosis, the vortices are formed upstream and downstream of the stenosed region for all stenotic shapes. The vorticity magnitude is found to be more than 21% higher for sinusoidal stenotic shape than round and elliptical ones. The effect of stenotic profile on the pressure drop characteristics shows that blood experiences maximum wall shear stress for the sinusoidal stenotic geometry, whereas the pressure drop is the maximum for the elliptical stenotic shape. The elliptical stenotic shape is more prone to further plaque formation than round and sinusoidal stenotic shapes. At lower Womersley number (Wo=2.76) corresponding to 60 beats per min heart beat rate, secondary vortices are formed downstream of the channel, causing higher OSI.

Numerical Investigations of a Two-Way Coupled Fluid–Structure Interaction Approach for Fast Transients in Fluid-Filled Flexible Piping Systems

Numerical Investigations of a Two-Way Coupled Fluid–Structure Interaction Approach for Fast Transients in Fluid-Filled Flexible Piping Systems

Structural optimisation of H-Rotor wind turbine blade based on one-way Fluid Structure Interaction approach

The paper aims to present structural topology optimisation of the structural topology of the H-Rotor wind turbine combined with the one-way Fluid Structure Interaction (FSI) approach. The developed methodology couples the set of Unsteady Reynolds Averaged Navier–Stokes Equations (URANS) with the steady state linear elasticity equations and density based topology optimisation method. This approach allowed us to reduce the mass of the blades with respect to the maximum stress and deformation limits. The results obtained provide detailed information on the unsteady flow fields around the operating wind turbine, as well as optimised topology of the blade interior without affecting the external aerodynamic profile. The load of pressure profile from the CFD simulation was implemented at the point of rotation at which the force values observed on the blade surface were the highest. The centrifugal force was considered in the structural model. Computations were carried out for wind speed equal to 30 m/s and rotational velocity equal to 90 rad/s. The turbine operating parameters were chosen to represent extreme operating conditions. The results obtained showed that the blade mass, made of the same material, can be easily reduced by around 60% with the reduction of the observed stress values and deformation.

Forest density is more effective than tree rigidity at reducing the onshore energy flux of tsunamis

Communities around the world are increasingly interested in nature-based solutions to the mitigation of coastal risks by coastal forests, but it remains unclear how much protective benefits vegetation provides, particularly in the limit of highly energetic flows after tsunami impact. The current study, using a three-dimensional incompressible computational fluid dynamics model with a fluid–structure interaction approach, aims to quantify how energy reflection and dissipation vary with different degrees of rigidity and vegetation density of a coastal forest. We represent tree trunks as cylinders and use the elastic modulus of hardwood trees such as pine or oak to characterize the rigidity of these cylinders. The numerical results show that energy reflection increases with rigidity only for a single cylinder. In the presence of multiple cylinders, the difference in energy reflection created by varying rigidity diminishes as the number of cylinders increases. Instead of rigidity, we find that the blockage area created by the presence of multiple tree trunks dominates energy reflection. As tree trunks are deformed by the hydrodynamic forces, they alter the flow field around them, causing turbulent kinetic energy generation in the wake region. As a consequence, trees dissipate flow energy, highlighting coastal forests reducing the onshore energy flux of tsunamis by means of both reflection and dissipation.

Deformation features and failure mechanism of subsea shield tunnels with different burial depths crossing fault-zone

Shield tunneling processes can result in accidents like water inrush and collapse due to instability in the tunnel face. This paper conducts a systematic investigation of the deformation features and failure mechanism of subsea shield tunnels with different burial depths crossing fault zone. The research entails combing theoretical analysis, numerical simulation, and data mining with the second subsea tunnel in Jiaozhou Bay, Qingdao. The impact of the shallow and deep burial depths on the damage modes of the shield tunnel face is analyzed from the perspective of the pressure arch effect. When the tunnel burial depth exceeds the limit, the overlying rock of the tunnel will form a pressure arch that supports the overlying rock. The study investigates the impact of various factors on the tunnel face’s failure modes using a fluid-structure interaction approach. Among them, as the thickness of the overlying bearing stratum increased, the displacement of the tunnel face initially decreased but later increased. The increase of cohesion and internal friction angle contributed to decreasing the displacement of the tunnel face, while seawater depth and tunnel diameter exhibited diverse effects. In the grey correlation analysis, the seawater depth exhibited the most sensitive impact on the pressure arch height.