Fluid-structure Interaction Simulations Research Articles

Towards blade-resolved ice accretion simulations of flexible blades using mixed-fidelity models

We present a modular framework for blade-resolved fluid-structure interaction simulations of rotors. The framework is based on the coupling of the multi-body solver MBDyn with the fluid solver SU2 using the library preCICE. A preliminary validation is carried out against other open-source codes analysing the UAE Phase VI experiment. In the future, the framework will be used to study ice accretion on flexible blades. Indeed, blade flexibility has not been considered yet in the numerical modelling of ice accretion due to lacking computational tools and high computational costs. A first application on ice accretion is shown by analyzing the temporal evolution of the modal frequencies of an isolated blade during an icing event simulated numerically with a quasi-3D approach. The results provide insight into the nature of ice throw from wind turbines.

Fluid-structure interaction simulation of dynamic response and wake of floating offshore wind turbine considering tower shadow effect

Fluid-structure interaction simulation of dynamic response and wake of floating offshore wind turbine considering tower shadow effect

Additional study of failure behavior for ship propeller: The comparative analysis of life evaluation approaches between BP neural network, LSTM and TV-HSMM

The study focuses on investigating the appropriate method for evaluating the residual life of a propeller and analyzing its fracture characteristics. A fluid-structure interaction simulation is conducted to qualitatively analyze the stress distribution of the propeller during operation, identifying regions prone to stress concentration. Subsequently, fracture surface characteristics are examined, revealing quasi-cleavage fracture features such as tongue striations, dimples, and tearing ridges in the initiation and propagation regions. Additionally, composite fracture characteristics and geometric mutations are observed on the edge of the fracture surface. Three residual life evaluation methods are then proposed, utilizing BP neural network, LSTM, and TV-HSMM. Comparative analysis suggests that the TV-HSMM method offers the advantage of minimal error in residual life evaluation. However, the models based on BP neural network and LSTM are preferred for their lower average error and distribution trends that closely resemble actual life data.

On the role of aortic valve architecture for physiological hemodynamics and valve replacement, Part I: Flow configuration and vortex dynamics

Aortic valve replacement has become an increasing concern due to the rising prevalence of aortic stenosis in an ageing population. Existing replacement options have limitations, necessitating the development of improved prosthetic aortic valves. In this study, flow characteristics during systole in a stenotic aortic valve case are compared with those downstream of two newly designed surgical bioprosthetic aortic valves (BioAVs). To do so, advanced three-dimensional fluid–structure interaction simulations are conducted and dedicated analysis methods to investigate jet flow configuration and vortex dynamics are developed. Our findings reveal that the stenotic case maintains a high jet flow eccentricity due to a fixed orifice geometry, resulting in flow separation and increased vortex stretching and tilting in the commissural low-flow regions. One BioAV design introduces non-axisymmetric leaflet motion, which reduces the maximum jet velocity and forms more vortical structures. The other BioAV design produces a fixed symmetric triangular jet shape due to non-moving leaflets and exhibits favourable vorticity attenuation, revealed by negative temporally and spatially averaged projected vortex stretching values, and significantly reduced drag. Therefore, this study highlights the benefits of custom-designed aortic valves in the context of their replacement through comprehensive and novel flow analyses. The results emphasise the importance of analysing jet flow, vortical structures, momentum balance and vorticity transport for thoroughly evaluating aortic valve performance.

On the role of aortic valve architecture for physiological hemodynamics and valve replacement, Part II: Spectral analysis and anisotropy

Severe aortic valve stenosis can lead to heart failure and aortic valve replacement (AVR) is the primary treatment. However, increasing prevalence of aortic stenosis cases reveal limitations in current replacement options, necessitating improved prosthetic aortic valves. We investigate flow disturbances downstream of severe aortic stenosis and two bioprosthetic aortic valve (BioAV) designs using advanced energy-based analyses. Three-dimensional high-fidelity fluid–structure interaction simulations have been conducted and a dedicated and novel spectral analysis has been developed to characterise the kinetic energy (KE) carried by eddies in the wavenumber space. In addition, new field quantities, i.e. modal KE anisotropy intensity as well as normalised helicity intensity, are introduced. Spectral analysis shows kinetic energy (KE) decay variations, with the stenotic case aligning with Kolmogorov’s theory, while BioAV cases differing. We explore the impact of flow helicity on KE transfer and decay in BioAVs. Probability distributions of modal KE anisotropy unveil flow asymmetries in the stenotic and one BioAV cases. Moreover, an inverse correlation between temporally averaged modal KE anisotropy and normalised instantaneous helicity intensity is noted, with the coefficient of determination varying among the valve configurations. Leaflet dynamics analysis highlights a stronger correlation between flow and biomechanical KE anisotropy in one BioAV due to higher leaflet displacement magnitude. These findings emphasise the role of valve architecture in aortic turbulence as well as its importance for BioAV performance and energy-based design enhancement.

An assessment of CFD-scale fluid–structure interaction simulations through comprehensive experimental data in cross-flow

Experimental tests accessing both fluid and structure behaviors are mandatory for a consistent and comprehensive assessment of fluid–structure interaction (FSI) numerical simulations. In this paper, recently published results of an experimental configuration of two in-line cantilever cylinders subjected to water cross-flow were considered: instantaneous fluid velocities measurements are available on several positions inside the test section, together with the cylinder vibrations. FSI simulations were performed by coupling Ansys Fluent (for the fluid domain) with Ansys Mechanical (for the solid domain). URANS simulations and Scale-Adaptive Simulations (SAS) were employed as CFD simulation approach. The structure displacements were taken into account through an Arbitrary Lagrangian–Eulerian approach. The fluid–structure coupling was 2-way explicit. Simulations were performed for two different water mass flow rates. For the highest one, vortex-induced resonance was observed experimentally. The numerical results show consistent agreement in terms of shedding frequency and velocity spectra behind the cylinders. The calculated vibration response is overall consistent, despite some underestimations, for the cylinders not featuring vortex-induced resonance; nevertheless, the experimentally observed vortex-induced resonance could not be reproduced by the numerical simulations.

Accelerating the force-coupling method for hydrodynamic interactions in periodic domains

The efficient simulation of fluid-structure interactions at zero Reynolds number requires the use of fast summation techniques in order to rapidly compute the long-ranged hydrodynamic interactions between the structures. One approach for periodic domains involves utilising a compact or exponentially decaying kernel function to spread the force on the structure to a regular grid where the resulting flow and interactions can be computed efficiently using an FFT-based solver. A limitation to this approach is that the grid spacing must be chosen to resolve the kernel and thus, these methods can become inefficient when the separation between the structures is large compared to the kernel width. In this paper, we address this issue for the force-coupling method (FCM) by introducing a modified kernel that can be resolved on a much coarser grid, and subsequently correcting the resulting interactions in a pairwise fashion. The modified kernel is constructed to ensure rapid convergence to the exact hydrodynamic interactions and a positive-splitting of the associated mobility matrix. We provide a detailed computational study of the methodology and establish the optimal choice of the modified kernel width, which we show plays a similar role to the splitting parameter in Ewald summation. Finally, we perform example simulations of rod sedimentation and active filament coordination to demonstrate the performance of fast FCM in application.

One-Way Fluid–Structure Interaction Simulation of the Added Damping of T(0,1) Modes of Two Partially Overlapped Identical Plates in Still Water

Abstract Structures with a partially overlapped status in water can be seen in some engineering applications, and the fluid-structure coupling vibration behavior of two partially overlapped identical plates has been studied previously through the experimental method. In this study, the added damping of the T(0,1) modes of two submerged partially overlapped plates is numerically investigated through the one-way fluid-structure interaction (FSI) method. The relevant numerical settings, like the mesh size, calculated periods, time-step size, and vibration amplitude, were tested first. Then, the numerical results were compared with experimental results, and good agreements were found. Finally, numerical results were analyzed. The vibration status of two plates in the joint abutting area or the overlapped area has an important influence on the added mass variation. When the added mass is higher, the phase difference between modal force and vibration displacement is also greater, which is the main reason for the higher added damping. The relationship between the phase difference and the frequency in water can be approximately fitted to a straight line, which can probably be used to predict the added damping variations caused by fluid boundary changes of submerged structures.

Investigation of oblique water entry of high-speed supercavitating projectiles using transient fluid-structure interaction simulation

This study employs a fluid-structure interaction simulation to examine the structural deformation of supercavitating projectiles during high-speed oblique water entry. A cavitation flow field simulation model is established using the VOF multiphase and cavitation models, while the structural dynamics model of the projectile is implemented using the Johnson-Cook constitutive equation. The transient fluid-structure interaction simulation is conducted through data transfer between the two models. Numerical simulations are then conducted with the assumptions of a rigid body, linear elasticity, and nonlinearity, respectively. The results indicate that the high-speed water entry induced substantial hydrodynamic force and plastic strain accumulation in the projectile, leading to considerable bending deformation. It was found that the deformation significantly influenced the hydrodynamic forces, highlighting the importance of fluid-structure interaction modeling. Moreover, the attack angle directly impacts the resulting deformation mode. Specifically, under high-speed conditions, three distinct deformation modes are observed at different attack angles: minimal deformation at low angles, an L-shaped bending at higher angles, and an S-shaped bending when the angle exceeds the threshold. This study highlights the importance of fluid-structure interaction methods in analyzing the deformation behavior of supercavitating projectiles during high-speed oblique water entry. The results show that the bidirectional dynamic interaction between fluid and structure significantly influences the structural stability of the supercavitating projectiles. The relevant conclusions can provide a reference for the optimal design of supercavitating projectiles to ensure the structural stability under high-speed conditions.

On the role of tissue mechanics in fluid–structure interaction simulations of patient‐specific aortic dissection

AbstractModeling an aortic dissection represents a particular challenge from a numerical perspective, especially when it comes to the interaction between solid (aortic wall) and liquid (blood flow). The complexity of patient‐specific simulations requires a variety of parameters, modeling assumptions and simplifications that currently hinder their routine use in clinical settings. We present a numerical framework that captures, among other things, the layer‐specific anisotropic properties of the aortic wall, the non‐Newtonian behavior of blood, patient‐specific geometry, and patient‐specific flow conditions. We compare hemodynamic indicators and stress measurements in simulations with increasingly complex material models for the vessel tissue ranging from rigid walls to anisotropic hyperelastic materials. We find that for the present geometry and boundary conditions, rigid wall simulations produce different results than fluid–structure interaction simulations. Considering anisotropic fiber contributions in the tissue model, stress measurements in the aortic wall differ, but shear stress‐based biomarkers are less affected. In summary, the increasing complexity of the tissue model enables capturing more details. However, an extensive parameter set is also required. Since the simulation results depend on these modeling choices, variations can lead to different recommendations in clinical applications.

Modified Overset Method in OpenFOAM for Simultaneous Motion and Deformation: A Case Study of a Flexible Flapping Foil

The influence of wing deformation in animal propulsion and motion is an intriguing phenomenon that has contributed to a growing interest in biomimetics in both academia and industry. The study focuses on developing a robust numerical tool for analyzing the motion of biological systems in fluid flows, crucial to understanding fluid-structure interaction (FSI) phenomena. The research introduces a new mesh deformation solver, dynamicOversetZoneFvMesh, designed to address limitations in OpenFOAM's conventional Overset method. This solver enables the application of both rigid body motion and deformation at specified patches within certain zones, allowing for more complex motion scenarios. Through meticulous verification, parallelization, validation, and mesh convergence studies, the modified Overset solver's reliability is confirmed. Two cases are examined: a rigid-body flapping foil and a flexible foil with leading-edge deformation. Results demonstrate the adapted solver's proficiency in managing complex fluid dynamics and accommodating simultaneous motion and deformation. Furthermore, the study reveals that leading-edge flexibility improves the power extraction efficiency of flapping foils, aligning with previous literature. Overall, the research lays the foundation for advancing FSI simulations and opens avenues for addressing previously challenging problems in the field.

Fast fluid–structure interaction simulation method based on deep learning flow field modeling

The rapid acquisition of high-fidelity flow field information is of great significance for engineering applications such as multi-field coupling. Current research in flow field modeling predominantly focuses on low Reynolds numbers and periodic flows, exhibiting weak generalization capabilities and notable issues with temporal inferring error accumulation. Therefore, we establish a reduced order model (ROM) based on Convolutional Auto-Encoder (CAE) and Long Short-Term Memory (LSTM) neural network and propose an unsteady flow field modeling method for the airfoil with a high Reynolds number and strong nonlinear characteristics. The attention mechanism and weak physical constraints are integrated into the model architecture to improve the modeling accuracy. A broadband excitation training strategy is proposed to overcome the error accumulation problem of long-term inferring. With only a small amount of latent codes, the relative error of the flow field reconstructed by CAE can be less than 5‰. By training LSTM with broadband excitation signals, stable dynamic evolution can be achieved in the time dimension. CAE-LSTM can accurately predict the forced response and complex limit cycle behavior of the airfoil in a wide range of amplitude and frequency under subsonic/transonic conditions. The relative errors of predicted variables and integral force are less than 1%. The fluid–structure interaction framework is built by coupling the ROM and motion equations of the structure. CAE-LSTM predicts the time series response of pitch displacement and moment coefficient at different reduced frequencies, which is in good agreement with computational fluid dynamics, and the simulation time savings exceed one order of magnitude.

Investigation of shear stress behavior on membrane surfaces in reciprocating membrane bioreactors using fluid-structure interaction simulations

Although reciprocating membrane bioreactors (rMBRs) have been developed as energy-efficient wastewater treatment technologies, fouling removal mechanisms are yet to be elucidated. This study conducted simulations of fluid-structure interaction and experiments using particle image velocimetry to investigate the behavior of shear stress on the membrane surface in rMBRs. Simulations were performed on the three-dimensional model, varying parameters such as average reciprocating velocity (ARV), moving distance (MD), and membrane slack ratio (MSR). The results found that fluid advection can hinder foulant adhesion to the membranes. Increasing the ARV from 4 to 12 cm/s under shear stress resulted in higher fluid resistance and inertial forces, leading to elevated shear stress. The effect of MD on shear stress was negligible within the range of 4–12 cm. An intermediate level of MSR from 0 to 3% promoted turbulence around the membrane, thereby increasing shear stress. The response surface suggested that it may be desirable to simultaneously increase the ARV and MSR to maximize shear stress. However, this increase can also lead to stability failure by increasing the principal stress in the membrane potting area. This study provides insights into fouling removal mechanisms to improve the efficiency of rMBRs.

An efficient partitioned framework to couple Arbitrary Lagrangian-Eulerian and meshless vector form intrinsic finite element methods for fluid-structure interaction problems with deformable structures

The Vector Form Intrinsic Finite Element (VFIFE) is an advanced meshless Lagrangian structural solver. This study introduces a partitioned framework to integrate Arbitrary Lagrangian-Eulerian and VFIFE methods into fluid-structure interaction (FSI) problems featured by large structural displacements and deformations. The VFIFE method enables the FSI simulations with deformable structures due to the nature of Lagrangian particles, which is challenging for traditional computational fluid mechanics techniques. In addition, the VFIFE can be effectively equipped with the partitioned coupling schemes because of the versatility of governing equations. In this algorithm system, several matrix problems in finite element methods and mapping and interpolation in the fluid-structure interface can be avoided. The significant potential of this algorithm in tackling complex FSI problems is demonstrated through various benchmarks, encompassing scenarios like Taylor-Green vortex, flow over a rigid cylinder (Re = 40 ∼ 200), forced vibration of a rigid cylinder, vortex-induced vibration (VIV) of an elastically mounted cylinder, and VIV of a rigid cylinder with an attached flexible cantilever beam. Finally, an experiment on flow past a curtain in a water tank is conducted to further evaluate the precision and stability of the current coupling strategy in three dimensions, demonstrating the feasibility of such methods in flexible structures.

Effect of pulsating flow on flow-induced vibrations of circular and square cylinders in the laminar regime

Through fluid-structure interaction simulations, this study assesses the dynamic response characteristics of elastically mounted circular and square cylinders subjected to pulsating inflow conditions, providing valuable insights into the analysis and optimization of these systems. The main focus of the present work is on analyzing the effects of two factors: (i) the ratio of the oscillatory velocity component to the steady velocity component in pulsating flow (flow ratio) and (ii) the ratio of the oscillation frequency of pulsating flow to the natural frequency of the structure (frequency ratio). The simulation results for different parameters of interest are analysed using Fourier analysis and Poincaré maps of time series data, and contour plots of vorticity. For the circular cylinder, it is found that cylinder loses synchronization in lock-in as the flow and frequency ratios are increased. Three distinct vibration patterns of vortex-induced vibration are observed for selected combinations of flow and frequency ratios at a Reynolds number of 110 for circular cylinder. For the galloping of square cylinder at a Reynolds number of 250, it is found that the instability and nonlinearity of vortex shedding become more pronounced as the flow ratio increases.

A phase field-immersed boundary-lattice Boltzmann coupling method for fluid-structure interaction analysis

Fluid-structure interaction (FSI) is critical in numerous engineering and scientific applications. This paper presents a numerical model for FSI simulations that integrates the lattice Boltzmann method (LBM) for fluid dynamics, the phase field method (PFM) for structural deformation and fracture, and the immersed boundary method (IBM) for fluid-solid coupling. The fluid employs an explicit time integration scheme with a smaller time step than the solid, which uses an implicit time integration scheme, satisfying the CFL condition. During the synchronization of fluid-solid interaction information, a multi-substep strategy is applied. The IBM is implemented using the iterative velocity method for boundary treatment. This approach facilitates the information transfer between domains with different time steps, enhancing the coupling accuracy through iterative computation. The effectiveness of the improved immersed boundary-lattice Boltzmann (IB-LBM) has been verified by the classical cylindrical flow case, which effectively prevents the particle penetration phenomenon and improves the calculation accuracy. The phase field-immersed boundary-lattice Boltzmann method (PF-IB-LBM) is subsequently deployed to simulate the deformation of an aluminum plate under hydrostatic pressure and the aftermath of a dam-break impact on an elastic plate. Verification of the coupling method substantiates its feasibility and efficacy. In addition, the proposed PF-IB-LBM is implemented to simulate the entire process of structural deformation and failure in response to fluid dynamics, which further confirms its capability to predict discontinuous mechanical behaviors such as damage and crack propagation under intricate fluid flows. The results emphasize the flexibility and efficiency of the proposed approach, providing a robust framework for solving more complex scenarios in engineering FSI problems.

Study on the optimal elastic modulus of flexible blades for right heart assist device supporting patients with single-ventricle physiologies.

Patients with single-ventricle physiologies continue to experience insufficient circulatory power after undergoing palliative surgeries. This paper proposed a right heart assist device equipped with flexible blades to provide circulatory assistance for these patients. The optimal elastic modulus of the flexible blades was investigated through numerical simulation. A one-way fluid-structure interaction (FSI) simulation was employed to study the deformation of flexible blades during rotation and its impact on device performance. The process began with a computational fluid dynamics (CFD) simulation to calculate the blood pressure rise and the pressure on the blades' surface. Subsequently, these pressure data were exported for finite element analysis (FEA) to compute the deformation of the blades. The fluid domain was then recreated based on the deformed blades' shape. Iterative CFD and FEA simulations were performed until both the blood pressure rise and the blades' shape stabilized. The blood pressure rise, hemolysis risk, and thrombosis risk corresponding to blades with different elastic moduli were exhaustively evaluated to determine the optimal elastic modulus. Except for the case at 8,000 rpm with a blade elastic modulus of 40 MPa, the pressure rise associated with flexible blades within the studied range (rotational speeds of 4,000 rpm and 8,000 rpm, elastic modulus between 10 MPa and 200 MPa) was lower than that of rigid blades. It was observed that the pressure rise corresponding to flexible blades increased as the elastic modulus increased. Additionally, no significant difference was found in the hemolysis risk and thrombus risk between flexible blades of various elastic moduli and rigid blades. Except for one specific case, deformation of the flexible blades within the studied range led to a decrease in the impeller's functionality. Notably, rotational speed had a more significant impact on hemolysis risk and thrombus risk compared to blade deformation. After a comprehensive analysis of blade compressibility, blood pressure rise, hemolysis risk, and thrombus risk, the optimal elastic modulus for the flexible blades was determined to be between 40 MPa and 50 MPa.

Spatiotemporal analysis of the effects of exercise on the hemodynamics of the aorta in hypertensive rats using fluid-structure interaction simulation

Abstract Background and Objective Hemodynamic changes that lead to increased blood pressure represent the main drivers of organ damage in hypertension. Prolonged increases to blood pressure can lead to vascular remodeling, which also affects vascular hemodynamics during the pathogenesis of hypertension. Exercise is beneficial for relieving hypertension, however the mechanistic link between exercise training and how it influences hemodynamics in the context of hypertension is not well understood. Methods n exercise model was developed using spontaneously hypertensive rats (SHR) subject to a 12-week treadmill training regime. The heart rates and blood pressures of rats were measured using the tail cuff method, while micro-computed tomography (CT) scanning was used to develop three-dimensional structures of rat aorta, and ultrasound was used to detect rat aortic blood flow and changes to vessel wall structures. Computational fluid dynamics (CFD) and fluid-structure interaction (FSI) models were used to simulate and measure hemodynamic parameters of the rat aortic vessels. In parallel, Masson staining was performed on fixed samples of blood vessels to investigate collagen volume fraction. Hypertensive rats in the sedentary and long-term exercise training groups were subjected to a single bout exercise training, and their aortic hemodynamic parameters were analyzed before, 5 min, 24 h, and 72 h after the single bout exercise. Results Of the two models, in comparison to actual ultrasonic measurement values recorded, we found that numerical simulation results from the FSI model could more accurately model blood flow in the ascending aorta of hypertensive rats, compared to the CFD model. Moreover, longterm exercise training improved local hemodynamic parameters of blood vessels, and led to improvements in adverse hemodynamic features documented, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT). Longterm exercise training of SHR also improved local vascular collagen deposition in the aorta, while improvements in vascular remodeling were also correlated with favorable hemodynamic parameters. Compared with sedentary SHR, signals for low TAWSS regions of the aortic arch in SHR on the long-term exercise regime shifted to the position of the ascending aorta after a single bout of exercise. Conclusions This study demonstrates that FSI is informative to study the spatiotemporal effects of long-term exercise training on hemodynamic changes within the aortas of hypertensive rats, and that long-term exercise is beneficial through its effects to modulate vascular hemodynamics in hypertension.

Smoothed particle hydrodynamics based FSI simulation of the native and mechanical heart valves in a patient-specific aortic model

The failure of the aortic heart valve is common, resulting in deterioration of the pumping function of the heart. For the end stage valve failure, bi-leaflet mechanical valve (most popular artificial valve) is implanted. However, due to its non-physiological behaviour, a significant alteration is observed in the normal haemodynamics of the aorta. While in-vivo experimentation of a human heart valve (native and artificial) is a formidable task, in-silico study using computational fluid dynamics (CFD) with fluid structure interaction (FSI) is an effective and economic tool for investigating the haemodynamics of natural and artificial heart valves. In the present work, a haemodynamic model of a natural and mechanical heart valve has been developed using meshless particle-based smoothed particle hydrodynamics (SPH). In order to further enhance its clinical relevance, this study employs a patient-specific vascular geometry and presents a successful validation against traditional finite volume method and 4D magnetic resonance imaging (MRI) data. The results have demonstrated that SPH is ideally suited to simulate the heart valve function due to its Lagrangian description of motion, which is a favourable feature for FSI. In addition, a novel methodology for the estimation of the wall shear stress (WSS) and other related haemodynamic parameters have been proposed from the SPH perspective. Finally, a detailed comparison of the haemodynamic parameters has been carried out for both native and mechanical aortic valve, with a particular emphasis on the clinical risks associated with the mechanical valve.

A Highly Sensitive Deep-Sea Hydrodynamic Pressure Sensor Inspired by Fish Lateral Line.

Hydrodynamic pressure sensors offer an auxiliary approach for ocean exploration by unmanned underwater vehicles (UUVs). However, existing hydrodynamic pressure sensors often lack the ability to monitor subtle hydrodynamic stimuli in deep-sea environments. In this study, we present the development of a deep-sea hydrodynamic pressure sensor (DSHPS) capable of operating over a wide range of water depths while maintaining exceptional hydrodynamic sensing performance. The DSHPS device was systematically optimized by considering factors such as piezoelectric polyvinylidene fluoride-trifluoroethylene/barium titanate [P(VDF-TrFE)/BTO] nanofibers, electrode configurations, sensing element dimensions, integrated circuits, and packaging strategies. The optimized DSHPS exhibited a remarkable pressure gradient response, achieving a minimum pressure difference detection capability of approximately 0.11 Pa. Additionally, the DSHPS demonstrated outstanding performance in the spatial positioning of dipole sources, which was elucidated through theoretical charge modeling and fluid-structure interaction (FSI) simulations. Furthermore, the integration of a high Young's modulus packaging strategy inspired by fish skull morphology ensured reliable sensing capabilities of the DSHPS even at depths of 1000 m in the deep sea. The DSHPS also exhibited consistent and reproducible positioning performance for subtle hydrodynamic stimulus sources across this wide range of water depths. We envision that the development of the DSHPS not only enhances our understanding of the evolutionary aspects of deep-sea canal lateral lines but also paves the way for the advancement of artificial hydrodynamic pressure sensors.