Time‐domain BEM solution of convection–diffusion‐type MHD equations

The two‐dimensional convection–diffusion‐type equations are solved by using the boundary element method (BEM) based on the time‐dependent fundamental solution. The emphasis is given on the solution of magnetohydrodynamic (MHD) duct flow problems with arbitrary wall conductivity. The boundary and time integrals in the BEM formulation are computed numerically assuming constant variations of the unknowns on both the boundary elements and the time intervals. Then, the solution is advanced to the steady‐state iteratively. Thus, it is possible to use quite large time increments and stability problems are not encountered. The time‐domain BEM solution procedure is tested on some convection–diffusion problems and the MHD duct flow problem with insulated walls to establish the validity of the approach. The numerical results for these sample problems compare very well to analytical results. Then, the BEM formulation of the MHD duct flow problem with arbitrary wall conductivity is obtained for the first time in such a way that the equations are solved together with the coupled boundary conditions. The use of time‐dependent fundamental solution enables us to obtain numerical solutions for this problem for the Hartmann number values up to 300 and for several values of conductivity parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

Unsteady magnetohydrodynamic stagnation point flow — closed-form analytical solutions

This paper investigates the unsteady stagnation point flow and heat transfer of magnetohydrodynamic(MHD) fluids over a moving permeable flat surface. The unsteady Navier-Stokes(NS) equations are transformed into a similarity nonlinear ordinary differential equation, and a closed form solution is obtained for the unsteadiness parameter of 2. The boundary layer energy equation is transformed into a similarity equation,and is solved for a constant wall temperature and a time-dependent uniform wall heat flux case. The solution domain, velocity, and temperature profiles are calculated for different combinations of parameters including the Prandtl number, mass transfer parameter, wall moving parameter, and magnetic parameter. Two solution branches are obtained for certain combinations of the controlling parameters, and a stability analysis demonstrates that the lower solution branch is not stable. The present solutions provide an exact solution to the entire unsteady MHD NS equations, which can be used for validating the numerical code of computational fluid dynamics.     

Application of a vortex method to free surface flows

Vortex methods have found wide applications in various practical problems. The use of vortex methods in free surface flow problems, however, is still very limited. This paper demonstrates a vortex method for practical computation of non-linear free surface flows produced by moving bodies. The method is a potential flow formulation which uses the exact non-linear free surface boundary condition at the exact location of the instantaneous free surface. The position of the free surface, on which vortices are distributed, is updated using a Lagrangian scheme following the fluid particles on the free surface. The vortex densities are updated by the non-linear dynamic boundary condition, derived from the Euler equations, with an iterative Lagrangian numerical scheme. The formulation is tested numerically for a submerged circular cylinder in unsteady translation. The iteration is shown to converge for all cases. The results of the unsteady simulations agree well with classical linearized solutions. The stability of the method is also discussed.     

Averaging of magnetohydrodynamic streams and applicability of the hydraulic approximation to the calculation of MHD flows in channels

At the present time the hydraulic approximation equations are used widely for calculating MHD flows in channels. Several years ago these solutions were considered as a method of expanding our ideas of the qualitative effect of various factors on the MHD flow in the channel of a MHD device. Today, however, the hydraulic analysis methods are beginning to be used for calculations on specific systems. In this case the selection of a particular design solution frequently is based on an analysis of the over-all characteristics (efficiency, power delivered to the external load, etc.) obtained from the hydraulic calculation, where a few percent rather than tens of percent are taken into account.On the other hand, it is known [1] that in gas dynamics the results of the hydraulic calculation for the same specific nonuniform stream may differ by an order of magnitude of tens of percent depending on the averaging method used, since the magnitude of this difference depends on the degree of nonuniformity of the actual stream.We may expect that the nonuniformity of the MHD streams will be far greater than for the gas dynamic flows as a result of the nonuniformities of the force and the thermal effect of the currents flowing in the stream. These nonuniformities may be associated, for example, with the nonuniform distribution of the currents in the channel cross section because of the nonuniform electrical conductivity, which may be significant in spite of the weak nonuniformity of the temperature distribution, or with the presence in the cross section of forces associated with the induced longitudinal component of the magnetic field, the presence of anisotropy of the electrical conductivity, etc.Moreover, in contrast with gas dynamics, in the design of various MHD devices several characteristics (power delivered to the external load, various efficiencies, etc.) which may be calculated in terms of the average value of the gas dynamic parameters are of great importance. Thus, it seems probable that the question of the applicability of the hydraulic approximation to the calculation of MHD flows in channels, the rational selection of the means for averaging the actual flows, the comparison of the results of the hydraulic calculations with the experimental data, and so on, may be far more significant than was the case for the study of gas dynamic flows.     

Solution of a variational problem for the flow in an MHD generator

The most complete study and construction of extremal plasma flow regimes in the channel of an MHD generator may be accomplished using the methods of variational calculus. The variational problem of conducting-gas motion in an MHD channel was first discussed in [1]. The general formulation of the problem for the MHD generator was considered in [2]. Solutions of variational problems for particular cases of extremal flows are given in [2–5].The present study obtains the solution of the variational problem of the flow of a variable conductivity plasma in an MHD generator which has maximal output power for given channel length or volume. An analysis of the solution is made, and a comparison of the extremal flows with optimized flow in a generator with constant values of the electrical efficiency and flow Mach number is carried out.     

Electrosolenoidal flows in a cylindrical cavity

The main difficulties in investigating three-dimensional magnetohydrodynamic (MHD) flows with vorticity arise, first, because it is necessary to solve an independent boundary-value problem in order to find the field of the electromagnetic forces and, second, because the regimes of these flows are strongly nonlinear for the majority of high-power technological MHD processes and a number of natural phenomena. Particular importance attaches to MHD flows generated by the interaction of an electric current applied to the fluid with the magnetic self-field. This class of MHD flows has become known as electrosolenoidal flows [1]. The presence of a definite symmetry in the distribution of the electromagnetic forces and the geometry of the region of the liquid conductor makes it possible to find a solution in self-similar form. The present paper is devoted to exact solutions of the nonlinear equations for axisymmetric electrosolenoidal flows of a conducting incompressible fluid in infinite cylindrical cavities.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 48–53, May–June, 1991.     

Extension of a combined analytical/numerical initial value problem solver for unsteady periodic flow

Here we describe analytical and numerical modifications that extend the Differential Reduced Ejector/ mixer Analysis (DREA), a combined analytical/numerical, multiple species ejector/mixing code developed for preliminary design applications, to apply to periodic unsteady flow. An unsteady periodic flow modelling capability opens a range of pertinent simulation problems including pulse detonation engines (PDE), internal combustion engine ICE applications, mixing enhancement and more fundamental fluid dynamic unsteadiness, e.g. fan instability/vortex shedding problems. Although mapping between steady and periodic forms for a scalar equation is a classical problem in applied mathematics, we will show that extension to systems of equations and, moreover, problems with complex initial conditions are more challenging. Additionally, the inherent large gradient initial condition singularities that are characteristic of mixing flows and that have greatly influenced the DREA code formulation, place considerable limitations on the use of numerical solution methods. Fortunately, using the combined analytical–numerical form of the DREA formulation, a successful formulation is developed and described. Comparison of this method with experimental measurements for jet flows with excitation shows reasonable agreement with the simulation. Other flow fields are presented to demonstrate the capabilities of the model. As such, we demonstrate that unsteady periodic effects can be included within the simple, efficient, coarse grid DREA implementation that has been the original intent of the DREA development effort, namely, to provide a viable tool where more complex and expensive models are inappropriate. Copyright © 2002 John Wiley & Sons, Ltd.     

Magnetohydrodynamic flow in electrodynamically coupled rectangular ducts

In Sezgin^1,2 the problems considered are the magnetohydrodynamic (MHD) flows in an electrodynamically conducting infinite channel and in a rectangular duct respectively, in the presence of an applied magnetic field. In the present paper we extend the solution procedure of these papers to two rectangular channels connected by a barrier which is partially conductor and partially insulator. The problem has been reduced to the solution of a pair of dual series equations and then to the solution of a Fredholm's integral equation of the second kind. The infinite series obtained were transformed to finite integrals containing Bessel Junctions of the second kind to avoid the computations of slowly converging infinite series and infinite integrals with oscillating integrands. The results obtained compared well with those of Butsenieks and Shcherbinin3 which were obtained for the perfectly conducting barrier separating the flows.     

Application of unsteady mixing equations to some aerodynamic problems

Tangential discontinuities [1] are introduced in solving several transient and steady-state problems of gas dynamics. These discontinuities are unstable [2] as a result of the effects of viscosity and thermal conductivity. Therefore it is advisable to replace the tangential discontinuity by a mixing region and account for its interaction with the inviscid flows, establishing on the boundaries of this region the conditions of vanishing friction stress and equality of the velocity and temperature components to the corresponding velocity and temperature components of the inviscid flows. This formulation improves the accuracy of the solution of such problems by posing them as problems with irregular reflection and intersection of shock waves [1].The consideration of the interaction of unsteady turbulent mixing regions with the inviscid flow also permits the formulation of several problems in which the effects of viscosity lead to complete rearrangement of the flow pattern (the lambda-configuration) with the interaction of the reflected shock wave with the boundary layer in the shock tube [3,4], the formation of zones of developed separation ahead of obstacles, etc.).In this connection, §1 presents an analysis of the self-similar solutions of the unsteady turbulent mixing equations (a corresponding analysis of the laminar mixing equations which coincide with the boundary layer equations is presented in [1]). It is shown that these self-similar solutions describe, along with the several problems noted above, the problems of the formation of steady jets and mixing zones in the base wake.As an example, §2 presents, within the framework of the proposed schematization, an approximate solution of the problem of the interaction of a shock wave reflected from a semi-infinite wall with the boundary layer on a horizontal plate behind the incident shock wave. The results obtained are applied to the analysis of reflection in a shock tube. Computational results are presented which are in qualitative agreement with experiment [3, 4].     

The (9,5) HOC formulation for the transient Navier–Stokes equations in primitive variable

We recently proposed an improved (9,5) higher order compact (HOC) scheme for the unsteady two‐dimensional (2‐D) convection–diffusion equations. Because of using only five points at the current time level in the discretization procedure, the scheme was seen to be computationally more efficient than its predecessors. It was also seen to capture very accurately the solution of the unsteady 2‐D Navier–Stokes (N–S) equations for incompressible viscous flows in the stream function–vorticity (ψ – ω) formulation. In this paper, we extend the scope of the scheme for solving the unsteady incompressible N–S equations based on primitive variable formulation on a collocated grid. The parabolic momentum equations are solved for the velocity field by a time‐marching strategy and the pressure is obtained by discretizing the elliptic pressure Poisson equation by the steady‐state form of the (9,5) scheme with the Neumann boundary conditions. In particular, for pressure, we adopt a strategy on the collocated grid in conjunction with ideas borrowed from the staggered grid approach in finite volume. We first apply this extension to a problem having analytical solution and then to the famous lid‐driven square cavity problem. We also apply our formulation to the backward‐facing step problem to see how the method performs for external flow problems. The results are presented and are compared with established numerical results. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2007 John Wiley & Sons, Ltd.     

Functional and generalized separable solutions to unsteady Navier–Stokes equations

The paper studies unsteady Navier–Stokes equations with two space variables. It shows that the non-linear fourth-order equation for the stream function with three independent variables admits functional separable solutions described by a system of three partial differential equations with two independent variables. The system is found to have a number of exact solutions, which generate new classes of exact solutions to the Navier–Stokes equations. All these solutions involve two or more arbitrary functions of a single argument as well as a few free parameters. Many of the solutions are expressed in terms of elementary functions, provided that the arbitrary functions are also elementary; such solutions, having relatively simple form and presenting significant arbitrariness, can be especially useful for solving certain model problems and testing numerical and approximate analytical hydrodynamic methods. The paper uses the obtained results to describe some model unsteady flows of viscous incompressible fluids, including flows through a strip with permeable walls, flows through a strip with extrusion at the boundaries, flows onto a shrinking plane, and others. Some blow-up modes, which correspond to singular solutions, are discussed.     

Unsteady flow of an elastico-viscous fluid in tubes of uniform cross-section

In this paper, the uniqueness of solution for internal bounded unsteady flows of a shortmemory fluid is first established. Closed-form solutions are then obtained for the equations characterizing flows of such fluids in circular and rectangular tubes of uniform cross-section under an arbitrary pressure gradient. Special cases including the oscillatory flow between two parallel plates are discussed.     

Optimization of unsteady incompressible Navier–Stokes flows using variational level set method

This paper presents the optimization of unsteady Navier–Stokes flows using the variational level set method. The solid–liquid interface is expressed by the level set function implicitly, and the fluid velocity is constrained to be zero in the solid domain. An optimization problem, which is constrained by the Navier–Stokes equations and a fluid volume constraint, is analyzed by the Lagrangian multiplier based adjoint approach. The corresponding continuous adjoint equations and the shape sensitivity are derived. The level set function is evolved by solving the Hamilton–Jacobian equation with the upwind finite difference method. The optimization method can be used to design channels for flows with or without body forces. The numerical examples demonstrate the feasibility and robustness of this optimization method for unsteady Navier–Stokes flows.Copyright © 2012 John Wiley & Sons, Ltd.     

Optimal control of unsteady compressible viscous flows

The control of complex, unsteady flows is a pacing technology for advances in fluid mechanics. Recently, optimal control theory has become popular as a means of predicting best case controls that can guide the design of practical flow control systems. However, most of the prior work in this area has focused on incompressible flow which precludes many of the important physical flow phenomena that must be controlled in practice including the coupling of fluid dynamics, acoustics, and heat transfer. This paper presents the formulation and numerical solution of a class of optimal boundary control problems governed by the unsteady two‐dimensional compressible Navier–Stokes equations. Fundamental issues including the choice of the control space and the associated regularization term in the objective function, as well as issues in the gradient computation via the adjoint equation method are discussed. Numerical results are presented for a model problem consisting of two counter‐rotating viscous vortices above an infinite wall which, due to the self‐induced velocity field, propagate downward and interact with the wall. The wall boundary control is the temporal and spatial distribution of wall‐normal velocity. Optimal controls for objective functions that target kinetic energy, heat transfer, and wall shear stress are presented along with the influence of control regularization for each case. Copyright © 2002 John Wiley & Sons, Ltd.     

The effect of the slip condition on unsteady flow due to non-coaxial rotations of disk and a fluid at infinity

An exact solution of the Navier-Stokes equation is constructed for the magnetohydrodynamic (MHD) flow. The flow is due to non-coaxially rotations of a porous disk with slip condition and a fluid at infinity. The solutions for steady and unsteady cases are obtained by Laplace transform method. The effects of magnetic field and slip parameters are shown and discussed.     

Numerical solution of the system of one-dimensional unsteady Navier-Stokes equations for a compressible gas

In many problems encountered in modern gasdynamics, the boundary layer approximations are inadequate to account for the dissipative factors-viscosity and thermal conductivity of the gas-and the solution of the complete system of Navier-Stokes equations is required. This includes, for example, flows with large longitudinal pressure gradients, which in order of magnitude are comparable with or exceed the transverse gradients (temperature jumps, sharp flow rotations, compression shocks, etc.). In many cases, for example in flows with low density, the scale of action of the longitudinal gradients becomes significant, which leads to the need for considering the flow structure in the vicinity of the large gradients. The formulation of certain problems of this type leads to a system of one-dimensional Navier-Stokes equations.We present a difference scheme for the solution of the system of one-dimensional stationary and nonstationary Navier-Stokes equations and give examples of the calculation of the structure of the stationary shock wave front, unsteady gas flow under the influence of sudden heating of one of the boundaries, and unsteady gas flow in the vicinity of the decay of an initial discontinuity. The solution of the stationary problems is accomplished as a result of stabilization as t .The author wishes to thank V. Ya. Likhushin and V. S. Avduevskii for interest in the study and for their valuable counsel during the investigation.     

On the flow of a viscous fluid driven along a channel by suction at porous walls

This paper is a theoretical treatment of the flow of a viscous incompressible fluid driven along a channel by steady uniform suction through porous parallel rigid walls. Many authors have found such flows when they are symmetric, steady and two-dimensional, by assuming a similarity form of solution due to Berman in order to reduce the Navier-Stokes equations to a nonlinear ordinary differential equation. We generalise their work by considering asymmetric flows, unsteady flows and three-dimensional perturbations. By use of numerical calculations, matched asymptotic expansions for large values of the Reynolds number, and the theory of dynamical systems, we find many more exact solutions of the Navier-Stokes equations, examine their stability, and interpret them. In particular, we show that most previously found steady solutions are unstable to antisymmetric two-dimensional disturbances. This leads to a pitchfork bifurcation, stable asymmetric steady solutions, a Hopf bifurcation, stable time-periodic solutions, stable quasi-periodic solutions, phase locking and chaos in succession as the Reynolds number increases.     

Self-similar solutions of non-Newtonian fluid boundary layer in MHD

The self-similar solutions of the boundary layer for a non-Newtonian fluid in MHD were considered in [1, 2] for a power-law velocity distribution along the outer edge of the layer and constant electrical conductivity through the entire flow. However, the MHD flows of many conducting media, which are solutions or molten metals, cannot be described by the MHD equations for non-Newtonian fluids.The self-similar solutions of the boundary layer for a non-Newtonian fluid without account for interaction with the electromagnetic field were studied in [3].In the following we present the self-similar solutions for the boundary layer of pseudoplastic and dilatant fluids with account for the interaction with an electromagnetic field for the case of a power-law velocity distribution along the outer edge of the layer, when the conductivity of the fluid is constant throughout the flow and the magnetic Reynolds number is small.Izv. AN SSSR. Mekhanika Zhidkosti i Gaza, Vol. 2, No. 6, pp. 77–82, 1967The author wishes to thank S. V. Fal'kovich for his interest in this study.     

UNSTEADY CONFINED VISCOUS FLOWS WITH OSCILLATING WALLS AND MULTIPLE SEPARATION REGIONS OVER A DOWNSTREAM-FACING STEP

This paper presents computational solutions for unsteady viscous flows in channels with a downstream-facing step, followed by an oscillating floor. These solutions of the unsteady Navier–Stokes equations are obtained with a time-integration method using artificial compressibility in a fixed computational domain, which is obtained via a time-dependent coordinate transformation from the fluid domain with moving boundaries. The computational method is first validated for steady viscous flows past a downstream-facing step by comparison with previous numerical solutions and experimental results. This method is then used to obtain solutions for unsteady viscous flows with multiple separation regions over a downstream-facing step with oscillating walls, for which there are no previously known solutions. Thus, the present results may be used as benchmark solutions for the unsteady viscous flows with multiple separation regions between fixed and oscillating walls.