Computation of 3-D magnetostatic fields using a reduced scalarpotential

Some improvements to the finite element computation of static magnetic fields in three dimensions using a reduced magnetic scalar potential are presented. Methods are described for obtaining an edge element representation of the rotational part of the magnetic field from a given source current distribution. When the current distribution is not known in advance, a boundary value problem is set up in terms of a current vector potential. An edge element representation of the solution can be directly used in the subsequent magnetostatic calculation. The magnetic field in a DC arc furnace is calculated by first determining the current distribution in terms of a current vector potential. A 3-D problem involving a permanent magnet as well as a coil is solved, and the magnetic field in some points is compared with measurement results     

Numerical solution of Stokes flow generated by vortices: part 2, inside an elliptical cylinder

In this second part paper, the two-dimensional flow inside an elliptical cylinder is studied in the presence of no-slip boundary conditions. For simplicity, line vortices are assumed to be parallel to the elliptical cylinder axis, all axes in the same plane. The interior boundary value problem is solved in terms of a stream function. Numerical solutions for the flow field are obtained by application of the boundary element method. The streamline patterns are sketched for a number of special cases where the elliptical cylinder is either stationary or rotating about its own axis. In particular, some interesting flow patterns are observed in the parameter space which may have potential significance in studies of various flows. We also investigate the change in streamline topologies as the parameters are varied. Eddies of various sizes and shapes appear depending on the primary vortices and their locations. The results presented may be relevant for a variety of applications including vortex mixing. The analytical closed-form expressions for the single vortex inside an elliptical cylinder and double vortices inside circular a cylinder are found.     

三边夹紧一边铰支正交各向异性矩形薄板的几何非线性分析

利用Galerkin方法分析了von-Karman型三边夹紧一边铰支正交各向异性矩形板。所设的位移函数为梁的主振型函数，它不仅能精确地满足边界条件，而且具有正交的特性，从而把复杂的非齐次非线性偏微分方程组化为一组非线性代数方程组，通过非线性方程组的线性化，用稳定化双共轭梯度法求解稀疏矩阵线性方程组以及可调节参数的修正迭代法求解非线性代数方程组。实践证明，梁的主振型函数收敛很快，只须取出级数的前几项即可满足精度要求。最后求出了不同复合材料的挠度和应力值。     

Multi-domain dual reciprocity for the solution of inelastic non-Newtonian flow problems at low Reynolds number

An efficient boundary element solution of the motion of inelastic non-Newtonian fluids at low Reynolds number is presented in this paper. For the numerical solution all the domain integrals of the boundary element formulation have been transformed into equivalent boundary integrals by means of the dual reciprocity method (DRM). To achieve an accurate approximation of the non-linear and non-Newtonian terms two major improvements have been made to the DRM, namely the use of augmented thin plate splines as interpolation functions, and the partition of the entire domain into smaller subregions or domain decomposition. In each subregion or domain element the DRM was applied together with some additional equations that ensure continuity on the interfaces between adjacent subdomains. After applying the boundary conditions the final systems of equations will be sparse and the approximation of the nonlinear terms will be more localised than in the traditional DRM. This new method known as multidomain dual reciprocity (MD-DRM) has been used to solve several non-Newtonian problems including the pressure driven flow of a power law fluid, the Couette flow and two simulations of industrial polymer mixers. Received 7 February 2001     

A coupled numerical scheme of dual reciprocity BEM with DQM for the transient elastodynamic problems

The two‐dimensional transient elastodynamic problems are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in spatial domain with the differential quadrature method (DQM) in time domain. The DRBEM with the fundamental solution of the Laplace equation transforms the domain integrals into the boundary integrals that contain the first‐ and the second‐order time derivative terms. Thus, the application of DRBEM to elastodynamic problems results in a system of second‐order ordinary differential equations in time. This system is then discretized by the polynomial‐based DQM with respect to time, which gives a system of linear algebraic equations after the imposition of both the boundary and the initial conditions. Therefore, the solution is obtained at any required time level at one stroke without the use of an iterative scheme and without the need of very small step size in time direction. The numerical results are visualized in terms of graphics. Copyright © 2008 John Wiley & Sons, Ltd.     

Dual boundary element analysis of cracked plates: singularity subtraction technique

The present paper is concerned with the formulation of the singularity subtraction technique in the dual boundary element analysis of the mixed-mode deformation of general homogeneous cracked plates.The equations of the dual boundary element method are the displacement and the traction boundary integral equations. When the displacement equation is applied on one of the crack surfaces and the traction equation is applied on the other, general mixed-mode crack problems can be solved in a single region boundary element formulation, with both crack surfaces discretized with discontinuous quadratic boundary elements.The singularity subtraction technique is a regularization procedure that uses a singular particular solution of the crack problem to introduce the stress intensity factors as additional problem unknowns. The single-region boundary element analysis of a general crack problem restricts the availability of singular particular solutions, valid in the global domain of the problem. A modelling strategy, that considers an automatic partition of the problem domain in near-tip and far-tip field regions, is proposed to overcome this difficulty. After the application of the singularity subtraction technique in the near-tip field regions, regularized locally with the singular term of the Williams' eigenexpansion, continuity is restored with equilibrium and compatibility conditions imposed along the interface boundaries. The accuracy and efficiency of the singularity subtraction technique make this formulation ideal for the study of crack growth problems under mixed-mode conditions.     

A new finite element technique for the analysis of steady viscous flow problems

A new finite element technique for two-dimensional viscous incompressible fluid flow problems is presented in this paper. The vorticity transport equation is integrated in a small control volume, which results in the conservation law of vorticity. The finite element technique is applied to this equation together with the continuity equation, where simple linear triangular elements with three nodes are used for the formulation. Resulting sets of algebraic equations are solved by the use of a kind of relaxation method. Numerical results for viscous flow past a cavity show good agreement with experimental results.     

DRBEM solution of the thermo-solutal buoyancy induced mixed convection flow problems

Mixed convective heat and solutal transport is important in engineering applications such as nuclear waste disposal, crystal growth and oceanography. Mixed convective in a lid-driven cavity and through channels with backward-facing step is studied by solving the equations of conservation of mass, momentum, energy and solutal concentration numerically. The governing equations are solved by using dual reciprocity boundary element method (DRBEM) with constant elements in terms of stream function, vorticity, temperature and concentration. Vorticity, energy and concentration equations are transformed to the form of modified Helmholtz equations by utilizing forward difference with relaxation parameters for the time derivatives, and also approximating Laplacian terms at two consecutive time levels. The DRBEM application is carried out with the fundamental solution of modified Helmholtz equation, resulting in linear algebraic systems for the time dependent unknowns. Inhomogeneities in modified Helmholtz equations are approximated with the thin plated radial basis functions. Computations are carried out for several values of Richardson number, buoyancy ratio and Reynolds number. When the temperature and solutal concentration boundary conditions are changed, the thermal and solutal buoyancy forces can oppose or aid each other. The effects of these parameters on the flow behavior and heat transfer are shown in terms of graphics.     

The boundary element method for the solution of Stokes equations in two-dimensional domains

A boundary element method for the solution of Stokes equations governing creeping flow or Stokes flow in the interior of an arbitrary two-dimensional domain is presented. A procedure for introducing pressure data on the boundary of the domain is also included and the integral coefficients of the resulting linear algebraic equations are evaluated analytically. Calculations are performed in a circular domain using a variety of different boundary conditions, including a combination of the fluid velocity and the pressure. Results are presented both on the boundary and inside the solution domain in order to illustrate that the boundary element method developed here provides an efficient technique, in terms of accuracy and convergence, to investigate Stokes flow numerically.     

Turbulent biomagnetic fluid flow in a rectangular channel under the action of a localized magnetic field

The fundamental problem of the turbulent flow of a biomagnetic fluid (blood) between two parallel plates under the action of a localized magnetic field is studied. The blood is considered to be an electrically conducting, incompressible and Newtonian fluid and its flow is steady, two-dimensional and turbulent. The turbulent flow is described by the Reynolds averaged Navier–Stokes (RANS) equations. For the numerical solution of the problem under consideration, which is described by a coupled and non-linear system of PDEs, with appropriate boundary conditions, the stream function–vorticity formulation is used. For the eddy-kinematic viscosity, the low Reynolds number k–ε turbulence model is adopted. The solution of the problem, for different values of the dimensionless parameter entering into it, is obtained by developing and applying an efficient numerical technique based on finite differences scheme. Results concerning the velocity and temperature field, skin friction and rate of heat transfer, indicate that the presence of the localized magnetic field, appreciable influences the turbulent flow field. A comparison is also made with the corresponding laminar flow, indicating that the influence of the magnetic field decreases in the presence of turbulence.     

Radiation and thermal diffusion effects on an unsteady MHD free convection mass-transfer flow past an infinite vertical porous plate with the hall current and a heat source

An analysis is performed to study the effects of radiation and thermal diffusion on an unsteady MHD free convection heat- and mass-transfer flow of an incompressible, electrically conducting, viscous fluid past an infinite vertical porous plate with the Hall current and a heat source. The flow is considered under the influence of a constant suction velocity and a uniform magnetic field applied normally to the flow. The dimensionless governing equations are solved numerically by the Galerkin finite element method. The effects of the flow parameters on the primary and secondary velocities, temperature, species concentration, shearing stresses, Nusselt number, and Sherwood number are calculated and presented in figures and tables. The results obtained show that a decrease in the temperature boundary layer thickness occurs when the Prandtl number and radiation parameter are increased and an increase in the Schmidt number leads to a decrease in the concentration boundary layer thickness.     

The dual boundary element method for thermoelastic crack problems

A boundary element formulation, which does not require domain discretization and allows a single region analysis, is presented for steady-state thermoelastic crack problems. The problems are solved by the dual boundary element method which uses displacement and temperature equations on one crack surface and traction and flux equations on the other crack surface. The domain integrals are transformed to boundary integrals using the Galerkin technique. Stress intensity factors are calculated using the path independent -integral. Several numerical problems are solved and the results are compared, where possible, with existing solutions.     

Magnetohydrodynamic heat transfer over a non-isothermal stretching sheet

Summary This paper presents solutions of the energy equation for the boundary layer flow of an electrically conducting fluid under the influence of a constant transverse magnetic field over a linearly stretching non-isothermal flat sheet. Effects due to dissipation, stress work and heat generation are considered. Analytical solutions of the resulting linear nonhomogeneous boundary value problems, expressed in terms of Kummer's functions, are presented for the case of prescribed surface temperature as well as the case of prescribed wall heat flux, both of which are assumed to be quadratic functions of distance. The boundary value problems are also solved by direct numerical integration yielding results in excellent agreement with the analytical solutions.     

Hydromagnetic couette flow of an Oldroyd-B fluid in a rotating system

The motion of electrically conducting, Oldroyd-B and incompressible fluid between two infinitely extended non-conducting parallel plates under a uniform transverse magnetic field, fixed relative to the fluid has been considered. The lower plate is at rest and the upper plate is oscillating in its own plane. The governing partial differential equation of this problem, subject to boundary conditions are solved analytically. The expressions for the steady and unsteady velocity fields for the conducting Oldroyd-B fluid are obtained. The graphs are plotted for different values of dimensionless parameters of the problem and the analysis of the results showed that the flow field is appreciably influenced by the applied magnetic field, the rotation and the material parameters of the fluid.     

On the flow of an electrically conducting nonlocal viscous fluid between parallel plates in the presence of a transverse magnetic field in magnetohydrodynamics

A continuum theory is constructed for the flow of an electrically conducting nonlocal viscous fluid between two nonconducting parallel plates. The flow is subject to the influence of a transverse magnetic field. The effects of long range or nonlocal interactions at a material point in the fluid arising from all material points in the rest of the fluid are taken into account by means of a nonlocal influence function. Equations of motion governing the nonlocal viscous flow are derived from localized forms of global balance laws and constitutive equations appropriate for electromagnetically active media. These field equations are analytically solved for the nonlocal velocity and the nonlocal stress fields. The effects of varying the magnetic field strength on the shear stress are investigated. The effects of such variations on the shear stress exerted on the walls of microscopic channels are also determined. Numerical computations are provided for these results.     

Transverse high gradient magnetic filter cell with bounded flowfield

The capture of fine paramagnetic particles from a fluid suspension in a magnetic filter element of a novel design is analyzed. Unlike the systems previously analyzed, in our model the flow is bounded by two-by-two parallel planar plates, and the ferromagnetic wires are installed outside these spaces, within planes parallel with the plates. Our analysis is based on the study of particle trajectories, considering the laminar flow of carrier fluid, From these we establish the conditions for the maximum recovery of the particles in suspension. This study is useful in designing magnetic filter batteries with corrosion-protected ferromagnetic wires     

驻波场中圆柱形换热管外声流特性的数值研究

基于非线性Navier-Stokes方程推导了一般性的声流控制方程,并利用有限元软件COMSOL对平面驻波声场中单根圆柱形换热管和双圆柱形换热管外包含二阶流场信息的Navier-Stokes方程进行全域求解,得到了一阶和二阶流场信息。研究了雷诺数Re和斯特劳哈尔数Sr对换热管外声流结构的影响规律。研究发现:随Re和Sr的增大,换热管边界层内的涡流区域变小,直至消失,而边界层外的涡流区域逐渐增大;且单换热管外涡流个数由8个减少到4个,双换热管外涡流个数从12个减少到8个。此外,边界层内和边界层外的涡流结构呈反向旋转;边界层外流场沿着振荡方向远离换热管,而在垂直方向流向换热管。     

A note on the Blasius and Sakiadis flow of a non-Newtonian power-law fluid in a constant transverse magnetic field

The Blasius and Sakiadis flows of a non-Newtonian power-law conducting fluid under the effect of a constant transverse magnetic field is considered. The boundary layer equations are transformed into a non-dimensional form and a new dimensionless magnetic parameter is introduced. The transformed boundary layer equations are solved with a finite difference method. Both Blasius and Sakiadis flows reach an asymptotic state and become one-dimensional with the increase of the coordinate in the streamwise direction. The flow behavior in the intermediate region is studied and exact analytical solutions have been found for the asymptotic state. The characteristics of physical and engineering interest are discussed in the paper in detail.     

Pinning of planar vortices in a three-dimensional Josephson medium and the applicability of the Bean model

Calculations are made of a system of planar vortices formed at the boundary of a sample with a monotonically increasing external magnetic field with allowance for pinning caused by the cellular nature of the medium. The spacing between the vortices remains approximately constant and does not decrease near the boundary. This behavior indicates that the magnetic field does not depend on the coordinate in the region where the vortices have penetrated. These observations contradict the conventional Bean model. Pis’ma Zh. Tekh. Fiz. 23, 66–69 (June 12, 1997)     

Solution of magnetohydrodynamic flow problems using the boundary element method

A boundary element solution is implemented for magnetohydrodynamic (MHD) flow problem in ducts with several geometrical cross-section with insulating walls when a uniform magnetic field is imposed perpendicular to the flow direction. The coupled velocity and induced magnetic field equations are first transformed into uncoupled inhomogeneous convection–diffusion type equations. After introducing particular solutions, only the homogeneous equations are solved by using boundary element method (BEM). The fundamental solutions of the uncoupled equations themselves (convection–diffusion type) involve the Hartmann number (M) through exponential and modified Bessel functions. Thus, it is possible to obtain results for large values of M (M≤300) using only the simplest constant boundary elements. It is found that as the Hartmann number increases, boundary layer formation starts near the walls and there is a flattening tendency for both the velocity and the induced magnetic field. Also, velocity becomes uniform at the center of the duct. These are the well-known behaviours of MHD flow. The velocity and the induced magnetic field contours are graphically visualized for several values of M and for different geometries of the duct cross-section.