Mixed boundary value problem of elasticity for a quarter space

We consider the static elasticity problem for a quarter space with zero displacements on one of its surfaces and with given stresses on the other. The method for solving this problem is based on the use of newunknown functions in the formof a linear combination of the desired displacements, which reduces the system of three Lamé equations to two equations to be solved simultaneously and one equation to be solved separately. The exact solution of this problem was obtained earlier by the same method [1]. But it was shown in [2] that such a solution is exact only under certain restrictions on the given functions. In the present paper, the solution of this problem is constructed without restrictions on the given functions, which necessitates solving a one-dimensional integro-differential equation; this can be done approximately by the orthogonal polynomial method. We present numerical results obtained on the basis of our solution.     

A code-independent technique for computational verification of fluid mechanics and heat transfer problems

The goal of this paper is to present a versatile framework for solution verification of PDE＇s. We first generalize the Richardson Extrapolation technique to an optimized extrapolation solution procedure that constructs the best consistent solution from a set of two or three coarse grid solution in the discrete norm of choice. This technique generalizes the Least Square Extrapolation method introduced by one of the author and W. Shyy. We second establish the conditioning number of the problem in a reduced space that approximates the main feature of the numerical solution thanks to a sensitivity analysis. Overall our method produces an a posteriori error estimation in this reduced space of approximation. The key feature of our method is that our construction does not require an internal knowledge of the software neither the source code that produces the solution to be verified. It can be applied in principle as a postprocessing procedure to off the shelf commercial code. We demonstrate the robustness of our method with two steady problems that are separately an incompressible back step flow test case and a heat transfer problem for a battery. Our error estimate might be ultimately verified with a near by manufactured solution. While our pro- cedure is systematic and requires numerous computation of residuals, one can take advantage of distributed computing to get quickly the error estimate.     

A formula of solution to the integral of rational functions

ＡＦＯＲＭＵＬＡＯＦＳＯＬＵＴＩＯＮＴＯＴＨＥＩＮＴＥＧＲＡＬＯＦＲＡＴＩＯＮＡＬＦＵＮＣＴＩＯＮＳＧｕｉＺｕ－ｈｕａ（桂祖华）（ＴｈｅＤｅｐａｒｔｍｅｎｔｏｆＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓ,ＳｈａｎｇｈａｉＪｉａｏ－ｔｏｎｇＵｎｉｖｅｒｓｉｔ...     

Diffusional growth of cloud particles: existence and uniqueness of solutions

Diffusional growth of cloud particles is commonly described by a coupled system of parabolic equations and ordinary differential equations. The Dirichlet boundary condition for the parabolic equation is obtained from the solution of the ordinary differential equations, but this solution itself depends on the solution of the parabolic equations. We first present the governing equations describing diffusional growth of cloud particles. In a second step, we consider a simplified model problem, motivated by the diffusional growth equations. The main difference between the simplified model problem and the diffusional growth equations consists in neglecting the dependence of the domain for the parabolic equations on the solution. For the model problem, we show unique solvability using a fixed point method. Finally, we discuss application of the main result for the model problem to the diffusional growth equations and illustrate these equations with the help of a numerical solution.     

Structure of polynomial solutions of elasticity equations in stresses

For an isotropic elasticity problem in stresses in three-dimensional space minus the origin, we study solutions that have the singularity 1/r ² and, after the multiplication by r ², polynomially depend on the direction cosines. In this polynomial class, for the equilibrium equation we write out the general solution that is a statically admissible (in the sense of Castigliano) solution of the Kelvin problem. We show that if one or several Beltrami equations are not satisfied, then the classical Kelvin solution becomes nonunique. A method for constructing nonunique solutions of this kind is given. The equivalence of various statements of the elasticity problem in stresses is discussed. For the problem on the action of a lumped force at the vertex of an arbitrary conical elastic body, we write out the exact solution in stresses for the case of an incompressible material. The solution for a compressible material is represented in the form of series in a parameter characterizing the deviation of the Poisson ratio from 1/2. We obtain iterative chains of problems in stresses and conditions for the finiteness of these chains. We also analyze the realizability of a linear-fractional dependence of the solution in stresses on the Poisson ratio.     

Multipoint problem for higher-order parabolic equations in a parallelepiped

We study the well-posedness of a problem for a Petrovskii-parabolic equation with coefficients depending on the space coordinates and with multipoint conditions with respect to the time variable. We establish conditions for the existence and uniqueness of the classical solution of the problem. For the proof of the existence of a solution of the problem, the method of divided differences is used. We prove a metric-type theorem on lower bounds for the small denominators that appear in the construction of the solution.     

Motion of a strong shock front in a two-dimensional gas layer with moving internal boundary

We examine the problem of planar one-dimensional motion of a strong shock wave with moving internal boundary in which the initial position of the front, its intensity, the mass of the gas involved in the motion, and the energy contained in this gas are known. The problem is not self-similar and its exact solution, which involves working with partial differential equations, presents serious difficulties. In the following we determine the law of shock-front motion in this problem via the method of [1], which makes it possible to find a system of ordinary differential equations for the problem. The method is based on an initial specification of the power-law coupling between the dimensionless Lagrangian and Eulerian variables and replacement of the energy equation by this coupling and the energy integral. The solution is sought in the first approximation.     

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.     

Acoustic scattering and radiation problems, and the null-field method

The best known methods for solving the scattering and radiation problems of acoustics are integral-equation methods. However, it is also known that the simplest of these methods yield equations which are not uniquely solvable at certain discrete sets of frequencies (the irregular frequencies). In this paper, we shall analyse an alternative method (the null-field method, or T-matrix method). We prove that the infinite system of null-field equations always has precisely one solution, i.e. the unphysical irregular frequencies do not occur with this method. Moreover, we also prove that the solution of the original boundary-value problem can always be determined (at any point exterior to the scatterer) from the solution of the null-field equations. We prove these results in two dimensions, for two radiation problems (the exterior Neumann problem and the exterior Dirichlet problem) and two scattering problems (scattering by a sound-hard body and scattering by a sound-soft body); similar results can be proved in three dimensions. We also prove some subsidiary results, concerning the solvability of certain boundary integral equations and the completeness of certain sets of radiating wave-functions, and give a discussion of related numerical techniques.     

Use of the energy integral in optimal control of the spacecraft spatial attitude

In this paper, we study the problem of spacecraft optimal rotation into the position with given attitude. The rotation takes a fixed time. The control program is optimized by minimizing the energy integral. We present the analytic solution of this problem and obtain formalized equations and computational expressions for constructing the optimal rotation program. We also present an example and the results of mathematical simulation of the spacecraft motion dynamics under optimal control, which demonstrate the practical realizability of our method for the spacecraft spatial attitude control.     

On the construction of asymptotic solutions of two-point boundary-value problems for degenerate singularly perturbed systems of differential equations

We consider the problem of the existence of a solution of a two-point boundary-value problem for degenerate singularly perturbed linear systems of differential equations. We obtain asymptotic formulas for this solution.     

Differential-algebraic equations of programmed motions of Lagrangian dynamical systems

We suggest a method for constructing the dynamic equations of manipulator systems in canonical variables. The system of differential dynamic equations has an integral manifold corresponding to the holonomic and nonholonomic constraint equations. The controls are determined so as to ensure the stability of this manifold. We state conditions for the exponential stability of the manifold and for constraint stabilization when solving the dynamic equations numerically by a simplest difference method. We also present the solution of the problem of control of a plane two-link manipulator.     

Axisymmetric strain of a thick-walled pipe made of a highly elastic material

We reduce the plane strain problem to a nonlinear elasticity problem for inhomogeneous bodies by choosing a new form of the elastic potential, whose parameters are determined from the data known in the literature. By using the geometric linearization method, we reduce that problem to a sequence of linear elasticity problems for inhomogeneous bodies. We obtain an analytic solution of the corresponding linear elasticity problem in the case of an arbitrary continuously differentiable dependence of the shear modulus on the radial coordinate. We determine the pipe stress-strain state and parameters in the case of finite and large strains for given sets of initial data and estimate the accuracy of the solution thus obtained.     

Application of the Liouville–Green transformation to free vibrations of non-homogeneous rectangular membranes

We use the Liouville-Green transformation to derive a new approximate solution to Liouville’s problem. This approximate solution is then applied to the eigenvalue problem of free vibrations of non-homogeneous rectangular membranes. The examples show that the fundamental natural frequencies obtained in this work are in good agreement with those of previous studies. The main advantage of our method is that it is simple, easy to implement.     

Low-rank solution of an optimal control problem constrained by random Navier-Stokes equations

We develop a low-rank tensor decomposition algorithm for the numerical solution of a distributed optimal control problem constrained by two-dimensional time-dependent Navier-Stokes equations with a stochastic inflow. The goal of optimization is to minimize the flow vorticity. The inflow boundary condition is assumed to be an infinite-dimensional random field, which is parametrized using a finite- (but high-) dimensional Fourier expansion and discretized using the stochastic Galerkin finite element method. This leads to a prohibitively large number of degrees of freedom in the discrete solution. Moreover, the optimality conditions in a time-dependent problem require solving a coupled saddle-point system of nonlinear equations on all time steps at once. For the resulting discrete problem, we approximate the solution by the tensor-train (TT) decomposition and propose a numerically efficient algorithm to solve the optimality equations directly in the TT representation. This algorithm is based on the alternating linear scheme (ALS), but in contrast to the basic ALS method, the new algorithm exploits and preserves the block structure of the optimality equations. We prove that this structure preservation renders the proposed block ALS method well posed, in the sense that each step requires the solution of a nonsingular reduced linear system, which might not be the case for the basic ALS. Finally, we present numerical experiments based on two benchmark problems of simulation of a flow around a von Kármán vortex and a backward step, each of which has uncertain inflow. The experiments demonstrate a significant complexity reduction achieved using the TT representation and the block ALS algorithm. Specifically, we observe that the high-dimensional stochastic time-dependent problem can be solved with the asymptotic complexity of the corresponding deterministic problem.     

A robust numerical method for flow through a pipe driven by an oscillating pressure gradient

The problem of periodic flow of an incompressible fluid through a pipe, which is driven by an oscillating pressure gradient (e.g. a reciprocating piston), is investigated in the case of a large Reynolds number. This process is described by a singularly perturbed parabolic equation with a periodic right‐hand side, where the singular perturbation parameter is the viscosity ν. The periodic solution of this problem is a solution of the Navier–Stokes equations with cylindrical symmetry. We are interested in constructing a parameter‐robust numerical method for this problem, i.e. a numerical method generating numerical approximations that converge uniformly with respect to the parameter ν and require a bounded time, independent of the value of ν, for their computation. Our method comprises a standard monotone discretization of the problem on non‐standard piecewise uniform meshes condensing in a neighbourhood of the boundary layer. The transition point between segments of the mesh with different step sizes is chosen in accordance with the behaviour of the analytic solution in the boundary layer region. In this paper we construct the numerical method and discuss the results of extensive numerical experiments, which show experimentally that the method is parameter‐robust. Copyright © 2006 John Wiley & Sons, Ltd.     

Solution of the Eshelby problem in gradient elasticity for multilayer spherical inclusions

We consider gradient models of elasticity which permit taking into account the characteristic scale parameters of the material. We prove the Papkovich–Neuber theorems, which determine the general form of the gradient solution and the structure of scale effects. We derive the Eshelby integral formula for the gradient moduli of elasticity, which plays the role of the closing equation in the self-consistent three-phase method. In the gradient theory of deformations, we consider the fundamental Eshelby–Christensen problem of determining the effective elastic properties of dispersed composites with spherical inclusions; the exact solution of this problem for classical models was obtained in 1976.     

On complex potentials for finite plane deformation of a transversely isotropic material

We consider the finite deformation of plane equilibrium problem for a transversely isotropic layer, using the complex variable approach. We give the general expression for the pertinent complex potentials and state the corresponding fundamental problems. We discuss in detail the boundary value problem for fundamental problem-one. As an application of the espoused method, an analytical solution of “Lame's problem” for an infinite layer is obtained. The nonlinear effect of this is highlighted in the obtained figure.     

Effect of the physical nonlinearity of a material on the stress state of a medium with an elastic spherical inclusion under uniform loading

We have used the perturbation method as the basis for obtaining an approximate solution of the three-dimensional problem for a physically nonlinear elastic medium with an elastic inclusion under uniform tension— compression. From this solution, we can obtain as a special case a solution for an elastic medium with a stress-free cavity and for an elastic medium with a rigid inclusion. We have plotted the normal and tangential stresses as a function of the radius and the ratio of shear moduli for the inclusion and the medium. We have investigated their behavior under different loading conditions. Translated from Prikladnaya Mekhanika, Vol. 34, No. 11, pp. 46–51, November, 1998.     

Optimal spacecraft terminal attitude control synthesis by the quaternion method

We consider an optimal control problem of spacecraft retargeting from an arbitrary initial position to a given terminal angular position. The objective is to minimize a path functional. Using the quaternion method, we obtain an analytic solution of this problem and present formalized equations and computational expressions for synthesizing an optimal manoeuvre program. We also present an example and the results of mathematical simulation of the spacecraft motion dynamics under an optimal control, which demonstrate the practical implementability of the control method developed in the present paper.