The existing research for unsteady flow field and the corresponding flow induced vibration analysis of centrifugal pump are mainly carried out respectively without considering the interaction between fluid and structu...The existing research for unsteady flow field and the corresponding flow induced vibration analysis of centrifugal pump are mainly carried out respectively without considering the interaction between fluid and structure. The ignorance of fluid structure interaction (FSI) means that the energy transfer between fluid and structure is neglected. To some extent, the accuracy and reliability of unsteady flow and rotor deflection analysis should be affected by this interaction mechanism. In this paper, a combined calculation between two executables for turbulent flow and vibrating structure was established using two-way coupling method to study the effect of FSI. Pressure distributions, radial forces, rotor deflection and equivalent stress are analyzed. The results show that the FSI effect to pressure distribution in flow field is complex. The pressure distribution is affected not only around impeller outlet where different variation trends of pressure values with and without FSI appear according to different relative positions between blade and cutwater, but also in the diffusion section of volute. Variation trends of peak values of radial force amplitude calculated with and without FSI are nearly same under high flow rate and designed conditions while the peak value with FSI is slightly smaller, and differently, the peak value with FSI is larger with low flow rate. In addition, the effect of FSI on the angle of radial force is quite complex, especially under 0.5Q condition. Fluctuation of radial deflection of the rotor has obvious four periods, of which the extent is relatively small under design condition and is relatively large under off-design condition. Finally, fluctuations of equivalent stress with time are obvious under different conditions, and stress value is small. The proposed research establishes the FSI calculation method for centrifugal pump analysis, and ensures the existing affect by fluid structure interaction.展开更多
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.展开更多
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.展开更多
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.展开更多
The interaction between fluid and a down-pumping pitched blade turbine fixed with a flexible shaft in the stirred vessel, as a typical fluid structure interaction phenomenon, was simulated by coupling the Computationa...The interaction between fluid and a down-pumping pitched blade turbine fixed with a flexible shaft in the stirred vessel, as a typical fluid structure interaction phenomenon, was simulated by coupling the Computational Fluid Dynamics and Computational Structural Dynamics. Based on the verification of the simulated impeller torque and dimensionless shaft bending moment with experimental result, the dimensionless shaft bending moment and various loads acting on impeller(including lateral force, axial force and bending moment) were discussed in detail. By separating and extracting the fluid and structural components from those loads, the results show that the shaft bending moment mainly results from the lateral force on impeller although the axial force on impeller is much larger. The impeller mass imbalance increases the shaft bending moment and the lateral force on impeller, but has little influence on the axial force and bending moment acting on impeller. The dominant frequencies of impeller forces are macro-frequency, speed frequency and blade passing frequency, and are associated with the impeller mass imbalance.展开更多
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.展开更多
A blockage of blood vessels resulting from thrombus or plaque deposit causes serious cardiovascular diseases. This study developed a computational model of blood flow and drug transport to investigate the effectivenes...A blockage of blood vessels resulting from thrombus or plaque deposit causes serious cardiovascular diseases. This study developed a computational model of blood flow and drug transport to investigate the effectiveness of drug delivery to the stenotic sites. A three-dimensional (3D) model of the curved stenotic right coronary artery (RCA) was reconstructed based on the clinical angiogram image. Then, blood flow and drug transport with the flexible RCA wall were simulated using the fluid structure interaction (FSI) analysis and compared with the rigid RCA wall. Results showed that the maximal total displacement and von Mises stress of the flexible RCA model are 2.14 mm and 92.06 kPa. In addition, the effective injecting time point for the best performance of drug delivery was found to be between 0 s and 0.15 s (i.e., the fluid acceleration region) for both rigid and flexible RCA models. However, there was no notable difference in the ratio of particle deposition to the stenotic areas between the rigid and flexible RCA models. This study will be significantly useful to the design of a drug delivery system for the treatment of the stenotic arteries by targeting drugs selectively to the stenotic sites.展开更多
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.展开更多
A loss-of-coolant accident(LOCA)is one of the basic design considerations for nuclear reactor safety analysis.A LOCA induces propagation of a depressurization wave in the coolant,exerting hydrodynamic forces on struct...A loss-of-coolant accident(LOCA)is one of the basic design considerations for nuclear reactor safety analysis.A LOCA induces propagation of a depressurization wave in the coolant,exerting hydrodynamic forces on structures viafiuid–structure interaction(FSI).The analysis of hydrodynamic forces on the core structures during a LOCA process is indispensable.We describe the implementation of a numerical strategy for prestressed structures.It consists of an initialization and a restarted transient analysis process,all implemented via the ANSYS Workbench by system coupling of ANSYS and Fluent.Our strategy is validated by making extensive comparisons of the pressures,displacements,and strains on various locations between the simulation and reported measurements.The approach is appealing for dynamic analysis of other prestressed structures,owing to the good popularity and acknowledgement of ANSYS and Fluent in both academia and industry.展开更多
In this paper, the second order steady drift forces on the ships and other floating offshore structures are calculated by the far field method. The amplitudes of diffracted waves and radiated waves at infinity are obt...In this paper, the second order steady drift forces on the ships and other floating offshore structures are calculated by the far field method. The amplitudes of diffracted waves and radiated waves at infinity are obtained by the two-dimensional source distribution and strip-theory method. For the twin hull structure, the hydrodynamic interaction between the two hulls is taken into account. The drift forces on cross sections of Lewis type as well as on the semi-submersibles are computed. The theoretical results obtained by the present method agree fairly well with experimental results.展开更多
Conventional coupled BE/FE (Boundary-Element/Finite-Element) method and modeling of structural-acoustic interaction has shown its promise and potential in the design and analysis of various structural-acoustic inter...Conventional coupled BE/FE (Boundary-Element/Finite-Element) method and modeling of structural-acoustic interaction has shown its promise and potential in the design and analysis of various structural-acoustic interaction applications. Unified combined acoustic and aerodynamic loading on the structure is synthesized using two approaches. Firstly, by linear superposition of the acoustic pressure disturbance to the aeroelastic problem, the effect of acoustic pressure disturbance to the aeroelastic structure is considered to consist of structural motion independent incident acoustic pressure and structural motion dependent acoustic pressure, which is known as the scattering pressure, referred here as the acoustic aerodynamic analogy. Secondly, by synthesizing the acoustic and aerodynamic effects on elastic structure using an elegant, effective and unified approach, both acoustic and aerodynamic effect on solid structural boundaries can be formulated as a boundary value problem governed by second order differential equations which lead to solutions expressible as surface integral equations. The unified formulation of the acousto-aeroelastic problem is amenable for simultaneous solution, although certain prevailing situations allow the solution of the equations independently. For this purpose, the unsteady aerodynamic problem which was earlier utilizes well-established lifting surface method is reformulated using Boundary Element (BE) approach. These schemes are outlined and worked out with examples.展开更多
This paper develops a numerical code for modelling liquid sloshing.The coupled boundary element-finite element method was used to solve the Laplace equation for inviscid fluid and nonlinear free surface boundary condi...This paper develops a numerical code for modelling liquid sloshing.The coupled boundary element-finite element method was used to solve the Laplace equation for inviscid fluid and nonlinear free surface boundary conditions.Using Nakayama and Washizu’s results,the code performance was validated.Using the developed numerical mode,we proposed artificial neural network(ANN)and genetic algorithm(GA)methods for evaluating sloshing loads and comparing them.To compare the efficiency of the suggested methods,the maximum free surface displacement and the maximum horizontal force exerted on a rectangular tank’s perimeter are examined.It can be seen from the results that both ANNs and GAs can accurately predict η_(max) and F_(max).展开更多
Certain insect species have been observed to exploit the resonance mechanism of their wings.In order to achieve resonance and optimize aerodynamic performance,the conventional approach is to set the flapping frequency...Certain insect species have been observed to exploit the resonance mechanism of their wings.In order to achieve resonance and optimize aerodynamic performance,the conventional approach is to set the flapping frequency of flexible wings based on the Traditional Structural Modal(TSM)analysis.However,there exists controversy among researchers regarding the relationship between frequency and aerodynamic performance.Recognizing that the structural response of wings can be influenced by the surrounding air vibrations,an analysis known as Acoustic Structure Interaction Modal(ASIM)is introduced to calculate the resonant frequency.In this study,Fluid Structure Interaction(FSI)simulations are employed to investigate the aerodynamic performance of flapping wings at modal frequencies derived from both TSM and ASIM analyses.The performance is evaluated for various mass ratios and frequency ratios,and the findings indicate that the deformation and changes in vortex structure exhibit similarities at mass ratios that yield the highest aerodynamic performance.Notably,the flapping frequency associated with the maximum time-averaged vertical force coefficient at each mass ratio closely aligns with the ASIM frequency,as does the frequency corresponding to maximum efficiency.Thus,the ASIM analysis can provide an effective means for predicting the optimal flapping frequency for flexible wings.Furthermore,it enables the prediction that flexible wings with varying mass ratios will exhibit similar deformation and vortex structure changes.This paper offers a fresh perspective on the ongoing debate concerning the resonance mechanism of Flexible Flapping Wings(FFWs)and proposes an effective methodology for predicting their aerodynamic performance.展开更多
The data information transfer and time marching strategies between computational fluid dynamics (CFD) and computational structural dynamics (CSD) play crucial roles on the aeroelastic analysis in a time domain. An...The data information transfer and time marching strategies between computational fluid dynamics (CFD) and computational structural dynamics (CSD) play crucial roles on the aeroelastic analysis in a time domain. An improved CFD/CSD coupled system is designed, including an interpolation method and an improved loosely coupled algorithm. The interpolation method based on boundary element method (BEM) is developed to transfer aerodynamic loads and structural displacements between CFD and CSD grid systems, it can be universally used in fluid structural interaction solution by keeping energy conservation. The improved loosely coupled algo-rithm is designed, thus it improves the computational accuracy and efficiency. The new interface is performed on the two-dimensional (2-D) extrapolation and the aeroelastie response of AGARD445.6 wing. Results show that the improved interface has a superior accuracy.展开更多
Experimental study and theoretical analysis show that the critical value of relative wave height (H / d)b given by Goda and the critical wave steepness (H / L)b given by Michell and Miche can be adopted as the spillin...Experimental study and theoretical analysis show that the critical value of relative wave height (H / d)b given by Goda and the critical wave steepness (H / L)b given by Michell and Miche can be adopted as the spilling breaking indices of regular waves. According to the same principle, a systematic theoretical analysis and experiment of irregular wave have been done by the authors in order to solve the breaking problem of irregular waves. It is indicated that the authors' method for determining wave breaking of regular waves can also be used for irregular waves.展开更多
Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model(moving particle semi-implicit(MPS)) is c...Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model(moving particle semi-implicit(MPS)) is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model(MPS or HMPS) in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model(Enhanced MPS), thus, the developed coupled solvers for fluid structure interaction(FSI) are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics(SPH)-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH(ISPH) fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.展开更多
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).展开更多
In this paper, vortex-induced vibrations of a cylinder are simulated by use of ANSYS CFX simulation code. The cylinder is treated as a rigid body and transverse displacements are obtained by use of a one degree of fre...In this paper, vortex-induced vibrations of a cylinder are simulated by use of ANSYS CFX simulation code. The cylinder is treated as a rigid body and transverse displacements are obtained by use of a one degree of freedom spring damper system. 2-D as well as 3-D analysis is performed using air as the fluid. Reynolds number is varied from 40 to 16000 approx., covering the laminar and turbulent regimes of flow. The experimental results of (Khalak and Williamson, 1997) and other researchers are used for validation purposes. The results obtained are comparable.展开更多
基金supported by National Outstanding Young Scientists Founds of China (Grant No. 50825902)Jiangsu Provincial Innovative Scholars "Climbing" Project of China (Grant No. BK 2009006)+1 种基金National Natural Science Foundation of China (Grant No. 50979034)Jiangsu Provincial Project for Innovative Postgraduates of China (Grant No. CX10B_262Z)
文摘The existing research for unsteady flow field and the corresponding flow induced vibration analysis of centrifugal pump are mainly carried out respectively without considering the interaction between fluid and structure. The ignorance of fluid structure interaction (FSI) means that the energy transfer between fluid and structure is neglected. To some extent, the accuracy and reliability of unsteady flow and rotor deflection analysis should be affected by this interaction mechanism. In this paper, a combined calculation between two executables for turbulent flow and vibrating structure was established using two-way coupling method to study the effect of FSI. Pressure distributions, radial forces, rotor deflection and equivalent stress are analyzed. The results show that the FSI effect to pressure distribution in flow field is complex. The pressure distribution is affected not only around impeller outlet where different variation trends of pressure values with and without FSI appear according to different relative positions between blade and cutwater, but also in the diffusion section of volute. Variation trends of peak values of radial force amplitude calculated with and without FSI are nearly same under high flow rate and designed conditions while the peak value with FSI is slightly smaller, and differently, the peak value with FSI is larger with low flow rate. In addition, the effect of FSI on the angle of radial force is quite complex, especially under 0.5Q condition. Fluctuation of radial deflection of the rotor has obvious four periods, of which the extent is relatively small under design condition and is relatively large under off-design condition. Finally, fluctuations of equivalent stress with time are obvious under different conditions, and stress value is small. The proposed research establishes the FSI calculation method for centrifugal pump analysis, and ensures the existing affect by fluid structure interaction.
文摘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.
文摘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.
基金supported by the National Basic Research Program of China (‘‘973’’ Program, No. 2013CB035906)the National Natural Science Foundation of China (No. 51439005)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (No. 51621092)
基金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.
基金Supported by the National Natural Science Foundation of China (21376016).
文摘The interaction between fluid and a down-pumping pitched blade turbine fixed with a flexible shaft in the stirred vessel, as a typical fluid structure interaction phenomenon, was simulated by coupling the Computational Fluid Dynamics and Computational Structural Dynamics. Based on the verification of the simulated impeller torque and dimensionless shaft bending moment with experimental result, the dimensionless shaft bending moment and various loads acting on impeller(including lateral force, axial force and bending moment) were discussed in detail. By separating and extracting the fluid and structural components from those loads, the results show that the shaft bending moment mainly results from the lateral force on impeller although the axial force on impeller is much larger. The impeller mass imbalance increases the shaft bending moment and the lateral force on impeller, but has little influence on the axial force and bending moment acting on impeller. The dominant frequencies of impeller forces are macro-frequency, speed frequency and blade passing frequency, and are associated with the impeller mass imbalance.
基金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.
文摘A blockage of blood vessels resulting from thrombus or plaque deposit causes serious cardiovascular diseases. This study developed a computational model of blood flow and drug transport to investigate the effectiveness of drug delivery to the stenotic sites. A three-dimensional (3D) model of the curved stenotic right coronary artery (RCA) was reconstructed based on the clinical angiogram image. Then, blood flow and drug transport with the flexible RCA wall were simulated using the fluid structure interaction (FSI) analysis and compared with the rigid RCA wall. Results showed that the maximal total displacement and von Mises stress of the flexible RCA model are 2.14 mm and 92.06 kPa. In addition, the effective injecting time point for the best performance of drug delivery was found to be between 0 s and 0.15 s (i.e., the fluid acceleration region) for both rigid and flexible RCA models. However, there was no notable difference in the ratio of particle deposition to the stenotic areas between the rigid and flexible RCA models. This study will be significantly useful to the design of a drug delivery system for the treatment of the stenotic arteries by targeting drugs selectively to the stenotic sites.
基金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 Natural Science Foundation of China(Nos.11972277 and 11872060)。
文摘A loss-of-coolant accident(LOCA)is one of the basic design considerations for nuclear reactor safety analysis.A LOCA induces propagation of a depressurization wave in the coolant,exerting hydrodynamic forces on structures viafiuid–structure interaction(FSI).The analysis of hydrodynamic forces on the core structures during a LOCA process is indispensable.We describe the implementation of a numerical strategy for prestressed structures.It consists of an initialization and a restarted transient analysis process,all implemented via the ANSYS Workbench by system coupling of ANSYS and Fluent.Our strategy is validated by making extensive comparisons of the pressures,displacements,and strains on various locations between the simulation and reported measurements.The approach is appealing for dynamic analysis of other prestressed structures,owing to the good popularity and acknowledgement of ANSYS and Fluent in both academia and industry.
文摘In this paper, the second order steady drift forces on the ships and other floating offshore structures are calculated by the far field method. The amplitudes of diffracted waves and radiated waves at infinity are obtained by the two-dimensional source distribution and strip-theory method. For the twin hull structure, the hydrodynamic interaction between the two hulls is taken into account. The drift forces on cross sections of Lewis type as well as on the semi-submersibles are computed. The theoretical results obtained by the present method agree fairly well with experimental results.
文摘Conventional coupled BE/FE (Boundary-Element/Finite-Element) method and modeling of structural-acoustic interaction has shown its promise and potential in the design and analysis of various structural-acoustic interaction applications. Unified combined acoustic and aerodynamic loading on the structure is synthesized using two approaches. Firstly, by linear superposition of the acoustic pressure disturbance to the aeroelastic problem, the effect of acoustic pressure disturbance to the aeroelastic structure is considered to consist of structural motion independent incident acoustic pressure and structural motion dependent acoustic pressure, which is known as the scattering pressure, referred here as the acoustic aerodynamic analogy. Secondly, by synthesizing the acoustic and aerodynamic effects on elastic structure using an elegant, effective and unified approach, both acoustic and aerodynamic effect on solid structural boundaries can be formulated as a boundary value problem governed by second order differential equations which lead to solutions expressible as surface integral equations. The unified formulation of the acousto-aeroelastic problem is amenable for simultaneous solution, although certain prevailing situations allow the solution of the equations independently. For this purpose, the unsteady aerodynamic problem which was earlier utilizes well-established lifting surface method is reformulated using Boundary Element (BE) approach. These schemes are outlined and worked out with examples.
文摘This paper develops a numerical code for modelling liquid sloshing.The coupled boundary element-finite element method was used to solve the Laplace equation for inviscid fluid and nonlinear free surface boundary conditions.Using Nakayama and Washizu’s results,the code performance was validated.Using the developed numerical mode,we proposed artificial neural network(ANN)and genetic algorithm(GA)methods for evaluating sloshing loads and comparing them.To compare the efficiency of the suggested methods,the maximum free surface displacement and the maximum horizontal force exerted on a rectangular tank’s perimeter are examined.It can be seen from the results that both ANNs and GAs can accurately predict η_(max) and F_(max).
基金This study was co-supported by the National Natural Science Foundation of China(No.52275293)the Guangdong Basic and Applied Basic Research Foundation,China(No.2023A1515010774)+1 种基金the Basic Research Program of Shenzhen,China(No.JCYJ 20190806142816524)the National Key Laboratory of Science and Technology on Aerodynamic Design and Research,China(No.61422010301).
文摘Certain insect species have been observed to exploit the resonance mechanism of their wings.In order to achieve resonance and optimize aerodynamic performance,the conventional approach is to set the flapping frequency of flexible wings based on the Traditional Structural Modal(TSM)analysis.However,there exists controversy among researchers regarding the relationship between frequency and aerodynamic performance.Recognizing that the structural response of wings can be influenced by the surrounding air vibrations,an analysis known as Acoustic Structure Interaction Modal(ASIM)is introduced to calculate the resonant frequency.In this study,Fluid Structure Interaction(FSI)simulations are employed to investigate the aerodynamic performance of flapping wings at modal frequencies derived from both TSM and ASIM analyses.The performance is evaluated for various mass ratios and frequency ratios,and the findings indicate that the deformation and changes in vortex structure exhibit similarities at mass ratios that yield the highest aerodynamic performance.Notably,the flapping frequency associated with the maximum time-averaged vertical force coefficient at each mass ratio closely aligns with the ASIM frequency,as does the frequency corresponding to maximum efficiency.Thus,the ASIM analysis can provide an effective means for predicting the optimal flapping frequency for flexible wings.Furthermore,it enables the prediction that flexible wings with varying mass ratios will exhibit similar deformation and vortex structure changes.This paper offers a fresh perspective on the ongoing debate concerning the resonance mechanism of Flexible Flapping Wings(FFWs)and proposes an effective methodology for predicting their aerodynamic performance.
基金Supported by the Ph.D.Program Foundation of Ministry of Education of China (20070699054)~~
文摘The data information transfer and time marching strategies between computational fluid dynamics (CFD) and computational structural dynamics (CSD) play crucial roles on the aeroelastic analysis in a time domain. An improved CFD/CSD coupled system is designed, including an interpolation method and an improved loosely coupled algorithm. The interpolation method based on boundary element method (BEM) is developed to transfer aerodynamic loads and structural displacements between CFD and CSD grid systems, it can be universally used in fluid structural interaction solution by keeping energy conservation. The improved loosely coupled algo-rithm is designed, thus it improves the computational accuracy and efficiency. The new interface is performed on the two-dimensional (2-D) extrapolation and the aeroelastie response of AGARD445.6 wing. Results show that the improved interface has a superior accuracy.
文摘Experimental study and theoretical analysis show that the critical value of relative wave height (H / d)b given by Goda and the critical wave steepness (H / L)b given by Michell and Miche can be adopted as the spilling breaking indices of regular waves. According to the same principle, a systematic theoretical analysis and experiment of irregular wave have been done by the authors in order to solve the breaking problem of irregular waves. It is indicated that the authors' method for determining wave breaking of regular waves can also be used for irregular waves.
文摘Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model(moving particle semi-implicit(MPS)) is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model(MPS or HMPS) in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model(Enhanced MPS), thus, the developed coupled solvers for fluid structure interaction(FSI) are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics(SPH)-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH(ISPH) fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.
基金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).
文摘In this paper, vortex-induced vibrations of a cylinder are simulated by use of ANSYS CFX simulation code. The cylinder is treated as a rigid body and transverse displacements are obtained by use of a one degree of freedom spring damper system. 2-D as well as 3-D analysis is performed using air as the fluid. Reynolds number is varied from 40 to 16000 approx., covering the laminar and turbulent regimes of flow. The experimental results of (Khalak and Williamson, 1997) and other researchers are used for validation purposes. The results obtained are comparable.