Numerical simulation of axisymmetric afterbody flows with jet exhaust

A fifth-order accurate compact difference scheme was used to compute the flow over an axisymmetric afterbody with jet exhaust. The solution was based on the mass-averaged Navier-Stokes equations combined with a two-parameter differential model of turbulence. The computations were performed on a specially generated mesh such that the flow in the exterior and interior of the nozzle could be described simultaneously. Numerical results are presented for various external flow conditions and various chamber pressures.     

On transport equations in slow flows of a polyatomic gas in an external magnetic field

The transport equations in slow flows of a polyatomic gas in an external magnetic field varying with time are obtained by quantum kinetic theory methods. It is shown that local macroscopic variables, changing slowly with time, obey a system of linked equations. In the Burnett approximation, the Wigner operator and stress tensor are obtained. Transport coefficients depending on the magnetic field frequency are determined. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 110, No. 3, pp. 459–469, March, 1997.     

Influence of rapid changes in a channel bottom on free-surface flows

Two-dimensional non-linear free-surface flows in a channel boundedbelow by an uneven bottom with rapid changes are considered.Numerical solutions are computed by a boundary integral equationmethod similar to that first introduced by King & Bloor(1987, J. Fluid Mech., 182, 193–208). Free-surface flowspast localized disturbances, steps and sluice gates are calculated.In addition, weakly non-linear solutions are discussed.     

Transition of MHD boundary layer flow past a stretching sheet

In this paper a study is carried out to understand the transition effect of boundary layer flow: (1) due to a suddenly imposed magnetic field over a viscous flow past a stretching sheet and (2) due to sudden withdrawal of magnetic field over a viscous flow past a stretching sheet under a magnetic field. In both the cases the sheet stretches linearly along the direction of the fluid flow. Governing equations have been non-dimensionalised and the non-dimensionalised equations have been solved using the implicit finite difference method of Crank–Nicholson type. Comparison between the steady state exact solutions and the steady state computed solutions has been carried out. Graphical representation of the dimensionless horizontal velocity, vertical velocity and local skin friction profiles of the steady state and unsteady state has been presented. Computation has been carried out for various values of the magnetic parameter M. The obtained results has been interpreted and discussed.     

Couette and Poiseuille flows of an Oldroyd 6-constant fluid with magnetic field

In this paper, we develop a set of differential equations describing the steady flow of an Oldroyd 6-constant magnetohydrodynamic fluid. The fluid is electrically conducting in the presence of a uniform transverse magnetic field. The developed non-linear differential equation takes into account the effect of the material constants and the applied magnetic field. We presented the solution for three types of steady flows, namely,

(i)

Couette flow

(ii)

Poiseuille flow and

(iii)

generalized Couette flow.

Homotopy analysis method (HAM) is used to solve the non-linear differential equation analytically. It is found from the present analysis that for steady flow the obtained solutions are strongly dependent on the material constants (non-Newtonian parameters) which is different from the model of Oldroyd 3-constant fluid. Numerical solutions are also given and compared with the solutions by HAM.     

The influence of an external magnetic field on the dynamic stress of an elastic conducting one-sided layer with a longitudinal shear crack

We study the interaction of a magnetoelastic shear wave with a curvilinear tunnel crack in an ideally conducting diamagnetic (resp. paramagnetic) one-sided (resp. two-sided) layer in the presence of an external static magnetic field. The bases of the one-sided layer are free of mechanical load, and the rim of the face is clamped or free. The corresponding linearized boundary-value problem of magnetoelasticity is reduced to a singular integrodifferential equation with subsequent implementation on a computer. We give numerical results that characterize the influence of the size of the preliminary magnetic field, the frequencies of the load, the curvature, and the orientation of the crack on the stress intensity factor. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 23, 1992, pp. 96–102.     

Instability of a planar liquid layer in an alternating longitudinal magnetic field with non-zero mean

The conducting liquid interface is found to undulate in an alternating magnetic field. It was shown earlier that ifM =B ₀ ²/μηω, B₀, ω, μ andη being the amplitude (complex) of the alternating longitudinal magnetic field imposed at the interface, the angular frequency of the field, the magnetic permeability and the viscosity respectively, and ifM _c was the critical value ofM then the planar layer was stable or unstable according asM < M _c orM > M _c. In this paper we have determined the stability criterion when in addition to the alternating longitudinal field there acts a uniform field in the same direction. After comparing our results with those obtained earlier, in the absence of the uniform field, we find that the additional uniform field has a significant destabilizing effect.     

Peristaltic motion of a magnetohydrodynamic generalized second‐order fluid in an asymmetric channel

We investigate the peristaltic motion of a magnetohydrodynamic (MHD) generalized second‐order fluid in an asymmetric channel. The governing equations are first modeled and then numerically solved under the long wavelength approximation. Attention has been focused to analyze the shear‐thinning and shear‐thickening effects of the investigated non‐Newtonian fluid, the influence of the magnetic force on the flow, especially the trapping, pumping characteristics caused by the peristalsis of the walls. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Direct numerical simulation of turbulent flows in eccentric pipes

A numerical algorithm was developed for solving the incompressible Navier-Stokes equations in curvilinear orthogonal coordinates. The algorithm is based on a central-difference discretization in space and on a third-order accurate semi-implicit Runge-Kutta scheme for time integration. The discrete equations inherit some properties of the original differential equations, in particular, the neutrality of the convective terms and the pressure gradient in the kinetic energy production. The method was applied to the direct numerical simulation of turbulent flows between two eccentric cylinders. Numerical computations were performed at Re = 4000 (where the Reynolds number Re was defined in terms of the mean velocity and the hydraulic diameter). It was found that two types of flow develop depending on the geometric parameters. In the flow of one type, turbulent fluctuations were observed over the entire cross section of the pipe, including the narrowest gap, where the local Reynolds number was only about 500. The flow of the other type was divided into turbulent and laminar regions (in the wide and narrow parts of the gap, respectively).     

On a possibility of generalization of the Hopf lemma to the case of the Navier-Stokes system with nonzero flows

We examine the Hopf lemma (Leray inequality) which is used in proving the existence of a solution to a nonhomogeneous boundary value problem for the stationary Navier-Stokes equations of an incompressible fluid in a bounded domain. We study a possibility of generalization of a weakened variant of the lemma to the case of nonzero flows through the connected components of the boundary of the domain.     

Combined effect of magnetic field and heat absorption on unsteady free convection and heat transfer flow in a micropolar fluid past a semi-infinite moving plate with viscous dissipation using element free Galerkin method

The fully developed electrically conducting micropolar fluid flow and heat transfer along a semi-infinite vertical porous moving plate is studied including the effect of viscous heating and in the presence of a magnetic field applied transversely to the direction of the flow. The Darcy-Brinkman-Forchheimer model which includes the effects of boundary and inertia forces is employed. The differential equations governing the problem have been transformed by a similarity transformation into a system of non-dimensional differential equations which are solved numerically by element free Galerkin method. Profiles for velocity, microrotation and temperature are presented for a wide range of plate velocity, viscosity ratio, Darcy number, Forchhimer number, magnetic field parameter, heat absorption parameter and the micropolar parameter. The skin friction and Nusselt numbers at the plates are also shown graphically. The present problem has significant applications in chemical engineering, materials processing, solar porous wafer absorber systems and metallurgy.     

Bifurcations of limit cycles in a Z6-equivariant planar vector field of degree 5

A concrete numerical example of Z₆-equivariant planar perturbed Hamiltonian polynomial vector fields of degree 5 having at least 24 limit cycles and the configurations of compound eyes are given by using the bifurcation theory of planar dynamical systems and the method of detection functions. There is reason to conjecture that the Hilbert number H(2k + 1) ⩾ (2k + I)² - 1 for the perturbed Hamiltonian systems.     

Structural elastic and creep models of a UD composite in longitudinal shear

The transient creep of a UD composite with a quadratic arrangement of elastic fibers of quadratic cross section is investigated. The deformational properties of the composite are determined from the known properties of its constituents. A structural model of the UD composite is developed, whose minimal elementary cell contains four elements. The stress-strain state of the elements is assumed homogeneous. Two types of basic and resolving governing equations of transient creep are deduced, which are based on static or kinematic assumptions. In each of the cases, a formula for the longitudinal elastic shear modulus of the composite is found. The stationary solutions of creep equations allow one to obtain formulas of the steady-state creep of the composite in a form similar to Norton’s law. Numerical calculations are also performed, and a comparison of the results with data given in the literature bears witness to the efficiency of the models developed and the solutions obtained. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 4, pp. 437–448, July–August, 2007.     

One approach to a magnetostatic problem for bodies with inclusions in an inhomogeneous external field

A direct magnetostatic problem for magnets with a finite-size inclusion is considered in an integrodifferential form. An approach is used that, under certain conditions, reduces the problem to a single integral equation on a two-dimensional manifold-the inclusion surface. As an important illustrative example, finite formulas are derived to compute the resulting field of a magnetic half-space with a spherical cavity in an arbitrary external field.     

Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field

This paper presents a mathematical analysis of MHD flow and heat transfer to a laminar liquid film from a horizontal stretching surface. The flow of a thin fluid film and subsequent heat transfer from the stretching surface is investigated with the aid of similarity transformation. The transformation enables to reduce the unsteady boundary layer equations to a system of non-linear ordinary differential equations. Numerical solution of resulting non-linear differential equations is found by using efficient shooting technique. Boundary layer thickness is explored numerically for some typical values of the unsteadiness parameter S and Prandtl number Pr, Eckert number Ec and Magnetic parameter Mn. Present analysis shows that the combined effect of magnetic field and viscous dissipation is to enhance the thermal boundary layer thickness.     

Asymptotics of the parabolic Green function for an elliptic operator on a manifold with conical points

We construct the asymptotics ast→0 of the trace of the operator exp(−tP) for an elliptic operatorP on a manifold with conical points. Translated fromMatematicheskie Zametki, Vol. 63, No. 1, pp. 28–36, January, 1998. This research was supported by the Belarus Foundation for Basic Research under grant No. 96-01 00790.     

Effects of pulsatile flow and straight length on low-density lipoprotein transport in a curved arterial wall

Abnormal accumulation of macromolecules such as low-density lipoproteins (LDLs) in the arterial wall causes narrowing and blockage of vessels, which leads to atherosclerosis. Effects of pulsatile nature of blood flows as well as the initial length on transport of the LDL species in the arterial boundary layer region are analyzed numerically in the present work. The set of governing equations consisting of continuity, Navier-Stokes, and species transport is solved using a projection method based on the second-order central difference discretization. The obtained results are in excellent agreement with the pertinent data. The computational results imply that the flow field and concentration distribution are time dependent but the variation of the filtration velocity can be ignored. The LDL concentration boundary layer thickness decreases in the outer part and increases in the inner part for both with or without straight length. Presence of initial straight length generates about 26% growth in the boundary layer thickness, although its effect on the LDL surface concentration (LSC) is negligible. The maximum LSC is related to the regions with minimum wall shear stress in the inner part of the curved artery, which have more potential for formation of atherosclerosis. A new numerical correlation between the LSC and boundary layer thickness is proposed and examined.     

Problem of quality for a certain linear differential game with critical properties

A solution to the pursuit problem for one linear differential game, critical in the sense that it lies on the boundary of solvability of the approach and evasion problems, is given. The result thus obtained is used to answer two question connected with Pontryagin’s methods. Translated fromMatematicheskie Zametki, Vol. 67, No. 4, pp. 484–488, April, 2000.     

Numerical study of magnetohydrodynamic viscous plasma flow in rotating porous media with Hall currents and inclined magnetic field influence

A numerical solution is developed for the viscous, incompressible, magnetohydrodynamic flow in a rotating channel comprising two infinite parallel plates and containing a Darcian porous medium, the plates lying in the x–z plane, under constant pressure gradient. The system is subjected to a strong, inclined magnetic field orientated to the positive direction of the y-axis (rotational axis, normal to the x–z plane). The Navier–Stokes flow equations for a general rotating hydromagnetic flow are reduced to a pair of linear, viscous partial differential equations neglecting convective acceleration terms, for primary velocity (u′) and secondary velocity (v′) where these velocities are directed along the x and y axes. Only viscous terms are retained in the momenta equations. The model is non-dimensionalized and shown to be controlled by a number of dimensionless parameters. The resulting dimensionless ordinary differential equations are solved using a robust numerical method, Network Simulation Methodology. Full details of the numerics are provided. The present solutions are also benchmarked against the analytical solutions presented recently by Ghosh and Pop [Ghosh SK, Pop I. An analytical approach to MHD plasma behaviour of a rotating environment in the presence of an inclined magnetic field as compared to excitation frequency. Int J Appl Mech Eng 2006;11(4):845–856] for the case of a purely fluid medium (infinite permeability). We study graphically the influence of Hartmann number (Ha, magnetic field parameter), Ekman number (Ek, rotation parameter), Hall current parameter (Nh), Darcy number (Da, permeability parameter), pressure gradient (Np) and also magnetic field inclination (θ) on primary and secondary velocity fields. Additionally we investigate the effects of these multiphysical parameters on the dimensionless shear stresses at the plates. Both primary and secondary velocity are seen to be increased with a rise in Darcy number, owing to a simultaneous reduction in Darcian drag force. Primary velocity is seen to decrease with an increase in Hall current parameter (Nh) but there is a decrease in secondary velocity. The study finds important applications in magnetic materials processing, hydromagnetic plasma energy generators, magneto-geophysics and planetary astrophysics.