Hamiltonian structure for a neutrally buoyant rigid body interacting with N vortex rings of arbitrary shape: the case of arbitrary smooth body shape

We present a (noncanonical) Hamiltonian model for the interaction of a neutrally buoyant, arbitrarily shaped smooth rigid body with N thin closed vortex filaments of arbitrary shape in an infinite ideal fluid in Euclidean three-space. The rings are modeled without cores and, as geometrical objects, viewed as N smooth closed curves in space. The velocity field associated with each ring in the absence of the body is given by the Biot–Savart law with the infinite self-induced velocity assumed to be regularized in some appropriate way. In the presence of the moving rigid body, the velocity field of each ring is modified by the addition of potential fields associated with the image vorticity and with the irrotational flow induced by the motion of the body. The equations of motion for this dynamically coupled body-rings model are obtained using conservation of linear and angular momenta. These equations are shown to possess a Hamiltonian structure when written on an appropriately defined Poisson product manifold equipped with a Poisson bracket which is the sum of the Lie–Poisson bracket from rigid body mechanics and the canonical bracket on the phase space of the vortex filaments. The Hamiltonian function is the total kinetic energy of the system with the self-induced kinetic energy regularized. The Hamiltonian structure is independent of the shape of the body, (and hence) the explicit form of the image field, and the method of regularization, provided the self-induced velocity and kinetic energy are regularized in way that satisfies certain reasonable consistency conditions.      

Determining the stress intensity factor by using the principle of minimum potential energy

Expressing the total potential energy of the system of a cracked body П by Williams’ infinite series solution of stress and displacement components containing coefficients A_n(n = 1,2,...), we obtain a set of simultaneous linear equations of unknown coefficients A_n by using the principle of minimum potential energy. When the set of equations is solved, the stress intensity factor K₁ can be easily determined. It is equal to √2πaA₁ Take a sample plate as an example. A single-edgc-cracked plate under tension, with the ratio of crack length to the width of the plate being 0.5 and the ratio of half plate height to the width of the plate being 2.0 and 2. 5, has been calculated. Only 20 - 30 coefficients are taken, and the errors in stress intensity factors are within 5%.     

Optimal control of MHD-flow deceleration

A problem of pulsed control for a three-dimensional magnetohydrodynamic (MHD) model is considered. It is demonstrated that singularities of the solution of MHD equations do not develop with time because they are suppressed by a magnetic field. The existence of an optimal control is proved. An optimality system with the solution regular in time as a whole is constructed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 3–10, September–October, 2008.     

Mass transfer effects on the non-Newtonian fluids past a vertical plate embedded in a porous medium with non-uniform surface heat flux

 The present study is devoted to investigate the influences of mass transfer on buoyancy induced flow over vertical flat plate embedded in a non-Newtonian fluid saturated porous medium. The Ostwald–de Waele power-law model is used to characterize the non-Newtonian fluid behavior. Similarity solution for the transformed governing equations is obtained with prescribed variable surface heat flux. Numerical results for the details of the velocity, temperature and concentration profiles are shown on graphs. Excess surface temperature as well as concentration gradient at the wall associated with heat flux distributions, which are entered in tables, have been presented for different values of the power-law index n, buoyancy ration B and the exponent λ as well as Lewis number Le. Received on 26 April 2000     

Experimental study and direct numerical simulation of the evolution of disturbances in a viscous shock layer on a flat plate

The evolution of disturbances in a hypersonic viscous shock layer on a flat plate excited by slow-mode acoustic waves is considered numerically and experimentally. The parameters measured in the experiments performed with a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵ are the transverse profiles of the mean density and Mach number, the spectra of density fluctuations, and growth rates of natural disturbances. Direct numerical simulation of propagation of disturbances is performed by solving the Navier-Stokes equations with a high-order shock-capturing scheme. The numerical and experimental data characterizing the mean flow field, intensity of density fluctuations, and their growth rates are found to be in good agreement. Possible mechanisms of disturbance generation and evolution in the shock layer at hypersonic velocities are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 3–15, September–October, 2006.     

Effects of Temperature-Dependent Viscosity with Soret and Dufour Numbers on Non-Darcy MHD Free Convective Heat and Mass Transfer Past a Vertical Surface Embedded in a Porous Medium

An analysis is presented to investigate the effects of temperature-dependent viscosity, thermal dispersion, Soret number and Dufour number on non-Darcy MHD free convective heat and mass transfer of a viscous, incompressible and electrically conducting fluid past a vertical isothermal surface embedded in a saturated porous medium. The governing partial differential equations are transferred into a system of ordinary differential equations, which are solved numerically using a fourth order Runge–Kutta scheme with the shooting method. Comparisons with previously published work by Hong and Tien [Hong, J. T. and Tien, C. L.: 1987, Int. J. Heat Mass Transfer 30, 143–150] and Sparrow et al. [Sparrow, E. M. et al.: 1964, AIAA J. 2 652–659] are performed and good agreement is obtained. Numerical results of the skin friction coefficient, the local Nusselt number and the local Sherwood number as well as the velocity, temperature and concentration profiles are presented for different physical parameters.     

Analysis and Modeling of the Turbulent Diffusion of Turbulent Kinetic Energy in Natural Convection

Buoyant flows often contain regions with unstable and stable thermal stratification from which counter gradient turbulent fluxes are resulting, e.g. fluxes of heat or of any turbulence quantity. Basing on investigations in meteorology an improvement in the standard gradient-diffusion model for turbulent diffusion of turbulent kinetic energy is discussed. The two closure terms of the turbulent diffusion, the velocity-fluctuation triple correlation and the velocity-pressure fluctuation correlation, are investigated based on Direct Numerical Simulation (DNS) data for an internally heated fluid layer and for Rayleigh–Bénard convection. As a result it is decided to extend the standard gradient-diffusion model for the turbulent energy diffusion by modeling its closure terms separately. Coupling of two models leads to an extended RANS model for the turbulent energy diffusion. The involved closure term, the turbulent diffusion of heat flux, is studied based on its transport equation. This results in a buoyancy-extended version of the Daly and Harlow model. The models for all closure terms and for the turbulent energy diffusion are validated with the help of DNS data for internally heated fluid layers with Prandtl number Pr = 7 and for Rayleigh–Bénard convection with Pr = 0.71. It is found that the buoyancy-extended diffusion model which involves also a transport equation for the variance of the vertical velocity fluctuation gives improved turbulent energy diffusion data for the combined case with local stable and unstable stratification and that it allows for the required counter gradient energy flux.     

Numerical simulation of the flow in a variable-section channel with pulsed-periodic energy supply

Results of numerical simulations of a quasi-one-dimensional unsteady flow in a channel considered as an element of an air-breathing engine are presented. The influence of parameters of energy supplied in the pulsed-periodic mode (power, pulse frequency, and distribution of energy sources along the channel) on the characteristics of the flow with Mach numbers M ₀ = 2.4–4.0 at the channel entrance is determined. A channel configuration that allows the energy supply distribution to be found from the condition of restriction of the maximum value of the gas temperature is proposed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 3–11, January–February, 2009.     

One Equation Model for Turbulent Channel Flow with Second Order Viscoelastic Corrections

A modified second order viscoelastic constitutive equation is used to derive a k–l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P _w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P _w could be either positive or negative. When P _w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P _w  < 0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k–l formulation rather than from the closure for P _w . An alternative closure for P _w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P _w  > 0 but has also a small region of negative P _w in the buffer layer, was also successfully tested. This second model for P _w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.     

Experimental investigation of vortex rings impinging on inclined surfaces

Vortex–ring interactions with oblique boundaries were studied experimentally to determine the effects of plate angle on the generation of secondary vorticity, the evolution of the primary vorticity and secondary vorticity as they interact near the boundary, and the associated energy dissipation. Vortex rings were generated using a mechanical piston-cylinder vortex ring generator at jet Reynolds numbers 2,000–4,000 and stroke length to piston diameter ratios (L/D) in the range 0.75–2.0. The plate angle relative to the initial axis of the vortex ring ranged from 3 to 60°. Flow analysis was performed using planar laser-induced fluorescence (PLIF), digital particle image velocimetry (DPIV), and defocusing digital particle tracking velocimetry (DDPTV). Results showed the generation of secondary vorticity at the plate and its subsequent ejection into the fluid. The trajectories of the centers of circulation showed a maximum ejection angle of the secondary vorticity occurring for an angle of incidence of 10°. At lower incidence angles (<20°), the lower portion of the ring, which interacted with the plate first, played an important role in generation of the secondary vorticity and is a key reason for the maximum ejection angle for the secondary vorticity occurring at an incidence angle of 10°. Higher Reynolds number vortex rings resulted in more rapid destabilization of the flow. The three-dimensional DDPTV results showed an arc of secondary vorticity and secondary flow along the sides of the primary vortex ring as it collided with the boundary. Computation of the moments and products of kinetic energy and vorticity magnitude about the centroid of each vortex ring showed increasing asymmetry in the flow as the vortex interaction with the boundary evolved and more rapid dissipation of kinetic energy for higher incidence angles.     

Fokker–Planck Equations for a Free Energy Functional or Markov Process on a Graph

The classical Fokker–Planck equation is a linear parabolic equation which describes the time evolution of the probability distribution of a stochastic process defined on a Euclidean space. Corresponding to a stochastic process, there often exists a free energy functional which is defined on the space of probability distributions and is a linear combination of a potential and an entropy. In recent years, it has been shown that the Fokker–Planck equation is the gradient flow of the free energy functional defined on the Riemannian manifold of probability distributions whose inner product is generated by a 2-Wasserstein distance. In this paper, we consider analogous matters for a free energy functional or Markov process defined on a graph with a finite number of vertices and edges. If N ≧ 2 is the number of vertices of the graph, we show that the corresponding Fokker–Planck equation is a system of N nonlinear ordinary differential equations defined on a Riemannian manifold of probability distributions. However, in contrast to stochastic processes defined on Euclidean spaces, the situation is more subtle for discrete spaces. We have different choices for inner products on the space of probability distributions resulting in different Fokker–Planck equations for the same process. It is shown that there is a strong connection but there are also substantial discrepancies between the systems of ordinary differential equations and the classical Fokker–Planck equation on Euclidean spaces. Furthermore, both systems of ordinary differential equations are gradient flows for the same free energy functional defined on the Riemannian manifolds of probability distributions with different metrics. Some examples are also discussed.     

Wave motion of an ideal fluid in a narrow open channel

This paper considers a nonlinear integrodifferential model constructed for the motion of an ideal incompressible fluid in an open channel of variable section using the long-wave approximation. A characteristic equation for describing the perturbation propagation velocity in the fluid is derived. Necessary and sufficient conditions of generalized hyperbolicity for the equations of motion are formulated, and the characteristic form of the system is calculated. In the case of a channel of constant width, the model reduces to the Riemann integral invariants which are conserved along the characteristics. It is found that, during the evolution of the flow, the type of the equations of motion can change, which corresponds to long-wave instability for a certain velocity distribution along the channel width. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 2, pp. 61–71, March–April, 2009.     

Global Weak Solutions to Equations of Motion for Magnetic Fluids

We study the differential system governing the flow of an incompressible ferrofluid under the action of a magnetic field. The system consists of the Navier–Stokes equations, the angular momentum equation, the magnetization equation, and the magnetostatic equations. We prove, by using the Galerkin method, a global in time existence of weak solutions with finite energy of an initial boundary-value problem and establish the long-time behavior of such solutions. The main difficulty is due to the singularity of the gradient magnetic force.      

Molecular dynamics simulation of nonodroplets with the modified Lennard-Jones potential function

An Argon droplet in contact with a Platinum surface was simulated by molecular dynamics method. Argon molecules were modeled with modified Lennard-Jones potential function and the Platinum surface was represented by three layers of molecules. The system temperature is fixed at saturation temperature from 72 to 108 K. An Argon droplet in contact with a Platinum surface is also simulated, at different solid–fluid combinations ) \left( {{\frac{{\varepsilon_{sf} }}{\varepsilon }}} \right) (at fixed s_(sf = 0.9s \sigma_{sf} = 0.9\sigma ) from 0.4 to 0.8, and at a temperature of 84 K. It is concluded that the contact angle decreases by increasing system temperature and increases when solid–fluid interaction energy parameter decreases. When the temperature is high enough, the contact angle drops to zero.     

Singular Limits of the Equations of Magnetohydrodynamics

This paper studies the asymptotic limit for solutions to the equations of magnetohydrodynamics, specifically, the Navier–Stokes–Fourier system describing the evolution of a compressible, viscous, and heat conducting fluid coupled with the Maxwell equations governing the behavior of the magnetic field, when Mach number and Alfvén number tends to zero. The introduced system is considered on a bounded spatial domain in \mathbbR³{\mathbb{R}^{3}}, supplemented with conservative boundary conditions. Convergence towards the incompressible system of the equations of magnetohydrodynamics is shown.     

Developed turbulent flow in a plane channel with simultaneous injection through one porous wall and suction through the other

Results of a numerical-theoretical study of a developed turbulent flow of an incompressible fluid in a plane channel with simultaneous injection of mass through one porous wall and suction of the same mass through the other are presented. The system of equations of averaged motion is closed using a turbulent-stress model. The calculated data of the mean and fluctuational characteristics are in reasonable agreement with experimental results for two values of the Reynolds number of the main flow (Re=10,400 and34,000). Al-Farabi Kazakh National State University, Almaty 480121. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 1, pp. 61–68, January–February, 1998.     

Numerical simulation of receptivity of a hypersonic boundary layer to acoustic disturbances

Direct numerical simulations of the evolution of disturbances in a viscous shock layer on a flat plate are performed for a free-stream Mach number M _∞ = 21 and Reynolds number Re _L = 1.44 · 10⁵. Unsteady Navier-Stokes equations are solved by a high-order shock-capturing scheme. Processes of receptivity and instability development in a shock layer excited by external acoustic waves are considered. Direct numerical simulations are demonstrated to agree well with results obtained by the locally parallel linear stability theory (with allowance for the shock-wave effect) and with experimental measurements in a hypersonic wind tunnel. Mechanisms of conversion of external disturbances to instability waves in a hypersonic shock layer are discussed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 84–91, May–June, 2007.     

Effect of Vertical Modulation on the Onset of Filtration Convection

The present paper examines the effect of vertical harmonic vibration on the onset of convection in an infinite horizontal layer of fluid saturating a porous medium. A constant temperature distribution is assigned on the rigid boundaries, so that there exists a vertical temperature gradient. The mathematical model is described by equations of filtration convection in the Darcy–Oberbeck–Boussinesq approximation. The linear stability analysis for the quasi-equilibrium solution is performed using Floquet theory. Employment of the method of continued fractions allows derivation of the dispersion equation for the Floquet exponent σ in an explicit form. The neutral curves of the Rayleigh number Ra versus horizontal wave number α for the synchronous and subharmonic resonant modes are constructed for different values of frequency Ω and amplitude A of vibration. Asymptotic formulas for these curves are derived for large values of Ω using the method of averaging, and, for small values of Ω, using the WKB method. It is shown that, at some finite frequencies of vibration, there exist regions of parametric instability. Investigations carried out in the paper demonstrate that, depending on the governing parameters of the problem, vertical vibration can significantly affect the stability of the system by increasing or decreasing its susceptibility to convection.      

Dynamics of a chain system of rigid bodies with gravity-friction seismic dampers: fixed supports

A design model for a chain system of N elastically linked rigid bodies with a spheroidal gravity-friction damper is proposed. The Lagrange–Painlevé equations of the first kind are used to construct nonlinear dynamical models of a mechanical system undergoing translational vibrations about the equilibrium position. The conditions under which the system moves in one plane are established. The double nonstationary phase–frequency resonance of a system with N = 2 is analyze. After the numerical integration of the systems of differential equations, the phase–frequency surfaces are plotted and examined for several combinations of system parameters under two-frequency loading     

Equations of nonisothermal filtration in fast processes in elastic porous media

The problem of the nonisothermal joint motion of an elastic porous body and the fluid filling the pores is considered for the case where the duration of the physical process is fractions of a second. A rigorous derivation of averaged equations (equations not containing fast oscillating coefficients) based on the Nguetseng two-scale convergence method is proposed. For various combinations of physical parameters of the problem, these equations include anisotropic nonisothermal Stokes equations for the velocity of the fluid component and the equations of nonisothermal acoustics for the displacements of the solid component or anisotropic nonisothermal Stokes equations for a single-velocity continuum. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 113–129, July–August, 2008.