A Eulerian/Lagrangian model to calculate the evolution of a water droplet spray

We introduce a Eulerian/Lagrangian model to compute the evolution of a spray of water droplets inside a complex geometry. To take into account the complex geometry we define a rectangular mesh and we relate each mesh node to a node function which depends on the location of the node. The time-dependent incompressible and turbulent Navier-Stokes equations are solved using a projection method. The droplets are regarded as individual entities and we use a Lagrangian approach to compute the evolution of the spray. We establish the exchange laws related to mass and heat transfer for a droplet by introducing a mass transfer coefficient and a heat transfer coefficient. The numerical results from our model are compared with those from the literature in the case of a falling droplet in the atmosphere and from experimental investigation in a wind tunnel in the case of a polydisperse spray. The comparison is fairly good. We present the computation of a water droplet spray inside a complex and realistic geometry and determine the characteristics of the spray in the vicinity of obstacles.     

Ideal gas flow behind a finite-amplitude shock wave

The equations of one-dimensional (with a plane of symmetry) adiabatic motion of an ideal gas are transformed to a form convenient for studying flows between a moving piston and a shock wave of variable intensity. The solution is found for the equations of a motion containing a shock wave which propagates through a quiescent gas with variable initial density and constant pressure. This solution contains four arbitrary constants and, in a particular case, gives an example of adiabatic shockless compression by a piston of a gas initially at rest.     

Transient responses of beams and plates subject to travelling load. Miscellaneous results

The steady-state response of a free and infinite Timoshenko beam is specified analytically in terms of non-dimensional displacements and stresses. The beam is supposed loaded by a travelling concentrated force or a moving step load. By a validated explicit numerical calculation, it is shown how a load travelling on a beam at constant velocity, from defined time and abscissa, generates a response which evolves towards the steady-state solution for a part, and towards a quantified transient solution for another part. Asymptotic values are given for the transient displacements and stresses according to the time and the speed of the loading. The solution is also found for a plate subject to a pressure, which spreads respecting the cylindrical symmetry. It is possible to identify in the response a part which follows the pressure front, and which is comparable with the steady-state response of a beam, and another transient part, which generates displacements and stresses with a much less oscillating character. An asymptotic solution is also presented for the plate.The whole series of the results makes it possible to better understand qualitatively the beginning of the transient response of a beam or of a plate to a moving load, and also makes it possible to estimate the stresses and displacements without needing specialised numerical codes.     

Levitation of Magnets and Paramagnetic Bodies in Vessels Filled with Magnetic Fluid

Formulas are obtained for the forces and moments acting on a spherical body made of a paramagnetic material in an uniform applied magnetic field and a magnet in a spherical vessel filled with magnetic fluid. An approximate formula is found for the force acting on bodies in ellipsoidal and cylindrical vessels or in a plane channel with a magnetic fluid in an uniform magnetic field. An analogy between the forces acting on a magnet and a paramagnetic body is demonstrated. The possibility of levitation of magnets and paramagnetic bodies in a vessel with a magnetic fluid is investigated.     

Deformability of thin metal films on elastomer substrates

Recent applications in flexible electronics require that thin metal films grown on elastomer substrates be deformable. However, how such laminates deform is poorly understood. While a freestanding metal film subject to tension will rupture at a small strain by undergoing a necking instability, we anticipate that a substrate will retard this instability to an extent that depends on the relative stiffness and thickness of the film and the substrate. Using a combination of a bifurcation analysis and finite element simulations, we identify three modes of tensile deformation. On a compliant elastomer, a metal film forms a neck and ruptures at a small strain close to that of a freestanding film. On a stiff elastomer, the metal film deforms uniformly to large strains. On an elastomer of intermediate compliance, the metal film forms multiple necks, deforms much beyond the initial bifurcation, and ruptures at a large strain. Our theoretical predictions call for new experiments.     

Large deformation and electrochemistry of polyelectrolyte gels

Immersed in an ionic solution, a network of polyelectrolytes imbibes the solution and swells, resulting in a polyelectrolyte gel. The swelling is reversible, and the amount of swelling is regulated by ionic concentrations, mechanical forces, and electric potentials. This paper develops a field theory to couple large deformation and electrochemistry. A specific material model is described, including the effects of stretching the network, mixing the polymers with the solvent and ions, and polarizing the gel. We show that the notion of osmotic pressure in a gel has no experimental significance in general, but acquires a physical interpretation within the specific material model. The theory is used to analyze several phenomena: a gel swells freely in an ionic solution, a gel swells under a constraint of a substrate, electric double layer at the interface between the gel and the external solution, and swelling of a gel of a small size.     

Interaction of a heavy fluid flow in a channel with a transverse jet barrier

An experimental and numerical analysis of the interaction between a plane horizontal water flow in a rectangular channel (free water current) and a plane thin water jet (water jet curtain) is presented; the jet flows out vertically from either a slot nozzle in the bottom of the channel or the crest of a rigid spillway at a velocity appreciably (several times) greater than the water velocity in the channel. Numerical calculations were carried out using the STAR-CD software package preliminarily tested against the experimental data obtained. The dependence of the water level in the channel at a certain distance ahead of the jet barrier on the main jet parameters and the water flow rate in the horizontal channel is studied. It is found that in the region of the interface between the flows both steady and unsteady (self-oscillatory) flow patterns can be realized. Steady stream/jet interaction patterns of the “ejection” and “ejection-spillway” types are distinguished and a criterion separating these regimes is obtained. The notion of a rigid spillway equivalent to a jet curtain is introduced and an approximate dependence of its height on the relevant parameters of the problem is derived. The possibility of effectively controlling the water level ahead of a rigid spillway with a sharp edge by means of a plane water jet flowing from its crest is investigated. The boundary of transition to self-oscillation interaction patterns in the region of the flow interface is determined. The structure of these flows and a possible mechanism of their generation are described. Within the framework of the inviscid incompressible fluid model in the approximate formulation for a “thin” jet, an analytical dependence of the greatest possible depth of a reservoir filled with a heavy fluid at rest and screened by a vertical jet barrier on the jet parameters is obtained.     

Cell model of nonisothermal gas absorption in gas-liquid bubbly media

A model for combined mass and heat transfer during nonisothermal gas absorption in a two-phase gasliquid bubbly medium with a high gas content and/or large times of gas-liquid contact is suggested. Diffusion and thermal interactions between bubbles is taken into account in the approximation of a cellular model of a bubbly medium whereby a bubbly medium is viewed as a periodic structure consisting of identical spherical cells with periodic boundary conditions at a cell boundary. Distribution of concentration of dissolved gas, temperature distribution in liquid and coefficients of mass and heat transfer during nonisothermal absorption of a soluble pure gas from a bubble by liquid are determined. In the limiting case of absorption without heat release the derived formulas recover the expressions for isothermal absorption.     

Unidimensional SPH simulations of reactive shock tubes in an astrophysical perspective

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method widely used for the modelling of a large variety of astrophysical fluid flows in more than one dimension. Simulations of thermonuclear explosions in stars require, besides the hydrodynamic equations, a realistic equation of state, an energy source term, and a set of nuclear kinetic equations to follow the composition changes of the gas during the explosion. The implementation of a realistic stellar equation of state, and the coupling of hydrodynamics and nuclear burning are investigated in the framework of the simple shock tube geometry. We present and discuss the results of a series of SPH simulations of a detonation in the presence of (1) a single exothermic nuclear reaction, and (2) a restricted network of nuclear reactions. Our results are compared to those of identical simulations performed by other authors using a different hydrodynamic method.     

Function of region

The purpose of this paper is to extend points function and interval functions theoretics to an arbitrary region. For this, the new theory, the contraction of a region, and the retraction of a region; the extension of a region, and the kernel-preserving extension of a region are established by the author. Starting from these concepts, the new definitions of a region function is given. And a kernel (i.e. fixed point) of a region function is connected with a stable centre of defining region of such a region function. Thereby, the region theoretics and algorithms are established.In applications, to find a stable centre of a region, the author has utilized the measure theoretics of matrice defined by Hartfiel⁽⁷⁾ and other authors. The measure problems of coefficient matrice of system of equations of linear algebra associated with some region are discussed.     

Displacement of yield-stress fluids in a fracture

We consider a displacement of several yield-stress fluids in a Hele-Shaw cell. The topic is relevant to the development of a model for the flow of multiple phases inside a narrow fracture with application to hydraulically fracturing a hydrocarbon-bearing underground formation. Existing models for fracturing flows include only pure power-law models without yield stress, and the present work is aimed at filling this gap. The fluids are assumed to be immiscible and incompressible. We consider fluid advection in a plane channel in the presence of density gradients. Gravity is taken into account, so that there can be slumping and gravitational convection. We use the lubrication approximation so that governing equations are reduced to a 2D width-averaged system formed by the quasi-linear elliptic equation for pressure and transport equations for volume concentrations of fluids. The numerical solution is obtained using a finite-difference method. The pressure equation is solved using an iterative algorithm and the Multigrid method, while the transport equations are solved using a second-order TVD flux-limiting scheme with the superbee limiter. This numerical model is validated against three different sets of experiments: (i) gravitational slumping of fluids in a closed Hele-Shaw cell, (ii) viscous fingering of fluids with a high viscosity contrast due to the Saffman–Taylor (S–T) instability in a Hele-Shaw cell at microgravity conditions, (iii) displacement of Bingham fluids in a Hele-Shaw cell with the development of fingers due to the S–T instability. Good agreement is observed between simulations and laboratory data. The model is then used to investigate the joint effect of fingering and slumping. Numerical simulations show that the slumping rate of yield-stress fluid is significantly less pronounced than that of a Newtonian fluid with the same density and viscosity. If a low-viscosity Newtonian fluid is injected after a yield-stress one, the S–T instability at the interface leads to the development of fingers. As a result, fingers penetrating into a fluid with a finite yield stress locally decrease the pressure gradient and unyielded zones develop as a consequence.     

Interaction of a rigid punch and a circular plate with a fixed side and a stress-free face

The axisymmetric contact problem of a rigid punch indentation into an elastic circular plate with a fixed side and a stress-free face is considered. The problem is solved by a method developed for finite bodies which is based on the properties of a biorthogonal system of vector functions. The problem is reduced to a Volterra integral equation (IE) of the first kind for the contract pressure function and to a system of two Volterra IE of the first kind for functions describing the derivative of the displacement of the plate upper surface outside the punch and the normal (or tangential) stress on the plate lower fixed surface. The last two functions are sought as the sum of a trigonometric series and a power-law function with a root singularity. The obtained ill-conditioned systems of linear algebraic equations are regularized by introducing small parameters and have a stable solution. A method for solving the Volterra IE is given. The contact pressure functions, the normal and tangential stresses on the plate fixed surface, and the dimensionless indentation force are found. Several examples of a plane punch computation are given.     

On the problem of vortex interaction with a plane

We consider the well-known problem of the interaction of a vortex filament with a perpendicular plane in a viscous incompressible fluid. In this study, the vortex filament is represented by a semi-infinite rotating needle. Different models are considered: a zero-radius needle and fixed and movable in the axial direction needles of a finite radius. The ranges of the existence of the solution are found, and the correspondence of the flow around a finite-radius needle to that around a zero-radius needle, as the needle radius decreases, is studied.     

A new facility for time-resolved PIV measurements in rotating channels

A new facility to measure the time evolution of 2D velocity fields in a rotating channel is presented, and the accuracy is discussed in detail. Measurements are made by means of a time-resolved PIV system composed of a continuous laser diode, coupled by a fiber optics cable to a laser plane optical module, and a CMOS high-speed camera. Both the PIV system and divergent channel are fixed on a 2.5 m rotating disk. This allows a direct measurement of the relative velocity of flows with Reynolds numbers between 3 × 10³ and 3 × 10⁴ and Rotation numbers between 0.0 and 0.52. These values correspond to the flow conditions in small radial impellers and can be independently adjusted by a change of the relative flow velocity and RPM. It is shown that this new facility allows high signal-to-noise ratios, and that the direct acquisition of the data in a rotating frame drastically reduces the measurement error. The accuracy and high spatial and temporal resolution of the measurements allow a detailed analysis of boundary layer characteristics in stationary and rotating conditions.     

Stability for manifolds of equilibrium states of fractional generalized Hamiltonian systems

For a fractional generalized Hamiltonian system, in terms of Riesz derivatives, stability theory for the manifolds of equilibrium states is presented. The gradient representation and second order gradient representation of a fractional generalized Hamiltonian system are studied, and the conditions under which the system can be considered as a gradient system and a second order gradient system are given, respectively. Then, equilibrium equations, disturbance equations, and first approximate equations of a fractional generalized Hamiltonian system are obtained. A theorem for the stability of the manifolds of equilibrium states of the general autonomous system is used to a fractional generalized Hamiltonian system, and three propositions on the stability of the manifolds of equilibrium states of the system are investigated. As the special cases of this article, the conditions which a fractional generalized Hamiltonian system can be reduced to a generalized Hamiltonian system, a fractional Hamiltonian system and a Hamiltonian system are given, respectively, and the stability theory for the manifolds of equilibrium states of these systems are obtained. Further, a fractional dynamical system and a fractional Volterra model of the three species groups are given to illustrate the method and results of the application. Finally, by using the method in this paper, we construct a new kind of fractional dynamical model, i.e. the fractional Hénon–Heiles model, and we study its stability of the manifolds of equilibrium states.     

The “decaying springs” model for irreversibility in polymer melts

In this work, entanglements in a polymer melt are modeled as a system of parallel springs which form and decay spontaneously. The springs are assumed to be nonlinear, and a certain fraction of them is torn apart by a certain strain.Based on these assumptions, a model of behavior in simple shear is developed. This model is shown to predict a behavior comprising that of a Wagner fluid, and is generalized to a tensorial model of single integral type. The integrand depends on a product of a material function, modeling reversible behavior, and a material functional which takes irreversible processes into account.Irreversibility of network disentanglement, which may occur when deformation changes or reverses direction, can be modeled in this way. It is shown that the two well-known Wagner constitutive equations with and without irreversibility assumptions are special cases of the model developed. In case of a deformation which does not change directions, the new material function and the material functional are multiplied to yield Wagner's damping function.When the rate of spring formation is a function of temperature, the developed model is shown to predict thermorheologically simple behavior. A constitutive equation for non-isothermal flow of polymers is developed with this assumption.     

Lattice of plates in an unsteady supersonic flow

The supersonic unsteady flow of a gas past a lattice of thin, slightly curved profiles has been investigated in several studies. The paper [1] is devoted to an evaluation of the effect of wind tunnel walls on the unsteady aerodynamic characteristics of a profile, and [2] investigates the effects of the boundaries of a free jet. These cases are equivalent respectively to the anti-phase and in-phase oscillations of the profiles of an unstaggered grid. In [3] consideration is given to a more general case of gas flow past a profile in a wind tunnel with perforated walls. Flow past a lattice of profiles with stagger is studied in [4], where the magnitude of the stagger angle is limited by the condition that the lattice leading edge is located in the undisturbed stream.In the present paper we present a method of calculation of the unsteady supersonic flow of a gas past a lattice of profiles with arbitrary stagger. As an example the results are presented of the calculation of the aerodynamic forces and moments acting on an oscillating profile in a wind tunnel with solid walls and in a free jet.     

Integral formulation for a circular cylindrical cavity in infinite solid and a finite length coaxial cylindrical crack compressed axially

A method for solving problems of fracture of an infinite solid with a circular cylindrical cavity and a coaxial cylindrical crack near the surface under an uniform axial compression is proposed using a non-classical criterial approach associated with a mechanism of a local stability loss near the defect. The theory of integral Fourier transforms and series expansions are used to reduce these problems to a system of paired integral equations and then to a system of linear algebraic equations with respect to the contraction parameter.     

Experimental PIV/V3V measurements of vortex-induced fluid–structure interaction in turbulent flow—A new benchmark FSI-PfS-2a

The investigation of the bidirectional coupling between a fluid flow and a structure motion is a growing branch of research in science and industry. Applications of the so-called fluid–structure interactions (FSI) are widespread. To improve coupled numerical FSI simulations, generic experimental benchmark studies of the fluid and the structure are necessary. In this work, the coupling of a vortex-induced periodic deformation of a flexible structure mounted behind a rigid cylinder and a fully turbulent water flow performed at a Reynolds number of Re=30 470 is experimentally investigated with a planar particle image velocimetry (PIV) and a volumetric three-component velocimetry (V3V) system. To determine the structure displacements a multiple-point laser triangulation sensor is used. The three-dimensional fluid velocity results show shedding vortices behind the structure, which reaches the second swiveling mode with a frequency of about 11.2 Hz corresponding to a Strouhal number of St=0.177. Providing phase-averaged flow and structure measurements precise experimental data for coupled computational fluid dynamics (CFD) and computational structure dynamics (CSD) validations are available for this new benchmark case denoted FSI-PfS-2a. The test case possesses four main advantages: (i) the geometry is rather simple; (ii) kinematically, the rotation of the front cylinder is avoided; (iii) the boundary conditions are well defined; (iv) nevertheless, the resulting flow features and structure displacements are challenging from the computational point of view. In addition to the flow field and displacement data a PIV-based force calculation method is used to estimate the lift and drag coefficients of the moving structure.     

On the wake of a circular cylinder with nodal and saddle attachment

Flow field of a cylinder with a mid-span curvature was experimentally investigated in a wind tunnel and a water tunnel. The azimuthal orientation of the cylinder was changed to obtain a nodal, saddle and a mixed nodal–saddle type of flow attachment. Surface flow topology suggested that the nature of the attachment strongly influenced the spanwise distributions of foci structures that play a significant role in introducing three-dimensionality in the immediate wake. Flow visualization in the water tunnel revealed that the length of a vortex formation region also followed the changes in the nature of the attachment. A symmetric shedding of vortices was observed with a saddle type of attachment. Wake mean velocity profiles showed that the velocity defect and therefore the drag of a curved cylinder was minimum for nodal, and maximum for saddle type of attachment. Nomenclature of the wake was compared with asymptotic profiles and equilibrium parameters. Approach to self-preservation, similarity and other features are discussed.