A Fourth Order Accurate Finite Difference Scheme for the Elastic Wave Equation in Second Order Formulation

We present a fourth order accurate finite difference method for the elastic wave equation in second order formulation, where the fourth order accuracy holds in both space and time. The key ingredient of the method is a boundary modified fourth order accurate discretization of the second derivative with variable coefficient, (??(x)u _x ) _x . This discretization satisfies a summation by parts identity that guarantees stability of the scheme. The boundary conditions are enforced through ghost points, thereby avoiding projections or penalty terms, which often are used with previous summation by parts operators. The temporal discretization is obtained by an explicit modified equation method. Numerical examples with free surface boundary conditions show that the scheme is stable for CFL-numbers up to 1.3, and demonstrate a significant improvement in efficiency over the second order accurate method. The new discretization of (??(x)u _x ) _x has general applicability, and will enable stable fourth order accurate approximations of other partial differential equations as well as the elastic wave equation.     

The formula of the solution for some classes of initial boundary value problems for the hyperbolic equation with two independent variables

A new representation is proved of the solutions of initial boundary value problems for the equation of the form u _xx (x, t) + r(x)u _x (x, t) ? q(x)u(x, t) = u _tt (x, t) + μ(x)u _t (x, t) in the section (under boundary conditions of the 1st, 2nd, or 3rd type in any combination). This representation has the form of the Riemann integral dependent on the x and t over the given section.     

Compact Numerical Quadrature Formulas for Hypersingular Integrals and Integral Equations

In the first part of this work, we derive compact numerical quadrature formulas for finite-range integrals $I[f]=\int^{b}_{a}f(x)\,dx$ , where f(x)=g(x)|x?t| ^?? , ?? being real. Depending on the value of ??, these integrals are defined either in the regular sense or in the sense of Hadamard finite part. Assuming that g??C ^??[a,b], or g??C ^??(a,b) but can have arbitrary algebraic singularities at x=a and/or x=b, and letting h=(b?a)/n, n an integer, we derive asymptotic expansions for ${T}^{*}_{n}[f]=h\sum_{1\leq j\leq n-1,\ x_{j}\neq t}f(x_{j})$ , where x _j =a+jh and t??{x ₁,??,x _n?1}. These asymptotic expansions are based on some recent generalizations of the Euler?CMaclaurin expansion due to the author (A.?Sidi, Euler?CMaclaurin expansions for integrals with arbitrary algebraic endpoint singularities, in Math. Comput., 2012), and are used to construct our quadrature formulas, whose accuracies are then increased at will by applying to them the Richardson extrapolation process. We pay particular attention to the case in which ??=?2 and f(x) is T-periodic with T=b?a and $f\in C^{\infty}(-\infty,\infty)\setminus\{t+kT\}^{\infty}_{k=-\infty}$ , which arises in the context of periodic hypersingular integral equations. For this case, we propose the remarkably simple and compact quadrature formula $\widehat{Q}_{n}[f]=h\sum^{n}_{j=1}f(t+jh-h/2)-\pi^{2} g(t)h^{-1}$ , and show that $\widehat{Q}_{n}[f]-I[f]=O(h^{\mu})$ as h??0 ???>0, and that it is exact for a class of singular integrals involving trigonometric polynomials of degree at most n. We show how $\widehat{Q}_{n}[f]$ can be used for solving hypersingular integral equations in an efficient manner. In the second part of this work, we derive the Euler?CMaclaurin expansion for integrals $I[f]=\int^{b}_{a} f(x)dx$ , where f(x)=g(x)(x?t) ^?? , with g(x) as before and ??=?1,?3,?5,??, from which suitable quadrature formulas can be obtained. We revisit the case of ??=?1, for which the known quadrature formula $\widetilde{Q}_{n}[f]=h\sum^{n}_{j=1}f(t+jh-h/2)$ satisfies $\widetilde{Q}_{n}[f]-I[f]=O(h^{\mu})$ as h??0 ???>0, when f(x) is T-periodic with T=b?a and $f\in C^{\infty}(-\infty,\infty)\setminus\{t+kT\}^{\infty}_{k=-\infty}$ . We show that this formula too is exact for a class of singular integrals involving trigonometric polynomials of degree at most n?1. We provide numerical examples involving periodic integrands that confirm the theoretical results.     

Exact solution of mixed problems for variable coefficient one-dimensional diffusion equation

This paper deals with diffusion problems modeled by the equation a(t)u_xx = u_t, x > 0, t > 0, u(x, 0) = c(x) together with the boundary condition u(0, t) = b(t) or u_x(0, t) = b(t). By using Fourier transforms, existence conditions and exact solutions of the above mixed problems are given.     

Consistent initialization and perturbation analysis for abstract differential-algebraic equations

In this paper, we consider linear and time-invariant differential-algebraic equations (DAEs) Eẋ(t) = Ax(t) + f(t), x(0) = x ₀, where x(·) and f(·) are functions with values in Hilbert spaces X and Z. is assumed to be a bounded operator, whereas A is closed and defined on some dense subspace D(A). A transformation to a decoupling form leads to a DAE including an abstract boundary control system. Methods of infinite-dimensional linear systems theory can then be used to formulate sufficient criteria for an initial value being consistent with the given inhomogeneity. We will further derive estimates for the trajectory x(·) in dependence of the initial state x ₀ and the inhomogeneity f(·). In the theory of differential-algebraic equations, this is commonly known as perturbation analysis.     

Some remarks on Bianchi type-II, VIII, and IX models

Within the scope of anisotropic non-diagonal Bianchi type-II, VIII, and IX spacetimes it is shown that the off-diagonal components of the Einstein equations impose severe restrictions on the components of the energy-momentum tensor (EMT) in general. We begin with a metric with three functions of time, a(t), b(t), and c(t), and two spatial ones, f(z) and h(z). It is shown that if the EMT is assumed to be diagonal, and f = f(z), in all cosmological models in question b ∝ c, and the matter distribution is isotropic, i.e., T ₁ ¹ = T ₂ ² = T ₃ ³ . If f = const, which is a special case of BII models, the matter distribution may be anisotropic, but only the z axis is distinguished, and in this case b(t) is not necessarily proportional to c(t).     

Automata for Reduction Properties Solving

We introduce a new class of tree automata, which we call Reduction Automata (RA), and we use it to prove the decidability of the whole first order of theory of reduction (the theory of reduction is the set of true sentences built up with unary predicates red_t(x) "x is reducible by t", i.e. "some instance of t is a subterm of x"). So we link rewriting, logic and formal languages.     

Cylindrical decomposition for systems transcendental in the first variable

We present a generalization of the Cylindrical Algebraic Decomposition (CAD) algorithm to systems of equations and inequalities in functions of the form p(x,f₁(x),…,f_m(x),y₁,…,y_n), where p∈Q[x,t₁,…,t_m,y₁,…,y_n] and f₁(x),…,f_m(x) are real univariate functions such that there exists a real root isolation algorithm for functions from the algebra Q[x,f₁(x),…,f_m(x)]. In particular, the algorithm applies when f₁(x),…,f_m(x) are real exp-log functions or tame elementary functions.     

Transit cosmological models with domain walls in <Emphasis Type="Italic">f</Emphasis>(<Emphasis Type="Italic">R,T</Emphasis>) gravity

We study the physical behavior of the transition of a 5D perfect fluid universe from an early decelerating phase to the current accelerating phase in the framework of f(R, T) theory of gravity in the presence of domain walls. The fifth dimension is not observed because it is compact. To determine the solution of the field equations, we use the concept of a time-dependent deceleration parameter which yields the scale factor a(t) = sinh^1/n(αt), where n and α are positive constants. For 0 < n ≤ 1, this generates a class of accelerating models, while for n > 1 the universe attains a phase transition from an early decelerating phase to the present accelerating phase, consistent with the recent observations. Some physical and geometric properties of the models are also discussed.     

On the qualitative behaviors of solutions of third order nonlinear differential equations

We present some new results about oscillation and asymptotic behavior of solutions of third order nonlinear differential equations of the form

(r₂(t)(r₁(t)y^′)^′)^′+p(t)y^′+q(t)f(y(g(t)))=0.     

Comparison of the lattice Boltzmann and pseudo-spectral methods for decaying turbulence: Low-order statistics

We conduct a detailed comparison of the lattice Boltzmann equation (LBE) and the pseudo-spectral (PS) methods for direct numerical simulations (DNS) of the decaying homogeneous isotropic turbulence in a three-dimensional periodic cube. We use a mesh size of N³=128³ and the Taylor micro-scale Reynolds number 24.35?Re_λ?72.37, and carry out all simulations to t≈30τ₀, where τ₀ is the turbulence turnover time. In the PS method, the second-order Adam-Bashforth scheme is used to numerically integrate the nonlinear term while the viscous term is treated exactly. We compare the following quantities computed by the LBE and PS methods: instantaneous velocity u and vorticity ω fields, and statistical quantities such as, the total energy K(t) and the energy spectrum E(k,t), the dissipation rate ε(t), the root-mean-squared (rms) pressure fluctuation δp(t) and the pressure spectrum P(k,t), and the skewness and flatness of the velocity derivative. Our results show that the LBE method performs very well when compared to the PS method in terms of accuracy and efficiency: the instantaneous flow fields, u and ω, and all the statistical quantities — except the rms pressure fluctuation δp(t) and the pressure spectrum P(k,t) — computed from the LBE and PS methods agree well with each other, provided that the initial flow field is adequately resolved by both methods. We note that δp(t) and P(k,t) computed from the two methods agree with each other in a period of time much shorter than that for other quantities, indicating that the pressure field p computed by using the LBE method is less accurate than other quantities. The skewness and flatness computed from the LBE method contain high-frequency oscillations due to acoustic waves in the system, which are absent in PS methods. Our results indicate that the resolution requirement for the LBE method is δx/η₀?1.0, approximately twice of the requirement for PS methods, where δx and η₀ are the grid spacing and the initial Kolmogorov length, respectively. Overall, the LBE method is shown to be a reliable and accurate method for the DNS of decaying turbulence.     

Approximate Solutions of Polynomial Equations

In this paper, we introduce “approximate solutions" to solve the following problem: given a polynomial F(x, y) over Q, where x represents an n -tuple of variables, can we find all the polynomials G(x) such that F(x, G(x)) is identically equal to a constant c in Q ? We have the following: let F(x, y) be a polynomial over Q and the degree of y in F(x, y) be n. Either there is a unique polynomial g(x)  ∈ Q [ x ], with its constant term equal to 0, such that F(x, y)  = ∑_j = 0ⁿc_j(y − g(x))^jfor some rational numbers c_j, hence, F(x, g(x)  + a)  ∈ Q for all a ∈ Q, or there are at most t distinct polynomials g₁(x),⋯ , g_t(x), t ≤ n, such that F(x, g_i(x))  ∈ Q for 1  ≤ i ≤ t. Suppose that F(x, y) is a polynomial of two variables. The polynomial g(x) for the first case, or g₁(x),⋯ , g_t(x) for the second case, are approximate solutions of F(x, y), respectively. There is also a polynomial time algorithm to find all of these approximate solutions. We then use Kronecker’s substitution to solve the case of F(x, y).     

Stability of θ-methods for advanced differential equations with piecewise continuous arguments

This paper deals with the stability analysis of numerical methods for the solution of advanced differential equations with piecewise continuous arguments. We focus on the behaviour of the one-leg θ-method and the linear θ-method in the solution of the equation x′(t) = ax(t + a₀x([t]) + a₁x([t+1]), with real a, a₀, a₁ and [·] designates the greatest-integer function. The stability regions of two θ-methods are determined. The conditions under which the analytic stability region is contained in the numerical stability region are obtained and some numerical experiments are given.     

Scalar Field Cosmology in <Emphasis Type="Italic">f</Emphasis>(<Emphasis Type="Italic">R</Emphasis>,<Emphasis Type="Italic">T</Emphasis>) Gravity with Λ

We study the behavior of cosmological parameters, massive and massless scalar fields (normal or phantom) with a scalar potential in f(R, T) theory of gravity for a flat Friedmann-Robertson-Walker (FRW) universe. To get exact solutions to the modified field equations, we use the f(R, T) = R + 2f(T) model by Harko et al. (T. Harko et al., Phys. Rev. D 84, 024020 (2011)), where R is the Ricci scalar and T is the trace of the energy momentum tensor. Our cosmological parameter solutions agree with the recent observational data. Finally, we discuss our results with various graphics.     

Emergence of anti-<Emphasis Type="Italic">F</Emphasis>(<Emphasis Type="Italic">R</Emphasis>) gravity in type-IV bouncing cosmology as due to <Emphasis Type="Italic">M</Emphasis><Subscript>0</Subscript>-brane

Recently, some authors considered the origin of a type-IV singular bounce in modified gravity and obtained the explicit form of F(R) which can produce this type of cosmology. In this paper, we show that during the contracting branch of type-IV bouncing cosmology, the sign of gravity changes, and antigravity emerges. In our model, M0 branes get together and shape a universe, an anti-universe, and a wormhole which connects them. As time passes, this wormhole is dissolved in the universes, F(R) gravity emerges, and the universe expands. When the brane universes become close to each other, the squared energy of their system becomes negative, and some tachyonic states are produced. To remove these states, universes are assumed to be compact, the sign of compacted gravity changes, and anti-F(R) gravity arises, which causes getting away of the universes from each other. In this theory, a Type-IV singularity occurs at t = t _s , which is the time of producing tachyons between expansion and contraction branches.     

Stabilization of distributed control systems with delay

In this paper the asymptotic stabilization of linear distributed parameter control systems with delay is considered. Specifically, we are concerned with the class of control systems described by the equation x^′(t)=Ax(t)+L(x_t)+Bu(t), where A is the infinitesimal generator of a strongly continuous semigroup on a Banach space X. Assuming appropriate conditions, we will show that the usual spectral controllability assumption implies the feedback stabilization of the system. Applications to systems described by partial differential equations with delay are given.     

On the non-existence of stationary diffusions which satisfy the Bene condition

The one-dimensional diffusion x_t satisfying dx_t = f(x_t)dt + dw_t, where w_t is a standard Brownian motion and f(x) satisfies the Bene condition f′(x) + f²(x) = ax² + bx + c for all real x, is considered. It is shown that this diffusion does not admit a stationary probability measure except for the linear case f(x) = αx + β, α < 0.     

Additive Spanners and Distance and Routing Labeling Schemes for Hyperbolic Graphs

δ-Hyperbolic metric spaces have been defined by M. Gromov in 1987 via a simple 4-point condition: for any four points u,v,w,x, the two larger of the distance sums d(u,v)+d(w,x),d(u,w)+d(v,x),d(u,x)+d(v,w) differ by at most?2δ. They play an important role in geometric group theory, geometry of negatively curved spaces, and have recently become of interest in several domains of computer science, including algorithms and networking. In this paper, we study unweighted δ-hyperbolic graphs. Using the Layering Partition technique, we show that every n-vertex δ-hyperbolic graph with δ≥1/2 has an additive O(δlog?n)-spanner with at most O(δn) edges and provide a simpler, in our opinion, and faster construction of distance approximating trees of δ-hyperbolic graphs with an additive error O(δlog?n). The construction of our tree takes only linear time in the size of the input graph. As a consequence, we show that the family of n-vertex δ-hyperbolic graphs with δ≥1/2 admits a routing labeling scheme with O(δlog?² n) bit labels, O(δlog?n) additive stretch and O(log?₂(4δ)) time routing protocol, and a distance labeling scheme with O(log?² n) bit labels, O(δlog?n) additive error and constant time distance decoder.     

Oscillation theorems for second-order nonlinear neutral differential equations

The purpose of this paper is to study the oscillation of the second-order neutral differential equations of the form

(a(t)[z^′(t)]^γ)^′+q(t)x^β(σ(t))=0,     

Quantizing gravity using physical states of a superstring

A symmetric zero-mass tensor of rank two is constructed using the superstring modes of excitation, which satisfies the physical state constraints of a superstring. These states have a one to one correspondence with the quantized field operators and are shown to be the absorption and emission quanta of the Minkowski space Lorentz tensor, using the quantum field theory method of quantization. The principle of equivalence makes the tensor identical to the metric tensor at any arbitrary space-time point. The propagator for the quantized field is deduced. The gravitational interaction is switched on by going over from ordinary derivatives to co-derivatives. The Riemann-Christoffel affine connections are calculated, and the weak field Ricci tensor R _μν ⁰ is shown to vanish. The interaction part R _μν ^int is found, and the exact R _μν of the theory of gravity is expressed in terms of the quantized metric. The quantum-mechanical self-energy of the gravitational field in vacuum is shown to vanish. By the use of a projection operator, it is shown that gravitons are quanta of the general relativity field which gives the Einstein equation G _μν = 0. It is suggested that quantum gravity may be renormalizable by the use of the massless ground state of this superstring theory for general relativity, and a tachyonic vacuum creates and annihilates quanta of quantized gravitational field.