Bifurcations of relative equilibria of an oblate gyrostat with a discrete damper

We investigate relative equilibria of an oblate gyrostat with a discrete damper. Linear and nonlinear methods yield stability conditions for simple spins about the nominal principal axes. We use analytical and numerical methods to explore other equilibria, including bifurcations that occur for varying rotor momentum and damper parameters. These bifurcations are complex structures that are perturbations of the zero rotor momentum case. We use Lyapunov–Schmidt reduction to determine an analytic relationship between parameters to determine conditions for which a jump phenomenon occurs. This paper is declared a work of the U.S. Government and is not subject to copyright protection in the United States     

An Eulerian Description of Fluids Containing Visco-Elastic Particles

Equations governing the flow of fluid containing visco-hyperelastic particles are developed in an Eulerian framework. The novel feature introduced here is to write an evolution equation for the strain. It is envisioned that this will simplify numerical codes which typically compute the strain on Lagrangian meshes moving through Eulerian meshes. Existence results for the flow of linear visco-hyperelastic particles in a Newtonian fluid are established using a Galerkin scheme.     

Hamiltonian structure and relative equilibria of an axisymmetric rigid body in a second degree and order gravity field: cylindrical and generalized hyperbolic equilibria

The full dynamics of an axisymmetric rigid body in a uniformly rotating second degree and order gravity field are investigated, where orbit and attitude motions of the body are coupled through the gravity. Compared with the classical orbital dynamics with the body considered as a point mass, the full dynamics is a higher-precision model in close proximity of the central body where the gravitational orbit–attitude coupling is significant, such as a spacecraft about a small asteroid or an irregular-shaped natural satellite about a planet. The full dynamics are modeled by using the non-canonical Hamiltonian structure, in terms of variables expressed in the frame fixed with the central body. A Poisson reduction is carried out by means of the axial symmetry of the body, and a reduced system with lower dimension, as well as its non-canonical Hamiltonian structure and equations of motion, is obtained through the reduction process. With the second-order potential, three types of relative equilibria are found to be possible: cylindrical equilibria, generalized hyperbolic equilibria, and conic equilibria, which are counterparts to cylindrical equilibria, hyperbolic equilibria, and conic equilibria of an axisymmetric rigid body in a spherical gravity field, respectively. The geometrical properties and existence of the cylindrical equilibria and generalized hyperbolic equilibria are investigated in detail. It has been found that compared with the classical results in a spherical gravity field, the relative equilibria in this study are more complicated and diverse. The most significant difference is that the non-spherical gravity field enables the existence of non-Lagrangian hyperbolic equilibria, called generalized hyperbolic, which cannot exist in a spherical gravity.     

Eulerian–Lagrangian simulations of unsteady gas–liquid flows in bubble columns

We studied the dynamics of gas–liquid flows in a rectangular bubble column using Eulerian–Lagrangian simulations. Three-dimensional, unsteady simulations were performed to simulate the dynamic characteristics of the oscillating bubble plume. The effect of superficial gas velocity and aerated liquid height-to-column width (H/W) ratio on the dynamic and time-averaged flow properties was studied and the simulated results were validated using wall pressure and voidage fluctuation measurements. The effect of lift force and numerical diffusion on the dynamic and time-averaged properties is discussed in detail. Further, the results obtained using the Eulerian–Lagrangian simulations were compared with the Eulerian–Eulerian simulations. The bubble scale information, which is otherwise lost in the Eulerian–Eulerian simulations, was validated using the voidage fluctuation measurements. Such experimentally validated Eulerian–Lagrangian models will be useful for the simulation of mass transfer and reactions in bubble columns.     

Hamiltonian structures of dynamics of a gyrostat in a gravitational field

The dynamics of a gyrostat in a gravitational field is a fundamental problem in celestial mechanics and space engineering. This paper investigates this problem in the framework of geometric mechanics. Based on the natural symplectic structure, non-canonical Hamiltonian structures of this problem are derived in different sets of coordinates of the phase space. These different coordinates are suitable for different applications. Corresponding Poisson tensors and Casimir functions, which govern the phase flow and phase space structures of the system, are obtained in a differential geometric method. Equations of motion, as well as expressions of the force and torque, are derived in terms of potential derivatives. We uncover the underlying Lie group framework of the problem, and we also provide a systemic approach for equations of motion. By assuming that the gravitational field is axis-symmetrical and central, SO(2) and SO(3) symmetries are introduced into the general problem respectively. Using these symmetries, we carry out two reduction processes and work out the Poisson tensors of the reduced systems. Our results in the central gravitational filed are in consistent with previous results. By these reductions, we show how the symmetry of the problem affects the phase space structures. The tools of geometric mechanics used here provide an access to several powerful techniques, such as the determination of relative equilibria on the reduced system, the energy-Casimir method for determining the stability of equilibria, the variational integrators for greater accuracy in the numerical simulation and the geometric control theory for control problems.     

Eulerian formulation of constrained elastica

The formulation of the constrained elastica problem proposed in this paper is predicated on two key concepts: first, the deformed elastica is described by means of the distance from the conduit axis; second, the problem is formulated in terms of the Eulerian curvilinear coordinate of the conduit rather than the natural curvilinear coordinate of the elastica. This approach is further implemented within a segmentation algorithm, which transforms the global constrained elastica problem into a sequence of analogous auxiliary problems that result from dividing the conduit and the elastica into segments limited by contacts. Each auxiliary segment entails solving a segment of elastica subject to isoperimetric constraints corresponding to the assumed positions of the segment ends along the conduit. This new formulation resolves in one stroke a series of issues that afflict the classical Lagrangian approach: (i) the contact detection is reduced to checking whether a threshold on the distance function is violated, (ii) the isoperimetric conditions are transformed into regular boundary conditions, instead of being treated as external integral constraints, (iii) the method yields a well-conditioned set of equations that does not degenerate with decreasing flexural rigidity of the elastica and/or decreasing clearance between the conduit and the elastica.     

Some basic issues in traditional Eulerian formulations of finite elastoplasticity

In traditional Eulerian formulations of finite elastoplasticity, there are some basic issues that still need clarification and further investigation. Among them are the characterization of initial anisotropy for scalar- and tensor-valued constitutive functions, plastic consistency conditions involving time differentiation of anisotropic yield functions, a unified loading criterion for hardening, softening, and perfectly plastic behaviour, etc. Sometimes, it is thought that a satisfactory and complete treatment for these issues could not be achieved within a framework of traditional Eulerian formulation. In this work, efforts are made towards explaining and clarifying these basic issues. Introducing the notion of an Eulerian type rotation-conjugate group of the initial material symmetry group, we show that consistent Eulerian formulations of scalar- and tensorvalued constitutive functions may be achieved for any given type of initial anisotropy. With such consistent Eulerian formulations, we derive plastic consistency conditions in a corotating frame associated with the foregoing rotation-conjugate group. As to the loading conditions for hardening, softening, and perfectly plastic behaviour, we recall and study the unified criterion proposed by Hill [J Mech Phys Solids 6 (1958) 236]. It is pointed out that the tangential elastic stiffness tensor in this unified criterion has to fulfill rather complicated integrability conditions for the elastic rate equation with an objective stress rate. It appears to be far from being a simple or even trivial matter to give an explicit form of the tangential elastic stiffness tensor meeting these conditions. It is indicated that this issue may be resolved by using the logarithmic stress rate. With this rate, an explicit form of Hill's unified loading criterion, together with an explicit form of tangential elastic stiffness tensor, is presented in terms of a general complementary hyperelastic potential.     

Accounting for Polydispersion in the Eulerian Large Eddy Simulation of the Two-Phase Flow in an Aeronautical-type Burner

A major issue for the simulation of two-phase flows in engines concerns the modeling of the liquid disperse phase, either in the Lagrangian or the Eulerian approach. In the perspective of massively parallel computing, the Eulerian approach seems to be more suitable, as it uses the same algorithms as the gaseous phase solver. However taking into account the whole physics of a turbulent spray, especially in terms of polydispersity, requires an additional modeling effort. The Mesoscopic Eulerian Formalism (MEF) [13] accounts for the effect of turbulence on the disperse phase, and was extended to the Large Eddy Simulation framework [41], but is limited to monodisperse flows. In [38], the influence of polydispersity on resolved and unresolved turbulent motions of the disperse phase was highlighted, and a first model was proposed, based on size-conditioned statistics. Starting from this idea, a coupling between the MEF and the Multifluid Approach (MA) [30] is proposed. The MA decomposes the Eulerian phase into several fluid classes called sections, and corresponding to size intervals. Each section uses then size-conditioned closures. The original idea of this work is to use the MEF closures independently in each section, taking into account the mean droplet size of this section. This new approach, called Multifluid Mesoscopic Eulerian Formalism (MMEF), is then able to capture polydispersion with associated size-conditioned turbulent dynamics. First, the importance of polydispersity and the ability of MMEF to capture it are highlighted with a 0D evaporation case and a 2D vortex case, showing its impact on dynamics in both size and physical spaces. Then, the MMEF is applied to the MERCATO configuration of ONERA [18]. Results are compared to monodisperse Eulerian, Lagrangian and experimental results.     

Thermodynamic laws and consistent Eulerian formulation of finite elastoplasticity with thermal effects

Recently it has been demonstrated that, on the basis of the separation D=D^e+D^p arising from the split of the stress power and two consistency criteria for objective Eulerian rate formulations, it is possible to establish a consistent Eulerian rate formulation of finite elastoplasticity in terms of the Kirchhoff stress and the stretching, without involving additional deformation-like variables labelled “elastic” or “plastic”. It has further been demonstrated that this consistent formulation leads to a simple essential structure implied by the work postulate, namely, both the normality rule for plastic flow D^p and the convexity of the yield surface in Kirchhoff stress space. Here, we attempt to place such an Eulerian formulation on the thermodynamic grounds by extending it to a general case with thermal effects, where the consistency requirements are treated in a twofold sense. First, we propose a general constitutive formulation based on the foregoing separation as well as the two consistency criteria. This is accomplished by employing the corotational logarithmic rate and by incorporating an exactly integrable Eulerian rate equation for D^e for thermo-elastic behaviour. Then, we study the consistency of the formulation with thermodynamic laws. Towards this goal, simple forms of restrictions are derived, and consequences are discussed. It is shown that the proposed Eulerian formulation is free in the sense of thermodynamic consistency. Namely, a Helmholtz free energy function in explicit form may be found such that the restrictions from the thermodynamic laws can be fulfilled with positive internal dissipation for arbitrary forms of constitutive functions included in the constitutive formulation. In particular, that is the case for the foregoing essential constitutive structure in the purely mechanical case. These results eventually lead to a complete, explicit constitutive theory for coupled fields of deformation, stress and temperature in thermo-elastoplastic solids at finite deformations.     

Analytical solutions for dynamics of dual-spin spacecraft and gyrostat-satellites under magnetic attitude control in omega-regimes

The attitude dynamics of a dual-spin spacecraft (a gyrostat with one rotor) with magnetic actuators attitude control is considered in the constant external magnetic field at the presence of the spacecraft’s own magnetic dipole moment, which is created proportionally to the angular velocity components (this motion regime can be called as “the omega-regime” or “the omega-maneuver”). The research of the dual-spin spacecraft angular motion under the action of the magnetic restoring torque is fulfilled in the generalized formulation close to the classical mechanics’ task of the heavy body/gyrostat motion in the Lagrange top. Analytical exact solutions of differential equations of the motion are obtained for all parameters in terms of elliptic integrals and the Jacobi functions. New obtained analytical solutions can be classified as results developing the classical fundamental problem of the rigid body and gyrostat motion around the fixed point. The technical application of the omega-regime to the angular reorientation of the spacecraft longitudinal axis along the angular momentum vector is considered.     

An Eulerian approach for simulating frictional heating in disc-pad systems

Thermal stresses as a result from frictional heating must be considered when designing disc brakes, clutches or other rotating machine components with sliding contact conditions. The rotational symmetry of the disc in these kind of applications makes it possible to model these systems using an Eulerian approach instead of a Lagrangian framework. In this paper such an approach is developed and implemented. The disc is formulated in an Eulerian frame where the convective terms are defined by the angular velocity. By utilizing the Eulerian framework, a node-to-node formulation of the contact interface is obtained, producing most accurate frictional heat power solutions. The energy balance of the interface is postulated by introducing an interfacial temperature. Both frictional power and contact conductances are included in this energy balance. The contact problem is solved by a non-smooth Newton method. By adopting the augmented Lagrangian approach, this is done by rewriting Signorini’s contact conditions to an equivalent semi-smooth equation. The heat transfer in the disc is discretized by a Petrov–Galerkin approach, i.e. the numerical difficulties due to the non-symmetric convective matrix appearing in a pure Galerkin discretization is treated by following the streamline-upwind approach. In such manner a stabilization is obtained by adding artificial conduction along the streamlines. For each time step the thermo-elastic contact problem is first solved for the temperature field from the previous time step. Then, the heat transfer problem is solved for the corresponding frictional power. In such manner a temperature history is obtained sequentially via the trapezoidal rule. In particular the parameter is set such that both the Crank–Nicolson and the Galerkin methods are utilized. The method seems very promising. This is demonstrated by solving a two-dimensional benchmark as well as a real disc brake system in three dimensions.     

A multifractal model for linking Lagrangian and Eulerian velocity structure functions

A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.     

Nonlinear Eulerian thermoelasticity for anisotropic crystals

A complete continuum thermoelastic theory for large deformation of crystals of arbitrary symmetry is developed. The theory incorporates as a fundamental state variable in the thermodynamic potentials what is termed an Eulerian strain tensor (in material coordinates) constructed from the inverse of the deformation gradient. Thermodynamic identities and relationships among Eulerian and the usual Lagrangian material coefficients are derived, significantly extending previous literature that focused on materials with cubic or hexagonal symmetry and hydrostatic loading conditions. Analytical solutions for homogeneous deformations of ideal cubic crystals are studied over a prescribed range of elastic coefficients; stress states and intrinsic stability measures are compared. For realistic coefficients, Eulerian theory is shown to predict more physically realistic behavior than Lagrangian theory under large compression and shear. Analytical solutions for shock compression of anisotropic single crystals are derived for internal energy functions quartic in Lagrangian or Eulerian strain and linear in entropy; results are analyzed for quartz, sapphire, and diamond. When elastic constants of up to order four are included, both Lagrangian and Eulerian theories are capable of matching Hugoniot data. When only the second-order elastic constant is known, an alternative theory incorporating a mixed Eulerian–Lagrangian strain tensor provides a reasonable approximation of experimental data.     

万有引力场中陀螺体的混沌运动

研究万有引力场中沿圆轨道运行的非对称陀螺体的姿态运动,引入Ｄｅｐｒｉｔ正则变量建立系统的Ｈａｍｉｌｔｏｎ结构,利用Ｍｅｌｎｉｋｏｖ方法证明在万有引力短作用的昆体产生混沌运动的可能性。对Ｐｏｉｎｃａｒｅ截面的数值计算表明提高陀螺体的转子转速可对混沌起抑制作用。     

Chaotic motion in a class of asymmetrical Kelvin type gyrostat satellite

This work is an investigation on the roots of chaotic attitudinal motion in a class of asymmetrical gyrostat satellites. The result shows that for a class of Kelvin type gyrostat satellite, there is an equivalent rigid spinning satellite with the same attitude dynamics. Finding some constants of motion and eliminating the cyclic coordinates, the rotational kinetic energy is changed to a quadratic form and using Jordan canonical form of the associated inertia tensor and transforming the coordinate system, the Hamiltonian has been changed to those of a rigid satellite. The Hamiltonian has been split into integrable and non-integrable parts. Using Deprit canonical transformation and Andoyer variables the integrable part has been reduced to a one-dimensional form. The reduced Hamiltonian shows that the regular dynamics of the satellite can be chaotic, under the influence of gravitational effects. To demonstrate various attitudinal dynamics of the satellite, a second-order Poincaré map is employed. This research shows firstly, that the attitudinal dynamics of Kelvin type gyrostat satellites and rigid satellites follow the same dynamical patterns, secondly, for non-linear analysis of dynamics of gyrostat satellite based on the perturbation methods, there is a preferable form for Hamiltonian of the system in the near-integrable fashion and thirdly the chaotic motion is originated from the gravitational field effects that can be suppressed by increasing the attitudinal energy of the satellite in comparison with the translational energy.     

Fast computation of the wall distance in unsteady Eulerian fluid-structure computations

Highly nonlinear, turbulent, dynamic, fluid-structure interaction problems characterized by large structural displacements and deformations, as well as self-contact and topological changes, are encountered in many applications. For such problems, the Eulerian computational framework, which is often equipped with an embedded (or immersed) boundary method for computational fluid dynamics, is often the most appropriate framework. In many circumstances, it requires the computation of the time-dependent distance from each active mesh vertex of the embedding mesh to the nearest embedded discrete surface. Such circumstances include, for example, modeling turbulence using the Spalart-Allmaras or detached eddy simulation turbulence models and performing adaptive mesh refinement in order to track the boundary layer. Evaluating at each time step the distance to the wall is computationally prohibitive, particularly in the context of explicit-explicit fluid-structure time-integration schemes. Hence, this paper presents two complementary approaches for reducing this computational cost. The first one recognizes that many quantities depending on the wall distance are relatively insensitive to its inaccurate evaluation in the far field. Therefore, it simplifies a state-of-the-art algorithm for computing the wall distance accordingly. The second approach relies on an effective wall distance error estimator to update the evaluation of the wall distance function only when otherwise, a quantity of interest that depends on it would become tainted by an unacceptable level of error. The potential of combining both approaches for dramatically accelerating the computation of the wall distance is demonstrated with the Eulerian simulation of the inflation of a disk-gap-band parachute system in a supersonic airstream.     

Hamilton principle for an inviscid compressible fluid model in Eulerian coordinates

It is shown that the well-known variational principles for the ideal compressible fluid model in Eulerian coordinates have the following deficiencies:

1. They are not related to the corresponding variational principles in Lagrangian coordinates;
2. The variation procedure in these variational problems does not lead to the equations of motion themselves in the Euler form; rather it leads to relations which correspond to definite classes of solutions of the Euler equations. Here allowance for the equations of the constraints imposed by the adiabaticity and continuity conditions limits the region of application of these variational principles to only potential flows;
3. More general results, involving flows other than potential, are achieved by artificial selection of certain additional constraint conditions imposed on the quantities being varied, and in this case additional clarification is required to ascertain whether any inviscid compressible fluid flow is the extremum of the corresponding variational problem.

A new formulation of the Hamilton principle for the inviscid compressible fluid in Eulerian coordinates is suggested which is free from these deficiencies.     

Comparison of the Eulerian and Lagrangian approaches in studying the flow pattern and heat transfer in a separated axisymmetric turbulent gas-droplet flow

The Eulerian and Lagrangian approaches are used to perform a numerical study of the disperse phase dynamics, turbulence, and heat transfer in a turbulent gas-droplet flow in a tube with sudden expansion with the following ranges of two-phase flow parameters: initial droplet size d ₁ = 0–200 µm and mass fraction of droplets M _L1 = 0–0.1. The main difference between the Eulerian and Lagrangian approaches is the difference in the predictions of the droplet mass fraction: the Eulerian approach predicts a smaller value of M _L both in the recirculation region and in the flow core (the difference reaches 15–20%). It is demonstrated that the disperse phase mass fraction calculated by the Lagrangian approach agrees better with measured data than the corresponding value predicted by the Eulerian approach.