Flow of a biomagnetic viscoelastic fluid: application to estimation of blood flow in arteries during electromagnetic hyperthermia,a therapeutic procedure for cancer treatment

The paper deals with the theoretical investigation of a fundamental problem of biomaguetic fluid flow through a porous medium subject to a magnetic field by using the principles of biomagnetic fluid dynamics （BFD）. The study pertains to a situation where magnetization of the fluid varies with temperature. The fluid is considered to be non-Newtonian, whose flow is governed by the equation of a second-grade viscoelastic fluid. The walls of the channel are assumed to be stretchable, where the surface velocity is proportional to the longitudinal distance from the origin of coordinates. The problem is first reduced to solving a system of coupled nonlinear differential equations involving seven parameters. Considering blood as a biomagnetic fluid and using the present analysis, an attempt is made to compute some parameters of the blood flow by developing a suitable numerical method and by devising an appropriate finite difference scheme. The computational results are presented in graphical form, and thereby some theoretical predictions are made with respect to the hemodynamical flow of the blood in a hyperthermal state under the action of a magnetic field. The results clearly indicate that the presence of a magnetic dipole bears the potential so as to affect the characteristics of the blood flow in arteries to a significant extent during the therapeutic procedure of electromagnetic hyperthermia. The study will attract the attention of clinicians, to whom the results would be useful in the treatment of cancer patients by the method of electromagnetic hyperthermia.     

HOPF BIFURCATION OF A NONLINEAR RESTRAINED CURVED PIPE CONVEYING FLUID BY DIFFERENTIAL QUADRATURE METHOD

This paper proposes a new method for investigating the Hopf bifurcation of a curved pipe conveying fluid with nonlinear spring support. The nonlinear equation of motion is derived by forces equilibrium on microelement of the system under consideration. The spatial coordinate of the system is discretized by the differential quadrature method and then the dynamic equation is solved by the Newton-Raphson method. The numerical solutions show that the inner fluid velocity of the Hopf bifurcation point of the curved pipe varies with different values of the parameter,nonlinear spring stiffness. Based on this, the cycle and divergent motions are both found to exist at specific fluid flow velocities with a given value of the nonlinear spring stiffness. The results are useful for further study of the nonlinear dynamic mechanism of the curved fluid conveying pipe.     

Differential transformation method for studying flow and heat transfer due to stretching sheet embedded in porous medium with variable thickness,variable thermal conductivity,and thermal radiation

This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.     

3D numerical simulation on fluid-structure interaction of structure subjected to underwater explosion with cavitation

In the underwater-shock environment,cavitation occurs near the structural surface.The dynamic response of fluid-structure interactions is influenced seriously by the cavitation effects.It is also the difficulty in the field of underwater explosion.With the traditional boundary element method and the finite element method(FEM),it is difficult to solve the nonlinear problem with cavitation effects subjected to the underwater explosion.To solve this problem,under the consideration of the cavitation effects and fluid compressibility,with fluid viscidity being neglected,a 3D numerical model of transient nonlinear fluid-structure interaction subjected to the underwater explosion is built.The fluid spectral element method(SEM) and the FEM are adopted to solve this model.After comparison with the FEM,it is shown that the SEM is more precise than the FEM,and the SEM results are in good coincidence with benchmark results and experiment results.Based on this,combined with ABAQUS,the transient fluid-structure interaction mechanism of the 3D submerged spherical shell and ship stiffened plates subjected to the underwater explosion is discussed,and the cavitation region and its influence on the structural dynamic responses are presented.The paper aims at providing references for relevant research on transient fluid-structure interaction of ship structures subjected to the underwater explosion.     

Effects of thermophoresis particle deposition and of the thermal conductivity in a porous plate with dissipative heat and mass transfer

Network simulation method(NSM) is used to solve the laminar heat and mass transfer of an electricallyconducting,heat generating/absorbing fluid past a perforated horizontal surface in the presence of viscous and Joule heating problem. The governing partial differential equations are non-dimensionalized and transformed into a system of nonlinear ordinary differential similarity equations,in a single independent variable,η. The resulting coupled,nonlinear equations are solved under appropriate transformed boundary conditions. Computations are performed for a wide range of the governing flow parameters,viz Prandtl number,thermophoretic coeffcient(a function of Knudsen number),thermal conductivity parameter,wall transpiration parameter and Schmidt number. The numerical details are discussed with relevant applications. The present problem finds applications in optical fiber fabrication,aerosol filter precipitators,particle deposition on hydronautical blades,semiconductor wafer design,thermo-electronics and problems including nuclear reactor safety.     

A NUMERICAL METHOD FOR SIMULATING NONLINEAR FLUID-RIGID STRUCTURE INTERACTION PROBLEMS

A numerical method for simulating nonlinear fluid-rigid structure interaction problems is developed. The structure is assumed to undergo large rigid body motions and the fluid flow is governed by nonlinear, viscous or non-viscous, field equations with nonlinear boundary conditions applied to the free surface and fluid-solid interaction interfaces. An Arbitrary-Lagrangian-Eulerian (ALE) mesh system is used to construct the numerical model. A multi-block numerical scheme of study is adopted allowing for the relative motion between moving overset grids, which are independent of one another. This provides a convenient method to overcome the difficulties in matching fluid meshes with large solid motions. Nonlinear numerical equations describing nonlinear fluid-solid interaction dynamics are derived through a numerical discretization scheme of study. A coupling iteration process is used to solve these numerical equations. Numerical examples are presented to demonstrate applications of the model developed.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

Three-dimensional elasticity solutions for bending of generally supported thick functionally graded plates

Three-dimensional elasticity solutions for static bending of thick functionally graded plates are presented using a hybrid semi-analytical approach-the state-space based differential quadrature method （SSDQM）. The plate is generally supported at four edges for which the two-way differential quadrature method is used to solve the in-plane variations of the stress and displacement fields numerically. An approximate laminate model （ALM） is exploited to reduce the inhomogeneous plate into a multi-layered laminate, thus applying the state space method to solve analytically in the thickness direction. Both the convergence properties of SSDQM and ALM are examined. The SSDQM is validated by comparing the numerical results with the exact solutions reported in the literature. As an example, the Mori-Tanaka model is used to predict the effective bulk and shear moduli. Effects of gradient index and aspect ratios on the bending behavior of functionally graded thick plates are investigated.     

Numerical stabilities of loosely coupled methods for robust modeling of lightweight and flexible structures in incompressible and viscous flows

The growing interest to examine the hydroelastic dynamics and stabilities of lightweight and flexible materials requires robust and accurate fluid–structure interaction(FSI)models. Classically, partitioned fluid and structure solvers are easier to implement compared to monolithic methods;however, partitioned FSI models are vulnerable to numerical("virtual added mass") instabilities for cases when the solid to fluid density ratio is low and if the flow is incompressible.As a partitioned method, the loosely hybrid coupled(LHC)method, which was introduced and validated in Young et al.(Acta Mech. Sin. 28:1030–1041, 2012), has been successfully used to efficiently and stably model lightweight and flexible structures. The LHC method achieves its numerical stability by, in addition to the viscous fluid forces, embedding potential flow approximations of the fluid induced forces to transform the partitioned FSI model into a semi-implicit scheme. The objective of this work is to derive and validate the numerical stability boundary of the LHC. The results show that the stability boundary of the LHC is much wider than traditional loosely coupled methods for a variety of numerical integration schemes. The results also show that inclusion of an estimate of the fluid inertial forces is the most critical to ensure the numerical stability when solving for fluid–structure interaction problems involving cases with a solid to fluid-added mass ratio less than one.     

Parametric variational solution of linear-quadratic optimal control problems with control inequality constraints

A parametric variational principle and the corresponding numerical algo- rithm are proposed to solve a linear-quadratic （LQ） optimal control problem with control inequality constraints. Based on the parametric variational principle, this control prob- lem is transformed into a set of Hamiltonian canonical equations coupled with the linear complementarity equations, which are solved by a linear complementarity solver in the discrete-time domain. The costate variable information is also evaluated by the proposed method. The parametric variational algorithm proposed in this paper is suitable for both time-invariant and time-varying systems. Two numerical examples are used to test the validity of the proposed method. The proposed algorithm is used to astrodynamics to solve a practical optimal control problem for rendezvousing spacecrafts with a finite low thrust. The numerical simulations show that the parametric variational algorithm is ef- fective for LQ optimal control problems with control inequality constraints.     

Upwind local differential quandrature method for solving coupled viscous flow and heat transfer equations

The differential quadrature method (DQM) has been applied successfully to solve numerically many problems in the fluid mechanics. But it is only limited to the flow problems in regular regions. At the same time, here is no upwind mechanism to deal with the convective property of the fluid flow in traditional DQ method. A local differential quadrature method owning upwind mechanism (ULDQM) was given to solve the coupled problem of incompressible viscous flow and heat transfer in an irregular region. For the problem of flow past a contraction channel whose boundary does not parallel to coordinate direction, the satisfactory numerical solutions were obtained by using ULDQM with a few grid points. The numerical results show that the ULDQM possesses advantages including well convergence, less computational workload and storage as compared with the low-order finite difference method. Foundation item: the Municipal Key Subject Programs of Shanghai Biographies: A. S. J. Al-Saif (1964 ~)     

A mixture differential quadrature method for solving two-dimensional incompressible Navier-Stokes equations

Ｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ（ＤＱＭ）ｐｒｏｐｏｓｅｄｂｙＲ．Ｂｅｌｌｍａｎ［１,２］ｈａｓｂｅｅｎｓｕｃｃｅｓｓｆｕｌｌｙｅｍｐｌｏｙｅｄｉｎｎｕｍｅｒｉｃａｌｃｏｍｐｕｔａｔｉｏｎｓｏｆｐｒｏｂｌｅｍｓｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｐｈｙｓｉｃａｌｓｃｉｅｎｃｅ．ＢｅｃａｕｓｅｔｈｅｉｎｆｏｒｍａｔｉｏｎｏｎａｌｌｇｒｉｄｐｏｉｎｔｓｉｓｕｓｅｄｔｏｆｉｔｔｈｅｄｅｒｉｖａｔｉｖｅｓａｔｇｒｉｄｐｏｉｎｔｓｉｎｔｈｅＤＱＭ,ｉｔｉｓｅｎｏｕｇｈｔｏｏｂｔａ…     

稀薄流非线性模型方程离散速度坐标法有限差分解

从一般非线性Bo ltzm ann方程出发,发展并实现了一套适于大范围K nudsen数稀薄流问题数值模拟的统一算法。采用BGK模型和Shakov模型近似碰撞项,进而引入两个二速度无量纲简化分布函数,通过关于分子速度第三分量取矩积分,将三速度单一模型方程变换为二速度微分方程组。基于G auss-H erm ite积分公式和正交多项式G auss积分公式,借助离散速度坐标法消除简化模型方程对分子速度空间的连续依赖性,从相空间到物理空间得到一组带源项双曲守恒离散方程,并给出其显式和隐式二阶迎风TVD有限差分解。以二维圆柱A r气体超声速绕流算例,验证了数值算法的有效性,比较分析了漫反射和镜面反射两种气体分子壁面反射模型的计算结果。     

A decoupling numerical method for fluid flow

A first-order non-conforming numerical methodology, Separation method, for fluid flow problems with a 3-point exponential interpolation scheme has been developed. The flow problem is decoupled into multiple one-dimensional subproblems and assembled to form the solutions. A fully staggered grid and a conservational domain centred at the node of interest make the decoupling scheme first-order-accurate. The discretization of each one-dimensional subproblem is based on a 3-point interpolation function and a conservational domain centred at the node of interest. The proposed scheme gives a guaranteed first-order accuracy. It is shown that the traditional upwind (or exponentially weighted upstream) scheme is less than first-order-accurate. The pressure is decoupled from the velocity field using the pressure correction method of SIMPLE. Thomas algorithm (tri-diagonal solver) is used to solve the algebraic equations iteratively. The numerical advantage of the proposed scheme is tested for laminar fluid flows in a torus and in a square-driven cavity. The convergence rates are compared with the traditional schemes for the square-driven cavity problem. Good behaviour of the proposed scheme is ascertained.     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—Part II. Numerical algorithm and example applications

In the present study, an algorithm is presented for the dual-porosity model formulated in Part I of this series. The resultant flow equation with the dual-porosity formulation is of an integro-(partial) differential equation involving differential terms for the Darcy flow in large fractures and integrals in time for diffusion within matrix blocks. The algorithm developed here to solve this equation involves a step-by-step finite difference procedure combined with a quadrature scheme. The quadrature scheme, used for the integral terms, is based on the trapezoidal method which is of second-order precision. This order of accuracy is consistent with the first- and second-order finite difference approximations used here to solve the differential terms in the derived flow equation. In an approach consistent with many petroleum reservoir and groundwater numerical flow models, the example formulation presented uses a first-order implicit algorithm. A two-dimensional example is also demonstrated, with the proposed model and numerical scheme being directly incorporated into the commercial gas reservoir simulator SIMED II that is based on a fully implicit finite difference approach. The solution procedure is applied to several problems to demonstrate its performance. Results from the derived dual-porosity formulation are also compared to the classic Warren–Root model. Whilst some of this work confirmed previous findings regarding Warren–Root inaccuracies at early times, it was also found that inaccuracy can re-enter the Warren–Root results whenever there are changes in boundary conditions leading to transient variation within the domain.     

A fully conservative high‐order upwind multi‐moment method using moments in both upwind and downwind cells

We propose a fully conservative high‐order upwind multi‐moment method for the conservation equation. The proposed method is based on a third‐order polynomial interpolation function and semi‐Lagrangian formulation and is a variant of the constrained interpolation profile conservative semi‐Lagrangian scheme with third‐order polynomial function method. The third‐order interpolation function is constructed based on three constraints in the upwind cell (two boundary values and a cell average) and a constraint in the downwind cell (a cell center value). The proposed method shows fourth‐order accuracy in a benchmark problem (sine wave propagation). We also propose a less oscillatory formulation of the proposed method. The less oscillatory formulation can minimize numerical oscillations. These methods were validated through scalar transport problems, and compressible flow problems (shock tube and 2D explosion problems). Copyright © 2016 John Wiley & Sons, Ltd.     

A mixed modal-differential quadrature method for free and forced vibration of beams in contact with fluid

A simple and accurate mixed modal-differential quadrature formulation is proposed to study the dynamic behavior of beams in contact with fluid. Both free and forced vibration problems are considered. The proposed mixed methodology uses the modal technique for the structural domain while it applies the differential quadrature method (DQM) to the fluid domain. Thus, the governing partial differential equations of the beam and fluid are reduced to a set of ordinary differential equations in time. In the case of forced vibration, the Newmark time integration scheme is employed to solve the resulting system of ordinary differential equations. The proposed formulation, in general, combines the simplicity of the modal method and high accuracy and efficiency of the DQM. Its application is shown by solving some beam-fluid interaction problems. Comparisons with analytical solutions show that the present method is very accurate and reliable. To demonstrate its efficiency, the test problems are also solved using the finite element method (FEM). It is found that the proposed method can produce better accuracy than the FEM using less computational time. The technique presented in this investigation is general and can be used to solve various fluid-structure interaction problems.     

Homotopy analysis and differential quadrature solution of the problem of free-convective magnetohydrodynamic flow over a stretching sheet with the Hall effect and mass transfer taken into account

This paper presents an analytical solution of the problem of free-convective magnetohydrodynamic flow over a stretched sheet with the Hall effect and mass transfer taken into account. A similarity transform reduces the Navier-Stokes, energy, Ohm law, and mass-transfer equations to a system of nonlinear ordinary differential equations. The governing equations are solved analytically using an analytical method for solving nonlinear problems, namely, the homotopy analysis method. The results are compared with the results of a promising numerical method of differential quadrature developed by the authors. It is shown that there is very good agreement between analytical results and those obtained by the differential quadrature method. The differential quadrature method was validated, and the effects of non-dimensional parameters on the velocity, temperature and concentration profiles were studied.     

拟谱-微分求积混合方法求解一类双曲电报方程

拟谱方法和微分求积法是两类重要的无网格法，二者都已在科学和工程计算中获得了广泛应用。采用拉格朗日插值多项式作为二者的试函数，且采用同一种网格点分布，指出了在空间域上，微分求积法是拟谱方法的一种特殊形式。在此基础上，结合二者各自的特点，提出了拟谱-微分求积混合方法用于求解一类双曲电报方程。理论分析和数值测试表明，新方法在空间域上具有谱精度收敛性，在时间域上是A-稳定的，比较适合于求解多维电报方程。     

配置点谱方法-人工压缩法(SCM-ACM)求解同心圆筒内流体流动

发展了配置点谱方法SCM(Spectral collocation method)和人工压缩法ACM(Artificial compressibility method)相结合的SCM-ACM数值方法,计算了柱坐标系下稳态不可压缩流动N-S方程组。选取典型的同心圆筒间旋转流动Taylor-Couette流作为测试对象,首先,采用人工压缩法获得人工压缩格式的非稳态可压缩流动控制方程;再将控制方程中的空间偏微分项用配置点谱方法进行离散,得到矩阵形式的代数方程;编写了SCM-ACM求解不可压缩流动问题的程序;最后,通过与公开发表的Taylor-Couette流的计算结果对比,验证了求解程序的有效性。结果证明,本文发展的SCM-ACM数值方法能够用于求解圆筒内不可压缩流体流动问题,该方法既保留了谱方法指数收敛的特性,也具有ACM形式简单和易于实施的特点。本文发展的SCM-ACM数值方法为求解柱坐标下不可压缩流体流动问题提供了一种新的选择。