Semi‐analytical solution of MHD flow through boundary integrals on the pipe wall

A mathematical model is given for the magnetohydrodynamic (MHD) pipe flow as an inner Dirichlet problem in a 2D circular cross section of the pipe, coupled with an outer Dirichlet or Neumann magnetic problem. Inner Dirichlet problem is given as the coupled convection‐diffusion equations for the velocity and the induced current of the fluid coupling also to the outer problem, which is defined with the Laplace equation for the induced magnetic field of the exterior region with either Dirichlet or Neumann boundary condition. Unique solution of inner Dirichlet problem is obtained theoretically reducing it into two boundary integral equations defined on the boundary by using the corresponding fundamental solutions. Exterior solution is also given theoretically on the pipe wall with Poisson integral, and it is unique with Dirichlet boundary condition but exists with an additive constant obtained through coupled boundary and solvability conditions in Neumann wall condition. The collocation method is used to discretize these boundary integrals on the pipe wall. Thus, the proposed procedure is an improved theoretical analysis for combining the solution methods for the interior and exterior regions, which are consolidated numerically showing the flow behavior. The solution is simulated for several values of problem parameters, and the well‐known MHD characteristics are observed inside the pipe for increasing values of Hartmann number maintaining the continuity of induced currents on the pipe wall.     

The method of lines for Poisson's equation with nonlinear or free boundary conditions

Summary The method of lines is used to solve Poisson's equation on an irregular domain with nonlinear or free boundary conditions. The partial differential equation is approximated by a system of second order ordinary differential equations subject to multi-point boundary conditions. The system is solved with an SOR iteration which employs invariant imbedding for each one dimensional problem. An application of the method to a boundary control problem and to a free surface problem arising in electrochemical machining is described. Finally, some theoretical convergence results are presented for a model problem with radiative boundary conditions on fixed boundaries.This work was supported by the U.S. Army Research Office under Grant DA-AG29-76-G-0261     

Lie group analysis of unsteady MHD three dimensional by natural convection from an inclined stretching surface saturated porous medium

Lie group method is investigated for solving the problem of heat transfer in an unsteady, three-dimensional, laminar, boundary-layer flow of a viscous, incompressible and electrically conducting fluid over inclined permeable surface embedded in porous medium in the presence of a uniform magnetic field and heat generation/absorption effects. A uniform magnetic field is applied in the y-direction and a generalized flow model is presented to include the effects of the macroscopic viscous term and the microscopic permeability of porous medium. The infinitesimal generators accepted by the equations are calculated and the extension of the Lie algebra for the problem is also presented. The restrictions imposed by the boundary conditions on the generators are calculated. The investigation of the three-independent-variable partial differential equations is converted into a two-independent-variable system by using one subgroup of the general group. The resulting equations are solved numerically with the perturbation solution for various times. Velocity, temperature and pressure profiles, surface shear stresses, and wall-heat transfer rate are discussed for various values of Prandtl number, Hartmann number, Darcy number, heat generation/absorption coefficient, and surface mass-transfer coefficient.     

On the energetic Galerkin boundary element method applied to interior wave propagation problems

We consider two-dimensional interior wave propagation problems with vanishing initial and mixed boundary conditions, reformulated as a system of two boundary integral equations with retarded potential. These latter are then set in a weak form, based on a natural energy identity satisfied by the solution of the differential problem, and discretized by the related energetic Galerkin boundary element method. Numerical results are presented and discussed.     

On an iterative algorithm with superquadratic convergence for solving nonlinear operator equations

We study an iterative method with order for solving nonlinear operator equations in Banach spaces. Algorithms for specific operator equations are built up. We present the received new results of the local and semilocal convergence, in case when the first-order divided differences of a nonlinear operator are Hölder continuous. Moreover a quadratic nonlinear majorant for a nonlinear operator, according to the conditions laid upon it, is built. A priori and a posteriori estimations of the method’s error are received. The method needs almost the same number of computations as the classical Secant method, but has a higher order of convergence. We apply our results to the numerical solving of a nonlinear boundary value problem of second-order and to the systems of nonlinear equations of large dimension.     

Error estimates for MAC-like approximations to the linear Navier-Stokes equations

Summary In this paper a priori error estimates are derived for the discretization error which results when the linear Navier-Stokes equations are solved by a method which closely resembles the MAC-method of Harlow and Welch. General boundary conditions are permitted and the estimates are in terms of the discreteL ² norm. A solvability result is given which also applies to a generalization of the method to the nonlinear case. This generalization is used in the last section to produce a numerical solution to the problem of flow around an obstacle.This work supported in part by Westinghouse Nuclear Energy Systems. Research Report #76-13     

On the penalty-projection method for the Navier–Stokes equations with the MAC mesh

We deal with the time-dependent Navier–Stokes equations (NSE) with Dirichlet boundary conditions on the whole domain or, on a part of the domain and open boundary conditions on the other part. It is shown numerically that combining the penalty-projection method with spatial discretization by the Marker And Cell scheme (MAC) yields reasonably good results for solving the above-mentioned problem. The scheme which has been introduced combines the backward difference formula of second-order (BDF2, namely Gear’s scheme) for the temporal approximation, the second-order Richardson extrapolation for the nonlinear term, and the penalty-projection to split the velocity and pressure unknowns. Similarly to the results obtained for other projection methods, we estimate the errors for the velocity and pressure in adequate norms via the energy method.     

Turbulent boundary layer equationsÉquations de la couche limite turbulente

We study a boundary layer problem for the Navier–Stokes-alpha model obtaining a generalization of the Prandtl equations which we conjecture to represent the averaged flow in a turbulent boundary layer. We study the equations for the semi-infinite plate, both theoretically and numerically. Solutions agree with some experimental data in a part of the turbulent boundary layer. To cite this article: A. Cheskidov, C. R. Acad. Sci. Paris, Ser. I 334 (2002) 423–427.     

Acoustic gravity waves: a computational approach

This paper discusses numerical solutions of a hyperbolic initial boundary value problem that arises from acoustic wave propagation in the atmosphere. Field equations are derived from the atmospheric fluid flow governed by the Euler equations. The resulting original problem is nonlinear. A first-order linearized version of the problem is used for computational purposes. The main difficulty in the problem as with any open boundary problem is in obtaining stable boundary conditions. Approximate boundary conditions are derived and shown to be stable. Numerical results are presented to verify the effectiveness of these boundary conditions.     

On the implementation of boundary conditions for the method of lines

The method of lines for difference approximations of hyperbolic first order systems of partial differential equations is analyzed. The approximations are based on strictly semibounded difference operators including high order ones. The formulation of the ODE-system requires that the implementation of the boundary conditions is done carefully. We shall illustrate how different ways of implementation give rise to different stability properties. In particular, we derive a way of implementation that leads to an approximation that is strongly stable. It has been an open problem, whether for semidiscrete approximations with this strong stability property, the timestep for the ODE-solver is governed by the Cauchy problem. We present a counterexample showing that it is not. The analysis presented in this paper also serves as an illustration of the significant difference between different stability concepts.     

A modified accelerated monotone iterative method for finite difference reaction-diffusion-convection equations

This paper is concerned with monotone algorithms for the finite difference solutions of a class of nonlinear reaction-diffusion-convection equations with nonlinear boundary conditions. A modified accelerated monotone iterative method is presented to solve the finite difference systems for both the time-dependent problem and its corresponding steady-state problem. This method leads to a simple and yet efficient linear iterative algorithm. It yields two sequences of iterations that converge monotonically from above and below, respectively, to a unique solution of the system. The monotone property of the iterations gives concurrently improving upper and lower bounds for the solution. It is shown that the rate of convergence for the sum of the two sequences is quadratic. Under an additional requirement, quadratic convergence is attained for one of these two sequences. In contrast with the existing accelerated monotone iterative methods, our new method avoids computing local maxima in the construction of these sequences. An application using a model problem gives numerical results that illustrate the effectiveness of the proposed method.     

The modeling and numerical simulation of gas flow networks

Summary. A network formulation is introduced for the modeling and numerical simulation of complex gas transmission systems like a multi-cylinder internal combustion engine. Several simulation levels are discussed which result in different network representations of a specific system. Basic elements of a network are chambers of finite volume, straight pipes and connections like valves or nozzles. The pipe flow is modeled by the unsteady, one-dimensional Euler equations of gas dynamics. Semi-empirical approaches for the chambers and the connections yield differential-algebraic equations (DAEs) in time. The numerical solution is based on a TVD scheme for the pipe equations and a predictor-corrector method for the DAE-system. Simulation results for an internal combustion engine demonstrate the practical interest of the new approach. Received May 12, 1994 / Revised version received August 26, 1994     

The splitting finite-difference time-domain methods for Maxwell's equations in two dimensions

In this paper, we consider splitting methods for Maxwell's equations in two dimensions. A new kind of splitting finite-difference time-domain methods on a staggered grid is developed. The corresponding schemes consist of only two stages for each time step, which are very simple in computation. The rigorous analysis of the schemes is given. By the energy method, it is proved that the scheme is unconditionally stable and convergent for the problems with perfectly conducting boundary conditions. Numerical dispersion analysis and numerical experiments are presented to show the efficient performance of the proposed methods. Furthermore, the methods are also applied to solve a scattering problem successfully.     

Hydromagnetic rotating flow of a fourth‐order fluid past a porous plate

The rotating flow in the presence of a magnetic field is a problem belonging to hydromagnetics and deserves to be more widely studied than it has been to date. In the non‐linear regime the literature is scarce. We develop the governing equations for the unsteady hydromagnetic rotating flow of a fourth‐order fluid past a porous plate. The steady flow is governed by a boundary value problem in which the order of differential equations is more than the number of available boundary conditions. It is shown that by augmenting the boundary conditions based on asymptotic structures at infinity it is possible to obtain numerical solutions of the nonlinear hydromagnetic equations. Effects of uniform suction or blowing past the porous plate, exerted magnetic field and rotation on the flow phenomena, especially on the boundary layer structure near the plate, are numerically analysed and discussed. The flow behaviours of the Newtonian fluid and second‐, third‐ and fourth‐order non‐Newtonian fluids are compared for the special flow problem, respectively. Copyright © 2004 John Wiley & Sons, Ltd.     

流体诱发水平悬臂输液管的内共振和模态转换(Ⅰ)

运用牛顿法导出水平悬臂刚性输液管的非线性动力学数学模型．为了对该模型进行理论分析,通过对各个相关实际物理量的量级定性分析,给出了模型中各个物理参数的量级．在此基础上,应用多尺度法首先得到输液管自由振动模态的特征函数,利用悬臂管的边界条件给出了特征值满足的特征方程, 发现管内流体速度可以诱发第一阶模态和第二阶模态3种形式的内共振分别是3∶1、2∶1和1∶1内共振, 从理论上解释了流速诱发水平悬臂输液管系统内共振的机理．由于3∶1内共振所对应的流速最小,因此这种形式的内共振是最先出现的,也是最重要的．     

Numerical study of the flow in a three-dimensional thermally driven cavity

Solutions for the fully compressible Navier–Stokes equations are presented for the flow and temperature fields in a cubic cavity with large horizontal temperature differences. The ideal-gas approximation for air is assumed and viscosity is computed using Sutherland's law. The three-dimensional case forms an extension of previous studies performed on a two-dimensional square cavity. The influence of imposed boundary conditions in the third dimension is investigated as a numerical experiment. Comparison is made between convergence rates in case of periodic and free-slip boundary conditions. Results with no-slip boundary conditions are presented as well. The effect of the Rayleigh number is studied.     

Solvability of a nonstationary problem of a rigid body motion in the flow of a viscous incompressible fluid in a pipe of arbitrary cross-section

The existence of a generalized weak solution is proved for the nonstationary problem of motion of a rigid body in the flow of a viscous incompressible fluid filling a cylindrical pipe of arbitrary cross-section. The fluid flow conforms to the Navier–Stokes equations and tends to the Poiseuille flow at infinity. The body moves in accordance with the laws of classical mechanics under the influence of the surrounding fluid and the gravity force directed along the cylinder. Collisions of the body with the boundary of the flow domain are not admitted and, by this reason, the problem is considered until the body approaches the boundary.     

On general boundary-value problems for nonlinear elliptic equations of second order in a multiply connected plane domain

This paper deals with some general irregular oblique derivative problems for nonlinear uniformly elliptic equations of second order in a multiply connected plane domain. Firstly, we state the well-posedness of a new set of modified boundary conditions. Secondly, we verify the existence of solutions of the modified boundary-value problem for harmonic functions, and then prove the solvability of the modified problem for nonlinear elliptic equations, which includes the original boundary-value problem (i.e. boundary conditions without involving undertermined functions data). Here, mainly, the location of the zeros of analytic functions, a priori estimates for solutions and the continuity method are used in deriving all these results. Furthermore, the present approach and setting seems to be new and different from what has been employed before.The research was partially supported by a UPGC Grant of Hong Kong.