THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the up-stream artificial boundary and two vertical sides are introduced as the downstream arti-ficial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.     

A mapping technique for numerical computations of bed evolutions

Free surface flow is one of the most difficult problems in engineering to be solved, since velocity and pressure fields depend on the free surface. On the other hand, the position of the free surface is unknown previously. Furthermore, the boundary condition on the free surface is expressed by a complicated equation. In an alluvial stream, where the boundaries of the domain are not fixed, addition of free surface at the bed will increase this difficulty. A domain mapping technique is developed in this paper to study the bed evolutions. The flow is considered 2D, choosing two coordinates in streamwise and upward directions. With a proper transformation, the hydrodynamics and sediment transport governing equations in irregular domain will be mapped into a simple rectangular one. The new domain can be discretize by finite elements. The transformed governing equations are solved to obtain desired variables in the mapped domain. With a proper transformation, there is no need of inverse mapping to obtain the free water surface profile and bedform evolution and migration in the actual domain. The model has been applied to streams with movable bed and the results show a good agreement with the experimental experiences.     

First-order accurate finite difference schemes for boundary vorticity approximations in curvilinear geometries

In this work we develop first-order accurate, forward finite difference schemes for the first derivative on both a uniform and a non-uniform grid. The schemes are applied to the calculation of vorticity on a solid wall of a curvilinear, two-dimensional channel. The von Mises coordinates are used to transform the governing equations, and map the irregular domain onto a rectangular computational domain. Vorticity on the solid boundary is expressed in terms of the first partial derivative of the square of the speed of the flow in the computational domain, and the derived finite difference schemes are used to calculate the vorticity at the computational boundary grid points using combinations of up to five computational domain grid points. This work extends previous work (Awartani et al., 2005) [3] in which higher-order schemes were obtained for the first derivative using up to four computational domain grid points. The aim here is to shed further light onto the use of first-order accurate non-uniform finite difference schemes that are essential when the von Mises transformation is used. Results show that the best schemes are those that use a natural sequence of non-uniform grid points. It is further shown that for non-uniform grid with clustering near the boundary, solution deteriorates with increasing number of grid points used. By contrast, when a uniform grid is used, solution improves with increasing number of grid points used.     

Enhanced electroosmotic flow in a nano-channel patterned with curved hydrophobic strips

We consider the electroosmotic flow (EOF) in a nano-channel in which the channel walls are modulated with a periodic array of curved hydrophobic patches. The objective is to achieve an enhanced flow compared to a slit nano-channel. The shape of the hydrophobic strips are considered to be of sinusoidal form, which resembles the situation in which the channel indentations are filled with immiscible nonconducting fluid over which the electrolyte is considered to be in metastable Cassie state. The homogeneous no-slip portions of the channel walls are considered to posses a constant surface-potential (zeta-potential) or constant surface charge density, while the hydrophobic regions are uncharged. A mathematical model based on the Nernst–Planck–Navier-Stokes equations are considered to analyze the present EOF. A coordinate transformation is adopted to map the irregular physical domain to a regular computational domain. A pressure-correction based control volume approach is adopted to solve the governing equations. We have studied the EOF by varying the amplitude of the hydrophobic region. Our results show that an enhancement in EOF compared to a slit-channel is possible when the Debye length is in the order of the channel height. The EOF in the patterned channel varies with the planform length of the hydrophobic region as well as the relative span of the slip and no-slip regions. A comparison with the pressure-driven flow is also presented to analyze the hindrance created by the electric body force of the unbalanced ions.     

Shape-explicit constants for some boundary integral operators

Among the well-known constants in the theory of boundary integral equations are the coercivity constants of the single-layer potential and the hypersingular boundary integral operator, and the contraction constant of the double-layer potential. Whereas there have been rigorous studies how these constants depend on the size and aspect ratio of the underlying domain, only little is known on their dependency on the shape of the boundary. In this article, we consider the homogeneous Laplace equation and derive explicit estimates for the above-mentioned constants in three dimensions. Using an alternative trace norm, we make the dependency explicit in two geometric parameters, the so-called Jones parameter and the constant in Poincaré's inequality. The latter one can be tracked back to the constant in an isoperimetric inequality. There are many domains with quite irregular boundaries, where these parameters stay bounded. Our results provide a new tool in the analysis of numerical methods for boundary integral equations and in particular for boundary element based domain decomposition methods.     

Nonreflecting boundary conditions and numerical simulation of external flows

The specification of conditions on artificial boundaries of the computational domain in the simulation of subsonic viscous gas flows is considered. The steps in the construction and implementation of nonreflecting boundary conditions on the path from one-dimensional linearized Euler equations to real-life problems are described. The technique is intended for flow simulation at low Mach numbers. Numerical results for the essentially subsonic flow over a flat plate are presented     

Clifford algebra valued boundary integral equations for three-dimensional elasticity

Applications of Clifford analysis to three-dimensional elasticity are addressed in the present paper. The governing equation for the displacement field is formulated in terms of the Dirac operator and Clifford algebra valued functions so that a general solution is obtained analytically in terms of one monogenic function and one multiple-component spatial harmonic function together with its derivative. In order to solve numerically the three-dimensional problems of elasticity for an arbitrary domain with complicated boundary conditions, Clifford algebra valued boundary integral equations (BIEs) for multiple-component spatial harmonic functions at an observation point, either inside the domain, on the boundary, or outside the domain, are constructed. Both smooth and non-smooth boundaries are considered in the construction. Moreover, the singularities of the integrals are evaluated exactly so that in the end singularity-free BIEs for the observation point on the boundary taking values on Clifford numbers can be obtained. A Clifford algebra valued boundary element method (BEM) based on the singularity-free BIEs is then developed for solving three-dimensional problems of elasticity. The accuracy of the Clifford algebra valued BEM is demonstrated numerically.     

Expansion of a wedge of non-ideal gas into vacuum

We study the problem of expansion of a wedge of non-ideal gas into vacuum in a two-dimensional bounded domain. The non-ideal gas is characterized by a van der Waals type equation of state. The problem is modeled by standard Euler equations of compressible flow, which are simplified by a transformation to similarity variables and then to hodograph transformation to arrive at a second order quasilinear partial differential equation in phase space; this, using Riemann variants, can be expressed as a non-homogeneous linearly degenerate system provided that the flow is supersonic. For the solution of the governing system, we study the interaction of two-dimensional planar rarefaction waves, which is a two-dimensional Riemann problem with piecewise constant data in the self-similar plane. The real gas effects, which significantly influence the flow regions and boundaries and which do not show-up in the ideal gas model, are elucidated; this aspect of the problem has not been considered until now.     

Long Eddies in Sheared Flows

The characteristic feature of the wide variety of hydraulic shear flows analyzed in this study is that they all contain a critical level where some of the fluid is turned relative to the ambient flow. One example is the flow produced in a thin layer of fluid, contained between lateral boundaries, during the passage of a long eddy. The boundaries of the layer may be rigid, or flexible, or free; the fluid may be either compressible or incompressible. A further example is the flow produced when a shear layer separates from a rigid boundary producing a region of recirculating flow. The equations used in this study are those governing inviscid hydraulic shear flows. They are similar in form to the classical boundary layer equations with the viscous term omitted. The main result of the study is to show that when the hydraulic flow is steady and contained between lateral boundaries, the variation of vorticity ω(ψ) cannot be prescribed at any streamline which crosses the critical level. This variation is, in fact, determined by (1) the vorticity distribution at all streamlines which do not cross the critical level, by (2) the auxiliary conditions which must be satisfied at the boundaries of the fluid layer, and by (3) the dimensions of the region containing the turned flow. If at some instant the vorticity distribution is specified arbitrarily at all streamlines, generally the subsequent flow will be unsteady. In order to emphasize this point, a class of exact solutions describing unsteady hydraulic flows are derived. These are used to describe the flow produced by the passage of a long eddy which distorts as it is convected with the ambient flow. They are also used to describe the unsteady flow that is produced when a shear layer separates from a boundary. Examples are given both of flows in which the shear layer reattaches after separation and of flows in which the shear layer does not reattach. When the shear layer vorticity distribution has the form ωαyⁿ, where y is a distance measure across the layer, the steady flows are of Falkner-Skan type inside, and adjacent to, the separation region. The unsteady flows described in this paper are natural generalizations of these Falkner-Skan flows. One important result of the analysis is to show that if the unsteady flow inside the separation region is strongly sheared, then the boundary of the separation region moves upstream towards the point of separation, forming large transverse currents. Generally, the assumption of hydraulic flow becomes invalid in a finite time. On the other hand, if the flow inside the separation region is weakly sheared, this region is swept downstream and the flow becomes self-similar.     

An inverse problem of identifying the coefficient of first-order in a degenerate parabolic equation

This work studies an inverse problem of determining the first-order coefficient of degenerate parabolic equations using the measurement data specified at a fixed internal point. Being different from other ordinary parameter identification problems in parabolic equations, in our mathematical model there exists degeneracy on the lateral boundaries of the domain, which may cause the corresponding boundary conditions to go missing. By the contraction mapping principle, the uniqueness of the solution for the inverse problem is proved. A numerical algorithm on the basis of the predictor-corrector method is designed to obtain the numerical solution and some typical numerical experiments are also performed in the paper. The numerical results show that the proposed method is stable and the unknown function is recovered very well. The results obtained in the paper are interesting and useful, and can be extended to other more general inverse coefficient problems of degenerate PDEs.     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.     

Incompressible flows in elastic domains: an immersed boundary method approach

This study illustrates how the immersed boundary method may be applied to perform the numerical simulation of incompressible flows in two-dimensional domains bounded by elastic boundaries. It presents the basic intermediate steps involved in the derivation of a solution methodology, from a scientific motivation to the numerical results, which can be applied for both steady and transient problems, even when the boundaries have an arbitrary shape. Its motivation, briefly presented, was borne in a bioengineering problem: the numerical simulation of the performance of ventricular assist devices. The mathematical model is composed by the Navier–Stokes equations, where the forcing term contains singular forces which arise from the elastic stresses acting on the boundaries. The incompressibility constraint is modified to introduce the inflow and outflow conditions into the problem through the use of sources and sinks. The methodology is applied to simulate two problems: the steady flow between two parallel plates, for which the exact solution is known and can be used to validate the approach, and the periodic flow in a winding channel, a transient problem in a non-trivial domain.     

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.     

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.     

Linear and nonlinear dynamics of a circular cylindrical shell connected to a rigid disk

The dynamics of a circular cylindrical shell carrying a rigid disk on the top and clamped at the base is investigated. The Sanders–Koiter theory is considered to develop a nonlinear analytical model for moderately large shell vibration. A reduced order dynamical system is obtained using Lagrange equations: radial and in-plane displacement fields are expanded by using trial functions that respect the geometric boundary conditions.The theoretical model is compared with experiments and with a finite element model developed with commercial software: comparisons are carried out on linear dynamics.The dynamic stability of the system is studied, when a periodic vertical motion of the base is imposed. Both a perturbation approach and a direct numerical technique are used. The perturbation method allows to obtain instability boundaries by means of elementary formulae; the numerical approach allows to perform a complete analysis of the linear and nonlinear response.     

构造容许位移的转换边界法

本文提出了在复杂边界条件下构造容许位移的转换边界法.所谓转换边界法,就是首先根据叠加原理将实际系统转换为基本系统及附加边界力和附加边界位移,然后应用变分原理于基本系统,最后应用级数转换的办法求得实际系统的容许位移.本文还给出了边界条件变化的混合能量原理.这边界条件变化的混合能量原理和边界条件变化的势能原理和余能原理都是转换边界法的主要理论基础.应用转换边界法我们构造了复杂边界条件平面应力问题和弯曲问题矩形板的容许位移.由于转换边界法构造容许位移是遵循着变分原理和确定的程序进行的,因而克服了Rayleigh-Ritz法猜测、拼凑容许位移的困难.     

非线性自然弯扭闭口薄壁复合梁的广义变分原理

对复合材料自然弯扭闭口薄壁细长梁在小应变、大位移和大转动的情况作了研究,建立了两端边界均为完全约束的该梁大变形弹性理论的非完全广义变分原理的泛函。由泛函驻值条件可以导出所给问题的平衡方程及全部边界条件。上述方法可方便地推广到其它各种不完全约束边界的情况。此外,利用上述结果还可以得到该梁在小位移理论中的基本方程和有关公式。     

Transonic nozzle flows and free boundary problems for the full Euler equations

We establish the existence and uniqueness of transonic flows with a transonic shock through a two-dimensional nozzle of slowly varying cross-sections. The transonic flow is governed by the steady, full Euler equations. Given an incoming smooth flow that is close to a constant supersonic state (i.e., smooth Cauchy data) at the entrance and the subsonic condition with nearly horizontal velocity at the exit of the nozzle, we prove that there exists a transonic flow whose downstream smooth subsonic region is separated by a smooth transonic shock from the upstream supersonic flow. This problem is approached by a one-phase free boundary problem in which the transonic shock is formulated as a free boundary. The full Euler equations are decomposed into an elliptic equation and a system of transport equations for the free boundary problem. An iteration scheme is developed and its fixed point is shown to exist, which is a solution of the free boundary problem, by combining some delicate estimates for the elliptic equation and the system of transport equations with the Schauder fixed point argument. The uniqueness of transonic nozzle flows is also established by employing the coordinate transformation of Euler-Lagrange type and detailed estimates of the solutions.     

Quadrature formulas for periodic functions and their application to the boundary element method

Two-dimensional and axisymmetric boundary value problems for the Laplace equation in a domain bounded by a closed smooth contour are considered. The problems are reduced to integral equations with a periodic singular kernel, where the period is equal to the length of the contour. Taking into account the periodicity property, high-order accurate quadrature formulas are applied to the integral operator. As a result, the integral equations are reduced to a system of linear algebraic equations. This substantially simplifies the numerical schemes for solving boundary value problems and considerably improves the accuracy of approximation of the integral operator. The boundaries are specified by analytic functions, and the remainder of the quadrature formulas decreases faster than any power of the integration step size. The examples include the two-dimensional potential inviscid circulation flow past a single blade or a grid of blades; the axisymmetric flow past a torus; and free-surface flow problems, such as wave breakdown, standing waves, and the development of Rayleigh-Taylor instability.     

直角平面内圆孔对稳态SH波的散射

利用复变函数方法和多极坐标移动技术，研究了直角平面内圆孔在直边分布有反平面稳态载荷时的sH波散射问题．首先构造出直角平面内不含有圆孔时满足边界应力条件的Green函数解；其次提出直角平面内存在圆孔时满足边界应力自由条件的散射波解，并利用叠加原理写出问题的位移总波场．借助于多极坐标移动技术和圆孔边界处应力自由条件，列出求解散射波解中未知系数的无穷代数方程组，在满足计算精度的前提下，通过有限项截断进行求解．作为算例，具体讨论了圆孔边界处的环向动应力随不同波数、圆孔位置及载荷分布位置和分布范围大小的变化情况，算例结果说明了算法的有效实用性．