Dynamic coupling between horizontal vessel motion and two-layer shallow-water sloshing

Numerical and analytical results are presented for fluid sloshing, of a two-layer inviscid, incompressible and immiscible fluid with thin layers and a rigid lid, coupled to a vessel which is free to undergo horizontal motion governed by a nonlinear spring. Exact analytical results are obtained for the linear problem, giving the natural frequencies and the resonance structure, particularly between the fluid and vessel. A numerical method for the linear and nonlinear equations is developed based on the high-resolution f-wave-propagation finite volume methods due to Bale et al. (2002) [SIAM Journal on Scientific Computing 24, 955–978], adapted to include the pressure gradient at the rigid-lid, and coupled to a Runge–Kutta solver for the vessel motion. The numerical simulations in the linear limit are compared with the exact analytical solutions. The coupled nonlinear numerical solutions with simulations near the internal 1:1 resonance are presented. Of particular interest is the partition of energy between the vessel and fluid motion.     

The pendulum-slosh problem: Simulation using a time-dependent conformal mapping

Suspending a rectangular vessel which is partially filled with fluid from a single rigid pivoting pole produces an interesting theoretical model with which to investigate the dynamic coupling between fluid motion and vessel rotation. The exact equations for this coupled system are derived with the fluid motion governed by the Euler equations relative to the moving frame of the vessel, and the vessel motion governed by a modified forced pendulum equation. The nonlinear equations of motion for the fluid are solved numerically via a time-dependent conformal mapping, which maps the physical domain to a rectangle in the computational domain with a time dependent conformal modulus. The numerical scheme expresses the implicit free-surface boundary conditions as two explicit partial differential equations which are then solved via a pseudo-spectral method in space. The coupled system is integrated in time with a fourth-order Runge–Kutta method. The starting point for the simulations is the linear neutral stability contour discovered by Turner et al. (2015, Journal of Fluid & Structures 52, 166–180). Near the contour the nonlinear results confirm the instability boundary, and far from the neutral curve (parameterized by longer pole lengths) nonlinearity is found to significantly alter the vessel response. Results are also presented for an initial condition given by a superposition of two sloshing modes with approximately the same frequency from the linear characteristic equation. In this case the fluid initial conditions generate large nonlinear vessel motions, which may have implications for systems designed to oscillate in a confined space or on the slosh-induced-rolling of a ship.     

Instability of sloshing motion in a vessel undergoing pivoted oscillations

Suspending a rectangular vessel partially filled with an inviscid fluid from a single rigid pivoting rod produces an interesting physical model for investigating the dynamic coupling between the fluid and vessel motion. The fluid motion is governed by the Euler equations relative to the moving frame of the vessel, and the vessel motion is given by a modified forced pendulum equation. The fully nonlinear, two-dimensional, equations of motion are derived and linearised for small-amplitude vessel and free-surface motions, and the natural frequencies of the system analysed. It is found that the linear problem exhibits an unstable solution if the rod length is shorter than a critical length which depends on the length of the vessel, the fluid height and the ratio of the fluid and vessel masses. In addition, we identify the existence of 1:1 resonances in the system where the symmetric sloshing modes oscillate with the same frequency as the coupled fluid/vessel motion. The implications of instability and resonance on the nonlinear problem are also briefly discussed.     

多场耦合求解非线性气动弹性的研究综述

非线性气动弹性分析中,涉及到非线性流体动力学和非线性结构动力学的耦合问题.阐述了国内外应用计算流体动力学(computational fluid dynamics,CFD)/计算结构动力学(computational structural dynamics,CSD)耦合求解技术来处理非线性气动弹性问题的研究现状,全面分析了流体域和结构域的非线性特征模拟技术、高效网格运动策略、耦合界面相容性条件、耦合计算效率及降阶技术、气动弹性实验及算法验证、多学科耦合求解等关键技术,总结了相关最新研究的方法和成果, 对比分析了其优缺点并提出了展望.      

About a class of nonlinear oscillators with amplitude-independent frequency

This work is concerned with nonlinear oscillators that have a fixed, amplitude-independent frequency. This characteristic, known as isochronicity/isochrony, is achieved by establishing the equivalence between the Lagrangian of the simple harmonic oscillator and the Lagrangian of conservative oscillators with a position-dependent coefficient of the kinetic energy, which can stem from their mass that changes with the displacement or the geometry of motion. Conditions under which such systems have an isochronous center in the origin are discussed. General expressions for the potential energy, equation of motion as well as solutions for a phase trajectory and time response are provided. A few illustrative examples accompanied with numerical verifications are also presented.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Relationship between multibody dynamics computation and nonlinear vibration theory

This paper presents the ground-work of implementing the multibody dynamics codes to analyzing nonlinear coupled oscillators. The recent developments of the multibody dynamics have resulted in several computer codes that can handle large systems of differential and algebraic equations (DAE). However, these codes cannot be used in their current format without appropriate modifications. According to multibody dynamics theory, the differential equations of motion are linear in the acceleration, and the constraints are appended into the equations of motion through Lagrange's multipliers. This formulation should be able to predict the nonlinear phenomena established by the nonlinear vibration theory. This can be achieved only if the constraint algebraic equations are modified to include all the system kinematic nonlinearities. This modification is accomplished by considering secondary nonlinear displacements which are ignored in all current codes. The resulting set of DAE are solved by the Gear stiff integrator. The study also introduced the concept of constrained flexibility and uses an instantaneous energy checking function to improve integration accuracy in the numerical scheme. The general energy balance is a single scalar equation containing all the energy component contributions. The DAE solution is then compared with the solution predicted by the nonlinear vibration theory. It also establishes new foundation for the use of multibody dynamics codes in nonlinear vibration problems. It is found that the simulation CPU time is much longer than the simulation of the original equations of the system.     

Finite element formulation of slender structures with shear deformation based on the Cosserat theory

This paper addresses the derivation of finite element modelling for nonlinear dynamics of Cosserat rods with general deformation of flexure, extension, torsion, and shear. A deformed configuration of the Cosserat rod is described by the displacement vector of the deformed centroid curve and an orthogonal moving frame, rigidly attached to the cross-section of the rod. The position of the moving frame relative to the inertial frame is specified by the rotation matrix, parameterised by a rotational vector. The shape functions with up to third order nonlinear terms of generic nodal displacements are obtained by solving the nonlinear partial differential equations of motion in a quasi-static sense. Based on the Lagrangian constructed by the Cosserat kinetic energy and strain energy expressions, the principle of virtual work is employed to derive the ordinary differential equations of motion with third order nonlinear generic nodal displacements. A cantilever is presented as a simple example to illustrate the use of the formulation developed here to obtain the lower order nonlinear ordinary differential equations of motion of a given structure. The corresponding nonlinear dynamical responses of the structures are presented through numerical simulations using the MATLAB software. In addition, a MicroElectroMechanical System (MEMS) device is presented. The developed equations of motion have furthermore been implemented in a VHDL-AMS beam model. Together with available models of the other components, a netlist of the device is formed and simulated within an electrical circuit simulator. Simulation results are verified against Finite Element Analysis (FEA) results for this device.     

Relaxed representations and improving stability in time-stepping analysis of three-dimensional structural nonlinear dynamics

In the present paper, we propose a representation of the discrete motion equations in structural nonlinear dynamics to obtain an improvement in the stability of time numerical integrations. A geometrically nonlinear total Lagrangian formulation for three-dimensional beam elements in the hypotheses of large rotations and small strains is presented. In this formulation, slopes are used instead of rotation parameters to compute the nonlinear representations of the strain measures in the inertial frame of reference. Such representations of the internal strains??rotations compatibility are then imposed in their time derivatives version. The results, related to Newmark approximations for the variations in the displacement and velocity vectors, show a significant increase in the range of stability of the time integration process and a reduction in the number of Newton iterations required in the time integration steps. The numerical tests, furthermore, show that the variation in the total energy in the time steps has bounded oscillations about the zero value.     

A NUMERICAL METHOD FOR SIMULATING NONLINEAR FLUID-RIGID STRUCTURE INTERACTION PROBLEMS

A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-Lagrangian-Eulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.     

Chaotic attitude and reorientation maneuver for completely liquid-filled spacecraft with flexible appendage

The present paper investigates the chaotic attitude dynamics and reorientation maneuver for completely viscous liquid-filled spacecraft with flexible appendage. All of the equations of motion are derived by using Lagrangian mechanics and then transformed into a form consisting of an unperturbed part plus perturbed terms so that the system＇s nonlinear characteristics can be exploited in phase space. Emphases are laid on the chaotic attitude dynamics produced from certain sets of physical parameter values of the spacecraft when energy dissipation acts to derive the body from minor to major axis spin. Numerical solutions of these equations show that the attitude dynamics of liquid-filled flexible spacecraft possesses characteristics common to random, non- periodic solutions and chaos, and it is demonstrated that the desired reorientation maneuver is guaranteed by using a pair of thruster impulses. The control strategy for reorientation maneuver is designed and the numerical simulation results are presented for both the uncontrolled and controlled spins transition.     

Dynamics of trains and train-like articulated systems travelling in confined fluid—Part 1: Modelling and basic dynamics

The dynamics and stability of a train of flexibly interconnected rigid cylinders travelling in a confined cylindrical “tunnel” subjected to fluid dynamic forces is studied theoretically. Each cylinder, which is coupled to other cylinders and supported by springs and dampers, has degrees of freedom in the lateral translational and rotational directions. The kinetic, dissipation, and potential energies of the system and the generalized forces associated with the fluid dynamic forces acting on the system, such as inviscid fluid dynamic forces, viscous frictional forces, and form drag, are obtained first. Then the equations of motion are derived in a Lagrangian framework. The principal aim of this study is to investigate the effect of the aerodynamic forces on the dynamics of a high-speed train running in a tunnel, or more generally of a train-like system travelling in a coaxial cylindrical tube. The results of this study show that the system loses stability by flutter and that viscous frictional drag has a considerable effect on stability of the system. In addition, the mechanism of instability of the system is clarified with the aid of a study of the modal shapes and energy considerations.     

A numerical scheme for Euler–Lagrange simulation of bubbly flows in complex systems

An Eulerian–Lagrangian approach is developed for the simulation of turbulent bubbly flows in complex systems. The liquid phase is treated as a continuum and the Navier–Stokes equations are solved in an unstructured grid, finite volume framework for turbulent flows. The dynamics of the disperse phase is modeled in a Lagrangian frame and includes models for the motion of each individual bubble, bubble size variations due to the local pressure changes, and interactions among the bubbles and with boundaries. The bubble growth/collapse is modeled by the Rayleigh–Plesset (RP) equation. Three modeling approaches are considered: (a) one‐way coupling, where the influence of the bubble on the fluid flow is neglected, (b) two‐way coupling, where the momentum‐exchange between the fluid and the bubbles is modeled, and (c) volumetric coupling, where the volumetric displacement of the fluid by the bubble motion and the momentum‐exchange are modeled. A novel adaptive time‐stepping scheme based on stability‐analysis of the non‐linear bubble dynamics equations is developed. The numerical approach is verified for various single bubble test cases to show second‐order accuracy. Interactions of multiple bubbles with vortical flows are simulated to study the effectiveness of the volumetric coupling approach in predicting the flow features observed experimentally. Finally, the numerical approach is used to perform a large‐eddy simulation in two configurations: (i) flow over a cavity to predict small‐scale cavitation and inception and (ii) a rising dense bubble plume in a stationary water column. The results show good predictive capability of the numerical algorithm in capturing complex flow features. Copyright © 2010 John Wiley & Sons, Ltd.     

An ALE approach to mechano‐chemical processes in fluid–structure interactions

Mathematical modeling and simulation of fluid–structure interaction problems are in the focus of research already for a longer period. However, taking into account also chemical reactions, leading to structural changes, including changes of mechanical properties of the solid phase, is rather new but for many applications is highly important area. This paper formulates a model system for reactive flow and transport in a vessel system, the penetration of chemical substances into the solid wall. Inside the wall, reactions take place that lead to changes of volume and of the mechanical properties of the wall. Numerical algorithms are developed and used to simulate the dynamics of such a mechano‐chemical fluid–structure interaction problem. As a proof of concept scenario, plaque formation in blood vessels is chosen. The arbitrary Lagrangian Eulerian approach (ALE) is chosen to solve the systems numerically. Temporal discretization of the fully coupled monolithic model is accomplished by backward Euler scheme and spatial discretization by stabilized finite elements. The numerical approach is verified by numerical tests, and effective methods to maintain mesh qualities under large deformations are described. For realistic system parameters, the simulations show that the plaque formation in blood vessel is a long‐time effect. The time scale of the formation is in the simulation of comparable order as in reality. Copyright © 2016 John Wiley & Sons, Ltd.     

Direct numerical simulation of confined swirling jets

This study investigates the Lagrangian acceleration and velocity of fluid particles in swirling flows via direct numerical simulation. The intermittency characteristics of acceleration and velocity of fluid particles are investigated at different swirl numbers and Reynolds numbers. The flatness factor and trajectory curvature are used to analyse the effect of Lagrangian intermittency. The joint probability density function of Lagrangian acceleration and turbulence intensity is shown to explain the augmentation effect of Lagrangian intermittency by the strongly swirling levels under the relatively low intensity of turbulence. In addition, the correlation between the Lagrangian acceleration and the turbulence intensity is enhanced as the swirl level increases. It shows the important effect of swirl on the motion behaviour of fluid particles in the strongly swirling flows.     

Monolithic approach to thermal fluid-structure interaction with nonconforming interfaces

This paper presents a monolithic approach to the thermal fluid-structureinteraction(FSI) with nonconforming interfaces.The thermal viscous flow is governedby the Boussinesq approximation and the incompressible Navier-Stokes equations.Themotion of the fluid domain is accounted for by an arbitrary Lagrangian-Eulerian(ALE)strategy.A pseudo-solid formulation is used to manage the deformation of the fluid do-main.The structure is described by the geometrically nonlinear thermoelastic dynamics.An efficient data transfer strategy based on the Gauss points is proposed to guarantee theequilibrium of the stresses and heat along the interface.The resulting strongly coupledset of nonlinear equations for the fluid,structure,and heat is solved by a monolithicsolution procedure.A numerical example is presented to demonstrate the robustness andefficiency of the methodology.