As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises s...As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.展开更多
Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simpli...Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics （CFD）-fluid structure interaction （FSI） techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called ＂dynamic mesh＂ technique. A new mesh movement approacb is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the joumal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function （UDF）. The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.展开更多
A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a t...A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual（GMRES） method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor（HHT） time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC（Harten–Lax–van Leer contact） scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel（LUSGS） approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.展开更多
This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles a...This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles are analyzed synchronously using FLUENT software, and the effects of fluid-structure interaction on the performance of the labyrinth seal are revealed. The results indicate that with the growth of tooth profile angle, the leakage rate of labyrinth seal tends to decrease first, and then increase. The results of fluid-structure interaction analysis are close to those of actual engineering. The effect of fluid-structure interaction makes tiny deformation in calculation mesh of piston and cylinder structure, and the coupling interaction affects the performance of the labyrinth seal.展开更多
The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstr...The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.展开更多
Vascular diseases such as aneurysm,hemadostenosis,aortic dissection are the primary causes of people’s death around world.As a result,it is significant to improve our knowledge about them,which can help to treat the ...Vascular diseases such as aneurysm,hemadostenosis,aortic dissection are the primary causes of people’s death around world.As a result,it is significant to improve our knowledge about them,which can help to treat the disease.Measuring the hemodynamic factor like the blood pressure,the wall shear stress(WSS)and the oscillatory shear index(OSI)is,however,still beyond the capabilities of in-vivo measurement techniques.So the use of mathematical models and numerical simulations for the studies of the blood flow in arteries and,in general,of the cardiovascular system,both in physiological and pathological conditions,has received an increasing attention in the biomedical community during the last two decades.Indeed,such studies aims at enhancing the current knowledge of the physiology of the cardiovascular system,as well as providing reliable tools for the medical doctors to predict the natural course of pathologies and,possibly,the occurrence of cardiovascular accidents.The computational vascular fluid-structure interaction(FSI)methodology is a numerical simulation method which is used to explain the hemodynamic factors.The WSS on the luminal wall and the mechanical stress in the vascular wall are directly related to the location of the lesion,and the blood flow strongly interacts with the vascular wall motion.The arterial wall continually adapts to the charge of its mechanical environment(due to,for example,growth,atrophy,remodelling,repair,ageing,and disease)and consequently undergoes several irreversible processes.Primary acute mechanisms of vascularFSI numerical simulation seem to be associated with(1)the arterial histology and the patient-specific complex geometry,(2)the typical mechanical properties of the layer,(3)properties of the blood is assumed as Newtonian fluid or non-Newtonian fluid based on the scale ofthe diameter of a vessel,(4)residual stress in the zero-pressure configuration.The arterial system naturally function under permanent physiological loading conditions.Fung defined the residual stress and measured the opening angle which varies greatly along the aortic tree.Consequently,most of these systems never experience a stress-free state in their’service life’,so a stress and strain fields are present in any in vivo obtained patientspecific cardiovascular geometry.The residual stress always be ignored in FSI simulation or be assumed to equal zero,and the vivo patient-specific artery geometry is assumed as zero-pressure configuration.To define the in vivo stress state of artery,an inverse problem needs to be solved:the undeformed shape of a body or its stress state in its deformed state needs to be determined given the deformed configuration and the loads causing this deformation.The modular inverse elastostatics method is used to resolve the pressure-induced stress state for in vivo imaging based on cardiovascular modeling proposed by Peirlinck.Here,we build a living vessel FSI model based on 4 key factors.In order to get the universal simulation results,we focus on idealized geometries of the vessel that represent healthy(physiological)conditions of the cerebral vasculature.Blood can be assumed as the Newtonian fluid at this scale.The anisotropic hyperelastic constitutive law(Gasser-Holzapfel-Ogden)is used in zero-pressure configuration.Afterwards,we propose the material parameters for the different constitutive models and the computational configurations.We demonstrate the importance of introducing the residual stress into vascular blood flow modeling by performing a comparing zero-pressure configuration and no-resistance configuration.We get the conclusion that the zero-pressure status model has smaller displacement and larger stress distribution compared with no-resistance stress model.Hence,the methodology presented here will be particularly useful to study the mechanobiological processes in the healthy and diseased vascular wall.展开更多
Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for nume...Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for numerical simulation based on computational fluid dynamics and nonlinear finite element analysis methods. The typical results of the numerical simulation and wind tunnel experiment, including the overall lift and deformation, are in good agreement with each other. The results obtained indicate that the effect of fluid-structure interaction is noticeable and should be considered for non-rigid airships. Flow- induced deformation can further intensify the upward lift force and pitching moment, which can lead to a large deformation. Under a wind speed of 15 m/s, the lift force of the non-rigid model is increased to approximatelv 60% compared with that of the rigid model under a high angle of attack.展开更多
Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The b...Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.展开更多
Tuned Mass Damper （TMD） was applied to an offshore structure to control ocean wave-induced vibration, In the analysis of the dynamic response of the offshore structure, fluid-structure interaction is considered and ...Tuned Mass Damper （TMD） was applied to an offshore structure to control ocean wave-induced vibration, In the analysis of the dynamic response of the offshore structure, fluid-structure interaction is considered and the errors, which occur in the linearization of the interaction, are investigated. For the investigation of the performance of TMD in controlling the vibration, both regular waves with different periods and irregular waves with different significant wave heights are used. Based on the numerical analysis it is concluded that the fluid-structure interaction should be considered in the evaluation of the capability of TMD in vibration control of offshore structures.展开更多
Tension leg platform （TLP） for offshore wind turbine support is a new type structure in wind energy utilization. The strong-interaction method is used in analyzing the coupled model, and the dynamic characteristics ...Tension leg platform （TLP） for offshore wind turbine support is a new type structure in wind energy utilization. The strong-interaction method is used in analyzing the coupled model, and the dynamic characteristics of the TLP for offshore wind turbine support are recognized. As shown by the calculated results： for the lower modes, the shapes are water＇s vibration, and the vibration of water induces the structure＇s swing; the mode shapes of the structure are complex, and can largely change among different members; the mode shapes of the platform are related to the tower＇s. The frequencies of the structure do not change much after adjusting the length of the tension cables and the depth of the platform; the TLP has good adaptability for the water depths and the environment loads. The change of the size and parameters of TLP can improve the dynamic characteristics, which can reduce the vibration of the TLP caused by the loads. Through the vibration analysis, the natural vibration frequencies of TLP can be distinguished from the frequencies of condition loads, and thus the resonance vibration can be avoided, therefore the offshore wind turbine can work normally in the complex conditions.展开更多
Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing me...Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing mechanism of fluid structure interaction between circulation valve and liquid of hydraulic shock absorber. The fluid mesh distortion was controlled by the CEL language, and the fluid struc^tre interaction mathematical model was established. The finite element model was established by ANSYS CFX software and was analyzed by dynamic mesh technique. The local sensitive computational area was meshed by prismatic grid, which could reduce the negative volume problem during the simulation. The circulation valve and liquid of hydraulic shock absorber were simulated and analyzed under the condition of sinusoidal inlet velocity loads. Flow characteristic and dynamics characteristic were obtained. The pressure distribution and the displacement of circulation value were obtained, and the acceleration curve of circulation valve was simulated and analyzed. The conformity of the final simulation results with the experimental datum indicates that this method is accurate and reliable to analyze the dynamics characteristic between circulation valve and liquid of hydraulic shock absorber, which can provide a theoretical foundation for optimizing hydraulic shock absorber in the future.展开更多
An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a fin...An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished.In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact,a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces.In addition,the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method.Based on this method,fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained.Numerical results such as maximum Root Mean Square(RMS)drag,general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results.This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.展开更多
The objective of this study was to develop, as well as validate the strongly coupled method (two-way fluid structural interaction (FSI)) used to simulate the transient FSI response of the vertical axis tidal turbine (...The objective of this study was to develop, as well as validate the strongly coupled method (two-way fluid structural interaction (FSI)) used to simulate the transient FSI response of the vertical axis tidal turbine (VATT) rotor, subjected to spatially varying inflow. Moreover, this study examined strategies on improving techniques used for mesh deformation that account for large displacement or deformation calculations. The blade's deformation for each new time step is considered in transient two-way FSI analysis, to make the design more reliable. Usually this is not considered in routine one-way FSI simulations. A rotor with four blades and 4-m diameter was modeled and numerically analyzed. We observed that two-way FSI, utilizing the strongly coupled method, was impossible for a complex model; and thereby using ANSYS-CFX and ANSYS-MECHANICAL in work bench, as given in ANSYS-WORKBENCH, helped case examples 22 and 23, by giving an error when the solution was run. To make the method possible and reduce the computational power, a novel technique was used to transfer the file in ANSYS-APDL to obtain the solution and results. Consequently, the results indicating a two-way transient FSI analysis is a time- and resource-consuming job, but with our proposed technique we can reduce the computational time. The ANSYS STRUCTURAL results also uncover that stresses and deformations have higher values for two-way FSI as compared to one-way FSI. Similarly, fluid flow CFX results for two-way FSI are closer to experimental results as compared to one-way simulation results. Additionally, this study shows that, using the proposed method we can perform coupled simulation with simple multi-node PCs (core i5).展开更多
Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting...Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.展开更多
Deep-seated potassium-rich brines were identified in the Jiangling Basin,South China.Although magmatichydrothermal sources have been proposed,the relationship between brine-type potash mineralization and volcanism rem...Deep-seated potassium-rich brines were identified in the Jiangling Basin,South China.Although magmatichydrothermal sources have been proposed,the relationship between brine-type potash mineralization and volcanism remains unclear.In this study,U-Pb geochronology,geochemistry,fluid inclusion and C-O isotopic compositions of hydrothermal vein minerals in the Jiangling Basin are examined.Laser ablation U-Pb dating of calcite veins indicates that the ages are slightly younger than the formation age of the Balingshan basalt.Fluid inclusions in hydrothermal minerals show medium–low homogenization temperatures(160–220℃)and low salinities(0.14 to 4.9 wt%NaCl eqv.)and densities(0.882–0.944 g/cm^(3)).The liquid compositions of fluid inclusions in calcite veins from sedimentary strata have higher contents of potassium,compared with those from basalt.The coupled negativeδ^(13)CPDB(-10.3‰to-8.0‰)and positiveδ^(18)OSMOW(17.4‰to 20.7‰)values imply that calcite precipitation resulted from CO_(2)degassing of the basaltic magmatic fluids,as indicated by the gas composition of these inclusions in hydrothermal minerals.Rare earth element patterns indicate that water-rock interaction between hydrothermal fluids and sedimentary wall rocks contributed to the calcite precipitation in sedimentary strata.It is proposed that high-temperature water-rock interaction between magmatic fluids and sedimentary strata resulted in the potassium enrichment in fluids,interpreted as one of the sources of potassium-rich brines in the Jiangling Basin.展开更多
Numerical simulations of evolution characteristics of slug flow across a 90°pipe bend have been carried out to study the fluid−structure interaction response induced by internal slug flow.The two-phase flow patte...Numerical simulations of evolution characteristics of slug flow across a 90°pipe bend have been carried out to study the fluid−structure interaction response induced by internal slug flow.The two-phase flow patterns and turbulence were modelled by using the volume of fluid(VOF)model and the Realizable k−εturbulence model respectively.Firstly,validation of the CFD model was carried out and the desirable results were obtained.The different flow patterns and the time-average mean void fraction was coincident with the reported experimental data.Simulations of different cases of slug flow have been carried out to show the effects of superficial gas and liquid velocity on the evolution characteristics of slug flow.Then,a one-way coupled fluid-structure interaction framework was established to investigate the slug flow interaction with a 90°pipe bend under various superficial liquid and gas velocities.It was found that the maximum total deformation and equivalent stress increased with the increasing superficial gas velocity,while decreased with the increasing superficial liquid velocity.In addition,the total deformation and equivalent stress has obvious periodic fluctuation.Furthermore,the distribution position of maximum deformation and stress was related to the evolution of slug flow.With the increasing superficial gas velocity,the maximum total deformation was mainly located at the 90°pipe bend.But as the superficial liquid velocity increases,the maximum total deformation was mainly located in the horizontal pipe section.Consequently,the slug flow with higher superficial gas velocity will induce more serious cyclical impact on the 90°pipe bend.展开更多
As supercavitating projectiles move at high speed, the periodic impacts ("tail-slap") on the interior surface of the cavity generally occur due to disturbances. The interactions between the projectile and th...As supercavitating projectiles move at high speed, the periodic impacts ("tail-slap") on the interior surface of the cavity generally occur due to disturbances. The interactions between the projectile and the water/cavity interface are the sources of structural vibrations, which affect the guidance of the vehicle and undermine the structural reliability. The Fluid/Structure Interaction calculation procedure of the tail-slaps of supercavitating projectile is established, and the dynamic behaviours of the projectile operating in tail-slap conditions with and without considering Fluid/Structure Interaction are obtained and compared. The responses of the projectile riding a reducing cavity are studied, and the effect of Fluid/Structure Interaction is also analyzed. The results show that the angular velocity of projectile increases as the body slowing down, and the amplitude of the elastic displacement response decreases at the beginning and increases when the cavity size is close to the diameter of the tail of projectile. The effect of Fluid/Structure Interaction reduces the amplitudes and frequencies of the impact loads and the vibration responses of the body, and when the speed is higher, the effect is more apparent.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
The material damage of parachute may occur in parachutes at high speeds,and the growth of tearing may finally lead to failure of aerospace mission.In order to study the damage mechanism of parachute,a material failure...The material damage of parachute may occur in parachutes at high speeds,and the growth of tearing may finally lead to failure of aerospace mission.In order to study the damage mechanism of parachute,a material failure model is proposed to simulate the failure of canopy fabric.The inflation process of supersonic parachute is studied numerically based on Arbitrary Lagrange Euler(ALE)method.The ALE method with material failure can predict the transient parachute shape with damage propagation as well as the flow characteristics in the parachute inflation process,and the simulated dynamic opening load is consistent with the flight test.The damage propagation mechanism of parachute is then investigated,and the effect of parachute velocity on the damage process is discussed.The results show that the canopy tears apart by the fast flow from the initial damaged area and the damaged canopy shape leads to the asymmetric change of the flow structure.With the increase of Mach number,the canopy tearing speed increases,and the tearing directions become uncertain at high Mach numbers.The dynamic load when damage occurs increases with the Mach number,and is proportional to the dynamic pressure above the critical Mach number.展开更多
基金supported by the National Natural Science Foundation of China(No.12372233)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University,China(No.25GH01020005)the“111 Project”of China(No.B17037)。
文摘As a multidisciplinary phenomenon,panel aeroelasticity in shock-dominated flow is featured by two primary interactions:Fluid-Structure Interactions(FSIs)and Shock-Boundary Layer Interactions(SBLIs).The former raises structural concerns,and the latter is of aerodynamic interest.Thus,panel aeroelasticity in shock-dominated flow represents a vital topic for the development and optimization of supersonic vehicles and propulsion systems.This review systematically summarizes recent advances in the methodologies applied to capture structural and fluid dynamics,including theoretical models,numerical simulations,and wind tunnel experiments.The application of data-driven modal decomposition,an advanced technique to extract physically crucial features,on the topic is introduced.From the perspective of FSIs,the distinctive aeroelastic behaviors in shock-dominated flow,including hysteresis phenomena and nonlinear responses,are highlighted.From the perspective of SBLIs,the modifications in their spatial and temporal characteristics imposed by the aeroelastic responses are emphasized.Motivated by the interaction between the shock waves and structural response,different strategies have been proposed to implement aeroelastic suppression and shock control,which have the potential to enhance structural safety and aerodynamic performance in the next generation of high-speed flight vehicles.
基金supported by National Hi-tech Research and Development Program of China (863 Program, Grant No. 2009AA04Z413)Zhejiang Provincial Natural Science Foundation of China (Grant No. Y1110109)
文摘Journal bearings are important parts to keep the high dynamic performance of rotor machinery. Some methods have already been proposed to analysis the flow field of journal bearings, and in most of these methods simplified physical model and classic Reynolds equation are always applied. While the application of the general computational fluid dynamics （CFD）-fluid structure interaction （FSI） techniques is more beneficial for analysis of the fluid field in a journal bearing when more detailed solutions are needed. This paper deals with the quasi-coupling calculation of transient fluid dynamics of oil film in journal bearings and rotor dynamics with CFD-FSI techniques. The fluid dynamics of oil film is calculated by applying the so-called ＂dynamic mesh＂ technique. A new mesh movement approacb is presented while the dynamic mesh models provided by FLUENT are not suitable for the transient oil flow in journal bearings. The proposed mesh movement approach is based on the structured mesh. When the joumal moves, the movement distance of every grid in the flow field of bearing can be calculated, and then the update of the volume mesh can be handled automatically by user defined function （UDF）. The journal displacement at each time step is obtained by solving the moving equations of the rotor-bearing system under the known oil film force condition. A case study is carried out to calculate the locus of the journal center and pressure distribution of the journal in order to prove the feasibility of this method. The calculating results indicate that the proposed method can predict the transient flow field of a journal bearing in a rotor-bearing system where more realistic models are involved. The presented calculation method provides a basis for studying the nonlinear dynamic behavior of a general rotor-bearing system.
文摘A fluid–structure interaction method combining a nonlinear finite element algorithm with a preconditioning finite volume method is proposed in this paper to simulate parachute transient dynamics. This method uses a three-dimensional membrane–cable fabric model to represent a parachute system at a highly folded configuration. The large shape change during parachute inflation is computed by the nonlinear Newton–Raphson iteration and the linear system equation is solved by the generalized minimal residual（GMRES） method. A membrane wrinkling algorithm is also utilized to evaluate the special uniaxial tension state of membrane elements on the parachute canopy. In order to avoid large time expenses during structural nonlinear iteration, the implicit Hilber–Hughes–Taylor（HHT） time integration method is employed. For the fluid dynamic simulations, the Roe and HLLC（Harten–Lax–van Leer contact） scheme has been modified and extended to compute flow problems at all speeds. The lower–upper symmetric Gauss–Seidel（LUSGS） approximate factorization is applied to accelerate the numerical convergence speed. Finally,the test model of a highly folded C-9 parachute is simulated at a prescribed speed and the results show similar characteristics compared with experimental results and previous literature.
基金the Science and Technology Projects of Liaoning Province(No.2012219020)the China Postdoctoral Science Foundation(No.2013M541249)
文摘This paper addresses the issue of reciprocating compressors staggered labyrinth seal structure. The internal flow field of sealed structure, the displacement of cylinder and piston for different tooth profile angles are analyzed synchronously using FLUENT software, and the effects of fluid-structure interaction on the performance of the labyrinth seal are revealed. The results indicate that with the growth of tooth profile angle, the leakage rate of labyrinth seal tends to decrease first, and then increase. The results of fluid-structure interaction analysis are close to those of actual engineering. The effect of fluid-structure interaction makes tiny deformation in calculation mesh of piston and cylinder structure, and the coupling interaction affects the performance of the labyrinth seal.
文摘The effects of the renal artery stenosis （RAS） on the blood flow and vessel walls are investigated. The pulsatile blood flow through an anatomically realistic model of the abdominal aorta and renal arteries reconstructed from CT-scan images is simulated, which incorporates the fluid-structure interaction （FSI）. In addition to the investigation of the RAS effects on the wall shear stress and the displacement of the vessel wall, it is determined that the RAS leads to decrease in the renal mass flow. This may cause the activation of the renin-angiotension system and results in severe hypertension.
基金supported by the National Natural Science Foundation of China ( 11732001)
文摘Vascular diseases such as aneurysm,hemadostenosis,aortic dissection are the primary causes of people’s death around world.As a result,it is significant to improve our knowledge about them,which can help to treat the disease.Measuring the hemodynamic factor like the blood pressure,the wall shear stress(WSS)and the oscillatory shear index(OSI)is,however,still beyond the capabilities of in-vivo measurement techniques.So the use of mathematical models and numerical simulations for the studies of the blood flow in arteries and,in general,of the cardiovascular system,both in physiological and pathological conditions,has received an increasing attention in the biomedical community during the last two decades.Indeed,such studies aims at enhancing the current knowledge of the physiology of the cardiovascular system,as well as providing reliable tools for the medical doctors to predict the natural course of pathologies and,possibly,the occurrence of cardiovascular accidents.The computational vascular fluid-structure interaction(FSI)methodology is a numerical simulation method which is used to explain the hemodynamic factors.The WSS on the luminal wall and the mechanical stress in the vascular wall are directly related to the location of the lesion,and the blood flow strongly interacts with the vascular wall motion.The arterial wall continually adapts to the charge of its mechanical environment(due to,for example,growth,atrophy,remodelling,repair,ageing,and disease)and consequently undergoes several irreversible processes.Primary acute mechanisms of vascularFSI numerical simulation seem to be associated with(1)the arterial histology and the patient-specific complex geometry,(2)the typical mechanical properties of the layer,(3)properties of the blood is assumed as Newtonian fluid or non-Newtonian fluid based on the scale ofthe diameter of a vessel,(4)residual stress in the zero-pressure configuration.The arterial system naturally function under permanent physiological loading conditions.Fung defined the residual stress and measured the opening angle which varies greatly along the aortic tree.Consequently,most of these systems never experience a stress-free state in their’service life’,so a stress and strain fields are present in any in vivo obtained patientspecific cardiovascular geometry.The residual stress always be ignored in FSI simulation or be assumed to equal zero,and the vivo patient-specific artery geometry is assumed as zero-pressure configuration.To define the in vivo stress state of artery,an inverse problem needs to be solved:the undeformed shape of a body or its stress state in its deformed state needs to be determined given the deformed configuration and the loads causing this deformation.The modular inverse elastostatics method is used to resolve the pressure-induced stress state for in vivo imaging based on cardiovascular modeling proposed by Peirlinck.Here,we build a living vessel FSI model based on 4 key factors.In order to get the universal simulation results,we focus on idealized geometries of the vessel that represent healthy(physiological)conditions of the cerebral vasculature.Blood can be assumed as the Newtonian fluid at this scale.The anisotropic hyperelastic constitutive law(Gasser-Holzapfel-Ogden)is used in zero-pressure configuration.Afterwards,we propose the material parameters for the different constitutive models and the computational configurations.We demonstrate the importance of introducing the residual stress into vascular blood flow modeling by performing a comparing zero-pressure configuration and no-resistance configuration.We get the conclusion that the zero-pressure status model has smaller displacement and larger stress distribution compared with no-resistance stress model.Hence,the methodology presented here will be particularly useful to study the mechanobiological processes in the healthy and diseased vascular wall.
基金the National Natural Science Foundation of China (11202215 and 11332011)the Youth Innovation Promotion Association of CAS (2015015)
文摘Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for numerical simulation based on computational fluid dynamics and nonlinear finite element analysis methods. The typical results of the numerical simulation and wind tunnel experiment, including the overall lift and deformation, are in good agreement with each other. The results obtained indicate that the effect of fluid-structure interaction is noticeable and should be considered for non-rigid airships. Flow- induced deformation can further intensify the upward lift force and pitching moment, which can lead to a large deformation. Under a wind speed of 15 m/s, the lift force of the non-rigid model is increased to approximatelv 60% compared with that of the rigid model under a high angle of attack.
文摘Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.
文摘Tuned Mass Damper （TMD） was applied to an offshore structure to control ocean wave-induced vibration, In the analysis of the dynamic response of the offshore structure, fluid-structure interaction is considered and the errors, which occur in the linearization of the interaction, are investigated. For the investigation of the performance of TMD in controlling the vibration, both regular waves with different periods and irregular waves with different significant wave heights are used. Based on the numerical analysis it is concluded that the fluid-structure interaction should be considered in the evaluation of the capability of TMD in vibration control of offshore structures.
文摘Tension leg platform （TLP） for offshore wind turbine support is a new type structure in wind energy utilization. The strong-interaction method is used in analyzing the coupled model, and the dynamic characteristics of the TLP for offshore wind turbine support are recognized. As shown by the calculated results： for the lower modes, the shapes are water＇s vibration, and the vibration of water induces the structure＇s swing; the mode shapes of the structure are complex, and can largely change among different members; the mode shapes of the platform are related to the tower＇s. The frequencies of the structure do not change much after adjusting the length of the tension cables and the depth of the platform; the TLP has good adaptability for the water depths and the environment loads. The change of the size and parameters of TLP can improve the dynamic characteristics, which can reduce the vibration of the TLP caused by the loads. Through the vibration analysis, the natural vibration frequencies of TLP can be distinguished from the frequencies of condition loads, and thus the resonance vibration can be avoided, therefore the offshore wind turbine can work normally in the complex conditions.
基金Project(51275542) supported by the National Natural Science Foundation of Chinaproject(CDJXS12110010) supported by the Fundamental Research Funds for the Central Universities of China
文摘Based on the working principle and the damping characteristic of hydraulic shock absorber, a fluid structure interaction method was presented, which was used to analyze the microcosmic and high-frequency processing mechanism of fluid structure interaction between circulation valve and liquid of hydraulic shock absorber. The fluid mesh distortion was controlled by the CEL language, and the fluid struc^tre interaction mathematical model was established. The finite element model was established by ANSYS CFX software and was analyzed by dynamic mesh technique. The local sensitive computational area was meshed by prismatic grid, which could reduce the negative volume problem during the simulation. The circulation valve and liquid of hydraulic shock absorber were simulated and analyzed under the condition of sinusoidal inlet velocity loads. Flow characteristic and dynamics characteristic were obtained. The pressure distribution and the displacement of circulation value were obtained, and the acceleration curve of circulation valve was simulated and analyzed. The conformity of the final simulation results with the experimental datum indicates that this method is accurate and reliable to analyze the dynamics characteristic between circulation valve and liquid of hydraulic shock absorber, which can provide a theoretical foundation for optimizing hydraulic shock absorber in the future.
基金co-supported by National Nature Sciences Foundation of China(Nos 11972192,11172137)A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)。
文摘An Arbitrary Lagrangian-Eulerian(ALE)approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation.A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished.In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact,a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces.In addition,the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method.Based on this method,fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained.Numerical results such as maximum Root Mean Square(RMS)drag,general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results.This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.
基金supported by the National Natural Science Foundation of China (Nos. 51209060 and 51106034)the ‘111’ Project Foundation from Ministry of Education and State Administration of Foreign Experts Affairs (No. B07019), Chinathe National Special Foundation for Ocean Energy (No. GHME2010CY01)
文摘The objective of this study was to develop, as well as validate the strongly coupled method (two-way fluid structural interaction (FSI)) used to simulate the transient FSI response of the vertical axis tidal turbine (VATT) rotor, subjected to spatially varying inflow. Moreover, this study examined strategies on improving techniques used for mesh deformation that account for large displacement or deformation calculations. The blade's deformation for each new time step is considered in transient two-way FSI analysis, to make the design more reliable. Usually this is not considered in routine one-way FSI simulations. A rotor with four blades and 4-m diameter was modeled and numerically analyzed. We observed that two-way FSI, utilizing the strongly coupled method, was impossible for a complex model; and thereby using ANSYS-CFX and ANSYS-MECHANICAL in work bench, as given in ANSYS-WORKBENCH, helped case examples 22 and 23, by giving an error when the solution was run. To make the method possible and reduce the computational power, a novel technique was used to transfer the file in ANSYS-APDL to obtain the solution and results. Consequently, the results indicating a two-way transient FSI analysis is a time- and resource-consuming job, but with our proposed technique we can reduce the computational time. The ANSYS STRUCTURAL results also uncover that stresses and deformations have higher values for two-way FSI as compared to one-way FSI. Similarly, fluid flow CFX results for two-way FSI are closer to experimental results as compared to one-way simulation results. Additionally, this study shows that, using the proposed method we can perform coupled simulation with simple multi-node PCs (core i5).
基金supported by the National Natural Science Foundation of China(51809169,51879159)Chang Jiang Scholars Program(T2014099)+2 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)National Key Research and Development Program of China(2019YFB1704203,2019YFC0312400).
文摘Complex flow around floating structures is a highly nonlinear problem,and it is a typical feature in ship and ocean engineering.Traditional experimental methods and potential flow theory have limitations in predicting complex viscous flows.With the improvement of high-performance computing and the development of numerical techniques,computational fluid dynamics(CFD)has become increasingly powerful in predicting the complex viscous flow around floating structures.This paper reviews the recent progress in CFD techniques for numerical solutions of typical complex viscous flows in ship and ocean engineering.Applications to free-surface flows,breaking bow waves of high-speed ship,ship hull-propeller-rudder interaction,vortexinduced vibration of risers,vortex-induced motions of deep-draft platforms,and floating offshore wind turbines are discussed.Typical techniques,including volume of fluid for sharp interface,dynamic overset grid,detached eddy simulation,and fluid-structure coupling,are reviewed along with their applications.Some novel techniques,such as high-efficiency Cartesian grid method and GPU acceleration technique,are discussed in the last part as the future perspective for further enhancement of accuracy and efficiency for CFD simulations of complex flow in ship and ocean engineering.
基金supported by the Central Public Welfare Scientific Research Basic Scientific Research Business Expenses(Grant Nos.KK2005,KY1603)National Natural Science Foundation of China(Grant No.U20A2092)+1 种基金the National Basic Research Program of China(973 Program)(Grant No.2011CB403007)the China Geological Survey(Grant No.DD20190606)。
文摘Deep-seated potassium-rich brines were identified in the Jiangling Basin,South China.Although magmatichydrothermal sources have been proposed,the relationship between brine-type potash mineralization and volcanism remains unclear.In this study,U-Pb geochronology,geochemistry,fluid inclusion and C-O isotopic compositions of hydrothermal vein minerals in the Jiangling Basin are examined.Laser ablation U-Pb dating of calcite veins indicates that the ages are slightly younger than the formation age of the Balingshan basalt.Fluid inclusions in hydrothermal minerals show medium–low homogenization temperatures(160–220℃)and low salinities(0.14 to 4.9 wt%NaCl eqv.)and densities(0.882–0.944 g/cm^(3)).The liquid compositions of fluid inclusions in calcite veins from sedimentary strata have higher contents of potassium,compared with those from basalt.The coupled negativeδ^(13)CPDB(-10.3‰to-8.0‰)and positiveδ^(18)OSMOW(17.4‰to 20.7‰)values imply that calcite precipitation resulted from CO_(2)degassing of the basaltic magmatic fluids,as indicated by the gas composition of these inclusions in hydrothermal minerals.Rare earth element patterns indicate that water-rock interaction between hydrothermal fluids and sedimentary wall rocks contributed to the calcite precipitation in sedimentary strata.It is proposed that high-temperature water-rock interaction between magmatic fluids and sedimentary strata resulted in the potassium enrichment in fluids,interpreted as one of the sources of potassium-rich brines in the Jiangling Basin.
基金sponsored by the National Natural Science Foundation of China(Grant No.51779143)the Oceanic Interdisciplinary Program of Shanghai Jiao Tong University(Grant No.SL2020ZD101)the Cultivation of Scientific Research Ability of Young Talents of Shanghai Jiao Tong University(Grant No.19X100040072).
文摘Numerical simulations of evolution characteristics of slug flow across a 90°pipe bend have been carried out to study the fluid−structure interaction response induced by internal slug flow.The two-phase flow patterns and turbulence were modelled by using the volume of fluid(VOF)model and the Realizable k−εturbulence model respectively.Firstly,validation of the CFD model was carried out and the desirable results were obtained.The different flow patterns and the time-average mean void fraction was coincident with the reported experimental data.Simulations of different cases of slug flow have been carried out to show the effects of superficial gas and liquid velocity on the evolution characteristics of slug flow.Then,a one-way coupled fluid-structure interaction framework was established to investigate the slug flow interaction with a 90°pipe bend under various superficial liquid and gas velocities.It was found that the maximum total deformation and equivalent stress increased with the increasing superficial gas velocity,while decreased with the increasing superficial liquid velocity.In addition,the total deformation and equivalent stress has obvious periodic fluctuation.Furthermore,the distribution position of maximum deformation and stress was related to the evolution of slug flow.With the increasing superficial gas velocity,the maximum total deformation was mainly located at the 90°pipe bend.But as the superficial liquid velocity increases,the maximum total deformation was mainly located in the horizontal pipe section.Consequently,the slug flow with higher superficial gas velocity will induce more serious cyclical impact on the 90°pipe bend.
基金Sponsored by the Fundamental Research Funds for the Central Universities(Grant No.HIT.NSRIF.201159)the National Natural Science Foundation of China(Grant No.51149003)
文摘As supercavitating projectiles move at high speed, the periodic impacts ("tail-slap") on the interior surface of the cavity generally occur due to disturbances. The interactions between the projectile and the water/cavity interface are the sources of structural vibrations, which affect the guidance of the vehicle and undermine the structural reliability. The Fluid/Structure Interaction calculation procedure of the tail-slaps of supercavitating projectile is established, and the dynamic behaviours of the projectile operating in tail-slap conditions with and without considering Fluid/Structure Interaction are obtained and compared. The responses of the projectile riding a reducing cavity are studied, and the effect of Fluid/Structure Interaction is also analyzed. The results show that the angular velocity of projectile increases as the body slowing down, and the amplitude of the elastic displacement response decreases at the beginning and increases when the cavity size is close to the diameter of the tail of projectile. The effect of Fluid/Structure Interaction reduces the amplitudes and frequencies of the impact loads and the vibration responses of the body, and when the speed is higher, the effect is more apparent.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
基金the National Natural Science Foundation of China(No.11972192).
文摘The material damage of parachute may occur in parachutes at high speeds,and the growth of tearing may finally lead to failure of aerospace mission.In order to study the damage mechanism of parachute,a material failure model is proposed to simulate the failure of canopy fabric.The inflation process of supersonic parachute is studied numerically based on Arbitrary Lagrange Euler(ALE)method.The ALE method with material failure can predict the transient parachute shape with damage propagation as well as the flow characteristics in the parachute inflation process,and the simulated dynamic opening load is consistent with the flight test.The damage propagation mechanism of parachute is then investigated,and the effect of parachute velocity on the damage process is discussed.The results show that the canopy tears apart by the fast flow from the initial damaged area and the damaged canopy shape leads to the asymmetric change of the flow structure.With the increase of Mach number,the canopy tearing speed increases,and the tearing directions become uncertain at high Mach numbers.The dynamic load when damage occurs increases with the Mach number,and is proportional to the dynamic pressure above the critical Mach number.
