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.     

A computational model for the indentation and phase transformation of a martensitic thin film

We propose a computational model for a stress-induced martensitic phase transformation of a single-crystal thin film by indentation and its reverse transformation to austenite by heating. Our model utilizes a surface energy that allows sharp interfaces with finite energy and a penalty that forces the film to lie above the indenter and undergo a stress-induced austenite-to-martensite phase transformation. We introduce a method to nucleate the martensite-to-austenite phase transformation since in our model the film would otherwise remain in the martensitic phase in a local minimum of the energy.     

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 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.     

On moving heat sources

The two-dimensional thermal problem due to relative motion of a medium and a suddenly activated circular heat source is solved for several boundary conditions. The solutions can be interpreted as for a moving heat source in a stationary medium or a medium moving past a stationary heat source. Uniform and non-uniform temperature, and uniform and non-uniform heat flux boundary conditions are considered. The effect of velocity and radial direction on the temperature distribution is examined. Average, steady-state Nusselt numbers are derived. The transient response of a continuous line source is obtained as a limiting case of the prescribed heat flux solution. Received on 24 September 1996     

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.     

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.     

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.     

基于MEMS压阻传感器的低功耗高过载测试系统设计

为了检验弹体在高冲击环境下的工作情况,提出了一种基于压阻传感器的高过载、低功耗的测试系统设计方案。该系统可承受2×105g的过载冲击测试,且具有采样率高、体积小、功耗低的特点。通过打靶实验验证:在过载测试过程中该系统具有承受高冲击的能力,并能精确地采集到信号微弱变化,保证了弹体在飞行中记录数据的准确性。     

Progressive Frictional Delamination of an Infinite Elastic Film on a Rigid Substrate Due to In-Plane Point Loading

The paper presents a solution to a delamination problem of an infinite elastic film resting on a rigid substrate and loaded by a monotonically increasing in-plane point force. A?rigid-slip contact is assumed between the film and the substrate, leading to the development of two regions at the interface: a damaged zone with a relative slip between the materials, and a region where the interface remains intact. Both film natural and essential boundary conditions are zero on the boundary between these two interfacial zones with the shape of the boundary being a part of the solution. Problem??s self-similarity enables us to obtain an approximate distribution of interfacial traction within the delaminated zone and a shape of the zone itself. For film??s Poisson??s ratio ??=?1 the approximate solution becomes exact. It is argued that this can be treated as a special case of a rigid film sliding on a rigid substrate. The presented approach can be used to obtain approximate closed-form solutions to similar delamination problems.     

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.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

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.     

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.