It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on b...It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on bubble motion and gas-liquid mass transfer process in a bubble column.Particle image velocimetry(PIV)technology was employed to analyze the flow field and bubble motion behavior,and dynamic oxygen absorption technology was used to measure the gas-liquid volumetric mass transfer coefficient(kLa).The effect of embedding wire mesh,adding solid particles,and wire mesh coupling solid particles on the flow characteristic and kLa were analyzed and compared.The results show that the gas-liquid interface area increases by 33%-72%when using the wire mesh coupling solid particles strategy compared to the gas-liquid two-phase flow,which is superior to the other two strengthening methods.Compared with the system without reinforcement,kLa in the bubble column increased by 0.5-1.8 times with wire mesh coupling solid particles method,which is higher than the sum of kLa increases with inserting wire mesh and adding particles,and the coupling reinforcement mechanism for affecting gas-liquid mass transfer process was discussed to provide a new idea for enhancing gas-liquid mass transfer.展开更多
Hydrostatic mechanical face seals for reactor coolant pumps are very important for the safety and reliability of pressurized-water reactor power plants.More accurate models on the operating mechanism of the seals are ...Hydrostatic mechanical face seals for reactor coolant pumps are very important for the safety and reliability of pressurized-water reactor power plants.More accurate models on the operating mechanism of the seals are needed to help improve their performance.The thermal fluid–solid interaction（TFSI）mechanism of the hydrostatic seal is investigated in this study.Numerical models of the flow field and seal assembly are developed.Based on the mechanism for the continuity condition of the physical quantities at the fluid–solid interface,an on-line numerical TFSI model for the hydrostatic mechanical seal is proposed using an iterative coupling method.Dynamic mesh technology is adopted to adapt to the changing boundary shape.Experiments were performed on a test rig using a full-size test seal to obtain the leakage rate as a function of the differential pressure.The effectiveness and accuracy of the TFSI model were verified by comparing the simulation results and experimental data.Using the TFSI model,the behavior of the seal is presented,including mechanical and thermal deformation,and the temperature field.The influences of the rotating speed and differential pressure of the sealing device on the temperature field,which occur widely in the actual use of the seal,are studied.This research proposes an on-line and assembly-based TFSI model for hydrostatic mechanical face seals,and the model is validated by full-sized experiments.展开更多
The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contami...The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.展开更多
The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo val...The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo valve would reach 120℃ and the valve core and valve sleeve deform in a short amount of time. So the control precision of servo valve significantly decreases and the clamping stagnation phenomenon of valve core appears. In order to solve the problem of degraded control accuracy and clamping stagnation of servo valve under large temperature difference circumstance, the numerical simulation of heat-fluid-solid coupling by using finite element method is done. The simulation result shows that zero position leakage of servo valve is basically impacted by oil temperature and change of fit clearance. The clamping stagnation is caused by warpage-deformation and fit clearance reduction of the valve core and valve sleeve. The distribution roles of the temperature and thermal-deformation of shell, valve core and valve sleeve and the pressure, velocity and temperature field of flow channel are also analyzed. Zero position leakage and electromagnet＇s current when valve core moves in full-stroke are tested using Electro-hydraulic Servo-valve Characteristic Test-bed of an aerospace sciences and technology corporation. The experimental results show that the change law of experimental current at different oil temperatures is roughly identical to simulation current. The current curve of the electromagnet is smooth when oil temperature is below 80℃, but the amplitude of current significantly increases and the hairy appears when oil temperature is above 80℃. The current becomes smooth again after the warped valve core and valve sleeve are reground. It indicates that clamping stagnation is caused by warpage-deformation and fit clearance reduction of valve core and valve sleeve. This paper simulates and tests the heat-fluid-solid coupling of double flapper-nozzle servo valve, and the obtained results provide the reference value for the design of double flapper-nozzle force feedback servo valve.展开更多
Fluid and solid simulation is to generate a realistic simulation of fluids and solids,in particular for the fluids such as water and smoke,with computation of Euler equations or Navier-Stokes equations conducted to go...Fluid and solid simulation is to generate a realistic simulation of fluids and solids,in particular for the fluids such as water and smoke,with computation of Euler equations or Navier-Stokes equations conducted to govern the real fluid physics.Fluid simulation is an important field by its wide applications in many fields and industries,such as film and game simulation,weather forecasting,natural disaster simulation and protection,simulation in maritime and aviation.There are basically two main categories of methods for fluid simulation,data-driven methods and physically-based methods.The data-driven models establish a direct mapping between variables and extract their relationship from historically measured data by the algorithms developed in the fields of statistics,computational intelligence,machine learning,and data mining.展开更多
Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. On...Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. Only elastic deformations of hydraulic reciprocating seals were discussed, and hydrodynamic effects were neglected in many studies. The physical process of the fluid-solid interaction effect did not be clearly presented in the existing fluid-solid interaction models for hydraulic reciprocating O-ring seals, and few of these models had been simultaneously validated through experiments. By exploring the physical process of the fluid-solid interaction effect of the hydraulic reciprocating O-ring seal, a numerical fluid-solid interaction model consisting of fluid lubrication, contact mechanics, asperity contact and elastic deformation analyses is constructed with an iterative procedure. With the SRV friction and wear tester, the experiments are performed to investigate the elastohydrodynamic lubrication characteristics of the O-ring seal. The regularity of the friction coefficient varying with the speed of reciprocating motion is obtained in the mixed lubrication condition. The experimental result is used to validate the fluid-solid interaction model. Based on the model, The elastohydrodynamic lubrication characteristics of the hydraulic reciprocating O-ring seal are presented respectively in the dry friction, mixed lubrication and full film lubrication conditions, including of the contact pressure, film thickness, friction coefficient, liquid film pressure and viscous shear stress in the sealing zone. The proposed numerical fluid-solid interaction model can be effectively used to analyze the operation characteristics of the hydraulic reciprocating O-ring seal, and can also be widely used to study other hydraulic reciprocating seals.展开更多
As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and p...As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.展开更多
On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Sti...On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Stiffness matrix are obtained. Using QR method , pipe’s nature frequencies are calculated. The curves of the first four orders of natural frequency-flow velocity of pipe waw given .The influence of flowing velocity ,pressure, solid-liquid coupling damp and solid-liquid coupling stiffness on natural frequency are discussed respectively.The dynamic respondence of the pipes for stepload with different flow velocity are calculated by Newmark method .It is found that,with the flow velocity increased, the nature frequency of the pipes reduced, increased,reduced again and so on.展开更多
An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal str...An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.展开更多
For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. ...For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.展开更多
This paper is devoted to the two-dimensional nonlinear modeling of the fluid-solid interaction (FSI) between fabric and air flow, which is based on the Automatic Incremental Dynamic Nonlinear Analysis (AIDNA)-FSI prog...This paper is devoted to the two-dimensional nonlinear modeling of the fluid-solid interaction (FSI) between fabric and air flow, which is based on the Automatic Incremental Dynamic Nonlinear Analysis (AIDNA)-FSI program in order to study the dynamic bending features of fabrics in a specific air flow filed. The computational fluid dynamics (CFD) model for flow and the finite element model (FEM) for fabric was set up to constitute an FSI model in which the geometric nonlinear behavior and the dynamic stress-strain variation of the relatively soft fabric material were taken into account. Several FSI cases with different time-dependent wind load and the model frequency analysis for fabric were carried out. The dynamic response of fabric and the distribution of fluid variables were investigated. The results of numerical simulation and experiments fit quite well. Hence, this work contributes to the research of modeling the dynamic bending behavior of fabrics in air field.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between ...The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.展开更多
Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by ...Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.展开更多
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.展开更多
Currently, when magnesium alloy sheet is rolled, the method of controlling roll temperature is simple and inaccurate. Furthermore, roll temperature has a large influence on the quality of magnesium alloy sheet; theref...Currently, when magnesium alloy sheet is rolled, the method of controlling roll temperature is simple and inaccurate. Furthermore, roll temperature has a large influence on the quality of magnesium alloy sheet; therefore, a new model using circular fluid flow control roll temperature has been designed. A fluid heat transfer structure was designed, the heat transfer process model of the fluid heating roll was simplified, and the finite di erence method was used to cal?culate the heat transfer process. Fluent software was used to simulate the fluid?solid coupling heat transfer, and both the trend and regularity of the temperature field in the heat transfer process were identified. The results show that the heating e ciency was much higher than traditional heating methods(when the fluid heat of the roll and tempera?ture distribution of the roll surface was more uniform). Moreover, there was a bigger temperature di erence between the input and the output, and after using reverse flow the temperature di erence decreased. The axial and circum?ferential temperature distributions along the sheet were uniform. Both theoretical calculation results and numerical simulation results of the heat transfer between fluid and roll were compared. The error was 1.8%–12.3%, showing that the theoretical model can both forecast and regulate the temperature of the roll(for magnesium alloy sheets) in the rolling process.展开更多
Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenome...Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.展开更多
The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional can...The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.展开更多
In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The el...In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The electric current density, temperature and velocity distributions in the torch are obtained through the coupled iterative calculation about the electromagnetic equations described in a magnetic vector potential format and the modified fluid dynamics equations. The fluid-solid coupled calculation method is applied to guarantee the continuity of the electric current and heat transfer at the interface between the electrodes and fluid. The predicted location of the anodic arc root attachment and the arc voltage of the torch are consistent with corresponding experimental results. Through a specific analysis of the influence of mass flow rates and electric current on the torch outlet parameters, the total thermal efficiency, thermal loss of each part, and the laws of the variation of outlet parameters with the variation of mass flow rates and electric current was obtained. It is found that operation under a laminar condition with a limited area of the anode could increase the total thermal efficiency of the torch.展开更多
基金supported by the Key Research and Development Plan of Shandong Province(the Major Scientific and Technological Innovation Projects,2021ZDSYS13)the Natural Science Foundation of Shandong Province(ZR2021MB135)Natural Science Foundation of Shandong Province(ZR2021ME224).
文摘It is of vital significance to investigate mass transfer enhancements for chemical engineering processes.This work focuses on investigating the coupling influence of embedding wire mesh and adding solid particles on bubble motion and gas-liquid mass transfer process in a bubble column.Particle image velocimetry(PIV)technology was employed to analyze the flow field and bubble motion behavior,and dynamic oxygen absorption technology was used to measure the gas-liquid volumetric mass transfer coefficient(kLa).The effect of embedding wire mesh,adding solid particles,and wire mesh coupling solid particles on the flow characteristic and kLa were analyzed and compared.The results show that the gas-liquid interface area increases by 33%-72%when using the wire mesh coupling solid particles strategy compared to the gas-liquid two-phase flow,which is superior to the other two strengthening methods.Compared with the system without reinforcement,kLa in the bubble column increased by 0.5-1.8 times with wire mesh coupling solid particles method,which is higher than the sum of kLa increases with inserting wire mesh and adding particles,and the coupling reinforcement mechanism for affecting gas-liquid mass transfer process was discussed to provide a new idea for enhancing gas-liquid mass transfer.
基金Supported by National Basic Research Program of China(973 Program,Grant No.2009CB724304)National Key Technology R&D Program(Grant No.2011BAF09B05)National Natural Science Foundation of China(Grant No.50975157)
文摘Hydrostatic mechanical face seals for reactor coolant pumps are very important for the safety and reliability of pressurized-water reactor power plants.More accurate models on the operating mechanism of the seals are needed to help improve their performance.The thermal fluid–solid interaction（TFSI）mechanism of the hydrostatic seal is investigated in this study.Numerical models of the flow field and seal assembly are developed.Based on the mechanism for the continuity condition of the physical quantities at the fluid–solid interface,an on-line numerical TFSI model for the hydrostatic mechanical seal is proposed using an iterative coupling method.Dynamic mesh technology is adopted to adapt to the changing boundary shape.Experiments were performed on a test rig using a full-size test seal to obtain the leakage rate as a function of the differential pressure.The effectiveness and accuracy of the TFSI model were verified by comparing the simulation results and experimental data.Using the TFSI model,the behavior of the seal is presented,including mechanical and thermal deformation,and the temperature field.The influences of the rotating speed and differential pressure of the sealing device on the temperature field,which occur widely in the actual use of the seal,are studied.This research proposes an on-line and assembly-based TFSI model for hydrostatic mechanical face seals,and the model is validated by full-sized experiments.
文摘The process of contaminant transport is a problem of multicomponent and multiphase flow in unsaturated zone. Under the presupposition that gas existence affects water transport, a coupled mathematical model of contaminant transport in unsaturated zone has been established based on fluid_solid interaction mechanics theory. The asymptotical solutions to the nonlinear coupling mathematical model were accomplished by the perturbation and integral transformation method. The distribution law of pore pressure, pore water velocity and contaminant concentration in unsaturated zone has been presented under the conditions of with coupling and without coupling gas phase. An example problem was used to provide a quantitative verification and validation of the model. The asymptotical solution was compared with Faust model solution. The comparison results show reasonable agreement between asymptotical solution and Faust solution, and the gas effect and media deformation has a large impact on the contaminant transport. The theoretical basis is provided for forecasting contaminant transport and the determination of the relationship among pressure_saturation_permeability in laboratory.
基金Supposed by National Natural Science Foundation of China(Grant No.51075348)Hebei Provincial Natural Science Foundation of China(Grant No.E2011203151)Research Fund for Doctoral Program of Higher Education of China(Grant No.20101333110002)
文摘The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo valve would reach 120℃ and the valve core and valve sleeve deform in a short amount of time. So the control precision of servo valve significantly decreases and the clamping stagnation phenomenon of valve core appears. In order to solve the problem of degraded control accuracy and clamping stagnation of servo valve under large temperature difference circumstance, the numerical simulation of heat-fluid-solid coupling by using finite element method is done. The simulation result shows that zero position leakage of servo valve is basically impacted by oil temperature and change of fit clearance. The clamping stagnation is caused by warpage-deformation and fit clearance reduction of the valve core and valve sleeve. The distribution roles of the temperature and thermal-deformation of shell, valve core and valve sleeve and the pressure, velocity and temperature field of flow channel are also analyzed. Zero position leakage and electromagnet＇s current when valve core moves in full-stroke are tested using Electro-hydraulic Servo-valve Characteristic Test-bed of an aerospace sciences and technology corporation. The experimental results show that the change law of experimental current at different oil temperatures is roughly identical to simulation current. The current curve of the electromagnet is smooth when oil temperature is below 80℃, but the amplitude of current significantly increases and the hairy appears when oil temperature is above 80℃. The current becomes smooth again after the warped valve core and valve sleeve are reground. It indicates that clamping stagnation is caused by warpage-deformation and fit clearance reduction of valve core and valve sleeve. This paper simulates and tests the heat-fluid-solid coupling of double flapper-nozzle servo valve, and the obtained results provide the reference value for the design of double flapper-nozzle force feedback servo valve.
文摘Fluid and solid simulation is to generate a realistic simulation of fluids and solids,in particular for the fluids such as water and smoke,with computation of Euler equations or Navier-Stokes equations conducted to govern the real fluid physics.Fluid simulation is an important field by its wide applications in many fields and industries,such as film and game simulation,weather forecasting,natural disaster simulation and protection,simulation in maritime and aviation.There are basically two main categories of methods for fluid simulation,data-driven methods and physically-based methods.The data-driven models establish a direct mapping between variables and extract their relationship from historically measured data by the algorithms developed in the fields of statistics,computational intelligence,machine learning,and data mining.
基金supported by National Basic Research Program of China(973 Program,Grant No. 2009CB724304)Key Research Program of the State Key Laboratory of Tribology of Tsinghua University,China (Grant No. SKLT08A06)National Natural Science Foundation of China(Grant No. 50975157)
文摘Elastohydrodynamic lubrication characteristics of hydraulic reciprocating seals have significant effects on sealing and tribology performances of hydraulic actuators, especially in high parameter hydraulic systems. Only elastic deformations of hydraulic reciprocating seals were discussed, and hydrodynamic effects were neglected in many studies. The physical process of the fluid-solid interaction effect did not be clearly presented in the existing fluid-solid interaction models for hydraulic reciprocating O-ring seals, and few of these models had been simultaneously validated through experiments. By exploring the physical process of the fluid-solid interaction effect of the hydraulic reciprocating O-ring seal, a numerical fluid-solid interaction model consisting of fluid lubrication, contact mechanics, asperity contact and elastic deformation analyses is constructed with an iterative procedure. With the SRV friction and wear tester, the experiments are performed to investigate the elastohydrodynamic lubrication characteristics of the O-ring seal. The regularity of the friction coefficient varying with the speed of reciprocating motion is obtained in the mixed lubrication condition. The experimental result is used to validate the fluid-solid interaction model. Based on the model, The elastohydrodynamic lubrication characteristics of the hydraulic reciprocating O-ring seal are presented respectively in the dry friction, mixed lubrication and full film lubrication conditions, including of the contact pressure, film thickness, friction coefficient, liquid film pressure and viscous shear stress in the sealing zone. The proposed numerical fluid-solid interaction model can be effectively used to analyze the operation characteristics of the hydraulic reciprocating O-ring seal, and can also be widely used to study other hydraulic reciprocating seals.
基金supported by the Major National Science and Technology Program(Nos.2008ZX05026-00411 and 2011ZX05026-004-08)the Program for Changjiang Scholars and Innovative Research Team in University(No.RT1086)
文摘As the oil or gas exploration and development activities in deep and ultra- deep waters become more and more, encountering gas hydrate bearing sediments （HBS） is almost inevitable. The variation in temperature and pressure can destabilize gas hydrate in nearby formation around the borehole, which may reduce the strength of the formation and result in wellbore instability. A non-isothermal, transient, two-phase, and fluid-solid coupling mathematical model is proposed to simulate the complex stability performance of a wellbore drilled in HBS. In the model, the phase transition of hydrate dissociation, the heat exchange between drilling fluid and formation, the change of mechanical and petrophysical properties, the gas-water two-phase seepage, and its interaction with rock deformation are considered. A finite element simulator is developed, and the impact of drilling mud on wellbore instability in HBS is simulated. Results indicate that the re- duction in pressure and the increase in temperature of the drilling fluid can accelerate hydrate decomposition and lead to mechanical properties getting worse tremendously. The cohesion decreases by 25% when the hydrate totally dissociates in HBS. This easily causes the wellbore instability accordingly. In the first two hours after the formation is drilled, the regions of hydrate dissociation and wellbore instability extend quickly. Then, with the soaking time of drilling fluid increasing, the regions enlarge little. Choosing the low temperature drilling fluid and increasing the drilling mud pressure appropriately can benefit the wellbore stability of HBS. The established model turns out to be an efficient tool in numerical studies of the hydrate dissociation behavior and wellbore stability of HBS.
文摘On the basis of Hamilton principle. the equation of sonlid-liquid coupling vibration of pipe conveying fluid is deduced. An asymmetrical sonlid-liquid coupling damp matrix and a symmetrical solid-liquid coupling Stiffness matrix are obtained. Using QR method , pipe’s nature frequencies are calculated. The curves of the first four orders of natural frequency-flow velocity of pipe waw given .The influence of flowing velocity ,pressure, solid-liquid coupling damp and solid-liquid coupling stiffness on natural frequency are discussed respectively.The dynamic respondence of the pipes for stepload with different flow velocity are calculated by Newmark method .It is found that,with the flow velocity increased, the nature frequency of the pipes reduced, increased,reduced again and so on.
基金the National Metal and Materials Technology Centerthe Thailand Research Fund+1 种基金the Office of Higher Education Commissionthe Chulalongkorn University for supporting the present research
文摘An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin （SUPG） method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid axe performed by the Galerkin method. The second-order semi- implicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are lineaxized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.
基金Project supported by the National Natural Science Foundation of China(Nos.11702029,11771054,U1730118,91852207,and 11801036)the China Postdoctoral Science Foundation(No.2016M600967)
文摘For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.
基金National Natural Science Foundations of China(No.50803010,No.60904056)
文摘This paper is devoted to the two-dimensional nonlinear modeling of the fluid-solid interaction (FSI) between fabric and air flow, which is based on the Automatic Incremental Dynamic Nonlinear Analysis (AIDNA)-FSI program in order to study the dynamic bending features of fabrics in a specific air flow filed. The computational fluid dynamics (CFD) model for flow and the finite element model (FEM) for fabric was set up to constitute an FSI model in which the geometric nonlinear behavior and the dynamic stress-strain variation of the relatively soft fabric material were taken into account. Several FSI cases with different time-dependent wind load and the model frequency analysis for fabric were carried out. The dynamic response of fabric and the distribution of fluid variables were investigated. The results of numerical simulation and experiments fit quite well. Hence, this work contributes to the research of modeling the dynamic bending behavior of fabrics in air field.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
文摘The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.
基金National Natural Science Foundation of China (40521002 and 40474013).
文摘Tsunami induced by earthquake is an interaction problem between liquid and solid.Shallow-water wave equation is often used to modeling the tsunami,and the boundary or initial condition of the problem is determined by the displacement or velocity field from the earthquake under sea floor,usually no interaction between them is consid-ered in pure liquid model.In this study,the potential flow theory and the finite element method with the interaction between liquid and solid are employed to model the dynamic processes of the earthquake and tsunami.For model-ing the earthquake,firstly the initial stress field to generate the earthquake is set up,and then the occurrence of the earthquake is simulated by suddenly reducing the elastic material parameters inside the earthquake fault.It is dif-ferent from seismic dislocation theory in which the relative slip on the fault is specified in advance.The modeling results reveal that P,SP and the surface wave can be found at the sea surface besides the tsunami wave.The surface wave arrives at the distance of 600 km from the epicenter earlier than the tsunami 48 minutes,and its maximum amplitude is 0.55 m,which is 2 times as large as that of the sea floor.Tsunami warning information can be taken from the surface wave on the sea surface,which is much earlier than that obtained from the seismograph stations on land.The tsunami speed on the open sea with 3 km depth is 175.8 m/s,which is a little greater than that pre-dicted by long wave theory,(gh)1/2=171.5 m,and its wavelength and amplitude in average are 32 km and 2 m,respectively.After the tsunami propagates to the continental shelf,its speed and wavelength is reduced,but its amplitude become greater,especially,it can elevate up to 10 m and run 55 m forward in vertical and horizontal directions at sea shore,respectively.The maximum vertical accelerations at the epicenter on the sea surface and on the earthquake fault are 5.9 m/s2 and 16.5 m/s2,respectively,the later is 2.8 times the former,and therefore,sea water is a good shock absorber.The acceleration at the sea shore is about 1/10 as large as at the epicenter.The maximum vertical velocity at the epicenter is 1.4 times that on the fault.The maximum vertical displacement at the fault is less than that at the epicenter.The difference between them is the amplitude of the tsunami at the epicenter.The time of the maximum displacement to occur on the fault is not at the beginning of the fault slipping but retards 23 s.
基金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.
基金National Natural Science Foundation of China(Grant No.U1510131)Key Research and Development Projects of Shanxi Province of China(Grant Nos.201603D121010,201603D111004)+3 种基金Science and Technology Project of Jin Cheng City of China(Grant No.20155010)Youth Program of National Natural Science Fund of China(Grant No.51604181)Project of Young Scholar of Shanxi ProvinceLeading Talent Project of Innovative Entrepreneurial Team of Jiangsu Province(Grant No.51501122)
文摘Currently, when magnesium alloy sheet is rolled, the method of controlling roll temperature is simple and inaccurate. Furthermore, roll temperature has a large influence on the quality of magnesium alloy sheet; therefore, a new model using circular fluid flow control roll temperature has been designed. A fluid heat transfer structure was designed, the heat transfer process model of the fluid heating roll was simplified, and the finite di erence method was used to cal?culate the heat transfer process. Fluent software was used to simulate the fluid?solid coupling heat transfer, and both the trend and regularity of the temperature field in the heat transfer process were identified. The results show that the heating e ciency was much higher than traditional heating methods(when the fluid heat of the roll and tempera?ture distribution of the roll surface was more uniform). Moreover, there was a bigger temperature di erence between the input and the output, and after using reverse flow the temperature di erence decreased. The axial and circum?ferential temperature distributions along the sheet were uniform. Both theoretical calculation results and numerical simulation results of the heat transfer between fluid and roll were compared. The error was 1.8%–12.3%, showing that the theoretical model can both forecast and regulate the temperature of the roll(for magnesium alloy sheets) in the rolling process.
文摘Fluid-structure interaction （FSI） is essentially a dynamic phenomenon and always exists in fluid-filled pipe system. The four-equation model, which has been proved to be effective to describe and predict the phenomenon of FSI due to friction coupling and Poisson coupling being taken into account, is utilized to describe the FSI of fluid-filled pipe system. Terse compatibility equations are educed by the method of characteristics （MOC） to describe the fluid-filled pipe system. To shorten computing time needed to get the solutions under the condition of keeping accuracy requirement, two steps are adopted, firstly the time step Δt and divided number of the straight pipe are optimized, sec-ondly the mesh spacing Δz close to boundary is subdivided in several submeshes automatically ac-cording to the speed gradient of fluid. The mathematical model and arithmetic are validated by com-parisons between simulation solutions of two straight pipe systems and experiment known from lit-erature.
基金Project supported by the National Natural Science Foundation of China(No.10832007)the Shanghai Leading Academic Discipline Project(No.B206)
文摘The closely coupled approach combined with the finite volume method （FVM） solver and the finite element method （FEM） solver is used to investigate the fluid-structure interaction （FSI） of a three-dimensional cantilevered hydrofoil in the water tunnel. The FVM solver and the coupled approach are verified and validated by compar- ing the numerical predictions with the experimental measurements, and good agreement is obtained concerning both the lift on the foil and the tip displacement. In the noncav- itating flow, the result indicates that the growth of the initial incidence angle and the Reynolds number improves the deformation of the foil, and the lift on the foil is increased by the twist deformation. The normalized twist angle and displacement along the span of the hydrofoil for different incidence angles and Reynolds numbers are almost uniform. For the cavitation flow, it is shown that the small amplitude vibration of the foil has limited influence on the developing process of the partial cavity, and the quasi two-dimensional cavity shedding does not change the deformation mode of the hydrofoil. However, the frequency spectrum of the lift on the foil contains the frequency which is associated with the first bend frequency of the hydrofoil.
基金supported by National Natural Science Foundation of China (No. 50476081)
文摘In this work, a magnetic fluid dynamics （MHD） model is used to simulate the electromagnetic field, heat transfer and fluid flow in a DC non-transferred arc plasma torch under laminar and turbulent conditions. The electric current density, temperature and velocity distributions in the torch are obtained through the coupled iterative calculation about the electromagnetic equations described in a magnetic vector potential format and the modified fluid dynamics equations. The fluid-solid coupled calculation method is applied to guarantee the continuity of the electric current and heat transfer at the interface between the electrodes and fluid. The predicted location of the anodic arc root attachment and the arc voltage of the torch are consistent with corresponding experimental results. Through a specific analysis of the influence of mass flow rates and electric current on the torch outlet parameters, the total thermal efficiency, thermal loss of each part, and the laws of the variation of outlet parameters with the variation of mass flow rates and electric current was obtained. It is found that operation under a laminar condition with a limited area of the anode could increase the total thermal efficiency of the torch.
