Magnetohydrodynamic flow in a rectangular duct

The magnetohydrodynamic flow of an incompressible, viscous, electrically conducting fluid in a rectangular duct, with an external magnetic field applied transverse to the flow, has been investigated. One of the duct's boundaries which is perpendicular to the magnetic field is taken partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 100. All the infinite series obtained are transformed to infinite integrals first and then to finite integrals which contain modified Bessel functions of the second kind. In this way, the difficulties associated with the computation of infinite integrals with oscillating integrands and slowly converging infinite series, the convergence of which is further affected for large values of M, have been avoided. It is found that, as M increases, boundary layers are formed near the non-conducting boundaries and in the interface region, and a stagnant region is developed in front of the conducting boundary for velocity field. The maximm value of magnetic field takes place on the conducting part. These behaviours are shown on some graphs.     

Magnetohydrodynamic flow in an infinite channel

The magnetohydrodynamic (MHD) flow of an incompressible, viscous, electrically conducting fluid in an infinite channel, under an applied magnetic field has been investigated. The MHD flow between two parallel walls is of considerable practical importance because of the utility of induction flowmeters. The walls of the channel are taken perpendicular to the magnetic field and one of them is insulated, the other is partly insulated, partly conducting. An analytical solution has been developed for the velocity field and magnetic field by reducing the problem to the solution of a Fredholm integral equation of the second kind, which has been solved numerically. Solutions have been obtained for Hartmann numbers M up to 200. All the infinite integrals obtained are transformed to finite integrals which contain modified Bessel functions of the second kind. So, the difficulties associated with the computation of infinite integrals with oscillating integrands which arise for large M have been avoided. It is found that, as M increases, boundary layers are formed near the nonconducting boundaries and in the interface region for both velocity and magnetic fields, and a stagnant region in front of the conducting boundary is developed for the velocity field. Selected graphs are given showing these behaviours.     

Line source diffraction by perfectly conducting successive steps

A rigorous investigation for the high frequency diffraction of the cylindrical waves by the perfectly conducting successive step discontinuities is considered by using the Fourier transform technique in conjunction with the mode-matching method. The corresponding boundary value problem is formulated into a modified scalar Wiener–Hopf equation of the third kind whose approximate solution, which is obtained by iterations of a pair of Fredholm integral equations of the second kind, contains five sets of infinite number of unknown constants satisfying five infinite systems of linear algebraic equations and two types of branch-cut integrals which can be evaluated approximately and numerically due to their types. The effects of the characteristic parameters of the problem such as width and height of the steps on the diffraction phenomenon are shown graphically.     

Surface wave scattering by a material filled rectangular impedance groove in a reactive plane

The problem of TE-polarized surface wave scattering from a rectangular impedance groove located on an infinite reactive plane which is filled with dielectric material is considered for a rather general case where the impedances of the horizontal and vertical sides of the groove have different values. The multiple interactions up to the second order between the edges of the groove are obtained to yield diffracted field. The diffraction problem is first reduced into a modified Wiener-Hopf equation and then solved approximately. The solution contains branch-cut integrals and two infinite sets of constants satisfying two infinite systems of linear algebraic equations. The approximate analytical evaluations of the corresponding integrals as well as the numerical solutions of the linear algebraic equation systems are obtained for various values of the parameters such as the surface reactance of the guiding plane, the vertical and horizontal wall impedances of the groove, the permittivity of the material loading, the width and the height of the groove which permit one to study the effect of these parameters on the diffraction phenomenon.     

Series solutions of unsteady MHD flows above a rotating disk

The series solutions of unsteady flows of a viscous incompressible electrically conducting fluid caused by an impulsively rotating infinite disk are given by means of an analytic technique, namely the homotopy analysis method. Using a set of new similarity transformations, we transfer the Navier–Stokes equations into a pair of nonlinear partial differential equations. The convergent series solutions are obtained, which are uniformly valid for all dimensionless time 0 ≤ τ < ∞ in the whole spatial region 0 ≤ η < ∞. To the best of our knowledge, such kind of series solutions have never been reported. The effect of magnetic number on the velocity is investigated.     

Stationary flow of a conducting liquid in a rectangular channel with two nonconducting walls and two walls having an arbitrary conductivity parallel to an external magnetic field

The flow of a conducting liquid in a channel of rectangular cross section with two walls (parallel to the external magnetic field) having an arbitrary conductivity, the other two being insulators, is considered. The solution of the problem is presented in the form of infinite series. The relationships obtained are used for numerical calculations of the velocity distribution and the distribution of the induced magnetic field over the cross section for several modes of flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkostt i Gaza, No. 5, pp. 46–52, September–October, 1970.     

New development in Poincaré's problem of irregular integrals

In connection with non-Fuchsian equations Poincaré has made an important conclusion: It is impossible to obtain explicit expressions of irregular integrals .To elucidate the essence of Poincaré's problem, we establish correspondence theorem. Irregular integrals are analytic functions of new kind, possessing tree structure; part of which can be represented by conventional recursive series, while its remaining part is expressed by the so-called tree series, not subjecting to any recursive relation at all.In contrast to the numerical solution calculated by infinite determinant of classical theory (Hill-Poincaré-von Koch), our method yields naturally exact analytic solution in explicit form. The method proposed may be used to construct a unifying theory for general equations with variable coefficients, having various kinds of singularities as singular lines.The significance of Poincaré conjecture is discussed, the tree series obtained belong to higher automorphic functions.     

Mixed boundary value problem of heat conduction for infinite slab

Summary An axisymmetric steady state heat conduction boundary value problem having mixed boundary conditions on both faces of an infinite slab, is reduced to a pair of Fredholm integral equations of the second kind. For large values of h, the slab thickness, a solution correct to O(h ^–6) is obtained by expanding the kernels in power series.Presently at Imperial College, London.     

Asymptotically Exact Solution of the Problem of Harmonic Vibrations of an Elastic Parallelepiped

An asymptotically exact solution of the classical problem of elasticity about the steadystate forced vibrations of an elastic rectangular parallelepiped is constructed. The general solution of the vibration equations is constructed in the form of double Fourier series with undetermined coefficients, and an infinite system of linear algebraic equations is obtained for determining these coefficients. An analysis of the infinite system permits determining the asymptotics of the unknowns which are used to convolve the double series in both equations of the infinite systems and the displacement and stress components. The efficiency of this approach is illustrated by numerical examples and comparison with known solutions. The spectrum of the parallelepiped symmetric vibrations is studied for various ratios of its sides.     

On unsteady unidirectional flows of a second grade fluid

Some properties of unsteady unidirectional flows of a fluid of second grade are considered for flows produced by the sudden application of a constant pressure gradient or by the impulsive motion of one or two boundaries. Exact analytical solutions for these flows are obtained and the results are compared with those of a Newtonian fluid. It is found that the stress at the initial time on the stationary boundary for flows generated by the impulsive motion of a boundary is infinite for a Newtonian fluid and is finite for a second grade fluid. Furthermore, it is shown that initially the stress on the stationary boundary, for flows started from rest by sudden application of a constant pressure gradient is zero for a Newtonian fluid and is not zero for a fluid of second grade. The required time to attain the asymptotic value of a second grade fluid is longer than that for a Newtonian fluid. It should be mentioned that the expressions for the flow properties, such as velocity, obtained by the Laplace transform method are exactly the same as the ones obtained for the Couette and Poiseuille flows and those which are constructed by the Fourier method. The solution of the governing equation for flows such as the flow over a plane wall and the Couette flow is in a series form which is slowly convergent for small values of time. To overcome the difficulty in the calculation of the value of the velocity for small values of time, a practical method is given. The other property of unsteady flows of a second grade fluid is that the no-slip boundary condition is sufficient for unsteady flows, but it is not sufficient for steady flows so that an additional condition is needed. In order to discuss the properties of unsteady unidirectional flows of a second grade fluid, some illustrative examples are given.     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

Remarks on the solution of extended Stokes' problems

The analytical solutions of first and second Stokes' problems are discussed, for infinite and finite-depth flows of a Newtonian fluid in planar geometries. Problems arising from the motion of the wall as a whole (one-dimensional flows) as well as of only one half of the wall (two-dimensional) are solved and the wall stresses are evaluated.The solutions are written in real form. In many cases, they improve the ones in literature, leading to simpler mathematical forms of velocities and stresses. The numerical computation of the solutions is performed by using recurrence relations and elementary integrals, in order to avoid the evaluation of integrals of rapidly oscillating functions.The main physical features of the solutions are also discussed. In particular, the steady-state solutions of the second Stokes' problems are analyzed by separating their “in phase” and “in quadrature” components, with respect to the wall motion. By using this approach, stagnation points have been found in infinite-depth flows.     

NEW THEORY FOR EQUATIONS OF NON-FUGHSIAN TYPE REPRESENTATION THEOREM OF TREE SERIES SOLUTION (Ⅰ)

In the analytic theory of differential equations the exact explicit analytic solution has not been obtained for equations of the non-Fuchsian type (Poincare's problem). The new theory proposed in this paper for the first time affords a general method of finding exact analytic expres-sion for irregular integrals.By discarding the assumption of formal solution of classical theory,our method consists in deriving a cor-respondence relation from the equation itself and providing the analytic structure of irregular integrals naturally by the residue theorem. Irregular integrals are made up of three parts: noncontracted part,represented by ordinary recursion series,all-and semi-contracted part by the so-called tree series. Tree series solutions belong to analytic function of the new kind with recursion series as the special case only.     

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.     

A novel non‐upwind,interconnected, multi‐grid,overlapping numerical procedure for problems involving fluid flow

A novel numerical procedure for heat, mass and momentum transfer in fluid flow is presented. The new scheme is passed on a non‐upwind, interconnected, multi‐grid, overlapping (NIMO) finite‐difference algorithm. In 2D flows, the NIMO algorithm solves finite‐difference equations for each dependent variable on four overlapping grids. The finite‐difference equations are formulated using the control‐volume approach, such that no interpolations are needed for computing the convective fluxes. For a particular dependent variable, four fields of values are produced. The NIMO numerical procedure is tested against the exact solution of two test problems. The first test problem is an oblique laminar 2D flow with a double step abrupt change in a passive scalar variable for infinite Peclet number. The second test problem is a rotating radial flow in an annular sector with a single step abrupt change in a passive scalar variable for infinite Peclet number. The NIMO scheme produced essentially the exact solution using different uniform and non‐uniform square and rectangular grids for 45 and 30° angle of inclination. All other schemes were unable to capture the exact solution, especially for the rectangular and non‐uniform grids. The NIMO scheme was also successful in predicting the exact solution for the rotating radial flow, using a uniform cylindrical‐polar coordinate grid. Copyright © 2008 John Wiley & Sons, Ltd.     

Radiation and scattering from bodies of revolution

The problem of electromagnetic radiation and scattering from perfectly conducting bodies of revolution of arbitrary shape is considered. The mathematical formulation is an integro-differential equation, obtained from the potential integrals plus boundary conditions at the body. A solution is effected by the method of moments, and the results are expressed in terms of generalized network parameters. The expansion functions chosen for the solution are harmonic in ø (azimuth angle) and subsectional in t (contour length variable). Because of rotational symmetry, the solution becomes a Fourier series in ø, each term of which is uncoupled to every other term.Illustrative computations are given for radiation from apertures and plane wave scattering from bodies of revolution. The impedance elements, currents, radiation patterns, and scattering patterns for a conducting sphere are computed both from the general solution and from the classical eigenfunction solution. The agreement obtained serves to check the general solution. Similar computations for a cone-sphere illustrate the application of the general solution to problems not solvable by classical methods.     

Fracture of a rectangular piezoelectromagnetic body

The singular stress, electric fields and magnetic fields in a rectangular piezoelectromagnetic body containing a center Griffith crack under longitudinal shear are obtained. Fourier transforms and Fourier sine series are used to reduce the mixed boundary value problems of the crack, which is assumed to be impermeable, to dual integral equations. The solution of the dual integral equations is then expressed in terms of Fredholm integral equations of the second kind. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Also obtained are the field intensity factors and the energy release rates. Numerical results obtained show that the geometry of the rectangular body have significant influence on the field intensity factors and the energy release rates.     

Slow steady flows of a conducting fluid at large Hartmann numbers

We consider slow steady flows of a conducting fluid at large values of the Hartmann number and small values of the magnetic Reynolds number in an inhomogeneous magnetic field. The general solution is obtained in explicit form for the basic portion (core) of the flow, where the inertia and viscous forces may be neglected. The boundary conditions which this solution must satisfy at the outer edges of the boundary layers which develop at the walls are considered. Possible types of discontinuity surfaces and other singularities in the flow core are examined. An exact solution is obtained for the problem of conducting fluid flow in a tube of arbitrary section in an inhomogeneous magnetic field.The content of this paper is a generalization of some results on flows in a homogeneous magnetic field, obtained in [1–8], to the case of arbitrary flows in an inhomogeneous magnetic field. The author's interest in the problems considered in this study was attracted by a report presented by Professor Shercliff at the Institute of Mechanics, Moscow State University, in May 1967, on the studies of English scientists on conducting fluid flows in a strong uniform magnetic field.     

A singular perturbation problem of laminar flow in a uniformly porous channel with large injection and with an applied transverse magnetic field

Summary The problem of the steady flow of an electrically conducting viscous fluid through porous walls of a channel in the presence of an applied transverse magnetic field is considered. A solution for the case of small M ²/R (where M = Hartmann number, R = suction Reynolds number) with large blowing at the walls has been given by Terrill and Shrestha [3]. Their solution, on differentiating three times, is found to become infinite at the centre of the channel. Physically this means that there must be a viscous layer at the centre of the channel and Terrill and Shrestha are neglecting the shear layer. In this paper the solution given by Terrill and Shrestha is extended by obtaining an extra term of the series of expansion and the method of inner and outer expansion is used to obtain the complete solution which includes the viscous layer. The resulting series solutions are confirmed by numerical results.     

Analytical and numerical study of double diffusive convection in a vertical enclosure

This paper presents an analytical and numerical study of natural convection of a double-diffusive fluid contained in a rectangular slot subject to uniform heat and mass fluxes along the vertical sides. Governing parameters of the problem under study are the thermal Rayleigh number, Ra _T ; buoyancy ratio, N; Lewis number, Le; Prandtl number, Pr and aspect ratio of the cavity, A. In the first part of the analytical study a scale analysis is applied to the two extreme cases of heat-transfer and mass-transfer-driven flows. In the second part, an analytical solution, based on the parallel flow approximation, is reported for tall enclosures (A?1). Solutions for the flow fields, temperature and concentration distributions and Nusselt and Sherwood numbers are obtained in terms of the governing parameters of the problem. In the limits of heat-driven and solute-driven flows a good agreement is obtained between the prediction of the scale analysis and those of the analytical solution. The numerical solutions are based on the complete governing equations for two-dimensional flows, and cover the range 1≤Ra _T ≤10⁷, 0≤N≤10⁵, 10^-3≤Le≤10³, 1≤A≤20 and Pr=7. A good agreement is found between the analytical predictions and the numerical simulation.