Three-dimensional elasticity based on quaternion-valued potentials

One of the most fruitful and elegant approach (known as Kolosov–Muskhelishvili formulas) for plane isotropic elastic problems is to use two complex-valued holomorphic potentials. In this paper, the algebra of real quaternions is used in order to propose in three dimensions, an extension of the classical Muskhelishvili formulas. The starting point is the classical harmonic potential representation due to Papkovich and Neuber. Alike the classical complex formulation, two monogenic functions very similar to holomorphic functions in 2D and conserving many of interesting properties, are used in this contribution. The completeness of the potential formulation is demonstrated rigorously. Moreover, body forces, residual stress and thermal strain are taken into account as a left side term. The obtained monogenic representation is compact and a straightforward calculation shows that classical complex representation for plane problems is embedded in the presented extended formulas. Finally the classical uniqueness problem of the Papkovich–Neuber solutions is overcome for polynomial solutions by fixing explicitly linear dependencies.     

The complex variable boundary element method for solving sandwich plate problems

The objective of this paper is to develop a new complex variable boundary element method for sandwich plates of Reissner's type and Hoff's type. The general solution of Helmhotz equation in complex field is given. Based on the Vekua's complex integral representation of the analytic function, the new boundary integral equations are formulated. The density function in the integral equation is determined directly by boundary element method. Some standard examples are presented, and the results of numerical solutions are accurate everywhere in the plate. The approach presented is only applicable for bounded simply connected regions. The project is supported by the National Science Foundation of China.     

正交异性体反平面问题另一形式自相似解的推导

正交异性体的反平面运动方程有几种表达方式，而以自相似解表达方式尤为简单。由于在弹性动力学问题的计算中可以获得解析解，使得考虑的问题相应地简化，并具有一定的普遍性，因此对自相似形式的解的推导具有重要的意义。本文针对这一问题进行研究，利用复变函数论的方法导出另一形式自相似解的一般表达式。应用该法可以很容易地将所讨论的问题转化为Riemann—Hilbert问题，而后一问题可以用同通常的Muskhelishvili方法求解，并且可以相当简单地得到问题的闭合解。这些解在断裂动力学以及弹性动力学问题当中具有重要的应用价值。     

Quasiperiodic solutions of the Navier-Stokes equations induced by the quasiperiodic shape of the boundaries of a two-dimensional flow domain

Steady quasiperiodic solutions of the Navier-Stokes equations in an infinite two-dimensional layer with quasiperiodic boundaries are obtained on the Reynolds number range 0 < Re* < 200. The calculations are performed using a spectral-difference method based on the representation of the quasiperiodic solutions in the form of convergent double Fourier series. The properties of these solutions and the distinctive features of their spectra are studied and their fundamental differences from periodic solutions are shown. The possibility of applying the quasiperiodic solutions for modeling flows in fractal layers is discussed.     

Improving the robustness of the control volume finite element method with application to multiphase porous media flow

Control volume finite element methods (CVFEMs) have been proposed to simulate flow in heterogeneous porous media because they are better able to capture complex geometries using unstructured meshes. However, producing good quality meshes in such models is nontrivial and may sometimes be impossible, especially when all or parts of the domains have very large aspect ratio. A novel CVFEM is proposed here that uses a control volume representation for pressure and yields significant improvements in the quality of the pressure matrix. The method is initially evaluated and then applied to a series of test cases using unstructured (triangular/tetrahedral) meshes, and numerical results are in good agreement with semianalytically obtained solutions. The convergence of the pressure matrix is then studied using complex, heterogeneous example problems. The results demonstrate that the new formulation yields a pressure matrix than can be solved efficiently even on highly distorted, tetrahedral meshes in models of heterogeneous porous media with large permeability contrasts. The new approach allows effective application of CVFEM in such models.     

Dimension Reduction in the Context of Structured Deformations

In the present paper the basic boundary value problems (BVPs) of the full coupled linear theory of elasticity for triple porosity materials are investigated by means of the potential method (boundary integral equation method) and some basic results of the classical theory of elasticity are generalized. In particular, the Green’s identities and the formula of Somigliana type integral representation of regular vector and regular (classical) solutions are presented. The representation of Galerkin type solution is obtained and the completeness of this solution is established. The uniqueness theorems for classical solutions of the internal and external BVPs are proved. The surface (single-layer and double-layer) and volume potentials are constructed and their basic properties are established. Finally, the existence theorems for classical solutions of the BVPs are proved by means of the potential method and the theory of singular integral equations.     

Analytic test cases for three-dimensional hydrodynamic models

Exact periodic solutions are generated for the 3-D hydrodynamic equations in linearized form. A linear slip condition is enforced at the bottom, based on the velocity at the bottom. It is shown that the bottom stress can be equivalently expressed in terms of the vertically averaged velocity, and expressions for this bottom stress coefficient are derived in terms of the primary parameters of the problem. As a result, the three-dimensional structure may be assembled from conventional solutions to (a) the 1-D vertical diffusion equation; and (b) the 2-D vertically averaged shallow water equations. In the latter, the bottom stress effects are shown to be complex and frequency-dependent, and an additional rotational term is required for their representation.     

The use of displacement potentials in second order elasticity

This paper is concerned with the use of a representation in terms of displacement potentials in second order elasticity for equilibrium problems of homogeneous and isotropic materials. After justifying the adoption of an existing representation for linear elasticity for the purpose at hand, appropriate representations for solutions of second order elasticity problems in terms of displacement potentials (for both compressible and incompressible materials) are discussed. The use of the representations in obtaining complete solutions for equilibrium boundary-value problems is then illustrated by application to two examples of plane strain problems of compressible materials.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

On the representation of elastic displacement fields in terms of three harmonic functions

In this paper the representation of displacement fields in linear elasticity in terms of harmonic functions is considered. In the original work of Papkovich and Neuber four harmonic functions were presented with a subsequent reduction to three on the grounds that only three are sufficient for the representation of displacements fields. This reduction is unsubstantiated and several authors have investigated the generality of the Papkovich-Neuber solutions. The paper derives by simple means the conditions under which it is possible to omit one of the four harmonic functions and considers the significance of the subsequent three function form.     

Diagonalization of a three-dimensional system of equations in terms of displacements of the linear theory of elasticity of transversely isotropic media

A dynamic three-dimensional system of linear equations in terms of displacements of the theory of elasticity of transversely isotropic media is given explicit expressions for phase velocities and polarization vectors of plane waves. All the longitudinal normals are found. For some values of the elasticity moduli, the system of equations is reduced to a diagonal shape. For static equations, all the conditions of the system ellipticity are determined. Two new representations of displacements through potential functions that satisfy three independent quasi-harmonic equations are given. Constraints on elasticity moludi, at which the corresponding coefficients in these representations are real, different, equal, or complex, are determined. It is shown that these representations are general and complete. Each representation corresponds to a recursion (symmetry) operator, i.e., a formula of production of new solutions.     

Periodic behavior of a nonlinear dynamical system

Nonlinear dynamical systems, being more of a realistic representation of nature, could exhibit a somewhat complex behavior. Their analysis requires a thorough investigation into the solution of the governing differential equations. In this paper, a class of third order nonlinear differential equations has been analyzed. An attempt has been made to obtain sufficient conditions in order to guarantee the existence of periodic solutions. The results obtained from this analysis are shown to be beneficial when studying the steady-state response of nonlinear dynamical systems. In order to obtain the periodic solutions for any form of third order differential equations, a computer program has been developed on the basis of the fourth order Runge-Kutta method together with the Newton-Raphson algorithm. Results obtained from the computer simulation model confirmed the validity of the mathematical approach presented for these sufficient conditions.     

Two-dimensional piezoelectricity. Part I: eigensolutions of nondegenerate and degenerate materials

General solutions of two-dimensional piezoelectricity, which yield all solutions of 2-D boundary values problems, are obtained by combining four complex conjugate pairs of independent eigensolutions, each containing an arbitrary analytic function. The forms of representation are fundamentally different for 14 different classes of nondegenerate and degenerate piezoelectric materials, as determined by the multiplicity and types of eigenvalues. Degenerate materials possess high-order eigensolutions, in which the eigenvectors of equal and lower orders are intrinsically coupled. Such coupling is nonexistent in nondegenerate cases including the well-known and analytically simple case with no multiple eigenvalues. The present analysis is drastically simplified by using the compliance-based formalism, instead of the stiffness-based, extended Eshelby–Stroh formalism. Explicit expressions are obtained for the eigensolutions, the pseudometrics, and the intrinsic tensors characterizing piezoelectric materials of every type.     

A marker-and-cell approach to viscoelastic free surface flows using the PTT model

This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.     

The fundamental solution in the theory of shallow shells

A series representation for the fundamental solution of the shallow shell equations is obtained by means of a plane-wave decomposition of the Dirac δ-function. From this solution we can produce the singular solutions which correspond to concentrated forces, couples and thermal hot spots applied to a shallow shell with an arbitrary quadratic middle surface. The solutions converge for the entire range of the Gaussian curvature. Numerical results are presented for the case of a concentrated normal force acting on infinite shells having positive, zero or negative Gaussian curvature.     

On use of the Dirac delta function in the vibration analysis of elastic structures

The Dirac delta function has often been employed to represent the amplitude of concentrated harmonic forces in the analysis of vibration of elastic structures such as beams and plates. It is known that this function, as represented by a truncated Fourier series, does not provide a true representation of a concentrated force, nevertheless, it is frequently employed and good convergence is usually, though not always, encountered in solutions thereby obtained. In this paper, the nature of the function is discussed and for illustrative purposes it is used to obtain series solutions for some selected beam and plate free vibration problems. In some cases problems are chosen for which exact solutions are already obtainable by analytical means. This permits powerful checks to be made on rates of convergence experienced when the series solutions are investigated. Rates of convergence are discussed in detail and it is explained why convergence is to be expected when analyzing certain families of problems when employing this function and a lack of convergence is to be expected in others.     

2D unstructured mesh generation for oceanic and coastal tidal models from a localized truncation error analysis with complex derivatives

A method for computing target element size for tidal, shallow water flow is developed and demonstrated. The method, Localized truncation error analysis with complex derivatives (LTEA-CD) utilizes localized truncation error estimates of the linearized shallow water momentum equations consisting of complex derivative terms. This application of complex derivatives is the chief way in which the method differs from a similar existing method, LTEA. It is shown that LTEA-CD produces results that are essentially equivalent to those of LTEA (which in turn has been demonstrated to be capable of producing practicable target element sizes) with reduced computational cost. Moreover, LTEA-CD is capable of computing truncation error and corresponding target element sizes at locations up to and including the boundary, whereas LTEA can be applied only on the interior of the model domain. We demonstrate the convergence of solutions over meshes generated with LTEA-CD using an idealized representation of the western North Atlantic Ocean, Caribbean Sea and Gulf of Mexico.     

Traveling Wave Solutions in a Generalized Theory for Macroscopic Capillarity

One-dimensional traveling wave solutions for imbibition processes into a homogeneous porous medium are found within a recent generalized theory of macroscopic capillarity. The generalized theory is based on the hydrodynamic differences between percolating and nonpercolating fluid parts. The traveling wave solutions are obtained using a dynamical systems approach. An exhaustive study of all smooth traveling wave solutions for primary and secondary imbibition processes is reported here. It is made possible by introducing two novel methods of reduced graphical representation. In the first method the integration constant of the dynamical system is related graphically to the boundary data and the wave velocity. In the second representation the wave velocity is plotted as a function of the boundary data. Each of these two graphical representations provides an exhaustive overview over all one-dimensional and smooth solutions of traveling wave type, that can arise in primary and secondary imbibition. Analogous representations are possible for other systems, solution classes, and processes.     

A three dimensional field formulation,and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

We present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains (Toupin Arch. Ration. Mech. Anal., 11 (1962) 385). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comput. Methods Appl. Mech. Eng., 278 (2014) 705). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.