Fluid-structure Interaction Research Articles

Numerical Prediction of Two-Dimensional Coupled Galloping and Vortex-Induced Vibration of Square Cylinders Under Symmetric/Asymmetric Flow Orientations

This study presents an advanced numerical model for predicting a two-dimensional coupled galloping and vortex-induced vibration (VIV) in cross-flow and in-line directions of square cylinders under symmetric and asymmetric flow orientations. The present model combines the quasi-steady theory for the galloping with the nonlinear structure-wake oscillators simulating VIV, capturing the time-varying drag and lift hydrodynamic forces with the time-averaged and fluctuating components. By placing a flexibly mounted square cylinder in uniform flow at an initial angle of incidence, the cylinder is subject to instantaneous changes in the dynamic angle of attack accounting for relative flow-structure velocities. Modelling of such features in cross-flow and in-line directions for low and high mass ratio systems extends previous studies which have mostly focused on cross-flow responses of square cylinders with high mass ratios at a zero angle of incidence. New sets of empirical coefficients governing the drag and lift fluid forces for both the quasi-steady and wake oscillator approaches are introduced by calibrating with available experimental data in the literature, applicable to predict several flow-induced vibration phenomena under arbitrary flow-structure orientations. Mathematical criteria for the onset of two- and one-dimensional galloping instability are presented, verifying the likelihood of galloping occurrence. Parametric investigations are carried out to highlight the important effects of flow incidence angle, mass-damping ratio (Scruton number) and in-line response on the prediction of galloping and VIV in comparison with experimental results. By varying the reduced velocity parameter, the present model captures key qualitative features of the dominant galloping, interfering galloping-VIV and dominant VIV through the response amplitudes, mean drift displacements, oscillation frequencies, fluid force components and motion trajectories. Contributions from in-line responses are found to be meaningful for the interfering galloping-VIV system with a low mass-damping ratio and for an asymmetric flow orientation. The present model could be further calibrated and applied to other fluid-structure interaction applications with non-circular cross-sectional geometries under omnidirectional flow directions.

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

Initialization Algorithms of a Rigid–Flexible Coupling Fluid–Structure Interaction System by Considering Length-and-Area Constraints

Rigid–flexible coupling fluid–structure interaction systems are expected to be future solutions for reducing energy lost in water. The dynamics of these systems is usually investigated via numerical simulations. However, in existing numerical works there is no accurate algorithm for the initialization of the flexible filament, which ensures both the length and area constraints, leading to inaccurate results or even severe numerical instabilities. We propose two alternative initialization algorithms, respectively, the “Trapezoidal arrangement” and the “Quartic curve arrangement”. The performances of both of these two algorithms are investigated in numerical simulations by using the immersed boundary method. The motion responses and force characteristics of the flexible filament are analyzed carefully, verifying the capability of the proposed algorithms. Specifically, “Quartic curve arrangement” is further recommended due to its good property of convergence.

Outer hair cells stir cochlear fluids.

We hypothesized that active outer hair cells drive cochlear fluid circulation. The hypothesis was tested by delivering the neurotoxin, kainic acid, to the intact round window of young gerbil cochleae while monitoring auditory responses in the cochlear nucleus. Sounds presented at a modest level significantly expedited kainic acid delivery. When outer-hair-cell motility was suppressed by salicylate, the facilitation effect was compromised. A low-frequency tone was more effective than broadband noise, especially for drug delivery to apical locations. Computational model simulations provided the physical basis for our observation, which incorporated solute diffusion, fluid advection, fluid-structure interaction, and outer-hair-cell motility. Active outer hair cells deformed the organ of Corti like a peristaltic tube to generate apically streaming flows along the tunnel of Corti and basally streaming flows along the scala tympani. Our measurements and simulations coherently suggest that active outer hair cells in the tail region of cochlear traveling waves drive cochlear fluid circulation.

Comparison and identification of human coronary plaques with/without erosion using patient-specific optical coherence tomography-based fluid-structure interaction models: a pilot study.

Plaque erosion (PE) with secondary thrombosis is one of the key mechanisms of acute coronary syndrome (ACS) which often leads to drastic cardiovascular events. Identification and prediction of PE are of fundamental significance for disease diagnosis, prevention and treatment. In vivo optical coherence tomography (OCT) data of eight eroded plaques and eight non-eroded plaques were acquired to construct three-dimensional fluid-structure interaction models and obtain plaque biomechanical conditions for investigation. Plaque stenosis severity, plaque burden, plaque wall stress (PWS) and strain (PWSn), flow shear stress (FSS), and ΔFSS (FSS variation in time) were extracted for comparison and prediction. A logistic regression model was used to predict plaque erosion. Our results indicated that the combination of mean PWS and mean ΔFSS gave best prediction (AUC = 0.866, 90% confidence interval (0.717, 1.0)). The best single predictor was max ΔFSS (AUC = 0.819, 90% confidence interval (0.624, 1.0)). The average of maximum FSS values from eroded plaques was 76% higher than that from the non-eroded plaques (127.96 vs. 72.69 dyn/cm2) while the average of mean FSS from erosion sites of the eight eroded plaques was 48.6% higher than that from sites without erosion (71.52 vs. 48.11 dyn/cm2). The average of mean PWS from plaques with erosion was 22.83% lower than that for plaques without erosion (83.2kPa vs. 107.8kPa). This pilot study suggested that combining plaque stress, strain and flow shear stress could help better identify patients with potential plaque erosion, enabling possible early intervention therapy. Further studies are needed to validate our findings.

A probabilistic coral rubble mechanical instability model applied with field observations from the Great Barrier reef

Unstable coral rubble hinders coral recruitment and recovery of coral reefs after damage from cyclones and bleaching events. If coral rubble remains unstable under typical everyday environmental conditions, areas of coral rubble will not be able to recover. Evaluating the probability of rubble instability over regional scale reef systems can assist the optimization of coral reef restoration efforts. Currently, robust and verified models for such applications do not exist. This paper presents a comprehensive assessment method to predict the probability of coral rubble instability, which combines a fluid-structural interaction approach with a statistical regional wave climate model. The hydrodynamic model employs non-linear wave theory to determine near-bed velocity, pressure gradients, and the corresponding drag and inertia forces acting on the coral rubble. The instability model assesses when overturning or sliding forces exceed resisting forces, considering thousands of combinations of different coral sizes and densities to calculate the proportion of instability under a given wave forcing. The model was calibrated and validated using prior laboratory experiments as reported by Kenyon et al. (2023b). The hydrodynamic and instability models use an extensive dataset of non-cyclonic wave climates (hindcast from over 30 years of wind measurements) specific to the region around Heron Reef, Great Barrier Reef, Australia, enabling a comprehensive evaluation of the probability of rubble instability in this area. Results indicate that the overall probability of rubble instability (Pr3) reaches 0.74 in water depths less than 2 m (typical of reef crests or reef flats), while it declines to below 0.21 at a depth of 12 m (typical deeper parts of the fore reef). Coral rubble on reef crests near Heron Reef, which are sheltered by surrounding formations, demonstrates low probability of instability. Thus, coral rubble instability is influenced by both its specific location within the reef and the position of the reef relative to other nearby reefs. By integrating the rubble instability model with non-cyclonic wave climate data, a map of the probability of rubble instability was generated for eight reefs in the Capricorn and Bunker Group (CBG). This map provides valuable guidance for coral reef restoration efforts, significantly reducing the need for extensive field-based data.

Three-Dimensional Flutter Numerical Simulation of Wings in Heavy Gas and Transonic Flutter Similarity Law Correction Method

Wind tunnel testing is a crucial method for studying aircraft flutter. Using heavy gas as the wind tunnel medium can mitigate the escalating issue of test models being overweight as advanced aircraft develop. This paper employs an analytical method for numerical calculations of three-dimensional (3D) wing flutter based on fluid–structure interaction (FSI). Flutter calculations for the Goland wing are conducted, and the results in the air medium are consistent with the literature. In contrast, significant differences in flutter behavior are observed in the heavy gas R134a medium. Compared to air, when the model reaches a critical state in R134a, the incoming flow velocity is lower, the incoming flow density is approximately 3 to 5 times air, and the incoming flow dynamic pressure is about 1.1 to 1.2 times that of air. The correction of heavy gas flutter data is crucial for wind tunnel testing. This paper proposes a correction method based on the unsteady transonic flow similarity law proposed by Bendiksen under quasi-steady conditions. Attempts are made to revise relevant published wind tunnel tests and heavy gas flutter calculation results. The transonic flutter similarity law effectively explains the flutter similarity of rigid models in both heavy gas and air media. Still, it fails in cases with highly reduced frequencies and low mass ratios, such as those encountered with flexible wings.

Impact of Geometric Attributes on Abdominal Aortic Aneurysm Rupture Risk: An InVivo FSI-Based Study.

Reported in this paper is a cutting-edge computational investigation into the influence of geometric characteristics on abdominal aortic aneurysm (AAA) rupture risk, beyond the traditional measure of maximum aneurysm diameter. A Comprehensive fluid-structure interaction (FSI) analysis was employed to assess risk factors in a range of patient scenarios, with the use of three-dimensional (3D) AAA models reconstructed from patient-specific aortic data and finite element method. Wall shear stress (WSS), and its derivatives such as time-averaged WSS (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) and transverse WSS (transWSS) offer insights into the force dynamics acting on the AAA wall. Emphasis is placed on these WSS-based metrics and seven key geometric indices. By correlating these geometric discrepancies with biomechanical phenomena, this study highlights the novel and profound impact of geometry on risk prediction. This study demonstrates the necessity of a multidimensional assessment approach, future efforts should complement these findings with experimental validations for an applicable approach for clinical use.

Performance analysis and improvement of a novel reed valve in compressors using fluid–structure interaction method

The reed valves used as suction valves in reciprocating compressors are more susceptible to damage due to the limitation of the lift limiter which can't completely block the entire valve plate. This article proposed a novel reed valve featuring a redesigned shape aimed to solve this issue. Compared to conventional reed valves, the proposed valve could effectively mitigate the highest stress occurring at the root of the suction valve. To thoroughly investigate the characteristics of this valve type and improve its performance, a three-dimensional fluid–structure interaction (FSI) model was established and validated using experiments. The leakage through the valve gap in the FSI model was solved using a combination of user-defined function (UDF) and Scheme file in Fluent. Based on this FSI model, the dynamic characteristics and stress distribution of the valve were mainly analyzed. The effect of the lift limiter on the performance of the valve was analyzed, and a suitable lift was determined. Moreover, the influence of thickness and shape parameters, was also analyzed and an efficient improved scheme was subsequently proposed. Through this scheme, the maximum stress of the valve was reduced by 23.77% compared with that of the original valves. This article provides important information for the design and optimization of complex-shaped reed valves.

Optimization Simulation of Mooring System of Floating Offshore Wind Turbine Platform Based on SPH Method

ABSTRACTIn this paper, the Smoothed Particle Hydrodynamics method is used to simulate the dynamic response of the floating offshore wind turbine mooring system in a complex marine environment by using its advantages of dealing with free surface flow and fluid–structure interaction. The single‐cable mooring system and the double‐cable mooring system are introduced into the numerical simulation model of fluid–solid coupling interaction. The first‐order irregular wave and the second‐order regular wave are used to simulate different wave conditions. The operating state of the platform in these two cases is analyzed, and the accurate data such as the velocity change of the geometric center of the platform and the change of six degrees of freedom are obtained. Using visual processing and data analysis, it is found that the optimized double‐cable mooring system structure improves the vibration reduction ability of the platform.

Investigating the role of thrombosis and false lumen orbital orientation in the hemodynamics of Type B aortic dissection

While much about the fundamental mechanisms behind the initiation and progression of Type B aortic dissection (TBAD) is still unknown, predictive models based on patient-specific fluid-structure interaction (FSI) simulations can help in risk stratification and optimal clinical decision-making. Aiming at the development of personalized treatment, FSI models can be leveraged to investigate the interplay between complex aortic flow patterns and anatomical features, while considering the deformation of the arterial wall and the dissection flap. In this study, the hemodynamics of false lumen thrombosis, a large fenestration, and the orbital orientation of the false lumen is studied through image-based FSI simulations on three TBAD patient-specific geometries. A new pipeline is developed leveraging the open-source software SimVascular and ParaView to analyze multiple patients simultaneously and to achieve large-scale parallelization in FSI results based on patients’ computed tomography (CT) images. The results of this study suggest that the internal orbital orientation of the false lumen contributes to maintaining a positive luminal pressure difference ΔPTL-FL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Delta P_{TL-FL}$$\\end{document} = PTL-PFL\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$P_{TL}-P_{FL}$$\\end{document} between the true lumen (TL) and the false lumen (FL), despite an impingement area in the false lumen near the entry tear. A positive and high luminal pressure difference is thought to promote TL expansion and FL compression. Moreover, it was also found that FL thrombosis at the entry tear region reduce the magnitude of the negative luminal pressure difference, which in turn may reduce FL expansion and the risk of unstable aortic growth. Finally, this FSI study suggests that the aortic wall and dissection flap stiffness determines the effects of a large fenestration in the descending thoracic aorta on the luminal pressure difference.

Theoretical model construction and static performance analysis of multi-leaf journal foil bearings

Abstract Multi-leaf journal foil bearing (MLJFB) structures are complex, making performance prediction challenging, and the accurate prediction of gas film thickness and pressure distributions is crucial. In this study, COMSOL finite-element software was used to establish an initial foil assembly and structural deformation simulation model for an MLJFB with bump foils. Using MATLAB, a compressible gas lubrication model was established and solved using the finite volume method (FVM). Real-time data transfer between COMSOL and MATLAB enabled a fully coupled fluid-structure interaction simulation. This model meticulously accounts for the pre-tension force in the multi-leaf bearing assembly, contact friction between adjacent foils, and the interaction between the top foil and gas film pressure. The effects of journal movement distance, direction, and top foil overlapping ratio on the bearing capacity as well as the effects of different external load directions on the journal equilibrium position were analyzed based on this model. The results indicate that increasing the top foil overlapping ratio enhances the bearing capacity of the MLJFB. Furthermore, compared to the direction of the top foil overlap position, external loads applied in the direction of the top foil midpoint are more beneficial for bearing operational stability.

On the Local Existence of Solutions to the Fluid–Structure Interaction Problem with a Free Interface

We address a system of equations modeling an incompressible fluid interacting with an elastic body. We prove the local existence when the initial velocity belongs to the space H1.5+ϵ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^{1.5+\\epsilon }$$\\end{document} and the initial structure velocity is in H1+ϵ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H^{1+\\epsilon }$$\\end{document}, where ϵ∈(0,1/20)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\epsilon \\in (0, 1/20)$$\\end{document}.

Computational fluid dynamics for naval hydrodynamics

This article describes key issues which have to be addressed to apply Computational Fluid Dynamics to Naval Hydrodynamics. The specific aspects of Naval Hydrodynamics are discussed and illustrated by recent simulations and comparisons with available experiments. Free-surface flows with or without waves and even violent phenomena such as ventilation or cavitation can be modelled with mixture-fluid surface capturing. Turbulence modelling of thick boundary layers and vortical flows requires anisotropic RANS models or hybrid RANS/LES in case of strongly separated flows. Moreover, fluid–structure interaction in the form of rigid or flexible body motion and multi-body systems is crucial to represent ship manoeuvring and propulsion. Finally, the paper underlines the central role played by anisotropic adaptive grid refinement in the accurate simulation of marine flows.

Computational aerodynamics of insect flight using volume penalization

The state-of-the-art of insect flight research using advanced computational fluid dynamics techniques on supercomputers is reviewed, focusing mostly on the work of the present authors. We present a brief historical overview, discuss numerical challenges and introduce the governing model equations. Two open source codes, one based on Fourier, the other based on wavelet representation, are succinctly presented and a mass-spring flexible wing model is described. Various illustrations of numerical simulations of flapping insects at low, intermediate and high Reynolds numbers are presented. The role of flexible wings, data-driven modeling and fluid–structure interaction issues are likewise discussed.

Spatiotemporal mode extraction for fluid–structure interaction using mode decomposition

A fluid–structure interaction (FSI) analysis model was developed based on the partitioned coupling approach. The study analyzed the fluid and structure behavior in a well-established validation case known as the Turek–Hron FSI benchmark test. This test involves strong, unsteady interaction between fluid flow and an elastic flap behind a rigid cylinder. Comparison of the elastic flap displacement frequencies and amplitudes from the present FSI model with reference data showed good agreement, validating the model. Dynamic mode decomposition (DMD) was employed to extract fundamental structures from the complex spatiotemporal data obtained from the FSI analysis. In addition, a mode–sensing technique based on a greedy algorithm was developed, and significant modes were extracted with elastic flap displacement as a feature. The crucial modes for the elastic flap displacement were the second bending modes at specific frequencies. However, these results differed significantly from the natural frequencies of the second bending modes obtained via eigenmode analysis. This discrepancy can be attributed to the close coupling between the fluid and structure, which alters elastic deformation behavior. The study demonstrates the potential for straightforward extraction of essential fluid–structure coupling using FSI-DMD.

Analysis of the influence of bottom flow holes in control rod guide tubes on flow field and control rod displacement

During the operation of pressurized water reactors, flow-induced vibration issues are commonly encountered, where control rod components vibrate due to the impact of coolant flow, leading to wear, thereby posing a threat to the safe operation of the reactor. This paper conducts a detailed numerical simulation of the flow field and force distribution inside a full-scale lower control rod guide tube based on the CFD method. The fluid–structure interaction analysis is performed to reveal the influence of cross-flow under various operating conditions on the vibration and displacement of control rods. The research subject consists of 24 control rods, one continuous guide section, and six control rod guide cards. The lower part of the guide tube has two flow holes on each surface, with one assumed to be a cross-flow inlet and the others as outlets. The results indicate that when considering the cross-flow effect at the flow holes, 90 % of the coolant flows out from other flow holes at the bottom of the guide tube, resulting in increased transverse velocity at the bottom and the presence of numerous vortices. The magnitude and location of the force exerted on the control rod by coolant impact are correlated. The control rod closest to the inlet of the flow hole experiences the greatest force, approximately 5.35 times the average at other positions. At a cross-flow velocity of 0.4 m/s and a low frequency of 10 Hz, the maximum displacement of the control rod is 0.3355 mm.

A High-Fidelity Computational Model for Predicting Blood Cell Trafficking and 3D Capillary Hemodynamics in Retinal Microvascular Networks.

To present a first principle-based, high-fidelity computational model for predicting full three-dimensional (3D) and time-resolved retinal microvascular hemodynamics taking into consideration the flow and deformation of individual blood cells. The computational model is a 3D fluid-structure interaction model based on combined finite volume/finite element/immersed-boundary methods. Three in silico microvascular networks are built from high-resolution in vivo motion contrast images of the superficial capillary plexus in the parafoveal region of the human retina. The maximum tissue area represented in the model is approximately 500 × 500 µm2, and vessel lumen diameters ranged from 5.5 to 25 µm covering capillaries, arterioles, and venules. Blood is modeled as a suspension of individual blood cells, namely, erythrocytes (RBC), leukocytes (WBC), and platelets in plasma. An accurate and detailed biophysical modeling of each blood cell and their flow-induced deformation is considered. A physiological, pulsatile boundary condition corresponding to an average cardiac cycle of 0.9 second is used. Detailed quantitative data and analysis of 3D retinal microvascular hemodynamics are presented, and their relationship to RBC flow dynamics is illustrated. Blood velocity is shown to have temporal oscillations superimposed on the background pulsatile variation, which arise because of the way RBCs partition at vascular junctions, causing repeated clogging and unclogging of vessels. Temporal variations in RBC velocity and hematocrit are anti-correlated in a given vessel, but their time-averaged distributions are positively correlated across the network. Whole blood velocity is 65% to 85% of RBC velocity, with the discrepancy related to the formation of an RBC-free region, adjacent to the vascular endothelium and typically 0.8 to 1.8 µm thick. The 3D velocity and RBC concentration profiles are shown to be oppositely skewed with respect to each other, because of the way that RBCs "hug" the apex of each bifurcation. RBC deformation is predicted to have biphasic behavior with respect to vessel diameter, with minimal cell length for vessels approximately 7 µm in diameter. The wall shear stress (WSS) exhibits a strongly 3D distribution with local regions of high value and gradient spanning a range of 10 to 80 dyn/cm2. WSS is highest where there is faster flow, greater curvature of the vessel wall, capillary bifurcations, and at locations of RBC crowding and associated thinning of the cell-free layer. This study highlights the usefulness of high-fidelity cell-resolved modeling to obtain accurate and detailed 3D, time-resolved retinal hemodynamic parameters that are not readily available through noninvasive imaging approaches. The results presented are expected to complement and enhance the interpretation of in vivo data, as well as open new avenues to study retinal hemodynamics in health and disease.

A coupled peridynamics–smoothed particle hydrodynamics model for fluid–structure interaction with large deformation

Fluid–structure interaction (FSI) is ubiquitous in various engineering disciplines, and effectively managing FSI often appears to be the key for successful failure analysis and safety-oriented design. Smoothed particle hydrodynamics (SPH) serves as a potent nonlocal meshfree method for fluid dynamics modeling, while peridynamics (PD) demonstrates exceptional capability in addressing structural dynamics involving large deformations and discontinuities. Thus, leveraging their respective strengths in a combined approach holds significant promise for tackling FSI challenges. In this work, we propose a new peridynamics–smoothed particle hydrodynamics (PD-SPH) coupling model for addressing FSI. A stable and efficient coupling algorithm for data transfer between PD and SPH is put forward. In this coupling strategy, a PD particle directly participates in solving the SPH governing equations when it is identified to be within the support domain of an SPH particle. This can be done since the SPH quantities including the density, velocity, and pressure of a PD particle are naturally attainable within the framework of non-ordinary state-based peridynamics theory. Concurrently, in solving PD governing equations, reaction forces from SPH particles act as external forces for PD particles, determined straightforwardly through Newton's third law. As such, the proposed PD-SPH coupling strategy is straightforward to implement and offers high computational efficiency. Validation examples demonstrate that the proposed PD-SPH coupling model is computationally robust and adept at capturing physical phenomena in diverse FSI scenarios involving breaking free surfaces of fluid and large structural deformations of solid. Moreover, the proposed PD-SPH coupling model is flexible introducing no constraint conditions for applications and can accommodate different particle resolutions for PD and SPH domains. These features enable a broad application range of the proposed PD-SPH coupling model including simulations of explosion-induced soil fragmentation, rock fracture, and concrete dam failure, which will be conducted by authors in the near future.

Accurate simulations of moving flexible objects with an improved immersed boundary-lattice Boltzmann method

Accurately capturing boundaries is crucial for simulating fluid–structure interaction (FSI) problems involving flexible objects undergoing large deformations. This paper presents a coupling of the immersed boundary-lattice Boltzmann method with a node-based partly smoothed point interpolation method (NPS-PIM) to enhance the accuracy of simulating moving flexible bodies in FSI problems. The proposed method integrates a multiple relaxation time scheme and employs a force correction technique to address boundary capturing inaccuracies. The effect of virtual fluid is accounted for through a Lagrangian point approximation, ensuring precise FSI force calculations for unsteady solid motions. NPS-PIM is utilized as the solid solver, constructing a moderately softened model stiffness by combining the finite element method (FEM) with the node-based smoothed PIM (NS-PIM). Simulations of flow fields near flexible objects with large deformations demonstrate that the proposed approach reduces numerical errors, improves computational efficiency compared to traditional FSI models using FEM and NS-PIM, and accurately captures the behavior of moving flexible bodies and detailed flow fields.