Sloshing dynamics-modulated fluid angular momentum and moment fluctuations driven by orbital gravity gradient and jitter accelerations in microgravity

The dynamical behavior of spacecraft propellant affected by the asymmetric combined gravity gradient and jitter acceleration are studied. In particular the effect of surface tension on partially-filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank has been investigated. Three different cases of orbital accelerations: (a) gravity gradient-dominated, (b) equally weighted between gravity gradients and jitter, and (c) gravity jitter-dominated accelerations are considered. Fluctuations of angular momentum, fluid moment and bubble mass center caused by slosh wave excitations driven by gravity gradient and jitter accelerations are also investigated.     

Asymmetric slosh wave excitation in a liquid-vapor interface under microgravity

The dynamical behavior of fluids affected by the asymmetric gravity jitter oscillations, in particular, the effect of surface tension on partially-filled rotating fluids in a Dewar tank imposed by time-dependent directions of background reduced gravity accelerations is investigated. Results show that the greater the components of background reduced gravity in radial and circumferential directions, the greater will be the tendency toward increasing amplitude and degrees of asymmetry of the liquid-vapor interface profiles.     

低重环境下俯仰运动圆柱贮箱内液体晃动

考虑低重环境下由于表面张力的影响使得圆柱贮箱内液体呈现弯曲自由液面的情况,以俯仰激励下液体晃动的占优模态振型函数为晃动速度势的基函数,利用傅里叶-贝塞尔级数对贮箱受俯仰激励时的自由液面处的运动边界条件进行展开,得到能描述晃动系统本质的广义状态方程,并分别给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特征的计算方法. 通过具体算例得到了俯仰激励下贮箱内液体晃动特征的动态响应,同时验证了文中方法的收敛性、可行性和正确性.      

ACTUATION OF SLOSHING MODULATED FORCE AND MOMENT ON LIQUID CONTAINER DRIVEN BY JITTER ACCELERATIONS ASSOCIATED WITH SIEW MOTION IN MICROGRAVITY

ＡＣＴＵＡＴＩＯＮＯＦＳＬＯＳＨＩＮＧＭＯＤＵＬＡＴＥＤＦＯＲＣＥＡＮＤＭＯＭＥＮＴＯＮＬＩＱＵＩＤＣＯＮＴＡＩＮＥＲＤＲＩＶＥＮＢＹＪＩＴＴＥＲＡＣＣＥＬＥＲＡＴＩＯＮＳＡＳＳＯＣＩＡＴＥＤＷＩＴＨＳＩＥＷＭＯＴＩＯＮＩＮＭＩＣＲＯＧＲＡＶＩＴＹＲ...     

Instability analysis of nonlinear surface waves in a circular cylindrical container subjected to a vertical excitation

Singular perturbation theory of two-time scale expansions was developed both in inviscid and weak viscous fluids to investigate the motion of single surface standing wave in a liquid-filled circular cylindrical vessel, which is subject to a vertical periodical oscillation. Firstly, it is assumed that the fluid in the circular cylindrical vessel is inviscid, incompressible and the motion is irrotational, a nonlinear evolution equation of slowly varying complex amplitude, which incorporates cubic nonlinear term, external excitation and the influence of surface tension, was derived from solvability condition of high-order approximation. It shows that when forced frequency is low, the effect of surface tension on mode selection of surface wave is not important. However, when forced frequency is high, the influence of surface tension is significant, and can not be neglected. This proved that the surface tension has the function, which causes free surface returning to equilibrium location. Theoretical results much close to experimental results when the surface tension is considered. In fact, the damping will appear in actual physical system due to dissipation of viscosity of fluid. Based upon weakly viscous fluids assumption, the fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates damping term and external excitation, was derived from linearized Navier–Stokes equation. The analytical expression of damping coefficient was determined and the relation between damping and other related parameters (such as viscosity, forced amplitude and depth of fluid) was presented. The nonlinear amplitude equation and a dispersion, which had been derived from the inviscid fluid approximation, were modified by adding linear damping. It was found that the modified results much reasonably close to experimental results. Moreover, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent. Finally, instability of the surface wave is analyzed and properties of the solutions of the modified amplitude equation are determined together with phase-plane trajectories. A necessary condition of forming stable surface wave is obtained and unstable regions are illustrated.     

An SPH multi‐fluid model based on quasi buoyancy for interface stabilization up to high density ratios and realistic wave speed ratios

We introduce a smoothed particle hydrodynamics (SPH) concept for the stabilization of the interface between 2 fluids. It is demonstrated that the change in the pressure gradient across the interface leads to a force imbalance. This force imbalance is attributed to the particle approximation implicit to SPH. To stabilize the interface, a pressure gradient correction is proposed. In this approach, the multi‐fluid pressure gradients are related to the (gravitational and fluid) accelerations. This leads to a quasi‐buoyancy correction for hydrostatic (stratified) flows, which is extended to nonhydrostatic flows. The result is a simple density correction that involves no parameters or coefficients. This correction is included as an extra term in the SPH momentum equation. The new concept for the stabilization of the interface is explored in 5 case studies and compared with other multi‐fluid models. The first case is the stagnant flow in a tank: The interface remains stable up to density ratios of 1:1000 (typical for water and air), in combination with artificial wave speed ratios up to 1:4. The second and third cases are the Rayleigh‐Taylor instability and the rising bubble, where a reasonable agreement between SPH and level‐set models is achieved. The fourth case is an air flow across a water surface up to density ratios of 1:100, artificial wave speed ratios of 1:4, and high air velocities. The fifth case is about the propagation of internal gravity waves up to density ratios of 1:100 and artificial wave speed ratios of 1:4. It is demonstrated that the quasi‐buoyancy model may be used to stabilize the interface between 2 fluids up to high density ratios, with real (low) viscosities and more realistic wave speed ratios than achieved by other weakly compressible SPH multi‐fluid models. Real wave speed ratios can be achieved as long as the fluid velocities are not very high. Although the wave speeds may be artificial in many cases, correct and realistic wave speed ratios are essential in the modelling of heat transfer between 2 fluids (eg, in engineering applications such as gas turbines).     

Liquid sloshing in partly-filled laterally-excited circular tanks equipped with baffles

Linear potential theory in conjunction with the conformal mapping technique are employed to develop rigorous mathematical models for two-dimensional transient sloshing in non-deformable baffled horizontal circular cylindrical vessels, filled with inviscid incompressible fluids to arbitrary depths, and subjected to arbitrary time-dependent lateral accelerations. Three common baffle configurations are considered, namely, a pair of free surface-touching horizontal side baffles, and a central surface-piercing or bottom-mounted vertical baffle of arbitrary extension. The first few normalized antisymmetric/symmetric sloshing frequencies of the partially-filled tanks are tabulated for selected baffle extension and fill depth ratios. Also, the effects of liquid fill depth or baffle length parameter on the impulsive, total and modal convective mass ratios are examined. A ramp-step function is used to replicate the lateral acceleration excitation encountered in an idealized turning maneuver. Durbin's numerical Laplace transform inversion scheme was applied to solve the resulting truncated linear sets of ordinary differential equations in the time-domain. The effects of excitation input time, fill level, and baffle configuration/extension on the force and moment amplification factors are illustrated through appropriate design charts. Furthermore, the transient hydrodynamic responses to a real seismic event are calculated and the effectiveness of baffle configuration/length on suppression of the induced destabilizing lateral forces are examined. Limiting cases are considered and rigorous verifications are made by comparison with the available data as well as with the numerical simulations performed by using a commercial CFD software package.     

Dynamics of rotating conveying mud drill string subjected to torque and longitudinal thrust

In the view of fluid-structure interactions and rotor dynamics, this paper models the lateral vibration of a vertical downward rotating elastic drill string conveying mud subjected to supporting stabilizers, bit torque and longitudinal thrust. The dynamic model involves the rotational inertia of the drill string tube cross section, the gyroscopic effect caused by rotation, the damping due to friction with the surrounding fluid, the gravity force and mud buoyancy. Damped natural frequency, stability and resonance of the drill string system are determined by quadratic eigenvalue problem and investigated at influences of the stabilizer, rotational angular speed, mud flowing velocity, bit torque and thrust. As a result, the drill string can lose stability both at simultaneous and separate influences of the mud conveying, bit torque and thrust, whereas the rotation, stabilizer and gravity of the drill string can improve system stability; the rotational angular speed causing system resonance decreases with the increase of the mud flowing velocity, bit torque and thrust.     

Sloshing dynamics modulated cryogenic helium fluids driven by gravity gradient or jitter accelerations associated with slew motion in microgravity

The mathematical formulation of the Dewar container sloshing dynamics for a partially filled liquid of cryogenic superfluid helium II driven by the gravity gradient or jitter accelerations associated with slew motion in a microgravity environment are studied. The numerical computation of sloshing dynamics is based on the non-inertia container bounded frame and the solution of time-dependent, three-dimensional partial differential equations subjected to the initial and boundary conditions. This study discloses that the capillary effect of sloshing dynamics governs the liquid-vapor interface fluctuations driven by the gravity gradient or jitter accelerations associated with slew motion in a microgravity environment. The peculiar behavior of superfluid helium in response to sloshing dynamics is also investigated.     

Evolution of distortions of the spherical shape of a cavitation bubble in acoustic supercompression

The evolution of small perturbations of the spherical shape of a vapor bubble in the process of its single strong expansion and compression in deuterated acetone is studied. In the mathematical model used the motion of vapor and liquid is broken down into the spherical component and its small nonspherical perturbation. The spherical component is described by the fluid dynamics equations with account for time-dependent heat conduction and evaporation and condensation on the liquid-vapor interface using equations of state constructed from experimental data. In describing the nonspherical component the liquid viscosity and the surface tension are taken into account, while the effect of the bubble content is disregarded. Certain simple analytical formulas are presented which describe the bubble radius at the moment of maximum expansion, its variation in the compression stage, and the evolution of the bubble sphericity distortion in compression.     

重力环境下液体大幅晃动运动脉动球模型及实验研究

现代航天器通常携带大量的液体推进剂, 在航天器的姿态发生变化的过程中, 由于惯性力和重力的作用, 可能会导致液体燃料发生剧烈晃动, 由此产生附加的晃动力会对航天器造成重要影响. 为了得到液体晃动的规律并满足星载计算机实时计算的要求, 本文研究并验证了一种用于等效液体大幅晃动的动力学模型. 首先将液体大幅晃动运动脉动球模型MPBM推广到重力环境中, 通过脉动球的牛顿?欧拉动力学方程和“呼吸运动”过程中能量关系式, 推导出晃动力法向分量的表达式. 同时, 引入不参与晃动的液体的等效模型, 使得液体质心位置的计算更加准确. 通过和文献中实验数据以及CFD软件的计算结果进行比较, 分别验证了推广的MPBM模型在大幅晃动、零动量机动工况下的有效性, 并基于该等效模型, 研究了脉冲激励的不同时序对航天器中液体晃动响应的影响. 最后, 设计并搭建了用于精确测量液体晃动力的实验平台, 验证了MPBM模型在等效非球形储箱的液体晃动时也同样可以很好地反应出晃动力的变化趋势. 本文的研究工作对进一步研究重力环境中充液航天器刚–液耦合动力学行为具有重要的参考价值.      

Direct numerical simulation of thermocapillary flow based on the Volume of Fluid method

A numerical method for direct simulation of thermal Marangoni effects at dynamically deformable interface of two-phase incompressible fluids is developed. The approach is based on the Volume of Fluid (VOF) method with special focus on the numerical treatment of the temperature surface gradient because of its decisive role as the driving force of the flow. The surface gradient calculation is split into computing its length and direction in order to satisfy the correct thermal boundary condition at the interface without losing mobility of the interface. The method is applied to three different types of thermocapillary flow, namely thermocapillary migration of a droplet in an ambient fluid with linear temperature gradient, thermocapillary convection in a liquid layer under linear temperature gradient along the interface, and Marangoni convection due to Bénard–Marangoni instability. In the first case, different aspects of the dynamics of the migration are considered for validation such as the terminal migration velocity, the initial acceleration and quantification of the wall effects. The simulation also considers high convective heat transfer and covers a wide range of Marangoni numbers up to 5000, where good agreement with both theoretical and experimental results is achieved. In the second case, the convection velocity in the liquid layer is compared with an analytical result. In the final application, pattern formation due to the Bénard–Marangoni instability in a liquid layer in square geometry of small aspect ratio is investigated for realistic Biot number and dynamically deformable fluid interface. The results show good agreement with experiments from literature, where our numerical simulation also predicts cell pattern for a particular aspect ratio which is outside the limitation of the above cited experimental work.     

Some properties of multi-fluid stokes flows

The solution of Stokes' equations for a rotating axisymmetric body which possesses reflection symmetry about a planar interface between two infinite immiscible quiescent viscous fluids is shown to be independent of the viscosities of the fluids and identical with the solution when the fluids have the same viscosity. The result is generalized to a rotating axisymmetric system of bodies which possesses reflection symmetry about each interface of a plane stratified system of fluids. An analogous result for two-fluid systems with a nonplanar static interface is also derived. The effect on torque reduction produced by the presence of a second fluid layer adjacent to a rotating axisymmetric body is considered and explicit calculations are given for the case of a sphere. A proof of uniqueness for unbounded multi-fluid Stokes' flow is given and the asymptotic far field structure of the velocity field is determined for axisymmetric flow caused by the rotation of axisymmetric bodies.     

The effect of vertical vibrations on the onset of convection in a planar rotating layer

The effect of vertical vibrations on the convection in a rotating planar fluid layer heated from below was studied. In this case a modulation parameter, the acceleration due to gravity, appears in the problem. The modulation of the parameter may have a significant effect on the onset of convective instability. Parameter modulation in nonrotating layers has been investigated in earlier work [1–3]. The presence of rotation significantly increases the complexity of the mathematical problem, introducing an additional dependence of the solution on the Taylor number Ta and the Prandtl number Pr. Furthermore, an oscillatory convection regime can occur at the stability limit in rotating fluids with Pr < 1. Parameter modulation in the rotating fluid may not only lead to a change in the stability limit and critical wavelength but also to a change in the eigenfrequency of the oscillatory convection. Rauscher and Kelly [4] examined the effect of parameter modulation on the convective stability of a rotating fluid only for the particular case of a sinusoidal variation in the temperature gradient with a small amplitude for Pr = 1, i.e., the effect of modulation was studied on only a steady convection regime.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 12–22, July–August, 1984.     

DAMPING OF VERTICALLY EXCITED SURFACE WAVE IN WEAKLY VISCOUS FLUID

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.     

An improved three-dimensional model for interface pressure calculations in free-surface flows

A three-dimensional method for the calculation of interface pressure in the computational modeling of free surfaces and interfaces is developed. The methodology is based on the calculation of the pressure force at the interfacial cell faces and is mainly designed for volume of fluid (VOF) interface capturing approach. The pressure forces at the interfacial cell faces are calculated according to the pressure imposed by each fluid on the portion of the cell face that is occupied by that fluid. Special formulations for the pressure in the interfacial cells are derived for different orientations of an interface. The present method, referred to as pressure calculation based on the interface location (PCIL), is applied to both static and dynamic cases. First, a three-dimensional motionless drop of liquid in an initially stagnant fluid with no gravity force is simulated as the static case and then two different small air bubbles in water are simulated as dynamic cases. A two-fluid, piecewise linear interface calculation VOF method is used for numerical simulation of the interfacial flow. For the static case, both the continuum surface force (CSF) and the continuum surface stress (CSS) methods are used for surface tension calculations. A wide range of Ohnesorge numbers and density and viscosity ratios of the two fluids are tested. It is shown that the presence of spurious currents (artificial velocities present in case of considerable capillary forces) is mainly due to the inaccurate calculation of pressure forces in the interfacial computational cells. The PCIL model reduces the spurious currents up to more than two orders of magnitude for the cases tested.

Also for the dynamic bubble rise case, it is shown that using the numerical solver employed here, without PCIL, the magnitude of spurious currents is so high that it is not possible to simulate this type of surface tension dominated flows, while using PCIL, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.     

Nonlinear modeling of tuned liquid dampers (TLDs) in rotating wind turbine blades for damping edgewise vibrations

Tuned liquid dampers (TLDs) utilize the sloshing motion of the fluid to suppress structural vibrations and become a natural candidate for damping vibrations in rotating wind turbine blades. The centrifugal acceleration at the tip of a wind turbine blade can reach a magnitude of 7–8g. This facilitates the use of a TLD with a relatively small fluid mass and with feasible geometric dimensions to mitigate the lightly-damped edgewise vibrations effectively. In the present paper, modal expansions are carried out directly on the velocity field and the free surface of the sloshing liquid in the rotating coordinate system. A formulation has been proposed leading to coupled nonlinear ordinary differential equations, which have been obtained through the Galerkin variational approach together with the modal expansion technique. Two models, with one sloshing mode and three sloshing modes, have been studied in the numerical simulation. It is shown that the one-mode model is able to predict the sloshing force and the damped structural response accurately, since the primary damping effect on the structure is achieved by the first sloshing mode of the fluid. Although it is unable to predict the fluid free-surface elevation equally well, the one-mode model can still be utilized for the design of TLD. Parametric optimization of the TLD is carried out based on the one-mode model, and the optimized damper effectively improves the dynamic response of wind turbine blades.     

On damping coefficient due to phase transformation

The damping coefficient of capillary waves due to the evaporation-condensation process at the interface of the two phases of a fluid is evaluated. To highlight the mechanism of the effect of heat and mass transfer across the interface between regions of liquid and vapor, potential flow of incompressible fluids are assumed. Thus other mechanisms of damping are neglected. To fascilitate the analysis, the method of multiple-scale is employed in the analysis, even though the problem is linear.     

Experimental study of the stability of long axisymmetric liquid bridges between solid supports at different temperatures

The stability of axisymmetric, long liquid bridges held captive between two coaxial, circular solid disks kept at different temperatures is considered. Because of the temperature difference between the supporting disks, a thermally-induced surface tension gradient and its associated flow (Marangoni convection) appear in the liquid column, modifying (decreasing) the capillary stability of the bridge. The influence on the stability limits of long, axisymmetric liquid bridges of the combined effect of gravity acceleration and thermally induced surface tension gradients was experimentally analyzed by using very small size liquid bridges (between disks 1 mm in diameter). Experimental results are compared with available analytical results.     

The Compressible Viscous Surface-Internal Wave Problem: Nonlinear Rayleigh–Taylor Instability

This paper concerns the dynamics of two layers of compressible, barotropic, viscous fluid lying atop one another. The lower fluid is bounded below by a rigid bottom, and t he upper fluid is bounded above by a trivial fluid of constant pressure. This is a free boundary problem: the interfaces between the fluids and above the upper fluid are free to move. The fluids are acted on by gravity in the bulk, and at the free interfaces we consider both the case of surface tension and the case of no surface forces.We are concerned with the Rayleigh–Taylor instability when the upper fluid is heavier than the lower fluid along the equilibrium interface. When the surface tension at the free internal interface is below the critical value, we prove that the problem is nonlinear unstable.