Three-dimensional thermo-electro-elastic field in one-dimensional hexagonal piezoelectric quasi-crystal weakened by an elliptical crack

In this paper, an analytical solution is developed for the problem of an infinite 1D hexagonal piezoelectric quasi-crystal medium weakened by an elliptical crack and subjected to mixed loads on the crack surfaces. The mixed loads comprise the phonon pressure, phason pressure, electric displacement, and temperature increment, and the crack surfaces can be electrically permeable or impermeable. Based on a general solution, combined with the generalized potential theory, the steady-state 3D thermo-electro-elastic field variables in the quasi-crystal are obtained in terms of elliptic integral functions and elementary functions. Several important physical quantities on the cracked plane, such as the generalized crack surface displacements, normal stresses, and stress intensity factors, are derived in closed forms. An illustrative numerical calculation verifies the presented analytical solution and shows the distribution of the 3D thermo-electro-elastic field. It is indicated that the influence of the phason field on the result is pronounced, especially for the electric field variables, and the electric permeability of crack surfaces has a significant effect on the electric displacement intensity factor at the crack tip.

Nonlinear oscillations in a constant gravitational field

Exploring non-linear oscillations is a challenging task since the related differential equations cannot be directly solved in terms of elementary functions. Thus, complicated mathematical or numerical methods are usually employed to find accurate or approximate expressions that describe the behavior of the system with respect to time. In this paper, the vertical oscillations of an object under the influence of its weight and an opposite force with magnitude F=cyn, where n>0 are being explored. Accurate and approximate simple solutions regarding the object’s position with respect to time are presented and the dependence of the oscillation’s period from the oscillation’s range of displacements and the exponent n is revealed. In addition, the special case in which n=3/2 (which describes the oscillation of a rigid sphere on an elastic half space) is also highlighted. Lastly, it is shown that similar cases (such as the case of a force with magnitude F=kx+λx2) can be also treated using the same approach.

Parameter N Analysis on the Rational-Talbot Algorithm for Numerical Inversion of Laplace Transform

This article investigates the numerical inversion of the Laplace Transform by the Rational-Talbot method and analyzes the influence on the variation of the free parameter N established by the technique when applied to certain functions. The set of elementary functions, for which the method is tested, has exponential and oscillatory characteristics. Based on the results obtained, it was concluded that the Rational-Talbot method is e cient for the inversion of decreasing exponential functions. At the same time, to perform the inversion process effectively for trigonometric forms, the algorithm requires a greater amount of terms in the sum. For higher values of N, the technique works well. In fact, this is observed in inverting the functions transform, that combine trigonometric and polynomial factors. The method numerical results have a good precision for the treatment of decreasing exponential functions when multiplied by trigonometric functions.

Calculation of integrals in MathPartner

We present the possibilities provided by the MathPartner service of calculating definite and indefinite integrals. MathPartner contains software implementation of the Risch algorithm and provides users with the ability to compute antiderivatives for elementary functions. Certain integrals, including improper integrals, can be calculated using numerical algorithms. In this case, every user has the ability to indicate the required accuracy with which he needs to know the numerical value of the integral. We highlight special functions allowing us to calculate complete elliptic integrals. These include functions for calculating the arithmetic-geometric mean and the geometric-harmonic mean, which allow us to calculate the complete elliptic integrals of the first kind. The set also includes the modified arithmetic-geometric mean, proposed by Semjon Adlaj, which allows us to calculate the complete elliptic integrals of the second kind as well as the circumference of an ellipse. The Lagutinski algorithm is of particular interest. For given differentiation in the field of bivariate rational functions, one can decide whether there exists a rational integral. The algorithm is based on calculating the Lagutinski determinant. This year we are celebrating 150th anniversary of Mikhail Lagutinski.

Subsurface stresses in an elliptical Hertzian contact

The traditional solution for the stresses below an elliptical Hertzian contact expresses the results in terms of incomplete Legendre elliptic integrals, so are necessarily based on the length of the semi-major axis a and the axis ratio k. The result is to produce completely different equations for the stresses in the x and y directions; and although these equations are now well-known, their derivation from the fundamental, symmetric, integrals is far from simple. When instead Carlson elliptic integrals are used, they immediately match the fundamental integrals, allowing the equations for the stresses to treat the two semi-axes equally, and so providing a single equation where two were needed before. The numerical evaluation of the Carlson integrals is simple and rapid, so the result is that more convenient answers are obtained more conveniently. A bonus is that the temptation to record the depth of the critical stresses as a fraction of the length of the semi-major axis is removed. Thomas and Hoersch’s method of finding all the stresses along the axis of symmetry has been extended to determine the full set of stresses in a principal plane. The stress patterns are displayed, and a comparison between the answers for the planes of the major and minor semiaxes is made. The results are unchanged from those found from equations given by Sackfield and Hills, but not previously evaluated. The present equations are simpler, not only in the simpler elliptic integrals, but also for the “tail” of elementary functions.

Non-Debye Relaxations: The Ups and Downs of the Stretched Exponential vs. Mittag–Leffler’s Matchings

Experimental data collected to provide us with information on the course of dielectric relaxation phenomena are obtained according to two distinct schemes: one can measure either the time decay of depolarization current or use methods of the broadband dielectric spectroscopy. Both sets of data are usually fitted by time or frequency dependent functions which, in turn, may be analytically transformed among themselves using the Laplace transform. This leads to the question on comparability of results obtained using just mentioned experimental procedures. If we would like to do that in the time domain we have to go beyond widely accepted Kohlrausch–Williams–Watts approximation and become acquainted with description using the Mittag–Leffler functions. To convince the reader that the latter is not difficult to understand we propose to look at the problem from the point of view of objects which appear in the stochastic processes approach to relaxation. These are the characteristic exponents which are read out from the standard non-Debye frequency dependent patterns. Characteristic functions appear to be expressed in terms of elementary functions whose asymptotics is simple. This opens new possibility to compare behavior of functions used to describe non-Debye relaxations. It turnes out that the use of Mittag-Leffler function proves very convenient for such a comparison.

Simple Equations Method and Non-Linear Differential Equations with Non-Polynomial Non-Linearity

We discuss the application of the Simple Equations Method (SEsM) for obtaining exact solutions of non-linear differential equations to several cases of equations containing non-polynomial non-linearity. The main idea of the study is to use an appropriate transformation at Step (1.) of SEsM. This transformation has to convert the non-polynomial non- linearity to polynomial non-linearity. Then, an appropriate solution is constructed. This solution is a composite function of solutions of more simple equations. The application of the solution reduces the differential equation to a system of non-linear algebraic equations. We list 10 possible appropriate transformations. Two examples for the application of the methodology are presented. In the first example, we obtain kink and anti- kink solutions of the solved equation. The second example illustrates another point of the study. The point is as follows. In some cases, the simple equations used in SEsM do not have solutions expressed by elementary functions or by the frequently used special functions. In such cases, we can use a special function, which is the solution of an appropriate ordinary differential equation, containing polynomial non-linearity. Specific cases of the use of this function are presented in the second example.

New solution generating algorithm for isotropic static Einstein-Gauss-Bonnet metrics

The static isotropic gravitational field equation, governing the geometry and dynamics of stellar structure, is considered in Einstein–Gauss–Bonnet (EGB) gravity. This is a nonlinear Abelian differential equation which generalizes the simpler general relativistic pressure isotropy condition. A gravitational potential decomposition is postulated in order to generate new exact solutions from known solutions. The conditions for a successful integration are examined. Remarkably we generate a new exact solution to the Abelian equation from the well known Schwarzschild interior seed metric. The metric potentials are given in terms of elementary functions. A physical analysis of the model is performed in five and six spacetime dimensions. It is shown that the six-dimensional case is physically more reasonable and is consistent with the conditions restricting the physics of realistic stars.

Option Pricing under Randomised GBM Models

By employing a randomisation procedure on the variance parameter of the standard geometric Brownian motion (GBM) model, we construct new families of analytically tractable asset pricing models. In particular, we develop two explicit families of processes that are respectively referred to as the randomised gamma (G) and randomised inverse gamma (IG) models, both characterised by a shape and scale parameter. Both models admit relatively simple closed-form analytical expressions for the transition density and the no-arbitrage prices of standard European-style options whose Black-Scholes implied volatilities exhibit symmetric smiles in the log-forward moneyness. Surprisingly, for integer-valued shape parameter and arbitrary positive real scale parameter, the analytical option pricing formulas involve only elementary functions and are even more straightforward than the standard (constant volatility) Black-Scholes (GBM) pricing formulas. Moreover, we show some interesting characteristics of the risk-neutral transition densities of the randomised G and IG models, both exhibiting fat tails. In fact, the randomised IG density only has finite moments of the order less than or equal to one. In contrast, the randomised G density has a finite first moment with finite higher moments depending on the time-to-maturity and its scale parameter. We show how the randomised G and IG models are efficiently and accurately calibrated to market equity option data, having pronounced implied volatility smiles across several strikes and maturities. We also calibrate the same option data to the wellknown SABR (Stochastic Alpha Beta Rho) model.