Simulating turbulent Dean flow in Cartesian coordinates

A simplified approach to simulate turbulent flows in curved channels is proposed. A set of governing equations of motion in Cartesian coordinates is derived from the full Navier–Stokes equations in cylindrical coordinates. Terms to first order in the dimensionless curvature parameter are retained, whereas higher‐order terms are neglected. The curvature terms are implemented in a conventional Navier–Stokes code using Cartesian coordinates. Direct numerical simulations (DNS) of turbulent flow in weakly curved channels are performed. The pronounced asymmetries in the mean flow and the turbulence statistics observed in earlier DNS studies are faithfully reproduced by the present simplified Navier–Stokes model. It is particularly rewarding that also distinct pairs of counter‐rotating streamwise‐oriented vortices are embedded in the simulated flow field. Copyright © 2008 John Wiley & Sons, Ltd.     

低重环境下俯仰运动圆柱贮箱中液体非线性晃动

在低重力环境下,用变分原理建立了液体晃动的压力体积分形式的Lagrange函数,并将速度势函数在自由液面处作波高函数的级数展开,从而导出自由液面运动学和动力学边界条件非线性方程组;最后用四阶Runge-Kutta法求解非线性方程组。计算结果表明,随俯仰激励频率的逐渐变化,由于面外主模态和次生模态同时失稳,致使整个系统各阶模态和波高函数由稳态运动过渡为不稳定运动。     

A σ‐coordinate three‐dimensional numerical model for surface wave propagation

A three‐dimensional numerical model based on the full Navier–Stokes equations (NSE) in σ‐coordinate is developed in this study. The σ‐coordinate transformation is first introduced to map the irregular physical domain with the wavy free surface and uneven bottom to the regular computational domain with the shape of a rectangular prism. Using the chain rule of partial differentiation, a new set of governing equations is derived in the σ‐coordinate from the original NSE defined in the Cartesian coordinate. The operator splitting method (Li and Yu, Int. J. Num. Meth. Fluids 1996; 23 : 485–501), which splits the solution procedure into the advection, diffusion, and propagation steps, is used to solve the modified NSE. The model is first tested for mass and energy conservation as well as mesh convergence by using an example of water sloshing in a confined tank. Excellent agreements between numerical results and analytical solutions are obtained. The model is then used to simulate two‐ and three‐dimensional solitary waves propagating in constant depth. Very good agreements between numerical results and analytical solutions are obtained for both free surface displacements and velocities. Finally, a more realistic case of periodic wave train passing through a submerged breakwater is simulated. Comparisons between numerical results and experimental data are promising. The model is proven to be an accurate tool for consequent studies of wave‐structure interaction. Copyright © 2002 John Wiley & Sons, Ltd.     

俯仰激励下三维液体大幅晃动问题研究

主要研究俯仰激励下三维液体大幅晃动问题,将任意Lagrange-Euler法 (arbitrary Lagrange-Euler, ALE)运动学描述引入到Navier-Stokes方程中,推导了俯仰激励下液体大幅晃 动数值模拟计算公式,并利用Galerkin加权余量法推导了有限元数值离散方程和分步有限 元计算格式,采用ALE分步有限元方法对圆筒形贮腔中的液体大幅晃动进行了数值模拟计 算. 得到了波高、晃动力及晃动力矩等晃动特性的时间变化历程,并对结果进行了分析,揭 示了俯仰激励下三维液体大幅晃动问题的非线性现象.     

Non-linear multimodal model for tuned liquid dampers of arbitrary tank geometry

Tuned liquid dampers (TLDs) utilize sloshing fluid to absorb and dissipate structural vibrational energy. Simple TLD tank geometries may not always be feasible due to space limitations. While the non-linear modelling of sloshing fluid is currently limited to tanks of simple geometries, this paper develops a non-linear multimodal model which describes the sloshing behaviour of a fluid in a flat-bottom tank of arbitrary geometry. The mode shapes of the sloshing fluid are found by solving the Helmholtz equation over the tank domain using the finite element method. The Bateman-Luke variational principle is used to develop a system of ordinary differential equations which account for the coupling of the sloshing modes through the non-linear free surface boundary conditions. Damping is incorporated into the model by considering the drag produced on a set of damping screens inserted in the fluid. The system of ordinary differential equations is solved using the Runge-Kutta-Gill Method to predict the wave heights and sloshing forces. In general, the mode shapes in an arbitrary tank will have components in two orthogonal (x- and y-) directions. This out-of-plane behaviour is an important consideration for TLD design. The model is validated with existing models for the special cases of rectangular and circular tanks. Lastly, new shake table tests are conducted on a tank of complex geometry.     

Modelling and simulation of fires in vehicle tunnels

Applying a low‐Mach asymptotic for the compressible Navier–Stokes equations, we derive a new fluid dynamics model,which should be capable to model large temperature differences in combination with the low‐Mach number limit. The model is used to simulate fires in vehicle tunnels, where the standard Boussinesq‐approximation for the incompressible Navier–Stokes seems to be inappropriate due to the high temperatures developing in the tunnel. The model is implemented using a modified finite‐difference approach for the incompressible Navier–Stokes equations and tested in some realistic fire events. Copyright © 2004 John Wiley & Sons, Ltd.     

Formation of a secondary droplet over the crown upon the liquid film rupture under the action of a laser beam

The free surface of a liquid film exposed to a laser beam is deformed and suffers a rupture. Depending on the thermal load intensity and the thermal properties of the liquid the rupture can be accompanied by the formation of secondary droplets over the film crown. This process is investigated using a mathematical model describing the motion of the thin layer of a viscous nonisothermal liquid. The model is based on the two-dimensional Navier–Stokes equations. The boundary conditions at the film-gas and film-liquid interfaces necessary for the solution of these equations are derived in the explicit form. The results of the solution of model problems are presented.     

A simulation of free surface waves for incompressible two‐phase flows using a curvilinear level set formulation

A level set formulation in a generalized curvilinear coordinate is developed to simulate the free surface waves generated by moving bodies or the sloshing of fluid in a container. The Reynolds‐averaged Navier–Stokes (RANS) equations are modified to account for variable density and viscosity in two‐phase (i.e. water–air) fluid flow systems. A local level set method is used to update the level set function and a least square technique adopted to re‐initialize it at each time step. To assess the developed algorithm and its versatility, a selection of different fluid–structure interaction problems are examined, i.e. an oscillating flow in a two‐dimensional square tank, a breaking dam involving different density fluids, sloshing in a two‐dimensional rectangular tank and a Wigley ship hull travelling in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental findings of the suppression of rotary sloshing on the dynamic response of a liquid storage tank

In liquid storage tanks, rotary sloshing occurs when the frequency of the lateral harmonic load is near the lowest frequency of the tank–liquid system. Rotary sloshing is a type of sloshing that modifies the tank response, which may cause instabilities of the tank wall. However, the consequences of rotary sloshing for the development of strain in the tank wall have not been elucidated. This paper presents an experimental determination of the effects of rotary sloshing on the development of strain and acceleration at various locations of a storage tank. A low-density-polyethylene tank containing water was tested using a shake table. Nine excitations with frequencies near the first free-vibration frequency of the tank–water system were employed. To suppress rotary sloshing, a high-density foam floating lid was utilised as a barrier. Results reveal rotary sloshing boosts not only the development of both hoop and axial strain but also the acceleration in the horizontal direction perpendicular to the excitation. The lid reduced the maximum hoop and axial strain by 500% and 400%, respectively compared to that when rotary sloshing occurs. Moreover, the lid suppressed the nonplanar sloshing by erasing the first three free-vibration frequencies of the tank–water system without the lid.     

圆筒形贮腔中微重力液体非线性晃动的数值模拟

本文讨论低重力液体在圆筒形贮腔中的非线性晃动问题。将ALE（任意的拉格朗日－欧拉）运动学描述关系引入到Navier-Stokes方程中,在时间域上采用分步离散方法中的速度修正格式,利用Galerkinlk加权余量方法推导了系统的有限元数值离散方程;推导了考虑表面张力效应时有限元边界条件的弱积分形式。推导了自由液面上法向矢量的计算公式。模拟了圆筒形贮腔中低重力液体的非线性晃动,并得到了自由液面、波高变化、压力响应等非线性动力特性。揭示了微重力液体非线性晃动的重要特征并将所得结论与现有的实验结果进行了比较。从而证实了本文方法的有效性与正确性。     

Motion stability dynamics for spacecraft coupled with partially filled liquid container

This paper presents a canonical Hamiltonian model of liquid sloshing for the container coupled with spacecraft. Elliptical shape of rigid body is considered as spacecraft structure. Hamiltonian system is an important form of mechanical system. It mostly used to stabilize the potential shaping of dynamical system. Free surface movement of liquid inside the container is called sloshing. If there is uncontrolled resonance between the motion of tank and liquid-frequency inside the tank then such sloshing can be a reason of attitude disturbance or structural damage of spacecraft. Equivalent mechanical model of simple pendulum or mass attached with spring for sloshing is used by many researchers. Mass attached with spring is used as an equivalent model of sloshing to derive the mathematical equations in terms of Hamiltonian model. Analytical method of Lyapunov function with Casimir energy function is used to find the stability for spacecraft dynamics. Vertical axial rotation is taken as the major axial steady rotation for the moving rigid body.     

Stable response of low-gravity liquid non-linear sloshing in a circle cylindrical tank

Under pitch excitation,the sloshing of liquid in circular cylindrical tank includes planar motion,rotary motion and rotary motion inside planar motion.The boundaries between stable motion and unstable motion depend on the radius of the tank,the liquid height,the gravitational intension,the surface tensor and the sloshing damping.In this article,the differential equations of nonlinear sloshing are built first. And by variational principle,the Lagrange function of liquid pressure is constructed in volume intergration form.Then the velocity potential function is expanded in series by wave height function at the free surface.The nonlinear equations with kinematics and dynamics free surface boundary conditions through variation are derived.At last,these equations are solved by multiple-scales method.The influence of Bond number on the global stable response of nonlinear liquid sloshing in circular cylinder tank is analyzed in detail.The result indicates that variation of amplitude frequency response characteristics of the system with Bond,jump,lag and other nonlinear phenomena of liquid sloshing are investigated.     

Nonlinear multimodal model for TLD of irregular tank geometry and small fluid depth

Tuned liquid dampers (TLDs) utilize sloshing fluid to absorb and dissipate structural vibrational energy. TLDs of irregular or complex tank geometry may be required in practice to avoid tank interference with fixed structural or mechanical components. The literature offers few analytical models to predict the response of this type of TLD, particularly when the fluid depth is small. In this paper, a multimodal model is developed utilizing a Boussinesq-type modal theory which is valid for small TLD fluid depths. The Bateman–Luke variational principle is employed to develop a system of coupled nonlinear ordinary differential equations which describe the fluid response when the tank is subjected to base excitation. Energy dissipation is incorporated into the model from the inclusion of damping screens. The fluid model is used to describe the response of a 2D structure–TLD system when the structure is subjected to external loading and the TLD tank geometry is irregular.Shake table experiments are conducted on a rectangular and chamfered tank subjected to unidirectional base excitation. Comparisons of the experimental and predicted sloshing forces and energy dissipation per cycle indicate that the model is able to predict the fluid response at fluid depth ratios greater than h/L=0.10. Next, structure–TLD system tests are conducted and it is found that the model can predict the structural and TLD responses. The simulated and experimental results show that the TLD tank transfers energy between orthogonal structural sway modes.     

Three-dimensional arbitrary Lagrangian–Eulerian numerical prediction method for non-linear free surface oscillation

A numerical prediction method has been proposed to predict non-linear free surface oscillation in an arbitrarily-shaped three-dimensional container. The liquid motions are described with Navier–Stokes equations rather than Laplace equations which are derived by assuming the velocity potential. The profile of a liquid surface is precisely represented with the three-dimensional curvilinear co-ordinates which are regenerated in each computational step on the basis of the arbitrary Lagrangian–Eulerian (ALE) formulation. In the transformed space, the governing equations are discretized on a Lagrangian scheme with sufficient numerical accuracy and the boundary conditions near the liquid surface are implemented in a complete manner. In order to confirm the applicability of the present computational technique, numerical simulations are conducted for the free oscillations of viscid and inviscid liquids and for highly non-linear oscillation. In addition, non-linear sloshing motions caused by horizontal and vertical excitations and a transition from non-linear sloshing to swirling are numerically predicted in three-dimensional cylindrical containers. Conclusively, it is shown that these sloshing motions associated with high non-linearity are reasonably predicted with the present numerical technique. © 1998 John Wiley & Sons, Ltd.     

Optimization of unsteady incompressible Navier–Stokes flows using variational level set method

This paper presents the optimization of unsteady Navier–Stokes flows using the variational level set method. The solid–liquid interface is expressed by the level set function implicitly, and the fluid velocity is constrained to be zero in the solid domain. An optimization problem, which is constrained by the Navier–Stokes equations and a fluid volume constraint, is analyzed by the Lagrangian multiplier based adjoint approach. The corresponding continuous adjoint equations and the shape sensitivity are derived. The level set function is evolved by solving the Hamilton–Jacobian equation with the upwind finite difference method. The optimization method can be used to design channels for flows with or without body forces. The numerical examples demonstrate the feasibility and robustness of this optimization method for unsteady Navier–Stokes flows.Copyright © 2012 John Wiley & Sons, Ltd.     

储液容器内液体晃荡的非线性动力学分析

基于非线性波动理论模型,求解储液容器内液体晃动的固有频率、模态及动力学响应问题。流体使用u_s-u_p状态方程,利用ABAQUS软件的自适应网格技术,建立储液容器液体晃动数学模型,通过施加水平简谐激励得到液体晃动的固有频率和模态,并与解析解对比,验证了该方法的准确性与可行性。然后,分析了矩形储液容器在多种激励作用下液体非线性晃动响应特性。     

EQUIVALENT DYNAMICS MODELING AND ANALYSIS OF LIQUID-FILLED SPACECRAFT WITH FUEL CONSUMPTION 1)

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.     

燃料消耗下充液航天器等效动力学建模与分析

在轨航天器贮腔内的液体可能表现出多种不同的运动模式, 主要包括液体相对于贮腔的整体性刚体运动、自由液面横向晃动、液体起旋后逐步发生明显的旋转晃动及液体自旋运动; 复合三自由度刚体摆晃动模型能够较为全面地描述这些液体运动模式, 同时为研究起旋阶段的液体晃动动力学问题提供了有效手段. 本文对非线性液体晃动刚体摆复合模型作进一步发展, 考虑模型等效参数随贮腔充液比的变化, 提出了变参数的刚体摆复合模型, 该模型适用于研究燃料消耗下非线性晃动类充液航天器大范围运动耦合动力学问题. 采用刚体摆复合模型对球形贮腔内的液体晃动进行等效后, 基于混合坐标意义下的拉格朗日方程推导了一类充液航天器轨道-姿态-晃动全耦合的动力学方程组, 并展开了充液航天器大角度三轴稳定姿态机动和零冲量轨道机动仿真以及航天器耦合动力学响应特性分析. 研究表明: 液体相对于贮腔的运动会造成航天器主刚体位置发生偏移, 当航天器在执行零冲量机动时, 燃料消耗会造成航天器的轨道平动速度无法收敛到零; 贮腔偏心布放时, 航天器在执行轨道机动过程中贮腔内液体易发生剧烈而且形式复杂的晃动行为, 进而可能造成航天器刚体运动的不稳定.      

晃荡液舱内水动冲击力的分段仿射模型

针对舱内晃荡液体与舱壁的相互作用,对舱内水动冲击力的等效力学模型进行了研究。基于混合系统理论,建立了强非线性液体晃荡的等效摆分段仿射模型,重点对矩形液舱的简化等效力学模型进行了分析。利用计算流体动力学软件Flow3D对矩形液舱内的强非线性液体晃荡进行了数值仿真。理论分析表明：分段仿射模型更符合刚性碰撞的假定,可以更有效地描述等效摆和舱壁碰撞时的速度跃变。仿真结果的对比表明：受到激励时,等效摆分段仿射模型所产生的力与Flow3D计算的结果比较接近,利用此模型可以恰当地描述强非线性液体晃荡。     

Nonlinear simulation of a tuned liquid damper with damping screens using a modal expansion technique

Tuned liquid dampers utilize sloshing fluid to control wind-induced structural motions. However, as a result of the nonlinear free surface boundary conditions of fluid sloshing in a two-dimensional rectangular container, a closed-form solution describing the response behaviour is unavailable. Modal expansions, which couple the sloshing modes, are carried out to the first, third and fifth order to construct a system of coupled nonlinear ordinary differential equations that are solved using the Runge–Kutta–Gill Method. Modal damping is incorporated to account for energy losses arising from the fluid viscosity and the inclusion of damping screens. The model is in general agreement with a previous third-order model that incorporated screen damping in the fundamental sloshing mode only. Sinusoidal shake table experiments are conducted to validate the proposed models. Response time histories and frequency response plots assess the model’s prediction of wave heights, sloshing forces, and screen forces. The first-order model accurately predicts the resonant sloshing forces, and forces on a mid-tank screen. The higher-order models better represent the wave heights and forces on an off-centre screen. Experimental results from structure–TLD system tests under random excitation are used to evaluate the performance of the proposed models. The first-order model is able to predict the variance of the structural response and the effective damping the TLD adds to the structure, but as a minimum, a third-order model should be employed to predict the fluid response. It is concluded that a first-order model can be utilized for preliminary TLD design, while a higher-order model should be used to determine the required tank freeboard and the loading on damping screens positioned at off-centre locations.