Infinite Element Approximation to Axial Symmetric Stokes Flow

We considered in [1] the finite element approximation to axial symmetric Stokes flow in a bounded domain. The problem for the flow passing an obstacle in an unbounded domain is also frequently encountered. In this paper, we are going to give approximate solutions for this problem by an approach stated in [2]. An iterative method is used to calculate the combined stiffness matrix.     

Sakiadis流动的一致有效级数解

通过引入一个变换式,克服了Sakiadis流动中半无限大流动区域以及无穷远处渐近边界条件所带来的数学处理上的困难.基于泛函分析中的不动点理论,采用不动点方法求解了变换后的非线性微分方程,获得了Sakiadis流动的近似解析解.该近似解析解用级数的形式来表达并在整个半无限大流动区域内一致有效.     

GL method for solving equations in math-physics and engineering

In this paper, we propose a GL method for solving the ordinary and the partial differential equation in mathematical physics and chemics and engineering. These equations govern the acustic, heat, electromagnetic, elastic, plastic, flow, and quantum etc. macro and micro wave field in time domain and frequency domain. The space domain of the differential equation is infinite domain which includes a finite inhomogeneous domain. The inhomogeneous domain is divided into finite sub domains. We present the solution of the differential equation as an explicit recursive sum of the integrals in the inhomogeneous sub domains. Actualy, we propose an explicit representation of the inhomogeneous parameter nonlinear inversion. The analytical solution of the equation in the infinite homogeneous domain is called as an initial global field. The global field is updated by local scattering field successively subdomaln by subdomain. Once all subdomains are scattered and the updating process is finished in all the sub domains, the solution of the equation is obtained. We call our method as Global and Local field method, in short , GL method. It is different from FEM method, the GL method directly assemble inverse matrix and gets solution. There is no big matrix equation needs to solve in the GL method. There is no needed artificial boundary and no absorption boundary condition for infinite domain in the GL method. We proved several theorems on relationships between the field solution and Green＇s function that is the theoretical base of our GL method. The numerical discretization of the GL method is presented. We proved that the numerical solution of the GL method convergence to the exact solution when the size of the sub domain is going to zero. The error estimation of the GL method for solving wave equation is presented. The simulations show that the GL method is accurate, fast, and stable for solving elliptic, parabolic, and hyperbolic equations. The GL method has advantages and wide applications in the 3D electromagnetic （EM）     

A variable boundary method for modelling two dimensional free surface flows with moving boundaries

A coordinate transformation method is proposed for modelling unsteady, depth-averaged shallow water equations for a open channel flow with moving lateral boundaries. The transformation technique which maps the changing domain onto a fixed domain and solves the governing equations in the mapped domain, facilitates the numerical treatment of an irregular boundary configuration. The transformed equations are solved on a staggered grid with a conditionally stable, explicit finite difference scheme. Several numerical experiments are carried out corresponding to different situations, viz., flow with constant discharge, flow with constant discharge and a closed boundary at the downstream, flow in a converging channel with constant discharge and finally flow with varying discharge. The experiments are used to verify the model ability to predict free surface elevation, circulatory pattern and displacement of the boundaries. The simulated results such as displaced area, depth, displacement–time and flow-field are used to evaluate the effects of excess discharge at the upstream on the movement of lateral boundaries.     

A Robin-Robin Domain Decomposition Method for a Stokes-Darcy Structure Interaction with a Locally Modified Mesh

A new numerical method based on locally modified Cartesian meshes is proposed for solving a coupled system of a fluid flow and a porous media flow. The fluid flow is modeled by the Stokes equations while the porous media flow is modeled by Darcy's law. The method is based on a Robin-Robin domain decomposition method with a Cartesian mesh with local modifications near the interface. Some computational examples are presented and discussed.     

Topology optimization methods for guided flow – comparison of optimality criteria vs. phase-field method

Optimization of guided flow problems is an important task for industrial applications especially those with high Reynolds numbers. There exist several optimization methods to increase the energy efficiency of these problems. Different optimization methods are shown bei Klimetzek [1], Hinterberger [2] and Pingen [3]. In recent years the phase-field method has been shown to be an applicable method for different kinds of topology optimization [4, 5]. We present results of topology optimization methods with optimality criterion and by using a phase-field model in the area of guided fluid flow problems. The two methods aim on the same main target reducing the pressure drop between the inlet and outlet of the flow domain. The first method is based on local optimality criterion, preventing the backflow in the flow domain [1, 6, 7]. The second method is based on a phase field model, which describes a minimization problem and uses a specially constructed driving force to minimize the total energy of the system [4, 5]. We investigate the capabilities and limits of both methods and present examples of different resulting geometries. The initial configurations are prepared in a way that the same optimization problem is solved with both methods. We discuss these results regarding the shape of the improved flow geometry. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rational Chebyshev collocation method for the similarity solution of two dimensional stagnation point flow

In this study, we propose an efficient and accurate numerical technique that is called the rational Chebyshev collocation (RCC) method to solve the two dimensional flow of a viscous fluid in the vicinity of a stagnation point named Hiemenz flow. The Navier-Stokes equations governing the flow, are reduced to a third-order ordinary differential equation of a boundary value problem with a semi-infinite domain by using similarity transformation. The rational Chebyshev method reduces this nonlinear ordinary differential equation to a system of algebraic equations. This technique is a powerful type of the collocation methods for solving the boundary value problems over a semi-infinite interval without truncating it to a finite domain. We also present the comparison of this work with others and show that the present method is more accurate and efficient.     

The influence of flow acceleration on tidal stream turbine wake dynamics: A numerical study using a coupled BEM–CFD model

Studies of tidal stream turbine performance and of wake development are often conducted in tow-tanks or in regulated flumes with uniform flows across the turbine. Whilst such studies can be very useful, it is questionable as to what extent the results would differ if the flows were more complex in nature, for instance if the flows were unsteady or non-uniform or even both. This study aims to explore whether the results would be affected once we move away from the uniform flow scenario. A numerical modelling study is presented in which tidal stream turbine performance and wake development in non-uniform flow conditions are assessed. The model implements the Blade Element Momentum method for characterising turbine rotor source terms which are used within a computational fluid dynamics model for predicting the interaction between the turbines and the surrounding flow. The model is applied to a rectangular domain and a range of slopes are implemented for the water surface to instigate an increase in flow velocity along the domain. Within an accelerated flow domain wake recovery occurred more rapidly although rotor performance was not affected.     

Modeling flow in porous media with fractures; Discrete fracture models with matrix-fracture exchange

This article is concerned with a numerical model for flow in a porous medium containing fractures. The fractures are modeled as (d − 1)-dimensional surfaces inside the d-dimensional matrix domain, and a mixed finite element method containing both d and (d − 1) dimensional elements is used. The method allows for fluid exchange between the fractures and the matrix. The method is defined for single-phase Darcy flow throughout the domain and for Forchheimer flow in the fractures. We also consider the case of two-phase flow in a domain in which the fractures and the matrix are of different rock type.     

二维无界自由衰减流的数值研究

无界区域上的流体运动是流体力学中的热点和难点问题.采用传统的扩大计算区域算法和新发展的基于无界区域的Hermite基函数算法对二维无界区域的自由衰减流动进行研究. 结果发现,对于只存在相同符号涡的初始流场而言,两种方法都可以得出正确的结果;而对于正负涡都存在的初始流场,传统方法即便利用非常大的计算区域也无法进行正确的长时间模拟,但是新方法却能高效求解.对算例的Hermite算法数值模拟验证了理论解Oseen涡的存在.     

Numerical Simulation of Red Blood Cell Suspensions Behind a Moving Interface in a Capillary

Computational modeling and simulation are presented on the motion of red blood cells behind a moving interface in a capillary. The methodology is based on an immersed boundary method and the skeleton structure of the red blood cell (RBC) membrane is modeled as a spring network. As by the nature of the problem, the computational domain is moving with either a designated RBC or an interface in an infinitely long two-dimensional channel with an undisturbed flow field in front of the computational domain. The tanking-treading and the inclination angle of a cell in a simple shear flow are briefly discussed for the validation purpose. We then present and discuss the results of the motion of red blood cells behind a moving interface in a capillary, which show that the RBCs with higher velocity than the interface speed form a concentrated slug behind the moving interface.     

Impedance eduction from energy dissipation measurements with flow

An optimization method for the calculation of the wall impedance of a perforate from measurements under grazing flow conditions is presented. It is in particular characterized by the simultaneous eduction of the parameter set defining the impedance as function of frequency, which allows interpolation and extrapolation. The Extended Helmholtz Resonator is applied as impedance model and the simultaneous eduction is performed using a time–domain CAA method. Furthermore, the method uses processed energy data from experiments to define the objective function. The results give an insight in the effect of bias flow on the wall impedance under grazing flow conditions. It proves to be an important acoustic design parameter. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence rate analysis of an asynchronous space decomposition method for convex minimization

We analyze the convergence rate of an asynchronous space decomposition method for constrained convex minimization in a reflexive Banach space. This method includes as special cases parallel domain decomposition methods and multigrid methods for solving elliptic partial differential equations. In particular, the method generalizes the additive Schwarz domain decomposition methods to allow for asynchronous updates. It also generalizes the BPX multigrid method to allow for use as solvers instead of as preconditioners, possibly with asynchronous updates, and is applicable to nonlinear problems. Applications to an overlapping domain decomposition for obstacle problems are also studied. The method of this work is also closely related to relaxation methods for nonlinear network flow. Accordingly, we specialize our convergence rate results to the above methods. The asynchronous method is implementable in a multiprocessor system, allowing for communication and computation delays among the processors.

     

Analysis of a FEM/BEM coupling method for transonic flow computations

A sensitive issue in numerical calculations for exterior flow problems, e.g.around airfoils, is the treatment of the far field boundary conditions on a computational domain which is bounded. In this paper we investigate this problem for two-dimensional transonic potential flows with subsonic far field flow around airfoil profiles. We take the artificial far field boundary in the subsonic flow region. In the far field we approximate the subsonic potential flow by the Prandtl-Glauert linearization. The latter leads via the Green representation theorem to a boundary integral equation on the far field boundary. This defines a nonlocal boundary condition for the interior ring domain. Our approach leads naturally to a coupled finite element/boundary element method for numerical calculations. It is compared with local boundary conditions. The error analysis for the method is given and we prove convergence provided the solution to the analytic transonic flow problem around the profile exists.

     

Finite element error estimates for 3D exterior incompressible flow with nonzero velocity at infinity

We consider a stationary incompressible Navier–Stokes flow in a 3D exterior domain, with nonzero velocity at infinity. In order to approximate this flow, we use the stabilized P1–P1 finite element method proposed by Rebollo (Numer Math 79:283–319, 1998). Following an approach by Guirguis and Gunzburger (Model Math Anal Numer 21:445–464, 1987), we apply this method to the Navier–Stokes system with Oseen term in a truncated exterior domain, under a pointwise boundary condition on the artificial boundary. This leads to a discrete problem whose solution approximates the exterior flow, as is shown by error estimates.     

Three dimensional fluid–structure interaction under pulsatile flow by using distributed Lagrange multiplier method

In this article we discuss an application of a Lagrange multiplier based fictitious domain method for the simulation of moving leaflets in an unsteady flow generated by pressure gradients during the systolic phase of the cardiac cycle in three dimensional geometry. The mathematical model includes the Navier–Stokes equations coupled with the Euler–Newton equations describing fluid–structure interaction for the generalized Neumann boundary conditions on upstream and downstream boundaries. The solution method includes the finite element method combined with an operator-splitting scheme, where the fictitious formulation allows the flow calculations to be in a fixed rectangular parallelepiped to predict the dynamic structural positions and flow field during the valve opening and closing phases.     

An advanced study on the solution of nanofluid flow problems via Adomian’s method

Nanofluid flow is one of the most important areas of research at the present time due to its wide and significant applications in industry and several scientific fields. The boundary layer flow of nanofluids is usually described by a system of nonlinear differential equations with boundary conditions at infinity. These boundary conditions at infinity cause difficulties for any of the series method, such as Adomian’s method, the variational iteration method and others.The objective of the present work is to introduce a reliable method to overcome such difficulties that arise due to an infinite domain. The proposed scheme, that we will introduce, is based on Adomian’s decomposition method, where we will solve a system of nonlinear differential equations describing the boundary layer flow of a nanofluid past a stretching sheet.     

A local parallel superconvergence method for the incompressible flow by coarsening projection

In this article, we investigate a local parallel superconvergence method by coarsening projection for the incompressible Stokes flow. The method is a combination of the local superconvergence technique and the given framework of local parallel method. For the smooth subdomains, the local superconvergence method is applied in a higher order finite dimensional space corresponding to an appropriate coarse mesh on interior domain. Moreover, a useful and flexible local parallel method is designed to obtain the local parallel superconvergence results of presented method, which offset theoretical limitation of the model without the smoothness of the exact solution and a priori regularity of the underlying problem over the whole domain. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1209–1223, 2015     

General method for producing a boundary-fitted orthogonal curvilinear mesh

A general method is described for computing an orthogonal mesh fitted to a two-dimensional physical domain with arbitrary closed boundary. The method allows optimum control of mesh spacing through the introduction of arbitrary (with weak constraints) ‘packing’ functions into the elliptic governing equations. Two particular aspects are addressed: first, the condition on a scaling factor which normalizes the mesh aspect ratio; second, the condition for avoiding run-out of the mesh beyond the boundaries of the physical domain.Conversion of the equations to finite difference form and appropriate iterative techniques are discussed. Finally applications of the method in the context of flow across a bundle of rods are presented.     

Longitudinal flow through an array of rectangular cylinders

The flow parallel to an array of rectangular cylinders is solvedby an efficient domain decomposition and matching method. Thepermeability constant is determined for a variety of geometries,including strips of zero thickness. Formulae for the drag ofa single rectangular cylinder and a single strip are derived.