An effective direct solution method for certain boundary element equations in 3D

The boundary element method for the Dirichlet problem in a three-dimensional rotational domain leads to a system of linear equations with a full dense matrix having a special block structure. A direct solution method for such systems is presented, which requires O(N^3/2 ln N) arithmetical operations only, using a Fast Fourier Transformation (FFT), where N denotes the number of unknowns on the boundary surface.     

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

A numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1], [2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution.     

解聯立方程的差分方程法

<正> 一. 引 聯立一次方程的求解早就不是一個理論問題,而是一個改進計算技術的問題.問題在如何組織計算使計算機械化從而節省工作量. 給定充分多始值後,線性差分方程是很容易解的.在本文內,我們把特種的及一般的聯立一次方程組的解看作線性差分方程滿足某種邊值的解,從而推求出求聯立一次方程組的準確解的一種簡單的機械的列表計算方法。     

Fredholm boundary value problems for perturbed systems of dynamic equations on time scales

This paper offers conditions ensuring the existence of solutions of linear boundary value problems for systems of dynamic equations on time scales. Utilizing a method of Moore–Penrose pseudo‐inverse matrices leads to an analytical form of a criterion for the existence of solutions in a relevant space and, moreover, to the construction of a family of linearly independent solutions of such problems in a general case with the number of boundary conditions (defined by a linear vector functional) not coinciding with the number of unknowns of a system of dynamic equations. As an example of an application of the presented results, the problem of bifurcation of solutions of boundary value problems for systems of dynamic equations on time scales with a small parameter is considered.     

Localized method of fundamental solutions for large-scale modeling of two-dimensional elasticity problems

The traditional method of fundamental solutions (MFS) based on the “global” boundary discretization leads to dense and non-symmetric coefficient matrices that, although smaller in sizes, require huge computational cost to compute the system of equations using direct solvers. In this study, a localized version of the MFS (LMFS) is proposed for the large-scale modeling of two-dimensional (2D) elasticity problems. In the LMFS, the whole analyzed domain can be divided into small subdomains with a simple geometry. To each of the subdomain, the traditional MFS formulation is applied and the unknown coefficients on the local geometric boundary can be calculated by the moving least square method. The new method yields a sparse and banded matrix system which makes the method very attractive for large-scale simulations. Numerical examples with up to 200,000 unknowns are solved successfully using the developed LMFS code.     

Pre-elimination approach to the modelling of segmented representations of cable shape with application to the positioning of a towed sensor

Discrete models for the determination of cable shapes and tensions require the solution of large systems of nonlinear algebraic equations. In these the unknowns are the orientations and coordinates of the segments as well as their tensions. Smaller systems of algebraic equations are obtained by the pre-elimination of the tensions and the coordinates for the lumped mass and rod models which leave the orientation variables as the only unknowns. Furthermore, this pre-elimination analysis yields conditions for the equivalence of the lumped mass and the rod models. This method is applied to the determination of the position of a towed underwater sensor used for seabed mapping or for oceanographic recordings.     

The difference method in a diffraction problem: The technique of interior boundary condition

The diffraction problem for a plane wave on a half-plane covered by thin layer with an interface is solved by the difference method. The system of difference equations is derived from the variational principle. A boundary solution at infinity must be imposed; this is a radiation condition, which is used in the form of the limit absorption principle. The arising infinite system of difference equations is reduced to a finite part of the boundary (the interface) by using the technique of so-called interior boundary conditions in the sense of Ryaben’kii. The real conditions are found by the Fourier method with respect to one spatial variable in the form of Fourier or Laurent series in the corresponding variable, which converge either inside, outside, or on the unit circle. Above the upper boundary of the layer, all unknowns are eliminated by using the so-called grid Green function, that is, the resolving function for the half-plane satisfying the radiation condition at infinity. For the unknowns on the upper boundary of the layer, an equation in terms of a function of a complex variable of Wiener-Hopf type is obtained, which is solved by factorization. Factorization is performed numerically: the logarithm of the function is expanded in a bi-infinite series, which is replaced by a discrete Fourier series. The closing system in a neighborhood of the interface has order proportional to the number of points on the interface. Solving this system yields all of the required characteristics of the solution.     

Compactness of the product topology and the solvability of infinite systems of infinite linear equations

If every finite subsystem of an infinite system of linear equations (say, over the field of real numbers) each with finitely many unknowns has a solution then the entire system has a solution. The situation is not so if the equations contain infinitely many unknowns. In this case, as shown below, the solvability of every finite subsystem implies the solva. bility of the entire system provided finite subsystems have solution with common upper and lower bounds and the coefficients of ever equation satisfy some boundedness or convergence conditions. The passage from the solvability of finite subsystem to the solvability of the entire system is achieved based on Tychnoff’s theorem stating that any product of compact topological spaces is compact in their product topology.     

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.     

A method for exact computation with rational numbers

An easily programmed method is presented for solving N linear equations in N unknowns exactly for the rational answers, given that all coefficients and constants appearing in the equations are rational numbers. The rational answers are deduced from floating point approximations to the answers obtained by any of the standard solution algorithms. Criteria are given for determining for a particular set of equations the floating point precision needed.     

Robust control problems for the multilayer quasi-geostrophic system

In this article, we study some robust control problems associated with the multilayer quasi-geostrophic equations of the ocean and related to data assimilation in oceanography. We consider higher norms (compared to [T. Tachim Medjo, Robust control problems associated with the multilayer quasi-geostrophic equations of the ocean, Appl. Math. Optim. 51(3) (2005) 333–360]) in the definition of the cost functionals. We prove the existence and uniqueness of solutions. The result relies on better a priori estimates on the solutions to the multilayer quasi-geostrophic system obtained using a new formulation that we introduce for the multilayer quasi-geostrophic equation of the ocean. The new formulation replaces the non-homogenous boundary conditions (and the non-local constraint) on the stream-function by a simple homogenous Dirichlet boundary condition.     

线载和弹性支承作用面内运动薄板磁固耦合双重共振

针对磁场环境中具有线载荷和弹性支承作用的面内运动薄板,给出了系统的势能、动能及电磁力表达式,应用Hamilton变分原理,推得面内运动条形板的磁固耦合非线性振动方程.考虑边界为夹支-铰支的约束条件,利用变量分离法和Galerkin积分法,得到了含简谐线载力和电磁阻尼力项的两自由度非线性振动微分方程组.应用多尺度法对主-...     

Multigrid solution of high order discretisation for three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind

In this paper, we propose two compact finite difference approximations for three-dimensional biharmonic equation with Dirichlet boundary conditions of second kind. In these methods there is no need to define special formulas near the boundaries and boundary conditions are incorporated with these techniques. The unknown solution and its second derivatives are carried as unknowns at grid points. We derive second-order and fourth-order approximations on a 27 point compact stencil. Classical iteration methods such as Gauss–Seidel and SOR for solving the linear system arising from the second-order and fourth-order discretisation suffer from slow convergence. In order to overcome this problem we use multigrid method which exhibit grid-independent convergence and solve the linear system of equations in small amount of computer time. The fourth-order finite difference approximations are used to solve several test problems and produce high accurate numerical solutions.     

Connection between the spectral problem for linear matrix pencils and some problems of algebra

We suggest a method by which the solution of systems of nonlinear algebraic equations in one and two variables can be reduced to the spectral problem for linear pencils of two matrices and for a system of two matrix pencils of two matrices, respectively. This method is substantially different from the traditional methods of solution; at the same time it is useful for the study of the solvability and the determination of the number of solutions of such systems. We propose a modification of the well-known elimination method for the solution of nonlinear algebraic systems of two equations in two unknowns which provides a new approach to studying and solving the problem.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 80, pp. 98–116, 1978.     

Multidomain direct method and its applicability to the local time-step concept

A new multidomain direct (noniterative) method for solving boundary-value problems is presented. Using this method, the solution is expressed by a linear combination of auxiliary functions and unknowns which pertain to the boundaries of a subdomain. This approach enables us to solve problems independently and exactly, without any iterations between subdomains. As a consequence, different types of equations and coordinate systems may be considered in different subdomains. Moreover, different boundary conditions and variable (in space) time steps may be imposed on the subdomains as well. Applications are given for initial boundary-value problems with known analytical solutions, including a highly nonlinear porosity equation. © 1994 John Wiley & Sons, Inc.     

中细柔性圆环壳整体弯曲的一般解及在波纹管计算中的应用(Ⅱ)——Ω型波纹管的计算

(Ⅱ)是(Ⅰ)的具体应用。计算了Ω型波纹管的角向刚度、横向刚度和应力分布,并将所得结果与有关的细环壳理论及实验进行了比较。结果表明,单独用(Ⅰ)的非齐次解能够计算Ω型波纹管的纯弯曲,而且比细环壳理论更接近实际;但在横向位移作用下,(Ⅰ)的非齐次解只能部分地满足边界条件,此时应同时考虑齐次解的作用,即完整的一般解(Ⅰ)才能满足所有的要求.     

A Boundary Integral Formulation for Poroelastic Materials

Wave propagation in porous media is an important topic, e.g. in geomechanics or the oil-industry. We formulate a linear system of coupled partial differential equations based on Biot's theory with the solid displacements and the pore pressure as the primary unknowns. To solve this system of coupled partial differential equations in a semi-infinite homogeneous domain the BEM (Boundary element method) is especially suitable. Starting from a representation formula a system of two boundary integral equations is derived. These boundary integral equations are used to solve related boundary value problems via a direct approach. Coercivity of the resulting bilinear form is shown, from which unique solvability of the variational formulation follows from injectivity. Using these results we derive the unique solvability of the related boundary integral equations. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical conformal mapping based on the generalised conjugation operator

An iterative procedure for numerical conformal mapping is presented which imposes no restriction on the boundary complexity. The formulation involves two analytically equivalent boundary integral equations established by applying the conjugation operator to the real and the imaginary parts of an analytical function. The conventional approach is to use only one and ignore the other equation. However, the discrete version of the operator using the boundary element method (BEM) leads to two non-equivalent sets of linear equations forming an over-determined system. The generalised conjugation operator is introduced so that both sets of equations can be utilised and their least-square solution determined without any additional computational cost, a strategy largely responsible for the stability and efficiency of the proposed method. Numerical tests on various samples including problems with cracked domains suggest global convergence, although this cannot be proved theoretically. The computational efficiency appears significantly higher than that reported earlier by other investigators.

     

Incremental Unknowns Method for Solving Three-Dimensional Convection-Diffusion Equations

We use the incremental unknowns method in conjunction with the iterative methods to approximate the solution of the nonsymmetric and positive-definite linear systems generated from a multilevel discretization of three-dimensional convection-diffusion equations. The condition numbers of incremental unknowns matrices associated with the convection-diffusion equations and the number of iterations needed to attain an acceptable accuracy are estimated. Numerical results are presented with two-level approximations, which demonstrate that the incremental unknowns method when combined with some iter- ative methods is very effcient.     

含参数线性方程组解的判别方法注记

当线性方程组中含有未知参数时,线性方程组解的情况往往需要进行讨论.本文给出了在非齐次线性方程组系数矩阵中含有未知参数且系数行列式等于零的情况下,判定对应参数值下方程组的解是无解还是有无穷多解的两个判定定理.和以前的方法比较,本文提出的讨论方法更直接.