Modal part analysis of slosh of the liquid in a cylindrical container in the case of small bond number

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｔｉｏｎｏｆｔｈｅｉｎｖｉｓｃｉｄ,ｉｒｒｏｔａｔｉｏｎａｌａｎｄｉｎｃｏｍｐｒｅｓｓｉｂｌｅｆｌｕｉｄｉｎａｃｙｌｉｎｄｒｉｃａｌｃｏｎｔａｉｎｅｒｗｉｔｈｒａｄｉｕｓａｉｓｇｏｖｅｒｎｅｄｂｙＦｉｇ．１Ｓｌｏｓｈｏｆａｌｉｑｕｉｄｉｎａｃｙｌｉｎｄｒｉｃａｌｃｏｎｔａｉｎｅｒ２＝０　　（ｉｎＶ）,（１）ｒ＝０　　（ｏｎｒ＝ａ）,（２）ｚ＝０　　（ｏｎｚ＝－ｈ）,（３）ηｔ＋·η＝ｚ　　（ｏｎｚ＝η）,（４）ｔ＋１２｜…     

Fundamental validation of the finite volume particle method for 3D sloshing dynamics

The use of the finite volume particle (FVP) method is validated for three‐dimensional sloshing dynamics with a free surface by comparing with results from experiments. In the first part, two typical sloshing experiments for a single liquid phase are simulated, and slosh characteristics that include the free surface behavior and hydrodynamic pressure are reported. Moreover, the influences of the circular wall geometry and spatial resolution in the simulation are studied in a sensitivity analysis. In the second part, two sloshing problems with solid bodies are simulated to preliminarily verify the applicability of the FVP method to three‐dimensional solid bodies' motion in liquid flow. Good agreement between simulations and corresponding experiments indicates that the present FVP method well reproduces three‐dimensional sloshing behavior. Copyright © 2010 John Wiley & Sons, Ltd.     

Slosh dynamics of liquid-filled composite containers—A two dimensional meshless local Petrov–Galerkin approach

This paper deals with the studies of sloshing of liquid in partially filled rectangular composite containers subjected to external excitation. The governing equation for inviscid fluid is written as pressure variable form. At each time step, the pressure is evaluated using the meshless local Petrov–Galerkin (MLPG) approach. A local symmetric weak form (LSWF) for linearized sloshing is developed, and a truly meshless method, based on LSWF and moving least squares (MLS) approximation, is presented for the solution of Laplace equation with the requisite boundary conditions. The effectiveness of the developed algorithm is demonstrated through few numerical examples. The comparison of results in terms of natural slosh frequencies, sloshing amplitudes and hydrodynamic pressures obtained in the present investigation are made with those available in the reported literature. To observe the change in the total liquid response due fluid–structure interaction effects, parametric studies are carried out for different cases by varying the fibre orientations and wall thicknesses in the laminated container wall. The present meshless method based on LSWF is found to be simple and attractive with a great potential in engineering applications.     

部分滑移边界条件下自由液面流动问题研究

数值方法进行相关问题的研究就以圆柱形储液罐为例, 考虑部分滑移边界条件, 对所得的高 度非线性微分方程进行了数值求解, 讨论了贮液腔体内液面接触线和液体高度线, 及部分滑 移条件对流体层微观半径的影响; 所得结论对自由液面晃动问题数值研究中的边界条件处理 有参考意义.     

Motion of a permeable shell in a spherical container filled with non-Newtonian fluid

This paper presents an analytical study of creeping motion of a permeable sphere in a spherical container filled with a micro-polar fluid. The drag experienced by the permeable sphere when it passes through the center of the spherical container is studied.Stream function solutions for the flow fields are obtained in terms of modified Bessel functions and Gegenbauer functions. The pressure fields, the micro-rotation components,the drag experienced by a permeable sphere, the wall correction factor, and the flow rate through the permeable surface are obtained for the frictionless impermeable spherical container and the zero shear stress at the impermeable spherical container. Variations of the drag force and the wall correction factor with respect to different fluid parameters are studied. It is observed that the drag force, the wall correction factor, and the flow rate are greater for the frictionless impermeable spherical container than the zero shear stress at the impermeable spherical container. Several cases of interest are deduced from the present analysis.     

流固耦合不可压物质点法及其在晃动问题中的应用

作为一种混合拉格朗日欧拉法,物质点法在流固耦合问题中具有重要的应用前景。对于自由液面的流动问题,基于物质点法框架已建立了弱可压物质点法和完全不可压物质点法,但在处理流固耦合问题时遇到了困难。弱可压物质点法由于采用可压缩状态方程,导致求解时间步长过小,压力振荡严重,产生了非物理的飞溅现象;完全不可压物质点法基于投影算法和不可压条件,消除了弱可压物质点法的压力振荡,提高了时间步长,但难以处理移动边界问题。基于变分形式的投影算法提出了一种新型流固耦合不可压物质点法,得到了体积加权的压力泊松方程PPE(Pressure Poisson Equation),解决了完全不可压物质点法无法处理不规则边界和移动边界的问题。采用流固耦合不可压物质点法研究了运动刚体容器中的液体晃动问题,并与已有实验和数值结果进行对比,验证了算法的正确性和精度。     

Stability of a spinning liquid-filled spacecraft

Summary The stability of a spinning liquid-filled spacecraft has been investigated in the present paper. Using Galerkin's method, the attitude dynamic equations have been given. The Liapunov direct method was employed to obtain a sufficient condition for stability. Three kinds of characteristic modals were investigated: free motion of inviscid fluid, slosh motion and non-slosh motion. All characteristic problems can be solved numerically by the Finite Element Method or the Boundary Element Method. It has been demonstrated that the viscosity of the fluid has a dissipative effect at large Reynolds number, while the slosh motion plays a destabilizing role. The non-slosh model of fluid does not affect the stability criterion. Accepted for publication 19 October 1996     

Experimental study of large amplitude vibrations of a thin plate in contact with sloshing liquids

The linear and geometrically nonlinear (large amplitude) dynamical response of a thin plate in contact with water on one or both sides has been experimentally studied, considering different filling levels. The free liquid surface is free to slosh and the water is delimited by practically rigid walls, except for the thin plate. An experimental method for monitoring and measuring the free surface waves of the fluid has been also used in order to analyze the behavior of the liquid free surface during the nonlinear vibration of the plate forced by harmonic excitation. The plate deflection due to hydrostatic pressure plays a significant role in changing the plate nonlinearity, but tests with liquid on both sides eliminating this effect have been also presented. For excitation in the frequency neighborhood of the fundamental mode of the plate, the oscillation of the free surface of the liquid is characterized by a very large 1/2-subharmonic component.     

Three‐dimensional sloshing: A consistent finite element approach

This paper presents a new computational methodology based on Legendre's polynomials to predict the slosh and acoustic motion in nearly incompressible fluids in both rigid and flexible structures with free surface. Here, we have used a finite element formulation based on Lagrangian frame of reference to model the fluid motion derived using Hamiltonian equation of the fluid system. We formulated three hexahedral finite elements based on strain fields expressed in terms of extended Legendre's polynomials. Sloshing and acoustic motion of liquid is investigated using these newly formulated elements and inf–sup test is performed on these new elements to check the performance of these elements in modeling sloshing under two severe constraints, namely incompressibility and irrotationality. Comparisons of slosh and acoustic frequencies, and mode shapes with exact solutions are given. Dynamic analysis with earthquake and harmonic kind of forcing function is carried out to validate the formulated hexahedral elements to analyze the sloshing response. Numerical results obtained with these new finite elements, and with the present finite element formulation of the mathematical model agree well with the exact solution and as well as with published experimental literature. Copyright © 2010 John Wiley & Sons, Ltd.     

A finite element method for free surface flows of incompressible fluids in three dimensions. Part I. Boundary fitted mesh motion

Computational fluid mechanics techniques for examining free surface problems in two‐dimensional form are now well established. Extending these methods to three dimensions requires a reconsideration of some of the difficult issues from two‐dimensional problems as well as developing new formulations to handle added geometric complexity. This paper presents a new finite element formulation for handling three‐dimensional free surface problems with a boundary‐fitted mesh and full Newton iteration, which solves for velocity, pressure, and mesh variables simultaneously. A boundary‐fitted, pseudo‐solid approach is used for moving the mesh, which treats the interior of the mesh as a fictitious elastic solid that deforms in response to boundary motion. To minimize mesh distortion near free boundary under large deformations, the mesh motion equations are rotated into normal and tangential components prior to applying boundary conditions. The Navier–Stokes equations are discretized using a Galerkin–least square/pressure stabilization formulation, which provides good convergence properties with iterative solvers. The result is a method that can track large deformations and rotations of free surface boundaries in three dimensions. The method is applied to two sample problems: solid body rotation of a fluid and extrusion from a nozzle with a rectangular cross‐section. The extrusion example exhibits a variety of free surface shapes that arise from changing processing conditions. Copyright © 2000 John Wiley & Sons, Ltd.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

A meshless method for numerical simulation of depth‐averaged turbulence flows using a k‐ϵ model

A three‐dimensional numerical model is developed to analyze free surface flows and water impact problems. The flow of an incompressible viscous fluid is solved using the unsteady Navier–Stokes equations. Pseudo‐time derivatives are introduced into the equations to improve computational efficiency. The interface between the two phases is tracked using a volume‐of‐fluid interface tracking algorithm developed in a generalized curvilinear coordinate system. The accuracy of the volume‐of‐fluid method is first evaluated by the multiple numerical benchmark tests, including two‐dimensional and three‐dimensional deformation cases on curvilinear grids. The performance and capability of the numerical model for water impact problems are demonstrated by simulations of water entries of the free‐falling hemisphere and cone, based on comparisons of water impact loadings, velocities, and penetrations of the body with experimental data. For further validation, computations of the dam‐break flows are presented, based on an analysis of the wave front propagation, water level, and the dynamic pressure impact of the waves on the downstream walls, on a specific container, and on a tall structure. Extensive comparisons between the obtained solutions, the experimental data, and the results of other numerical simulations in the literature are presented and show a good agreement. Copyright © 2015 John Wiley & Sons, Ltd.     

垂直参数激励表面波研究进展

受外激励的充液刚性容器中流体的波动问题有实际的工程应用背景.竖直方向的受周期性外激励的充液容器的自由表面波问题--Faraday波问题是流体力学三大不稳定性难题之一(另外两个不稳定性问题是Rayleigh-B\'enard对流和Taylor-Couette流).本文综述了在理想流体中和弱粘性流体中Faraday波的研究成果;介绍了作者在底部垂直激励的圆柱形容器中流体表面波图谱的实验研究和理论分析的结果.最后提出有待进一步研究的问题.图13,参74      

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.     

A fully hydrodynamic model for three‐dimensional,free‐surface flows

The hydrostatic pressure assumption has been widely used in studying water movements in rivers, lakes, estuaries, and oceans. While this assumption is valid in many cases and has been successfully used in numerous studies, there are many cases where this assumption is questionable. This paper presents a three‐dimensional, hydrodynamic model for free‐surface flows without using the hydrostatic pressure assumption. The model includes two predictor–corrector steps. In the first predictor–corrector step, the model uses hydrostatic pressure at the previous time step as an initial estimate of the total pressure field at the new time step. Based on the estimated pressure field, an intermediate velocity field is calculated, which is then corrected by adding the non‐hydrostatic component of the pressure to the estimated pressure field. A Poisson equation for non‐hydrostatic pressure is solved before the second intermediate velocity field is calculated. The final velocity field is found after the free surface at the new time step is computed by solving a free‐surface correction equation. The numerical method was validated with several analytical solutions and laboratory experiments. Model results agree reasonably well with analytical solutions and laboratory results. Model simulations suggest that the numerical method presented is suitable for fully hydrodynamic simulations of three‐dimensional, free‐surface flows. Copyright © 2003 John Wiley & Sons, Ltd.     

A boundary element method for small-amplitude viscous fluid sloshing in couple with structural vibration

A boundary element method is presented for the coupled motion analysis of structural vibration with small-amplitude fluid sloshing in two-dimensional space. The linearized Navier-Stokes equations are considered in frequency domain and transformed into boundary integral equations. An appropriate fundamental solution for the Helmholtz equation with pure imaginary constant is found. The condition of zero-stress is imposed on the free surface, and non-slip condition of fluid particles is imposed on the walls of the container. For rigid motion models, the expressions for added mass and added damping to the structural motion equations are obtained. Some typical numerical examples are presented.     

COUPLED VIBRATORY CHARACTERISTICS OF A RECTANGULAR CONTAINER BOTTOM PLATE

The natural frequencies of an elastic thin plate placed into a rectangular hole and connected to the rigid bottom slab of a rectangular container filled with fluid having a free surface are studied. The fluid is assumed to be incompressible, inviscid and irrotational, and the effect of surface waves is neglected. An analytical-Ritz method is developed to study the vibratory characteristics of the plate in contact with the fluid. First of all, the exact expression of the motion of the fluid is obtained, in which the unknown coefficients are determined by using the method of separation of variables and the method of Fourier series expansion. Then, the Ritz approach is used to obtain the frequency equation of the system. The vibrating beam functions are adopted as the admissible functions for the wet-mode expansion of the plate, and the added virtual mass incremental (AVMI) matrices are obtained for plates with arbitrary boundary conditions. Finally, a convergence study is carried out and some numerical results are given. The accuracy of AVMI factor solutions is discussed by comparing with the more accurate analytical-Ritz solutions presented in this paper. Furthermore, It is seen that the present method is also suitable for the vibration analysis of rectangular plates in contact with infinite fluid by taking the finite, but larger size fluid domain as an approximation in the computation.     

Stability of a numerical algorithm for gas bubble modelling

A free-surface-tracking algorithm based on the SOLA-VOF method is analysed for numerical stability when modelling gas bubble evolution in a fluid. It is shown that an instability can arise from the fact that the bubble pressure varies with its volume. A time step stability criterion is introduced which is a function of the natural oscillation period but does not depend on the mesh size. This dependence suggests that the instability is likely to arise in the case of a general motion of a bubble, especially if break-up occurs. The effect is shown using linear Fourier analysis of the discretized equation for radial bubble oscillation and demonstrated numerically using a CFD code FLOW-3D. One- and three-dimensional situations are considered: a bubble in a fluid bounded by two concentric surfaces and a bubble floating in a fluid chamber with and without gravity. In cases where no analytical solution is available, a numerical method for the stability time step limit calculation is suggested based on finding the natural oscillation frequency. The nature of the instability suggests that it can be a feature of any numerical algorithm which models transient fluid flow with a free surface.