A staggered grid‐based explicit jump immersed interface method for two‐dimensional Stokes flows

In this paper the explicit jump immersed interface method (EJIIM) is applied to stationary Stokes flows. The boundary value problem in a general, non‐grid aligned domain is reduced by the EJIIM to a sequence of problems in a rectangular domain, where staggered grid‐based finite differences for velocity and pressure variables are used. Each of these subproblems is solved by the fast Stokes solver, consisting of the pressure equation (known also as conjugate gradient Uzawa) method and a fast Fourier transform‐based Poisson solver. This results in an effective algorithm with second‐order convergence for the velocity and first order for the pressure. In contrast to the earlier versions of the EJIIM, the Dirichlét boundary value problem is solved very efficiently also in the case when the computational domain is not simply connected. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

An analytical method for a mixed boundary value problem of circular plates under arbitrary lateral loads

In this paper by using the concept of mixed boundary functions,an analytical method is proposed for a mixed boundary value problem of circular plates.The trial functions are constructed by using the series of particular solutions of the biharmonic equations in the polar coordinate system.Three examples are presented to show the stability and high convergence rate of the method.     

Three-Dimensional Analytical Solution for an Axisymmetric Biharmonic Problem

In this paper, we propose a method for the solution of the axisymmetric boundary value problem for a finite elastic cylinder with assigned stress and/or displacements acting on the ends and side. The technique utilizes the Love representation, which allows for reduction of the solution of the elastic problem to the search for a biharmonic function on a cylindrical domain. In the solution method suggested here, we write the Love function with a Bessel expansion and analyze in detail the conditions under which it is possible to differentiate the expansion term by term. We show that this is possible only for a restricted class of elastic solutions. In the general case, we introduce two new auxiliary functions of the z-coordinate. In this way, we obtain the general form of the axisymmetric biharmonic function, which is discussed in relation to certain specific boundary conditions applied on the side and ends of the cylinder. We obtain an exact explicit solution of practical interest for a cylinder with free ends and assigned displacements applied to the side.     

On Positivity for the Biharmonic Operator under Steklov Boundary Conditions

The positivity-preserving property for the inverse of the biharmonic operator under Steklov boundary conditions is studied. It is shown that this property is quite sensitive to the parameter involved in the boundary condition. Moreover, positivity of the Steklov boundary value problem is linked with positivity under boundary conditions of Navier and Dirichlet type.     

An exact Riemann solver for detonation products

Numerical methods based upon the Riemann Problem are considered for solving the general initial-value problem for the Euler equations applied to real gases. Most of such methods use an approximate solution of the Riemann problem when real gases are involved. These approximate Riemann solvers do not yield always a good resolution of the flow field, especially for contact surfaces and expansion waves. Moreover, approximate Riemann solvers cannot produce exact solutions for the boundary points. In order to overcome these shortcomings, an exact solution of the Riemann problem is developed, valid for real gases. The method is applied to detonation products obeying a 5th order virial equation of state, in the shock-tube test case. Comparisons between our solver, as implemented in Random Choice Method, and finite difference methods, which do not employ a Riemann solver, are given.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.     

A moving grid method applied to one-dimensional non-stationary flame propagation

The paper presents applications of a moving grid method to the combined problem of ignition and premixed flame propagation in a closed vessel. This method belongs to the general class of adaptive grid techniques for the numerical integration of evolutionary partial differential equations and is based on the method of lines with variable node position. In the present case the motion of the grid and the solution of the partial differential equations are strongly coupled by an implicit formulation. The problem is reduced to an initial value problem for a stiff differential-algebraic system. The continuously moving grid is determined by the equidistribution of a positive function which depends on the solution of the partial differential equations. A differential-algebraic system solver is used for the time integration of the initial value problem. The numerical results of the test problems demonstrate the computational efficiency and the capability of the method to resolve the main features of the solution.     

ODE conversion techniques and their applications in computational mechanics

In this paper, a number of ordinary differential equation (ODE) conversion techniques for transformation of nonstandard ODE boundary value problems into standard forms are summarised, together with their applications to a variety of boundary value problems in computational solid mechanics, such as eigenvalue problem, geometrical and material nonlinear problem, elastic contact problem and optimal design problems through some simple and representative examples. The advantage of such approach is that various ODE boundary value problems in computational mechanics can be solved effectively in a unified manner by invoking a standard ODE solver. The project is supported by National Natural Science Foundation of China.     

A necessary and sufficient set of boundary integral equations with indirect unknowns for plate bending problem

A boundary integral representation of plane biharmonic function is established rigorously by the method of unanalytical continuation in the present paper. In this representation there are two boundary functions and four constants which bear a one to one correspondence to biharmonic functions. Therefore the set of boundary integral equations with indirect unknowns based on this representation is equivalent to the original differential equation formulation.     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.     

Positivity and Almost Positivity of Biharmonic Green’s Functions under Dirichlet Boundary Conditions

In general, for higher order elliptic equations and boundary value problems like the biharmonic equation and the linear clamped plate boundary value problem, neither a maximum principle nor a comparison principle or—equivalently—a positivity preserving property is available. The problem is rather involved since the clamped boundary conditions prevent the boundary value problem from being reasonably written as a system of second order boundary value problems. It is shown that, on the other hand, for bounded smooth domains W ì \mathbbRⁿ{\Omega \subset\mathbb{R}^n} , the negative part of the corresponding Green’s function is “small” when compared with its singular positive part, provided n\geqq 3{n\geqq 3} . Moreover, the biharmonic Green’s function in balls B ì \mathbbRⁿ{B\subset\mathbb{R}^n} under Dirichlet (that is, clamped) boundary conditions is known explicitly and is positive. It has been known for some time that positivity is preserved under small regular perturbations of the domain, if n = 2. In the present paper, such a stability result is proved for n\geqq 3{n\geqq 3} .     

Solution of the elastic boundary value problems for a layer with tunnel stress raisers

A novel procedure for solving three-dimensional problems for elastic layer weakened by through-thickness tunnel cracks has been developed and is presented in this paper. This procedure reduces the given boundary value problem to an infinite system of one-dimensional singular integral equations and is based on a system of homogeneous solutions for a layer. Integral representations of single- and double-layer potentials are used for metaharmonic and harmonic functions entering in the singular integral equations. These representations provide a continuous extendibility of the stress vector while allowing a jump in the displacement vector in the transition through the cut.Expanding the potential and biharmonic solutions in the Fourier series over the thickness coordinate yields the integral representations of the displacement vector and stress tensor. The problem of reducing a denumerable set of the integral equations of the given boundary value problem to one-to-one correspondence with the set of unknown densities appearing in the Fourier’s coefficient representations has been settled efficiently. Numerical investigations show a rapid convergence of the proposed reduction procedure as applied to the solution of the infinite system of one-dimensional integral equations. Numerical examples illustrate the proposed method and demonstrate its advantages.     

Approximate method for solving a two-dimensional problem of elasticity theory

A numerical-analytical method based on approximation by harmonic or biharmonic functions is proposed for solving a mixed two-dimensional problem of elasticity theory. This method allows one to decrease the geometric dimensionality of the boundary-value problem by reducing it to minimization of the boundary residual. The resultant approximate analytical solution satisfies all equations of elasticity theory. Kazan’ State Technical University, Kazan’ 420111. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 4, pp. 179–185, July–August, 1999.     

The MFS for the Cauchy problem in two-dimensional steady-state linear thermoelasticity

We study the reconstruction of the missing thermal and mechanical fields on an inaccessible part of the boundary for two-dimensional linear isotropic thermoelastic materials from over-prescribed noisy (Cauchy) data on the remaining accessible boundary. This problem is solved with the method of fundamental solutions (MFS) together with the method of particular solutions (MPS) via the MFS-based particular solution for two-dimensional problems in uncoupled thermoelasticity developed in Marin and Karageorghis, 2012a, Marin and Karageorghis, 2013. The stabilisation/regularization of this inverse problem is achieved by using the Tikhonov regularization method (Tikhonov and Arsenin, 1986), whilst the optimal value of the regularization parameter is selected by employing Hansen’s L-curve method (Hansen, 1998).     

Representation of the solution of refined bending theory for isotropic plates

A solution of the bending problem for isotropic plates in a refined statement based on the system of six-order differential equations is proposed. A procedure for determining the general solutions of the corresponding biharmonic and metaharmonic equations is suggested. A method for satisfying the boundary conditions is given. The results of numerical studies of the stress state of an infinite plate with an elliptic cavity are given.     

On asymmetric bending of functionally graded solid circular plates

Based on the three-dimensional elasticity equations, this paper studies the elastic bending response of a transversely isotropic functionally graded solid circular plate subject to transverse biharmonic forces applied on its top surface. The material properties can continuously and arbitrarily vary along the thickness direction. By virtue of the generalized England’s method, the problem can be solved by determining the expressions of four analytic functions. Expanding the transverse load in Fourier series along the circumferential direction eases the theoretical construction of the four analytic functions for a series of important biharmonic loads. Certain boundary conditions are then used to determine the unknown constants that are involved in the four constructed analytic functions. Numerical examples are presented to validate the proposed method. Then, we scrutinize the asymmetric bending behavior of a transversely isotropic functionally graded solid circular plate with different geometric and material parameters.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

A method for solving the factorized vorticity-stream function equations by finite elements

A new finite element method for solving the time-dependent incompressible Navier-Stokes equations with general boundary conditions is presented. The two second-order partial differential equations for the vorticity and the stream function are factorized, apart from the non-linear advection term, by eliminating the coupling due to the double specification on the stream function at (a part of) the boundary. This is achieved by reducing the no-slip boundary conditions to projection integral conditions for the vorticity field and by evaluating the relevant quantities involved according to an extension of the method of Glowinski and Pironneau for the biharmonic problem. Time integration schemes and iterative algorithms are introduced which require the solution only of banded linear systems of symmetric type. The proposed finite element formulation is compared with its finite difference equivalent by means of a few numerical examples. The results obtained using 4-noded bilinear elements provide an illustration of the superiority of the finite element based spatial discretization.     

Waves generated by an inclined-plate wave generator

This paper describes the characteristics of small-amplitude waves generated by a sinusoidally oscillating, inclined paddle-type wavemaker operating in a constant-depth channel. Two-dimensional, linearized potential flow is assumed. A semi-analytical method, the boundary collocation method, is used to establish the relationship between wave amplitude and paddle stroke. The numerical results are compared with the numerical results of the boundary integral equation method. It is found that the boundary collocation method is simpler and more flexible to implement and faster to compute. In addition, the numerical results are in reasonably good agreement with the laboratory experimental data. For the vertical wavemaker, which is a special case of the inclined wavemaker, an analytical series solution can be found. By using the boundary collocation method and the boundary integral equation method to solve the vertical wavemaker problem and comparing the results with the analytical series solution, it is found that the boundary collocation method yields a solution which is much more accurate than that from the boundary integral equation method. Finally, the relationships between wave amplitude and paddle stroke are established for different inclinations of the paddle-type wavemaker, based on the boundary collocation method.