Smooth boundary integration method in a linear two-dimensional problem of incompressible fluid and solid

This paper presents a new application of a theoretical and computational method of smooth boundary integration which belongs to the methods of boundary integral equations. Smooth integration is not a method of approximation. In its final analytical form, a smooth-kernel integral equation is computerized easily and accurately.

Smooth integration is associated with a “pressure-vorticity” formulation which covers linear problems in elasticity and fluid mechanics. The solution presented herein is essentially the same as that reported in an earlier paper for regular elasticity. The constraint of incompressibility does not cause difficulties in the pressure-vorticity formulation.

The linear fluid mechanics problem formulated and solved in this paper covers Stokes' problem of a slow viscous flow, and has a wider interpretation. The translational inertia forces are incorporated in the linear problem, as in Euler's dynamic theory of inviscid flow. The centrifugal inertia forces are left for the non-linear problem. The linear problem is perceived as a step in solution of the non-linear problems.     

Contact Interaction of Faces of the Plane Crack under Harmonic Loading

The present work is devoted to the solution of the three-dimensional fracture mechanics problem for a linear elastic, homogeneous and isotropic solid with a stationary plane crack under normal time-harmonic loading. The problem has been solved by the method of boundary integral equations with the allowance for the contact interaction of the opposite faces of the crack. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Boundary integral equation formulation for generalized thermoelastic diffusion – Analytical aspects

The present work deals with the formulation of the boundary integral equations for the solution of equations under linear theory of generalized thermoelastic diffusion in a three-dimensional Euclidean space. A mixed initial-boundary value problem is considered in the present context and the fundamental solutions of the corresponding coupled differential equations are obtained in the Laplace transform domain by employing the treatment of scalar and vector potential theory. A reciprocal relation of Betti type is established. Then we formulate the boundary integral equations for generalized thermoelastic diffusion on the basis of these fundamental solutions and the reciprocal relation.     

Variational differential quadrature: A technique to simplify numerical analysis of structures

A new solution method in the area of computational mechanics is developed in this article, which is called variational differential quadrature (VDQ). The main idea of this method is based on the accurate and direct discretization of the energy functional in the structural mechanics. In the VDQ method, through developing an efficient matrix formulation and using an accurate integral operator, the discretized governing equations are derived directly from the weak form of the equations with no need for the analytical derivation of the strong form. This technique provides an alternative way to discretize the energy functional, which avoids the local interpolation and the assembly process of the methods of this kind. We first implement the VDQ method for the nonlinear elasticity theory considering the Green-St. Venant strain tensor; then we simplify the formulation further for the first-order shear deformable beam and plate theories. The final formulation of these cases demonstrates the simplicity of the implementation for the VDQ method in the numerical analysis of the structures, which is a major goal for this article. Using these examples, one can easily learn and apply this technique to other structures. To assess the performance of the VDQ method, we compare it with the generalized differential quadrature (GDQ) method and finite element method (FEM) in the case of bending analysis of Mindlin plates. It is indicated that computational cost of VDQ is less than that of GDQ, and the convergence rate of VDQ is faster than that of FEM.     

Application of an anisotropic growth formulation to computational structural design

Common structural optimisation problems consist of problem-specific objective functions which have to be minimised mathematically with respect to design and state variables taking into account particular constraints. In contrast to this, we adopt a conceptually different approach for the design of a structure which is not based on a topology-optimisation technique. Instead, we apply a one-dimensional energy-driven constitutive evolution equation for the referential density–originally proposed for the simulation of remodelling effects in bones–and embed this into the micro-sphere-concept to end up with a three-dimensional anisotropic growth formulation. The objective of this contribution is to investigate this approach with emphasis on its application to structural design problems by means of two three-dimensional benchmark-type boundary value problems. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Compression of an orthotropic prism

The formulation and mathematically correct solution of a three-dimensional boundary value problem are given. A numerical calculation is presented for a prism of Canadian spruce and the displacement and stress fields are determined.Presented at a scientific seminar on polymer mechanics held in the Mechanics-Mathematics Department of Moscow State University, February 7, 1968.M. V. Lomonosov Moscow State University. Translated from Mekhanika Polimerov, Vol. 4, No. 5, pp. 810–815, September–October, 1968.     

On boundary-volume element discretization of inhomogenous elastodynamic problems

This paper presents an integral equation formulation and its discretization scheme for the elastodynamic problem in which the material properties are prescribed as arbitrary, continuous and differentiable functions of the spatial coordinates. The formulation is made by using the Green's function for the corresponding problem in homogenous elasticity. From a weighted residual statement of the problem, the governing differential equation is transformed into a set of the integral equations in the inner domain as well as on the boundary. These integral equations are discretized by introducing a finite number of the boundary-volume-time elements, and the solution for the system of linear equations thus obtained is discussed.     

A finite element method for the double-layer potential solutions of the Neumann exterior problem

A variational formulation for the integral equation used for the double layer potential solution of the Neumann exterior problem in the Laplace equation was proposed in [4]. This formulation allows the use of a finite element method which we describe and experiment here.     

The interaction of punches on the face of an elastic wedge

The interaction of two punches, which are elliptic in plan, on the face of an elastic wedge is investigated in a three-dimensional formulation for different types of boundary conditions on the other face. The wedge material is assumed to be incompressible. An asymptotic solution is obtained for punches which are relatively distant from one another and from the edge of the wedge. For the case when the punches are arranged relatively close to the edge of the wedge (or reach the edge, the contact area is unknown) the numerical method of boundary integral equations is used. The mutual effect of the punches is estimated by means of calculations. The asymptotic solution of the generalized Galin problem, concerning the effect of a concentrated force applied on the edge of the three-dimensional wedge on the contact pressure distribution under a circular punch relatively far from the edge, is obtained.     

A three-dimensional analog of the Tricomi problem for a parabolic-hyperbolic equation

For a parabolic-hyperbolic equation, we study the three-dimensional analog of the Tricomi problem with a noncharacteritic plane on which the type of the equation changes. The uniqueness of the solution to the problem is proved by the method of a priori estimates, and the existence of a solution is reduced to the existence of a solution to a Volterra integral equation of the second kind.     

A volume integral equation method for periodic scattering problems for anisotropic Maxwell's equations

This paper presents a volume integral equation method for an electromagnetic scattering problem for three-dimensional Maxwell's equations in the presence of a biperiodic, anisotropic, and possibly discontinuous dielectric scatterer. Such scattering problem can be reformulated as a strongly singular volume integral equation (i.e., integral operators that fail to be weakly singular). In this paper, we firstly prove that the strongly singular volume integral equation satisfies a Gårding-type estimate in standard Sobolev spaces. Secondly, we rigorously analyze a spectral Galerkin method for solving the scattering problem. This method relies on the periodization technique of Gennadi Vainikko that allows us to efficiently evaluate the periodized integral operators on trigonometric polynomials using the fast Fourier transform (FFT). The main advantage of the method is its simple implementation that avoids for instance the need to compute quasiperiodic Green's functions. We prove that the numerical solution of the spectral Galerkin method applied to the periodized integral equation converges quasioptimally to the solution of the scattering problem. Some numerical examples are provided for examining the performance of the method.     

Boundary Integral Equation Solution of Non-linear Plane Potential Problems

This paper presents the boundary integral equation (BIE) formulation,and numerical solution procedure for two-dimensional problemsgoverned by Laplac's equation and subject to non-linear boundaryconditions. The introduction of non-linear terms constitutesa fundamental extension of the BIE method, as previous applicationshave been restricted entirely to linear problems. Furthermore,non-linearities necessitate the use of iterative solution techniqueswhich present the conceptual disadvantage that a solution isnot guaranteed. However, such difficulties were not encounteredwith the Newton—Raphson method employed in this study.The various features of the BIE technique are illustrated bythe application to a physical problem which is of significancein heat exchanger design.     

Existence of a Fixed Point of a Nonsmooth Function Arising in Numerical Mechanics

A recent work (Acary et al. 2010) introduces a formulation as a nonsmooth fixed-point problem of a basic problem in numerical mechanics (namely the dynamical Coulomb friction problem in finite dimension with discretized time). Using this new formulation, the existence of a solution to the problem and its numerical resolution are then guaranteed under a strong assumption on the data of this problem. In this paper, we show that the fixed point problem admits solution under a natural, weaker assumption. This existence proof uses a perturbation argument combined with continuity properties of a set-valued mapping associated with the constraints of the problem.     

The Infinite-Dimensional Optimal Filtering Problem with Mobile and Stationary Sensor Networks

In this article, we introduce a framework to address filtering and smoothing with mobile sensor networks for distributed parameter systems. The main problem is formulated as the minimization of a functional involving the trace of the solution of a Riccati integral equation with constraints given by the trajectory of the sensor network. We prove existence and develop approximation of the solution to the Riccati equation in certain trace-class spaces. We also consider the corresponding optimization problem. Finally, we employ a Galerkin approximation scheme and implement a descent algorithm to compute optimal trajectories of the sensor network. Numerical examples are given for both stationary and moving sensor networks.     

Pressure of a circular stamp on a composite layer

The method of dual integral equations is used to obtain a solution to the problem of a rigid circular stamp pressing on an elastic composite layer, with a cylindrical surface separating the materials. A large number of papers have already been published, dealing with the mechanics of multilayered media in which the surfaces separating the layers from each other do not intersect the outer boundary (see references in /1/). The formulation and methods of solution of the fundamental boundary value problems can be found for such media in the monographs /2,3/.

Considerably less attention has been given to the study of the boundary value problems for composite media in which the surfaces separating the layers do intersect the outer boundary. The authors of /4, 5/ call such media the regions with transverse (vertical) layer folding. Out of the publications dealing with the methods of solving contact problems for transversely layered regions, attention should be drawn to /4–11/.     

The dynamics of the weakly perturbed motion of a liquid-filled gyroscope and the control problem

The Cauchy problem for the motion of a dynamically symmetrical rigid body with a cavity, filled with an ideal liquid, which is perturbed from uniform rotation, is considered in a linear formulation. The problem of the simultanious solution of the equations of hydrodynamics and the mechanics of a rigid body is reduced to the solution of an eigenvalue problem which depends solely on the geometry of the cavity and the subsequent integration of a system of differential equations.     

A computer-assisted proof of the existence of solutions to a boundary value problem with an integral boundary condition

In this paper, we present a computer-assisted method that establishes the existence and local uniqueness of a stationary solution to the viscous Burgers’ equation. The problem formulation involves a left boundary condition and one integral boundary condition, which is a variation of a previous approach.     

A comparison of some inverse methods for estimating the initial condition of the heat equation

In this work we analyze two explicit methods for the solution of an inverse heat conduction problem and we confront them with the least-squares method, using for the solution of the associated direct problem a classical finite difference method and a method based on an integral formulation. Finally, the Tikhonov regularization connected to the least-squares criterion is examined. We show that the explicit approaches to this inverse heat conduction problem will present disastrous results unless some kind of regularization is used.     

Wavelet Solution of Plane Elasticity Problem in the Upper Half-plane

The plane elasticity problem includes plane strain problem and plane stress problem which are widely applied in mechanics and engineering. In this article, we first reduce the plane elasticity problem in the upper half-plane into natural boundary integral equation and then apply wavelet-Galerkin method to deal with the numerical solution of the natural boundary integral equation. The test and trial functions used are the scaling basis functions of Shannon wavelet. In our fast algorithm, the computational formulae of entries of the stiffness matrix yield simple close-form and only 3 K entries need to be computed for one 4 K ‐ 4 K stiffness matrix.     

Finding the coefficients in the new representation of the solution of the Riemann–Hilbert problem using the Lauricella function

The solution of the Riemann–Hilbert problem for an analytic function in a canonical domain for the case in which the data of the problem is piecewise constant can be expressed as a Christoffel–Schwartz integral. In this paper, we present an explicit expression for the parameters of this integral obtained by using a Jacobi-type formula for the Lauricella generalized hypergeometric function F _D ^(N). The results can be applied to a number of problems, including those in plasma physics and the mechanics of deformed solids.