Thermal Fluid-Solid Interaction Model and Experimental Validation for Hydrostatic Mechanical Face Seals

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.

Fluid-solid Interaction Model for Hydraulic Reciprocating O-ring Seals

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.

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

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.

Transport diffuse interface model for simulation of solid-fluid interaction

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.

Simulation on Dynamic Bending Features of Fabric Based on Fluid-Solid Interaction Technique

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.

Simulation and interaction of fluid and solid dynamics

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.

FLUID-SOLID COUPLING MATHEMATICAL MODEL OF CONTAMINANT TRANSPORT IN UNSATURATED ZONE AND ITS ASYMPTOTICAL SOLUTION

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.

A New Cadmium(Ⅱ) Coordination Polymer Extended through Hydrogen Bonds and π-π Stacking Interactions: Synthesis and Photoluminescence Property

A new coordination polymer, {[Cd(OPY)(tdc)(HO)]·H2 O}n(OPY = 4,4?-(oxybis(4,1-phenylene))dipyridine, H2 tdc = thiophene-2,5-dicarboxylic acid), has been synthesized hydrothermally based on a V-shaped ligand OPY. The structure was fully characterized by elemental analysis, FT-IR spectroscopy, and X-ray single-crystal diffraction analysis. In1, two OPY ligands and one water molecule acted as terminal ligands coordinating to Cdcation to form [Cd(OPY)HO]units, which are then linked by tdc2-ligands to generate a one-dimensional chain. Every two adjacent chains linked by extensive O–H···O hydrogen bonds constitute one-dimensional double-chains, and such chains are extended into two-dimensional layers via O–H···N hydrogen bonds. These layers are further connected to form a three-dimensional supramolecular architecture via π-π stacking interactions. In addition, the thermal stability and solid state fluorescence property of 1 were also investigated.

Modeling the dynamic process of tsunami earthquake by liquid-solid coupling model

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.

基于CFD-DEM的流-固耦合数值建模方法研究进展

随着土体渗流侵蚀研究的逐渐深入,对土颗粒流失和变形破坏机理的研究方法呈现出多尺度的特点。其中,计算流体力学-离散元耦合方法(CFD-DEM)为在细观尺度上研究流-固耦合相互作用对土体宏观力学特性的影响提供了一种行之有效的方法。针对CFD-DEM耦合方法在岩土工程领域应用现状,本文系统总结现有流-固耦合计算方法的优缺点,重点论述CFD-DEM耦合方法的建模策略,包括固相颗粒形状建模与粒间接触模型、流体相控制方程及参数计算方法,以及CFD-DEM耦合计算,并就相关问题进行深入探讨,最后提出了CFD-DEM耦合方法未来的发展方向。

水介质中冰-桥墩碰撞动力响应分析

为准确模拟冰-桥墩碰撞过程中墩体结构的动力响应,基于流固耦合(fluid-solid interaction, FSI)的计算方法,运用ANSYS/LS-DYNA软件数值模拟了水介质中冰排与桥墩在不同碰撞影响参数下的动力响应。结果表明:恒定附加质量(constant added mass, CAM)模型相比FSI模型冰载荷计算结果偏大,是因为CAM模型忽略了降低冰排速度的“水垫效应”,但其计算效率高,计算时间是FSI模型的1/10,更有利于桥墩防撞设计,而FSI模型能够更真实地模拟冰-桥墩碰撞场景。综合分析该研究所模拟的不同冰排工况,发现桥墩在冰排撞击作用下呈现明显的冰激结构振动特征,桥墩侧向位移云图随着应力波的传递呈现明显的层状分布,其顶部侧向位移幅度最大,因此在工程设计时应充分考虑其对桥梁整体安全的影响;随着冰排速度、冰排厚度和冰排压缩强度的增大,冰载荷平均值均呈现近似多项式函数关系。此外,研究发现当冰排压缩强度小于2.668 MPa(环境温度高于-15℃)时,冰载荷平均值明显减小,表明环境温度的变化对于桥墩受冰排撞击挤压作用所受的冰害问题具有显著的差异。其仿真研究成果拟为桥梁及桥墩防撞设施设计提供参考。

变含水量泥流本构及泥流-桥墩冲刷的流固耦合分析

泥石流是山区主要的地质灾害之一,对山区桥梁破坏严重.在流动过程中,泥石流受到含水量的影响较大,其本构模型会发生变化.本文对泥石流的一种特殊形式——泥流的质量含水量进行分析,将含水量引入泥浆的Herschel-Bulkley模型中,提出质量含水量及剪切速率同时变化下的泥流材料本构模型,并与成都黏土泥浆流变实验对比,验证该模型的正确性.而后基于两相流体的Navier-Stokes方程,将该本构模型应用于三维泥流-桥墩立柱冲击的流固耦合分析中,通过与相关规范及文献冲击压力经验公式计算结果对比,证实了模型的实用性.

CO-DIMENSION 2 BIFURCATIONS AND CHAOS IN CANTILEVERED PIPE CONVEYING TIME VARYING FLUID WITH THREE-TO-ONE INTERNAL RESONANCES

The nonlinear behavior of a cantilevered fluid conveying pipe subjected to principal parametric and internal resonances is investigated in this paper.The flow velocity is divided into constant and sinusoidal parts.The velocity value of the constant part is so adjusted such that the system exhibits 3:1 internal resonances for the first two modes.The method of multiple scales is employed to obtain the response of the system and a set of four first-order nonlinear ordinary- differential equations for governing the amplitude of the response.The eigenvalues of the Jacobian matrix are used to assess the stability of the equilibrium solutions with varying parameters.The co- dimension 2 derived from the double-zero eigenvalues is analyzed in detail.The results show that the response amplitude may undergo saddle-node,pitchfork,Hopf,homoclinic loop and period- doubling bifurcations depending on the frequency and amplitude of the sinusoidal flow.When the frequency of the sinusoidal flow equals exactly half of the first-mode frequency,the system has a route to chaos by period-doubling bifurcation and then returns to a periodic motion as the amplitude of the sinusoidal flow increases.

Numerical Study on Plume Interaction Above An AlternatingDiffuser in Stagnant Water

The plume interaction above an alternating diffuser in stagnant water is studied with 3D Reynolds-averaged NavierStokes equations （RANS） combined with a buoyancy-extended κ-ε model. The steady three-dimensional turbulent flow and temperature fields are computed by use of the finite volume method on a non-uniform high resolution orthogonal grid. The numerical predictions demonstrate a generic flow pattern for different turbulent heated jet discharges： the buoyant jets on each side of the diffuser first merge to form an essentially two-dimensional plume which bends back toward the diffuser centerline due to a low pressure cavity. In general, an under-pressure exists in the cavity until the plumes merge; the pressure increases to slightly positive afterwards. Two-dimensionality of the scalar and flow field is attained much later than the point of zero pressure. The position of merging point is governed by mainly four parameters - the discharge densimetfic Froude number, the port diameter and space, and the horizontal distance between alternating jet nozzles. A formula from numerical simulations is obtained through regression analysis and it is used to predict the position of plume merging point. The predicted temperature fields are comparable to previous experiments.

Flux vector splitting solutions for coupling hydraulic transient of gas-liquid-solid three-phase flow in pipelines

The gas-liquid-solid three-phase mixed flow is the most general in multiphase mixed transportation. It is significant to exactly solve the coupling hydraulic transient problems of this type of multiphase mixed flow in pipelines. Presently, the method of characteristics is widely used to solve classical hydraulic transient problems. However, when it is used to solve coupling hydraulic transient problems, excessive interpolation errors may be introduced into the results due to unavoidable multiwave interpolated calculations. To deal with the problem, a finite difference scheme based on the Steger- Warming flux vector splitting is proposed. A flux vector splitting scheme is established for the coupling hydraulic transient model of gas-liquid-solid three-phase mixed flow in the pipelines. The flux subvectors are then discretized by the Lax-Wendroff central difference scheme and the Warming-Beam upwind difference scheme with second-order precision in both time and space. Under the Rankine-Hugoniot conditions and the corresponding boundary conditions, an effective solution to those points located at the boundaries is developed, which can avoid the problem beyond the calculation region directly induced by the second-order discrete technique. Numerical and experimental verifications indicate that the proposed scheme has several desirable advantages including high calculation precision, excellent shock wave capture capability without false numerical oscillation, low sensitivity to the Courant number, and good stability.

基于SPH-FEM耦合方法的泥石流冲击输电塔基础的动力分析

泥石流是我国西南山区常见的地质灾害。架空输电杆塔在泥石流的冲击下往往发生基础破坏甚至会造成杆塔倒塌。首先采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法和有限元方法(finite element method,简称FEM)相耦合的三维数值方法模拟了泥石流对杆塔基础的冲击作用;在与相关模型试验结果验证的基础上,开展了不同泥石流密度、黏度系数及初始速度条件下对输电塔基础的冲击力作用的参数分析;研究结果表明:随着泥石流初始速度的增加,冲击力峰值会随之增大;前排基础的冲击力峰值均大于后排基础;泥石流冲击过程特性受到泥石流密度和黏度系数影响。与稀性泥石流相比:黏性泥石流冲击基础后,基础下游真空区相对要小;此外,将数值模拟结果与Kwan冲击力公式及铁二院推荐的冲击压力设计公式预测值进行对比分析可以发现:Kwan冲击力公式能较好地预测出基础所受泥石流冲击力的平均趋势,最大预测误差低于30%,铁二院公式预测的稀性和黏性泥石流的冲击压力平均偏低分别约17%和28%。相关研究结果有望为泥石流频发区域输电塔基础的设计和风险评估提供一定的参考依据。

Fluid Reaction Force on Solids

Steady flow past a circular arc wing produces a vertical lift force, which has been calculated before, but it also causes two horizontal reaction forces, equal and opposite, due to Newton’s third law: the accelerated and decelerated flows act back on the top surfaces of the wing. Reaction forces are computed here and compared with the lift force. The magnitude of the two perpendicular forces are equal when the ratio of the maximum thickness of the wing is about 1/5 the length of the wing’s flat bottom surface. An example of an asymmetric wing is discussed and the net reaction force is calculated, which is always directed horizontally away from the top surface of the wing that has the greater mean slope. Based on these results,it is predicted that a surface ship should go more easily through the water if the bow were blunt and the stern pointed, just the opposite configuration of what is commonly found traversing lakes and oceans. Surface gravity waves are argued not to change this conclusion in general.

Simulation of dynamic fluid-solid interactions with an improved direct-forcing immersed boundary method

Dynamic fluid-solid interactions are widely found in chemical engineering, such as in particle-laden flows, which usually contain complex moving boundaries. The immersed boundary method （IBM） is a convenient approach to handle fluid-solid interactions with complex geometries. In this work, Uhlmann＇s direct-forcing IBM is improved and implemented on a supercomputer with CPU-GPU hybrid architec- ture. The direct-forcing IBM is modified as follows： the Poisson＇s equation for pressure is solved before evaluation of the body force, and the force is only distributed to the Cartesian grids inside the immersed boundary. A multidirect forcing scheme is used to evaluate the body force. These modifications result in a divergence-free flow field in the fluid domain and the no-slip boundary condition at the immersed boundary simultaneously. This method is implemented in an explicit finite-difference fractional-step scheme, and validated by 2D simulations of lid-driven cavity flow, Couette flow between two concentric cylinders and flow over a circular cylinder. Finally, the method is used to simulate the sedimentation of two circular particles in a channel. The results agree very well with previous experimental and numerical data, and are more accurate than the conventional direct-forcing method, especially in the vicinity of a moving boundary.

Validation of fluid-particle interaction force relationships in binary-solid suspensions

In this work several relationships governing solid-fluid dynamic interaction forces were validated against experimental data for a single particle settling in a suspension of other smaller particles. It was observed that force relationships based on Lattice-Boltzmann simulations did not perform as well as other inter- action types tested. Nonetheless, it is apparent that, in the case of a suspension of different particle types, it is important that the correct choice is made as to how the contribution to the overall fluid-particle interaction force is split between buoyancy and drag. Experimental evidence clearly suggests that the ＂generalized＂ Archimedes＇ principle （where the foreign particle is considered to displace the whole suspension and not just the fluid） provides the best result.