Numerical simulation for two-dimensional Riesz space fractional diffusion equations with a nonlinear reaction term

Fractional differential equations have attracted considerable interest because of their ability to model anomalous transport phenomena. Space fractional diffusion equations with a nonlinear reaction term have been presented and used to model many problems of practical interest. In this paper, a two-dimensional Riesz space fractional diffusion equation with a nonlinear reaction term (2D-RSFDE-NRT) is considered. A novel alternating direction implicit method for the 2D-RSFDE-NRT with homogeneous Dirichlet boundary conditions is proposed. The stability and convergence of the alternating direction implicit method are discussed. These numerical techniques are used for simulating a two-dimensional Riesz space fractional Fitzhugh-Nagumo model. Finally, a numerical example of a two-dimensional Riesz space fractional diffusion equation with an exact solution is given. The numerical results demonstrate the effectiveness of the methods. These methods and techniques can be extended in a straightforward method to three spatial dimensions, which will be the topic of our future research.     

Weak-field approximation of effective gravitational theory with local Galilean invariance

We examine the weak-field approximation of locally Galilean invariant gravitational theories with general covariance in a (4+1)-dimensional Galilean framework. The additional degrees of freedom allow us to obtain Poisson, diffusion, and Schrödinger equations for the fluctuation field. An advantage of this approach over the usual (3+1)-dimensional General Relativity is that it allows us to choose an ansatz for the fluctuation field that can accommodate the field equations of the Lagrangian approach to MOdified Newtonian Dynamics (MOND) known as AQUAdratic Lagrangian (AQUAL). We investigate a wave solution for the Schrödinger equations.     

A Generally Covariant Wave Equation for Grand Unified Field Theory

A generally covariant wave equation is derived geometrically for grand unified field theory. The equation states most generally that the covariant d'Alembertian acting on the vielbein vanishes for the four fields which are thought to exist in nature: gravitation, electromagnetism, weak field and strong field. The various known field equations are derived from the wave equation when the vielbein is the eigenfunction. When the wave equation is applied to gravitation the wave equation is the eigenequation of wave mechanics corresponding to Einstein's field equation in classical mechanics, the vielbein eigenfunction playing the role of the quantized gravitational field. The three Newton laws, Newton's law of universal gravitation, and the Poisson equation are recovered in the classical and nonrelativistic, weak-field limits of the quantized gravitational field. The single particle wave-equation and Klein-Gordon equations are recovered in the relativistic, weak-field limit of the wave equation when scalar components are considered of the vielbein eigenfunction of the quantized gravitational field. The Schrödinger equation is recovered in the non-relativistec, weak-field limit of the Klein-Gordon equation). The Dirac equation is recovered in this weak-field limit of the quantized gravitational field (the nonrelativistic limit of the relativistic, quantezed gravitational field when the vielbein plays the role of the spinor. The wave and field equations of O(3) electrodynamics are recovered when the vielbein becomes the relativistic dreibein (triad) eigenfunction whose three orthonormal space indices become identified with the three complex circular indices (1), (2), (3), and whose four spacetime indices are the indices of non-Euclidean spacetime (the base manifold). This dreibein is the potential dreibein of the O(3) electromagnetic field (an electromagnetic potential four-vector for each index (1), (2), (3)). The wave equation of the parity violating weak field is recovered when the orthonormal space indices of the relativistic dreibein eigenfunction are identified with the indices of the three massive weak field bosons. The wave equation of the strong field is recovered when the orthonormal space indices of the relativistic vielbein eigenfunction become the eight indices defined by the group generators of the SU (3) group.     

Crank–Nicolson method for the fractional diffusion equation with the Riesz fractional derivative

We examine a numerical method to approximate to a fractional diffusion equation with the Riesz fractional derivative in a finite domain, which has second order accuracy in time and space level. In order to approximate the Riesz fractional derivative, we use the “fractional centered derivative” approach. We determine the error of the Riesz fractional derivative to the fractional centered difference. We apply the Crank–Nicolson method to a fractional diffusion equation which has the Riesz fractional derivative, and obtain that the method is unconditionally stable and convergent. Numerical results are given to demonstrate the accuracy of the Crank–Nicolson method for the fractional diffusion equation with using fractional centered difference approach.     

Kinetic Description of Particle Interaction with a Gravitational Wave

The interaction of charged particles, moving in a uniform magnetic field, with a plane polarized gravitational wave is considered using the Fokker-Planck-Kolmogorov (FPK) approach. By using a stochasticity criterion, we determine the exact locations in phase space, where resonance overlapping occurs. We investigate the diffusion of orbits around each primary resonance of order m by deriving general analytical expressions for an effective diffusion coefficient. A solution of the corresponding diffusion equation (Fokker-Planck equation) for the static case is found. Numerical integration of the full equations of motion and subsequent calculation of the diffusion coefficient verifies the analytical results.     

Linear gravitational waves and electrodynamic formalism in cosmology

The propagation of linear gravitational waves is studied in open and multiply connected Robertson-Walker cosmologies. In order for the group velocity of the gravitational wave packets to coincide with the speed of light, the linear wave equation must be conformally coupled. This opens the possibility of using the electromagnetic formalism. The gravitational analogue to the electromagnetic field tensor is introduced, and a tensorial counterpart to Maxwell's equations on the spacelike 3-slices is derived. The energy-momentum tensor for linear gravitational waves is constructed without averaging procedures, a strictly positive energy density is obtained, and it is shown that the overall energy of a gravitational pulse scales with the inverse of the expansion factor.     

General relativity in Newtonian form

It is shown how the use of coordinates where time is measured with clocks moving radially in a spherically symmetric gravitational field leads to general relativistic dynamical expressions that are exactly identical to corresponding expressions in Newtonian theory. The general formalism is developed for the case where the stress-energy tensor is that of a perfect fluid. Expressions like the Newtonian inverse square gravitational law, the Newtonian equation of continuity for fluid flow, Newtonian conservation of energy, etc., follow quite naturally from the fully-fledged exact general relativistic equations. Specific examples involving cosmology and gravitational collapse are given.     

Conformal Fluctuations of the Interior Schwarzschild Solution

The basic formalism for conformal fluctuations of the gravitational field is presented. After developing a master propagator for the interior Schwarzschild solution, the time development of the gravitational wave function is considered. The effect of the two classical singularities (resp. pseudo-singularities) of the Schwarzschild solution on the quantum wave function for the gravitational field is studied using a wave function initially localized on the classical solution. While the true singularity at r = 0 imparts consequences on the wave function that cannot be ignored, the pseudo-singularity at the event horizon does not seem to cause any divergences on the interior fluctuations of the Schwarzschild solution.     

Classical Gravitational Interactions and Gravitational Lorentz Force

In quantum gauge theory of gravity, the gravitational field is represented by gravitational gauge field.The field strength of gravitational gauge field has both gravitoelectric component and gravitomagnetic component. In classical level, gauge theory of gravity gives classical Newtonian gravitational interactions in a relativistic form. Besides,it gives gravitational Lorentz force, which is the gravitational force on a moving object in gravitomagnetic field The direction of gravitational Lorentz force is not the same as that of classical gravitational Newtonian force. Effects of gravitational Lorentz force should be detectable, and these effects can be used to discriminate gravitomagnetic field from ordinary electromagnetic magnetic field.     

Non-relativistic Limit of Dirac Equations in Gravitational Field and Quantum Effects of Gravity

Based on unified theory of electromagnetic interactions and gravitational interactions, the non-relativistic limit of the equation of motion of a charged Dirac particle in gravitational field is studied. From the Schrodinger equation obtained from this non-relativistic limit, we can see that the classical Newtonian gravitational potential appears as a part of the potential in the Schrodinger equation, which can explain the gravitational phase effects found in COW experiments.And because of this Newtonian gravitational potential, a quantum particle in the earth's gravitational field may form a gravitationally bound quantized state, which has already been detected in experiments. Three different kinds of phase effects related to gravitational interactions are studied in this paper, and these phase effects should be observable in some astrophysical processes. Besides, there exists direct coupling between gravitomagnetic field and quantum spin, and radiation caused by this coupling can be used to directly determine the gravitomagnetic field on the surface of a star.     

Classical Gravitational Interactions and Gravitational Lorentz Force

In quantum gauge theory of gravity, the gravitational field is represented by gravitational gauge field. The field strength of gravitational gauge field has both gravitoelectric component and gravitomagnetic component. In classical level, gauge theory of gravity gives classical Newtonian gravitational interactions in a relativistic form. Besides, it gives gravitational Lorentz force, which is the gravitational force on a moving object in gravitomagnetic field. The direction of gravitational Lorentz force is not the same as that of classical gravitational Newtonian force. Effects of gravitational Lorentz force should be detectable, and these effects can be used to discriminate gravitomagnetic field from ordinary electromagnetic magnetic field.     

Ideal fluid with short-range scalar interactions in the field of a plane gravitational wave

An exact solution of the self-consistent equations of relativistic hydrodynamics and the scalar field equation is obtained. The solution describes motion of a fluid with short-range scalar interactions in the field of a plane gravitational wave.     

A second-order accurate numerical method for the two-dimensional fractional diffusion equation

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are used in modeling practical superdiffusive problems in fluid flow, finance and others. In this paper, we present an accurate and efficient numerical method to solve a fractional superdiffusive differential equation. This numerical method combines the alternating directions implicit (ADI) approach with a Crank–Nicolson discretization and a Richardson extrapolation to obtain an unconditionally stable second-order accurate finite difference method. The stability and the consistency of the method are established. Numerical solutions for an example super-diffusion equation with a known analytic solution are obtained and the behavior of the errors are analyzed to demonstrate the order of convergence of the method.     

The analog of electric and magnetic fields in stationary gravitational systems

Newtonian and Machian aspects of the stationary gravitational field are brought into formal analogy with a stationary electromagnetic field. The electromagnetic vector potential equals (up to a factor) the timelike Killing vector field. The current density is given by the contraction of the Killing vector with the Ricci tensor. A coordinate-dependent split in electric and magnetic field vectors is given, and some results of classical electrodynamics are used to illustrate the analogy. In the linearized theory, the usual Maxwell equations are obtained. The analogy also holds from the point of view of particle motion. The geodesic equation is brought into a special form that exhibits an analog to the Lorentz force. Two examples (which have played an important role in the theoretical discovery of Machian effects) are considered.     

Analytical Solution of a Wave Equation in Cosmology

In this paper we derive a wave equation based on a cosmological model of the universe expansion in the four-dimensional space-velocity. An analytical solution of this equation is obtained using Nikiforov-Uvarov mathematical method. In this procedure we use the solutions of the Einstein gravitational field equations for the perfect fluid energy-momentum tensor.     

Time-dependent solutions with localized energy of the Einstein system of equations and a nonlinear scalar field

A soliton-like time-dependent solution in the form of a running wave is derived of a self-consistent system of the gravitational field equations of Einstein and Born-Infeld type of equations of a nonlinear scalar field in a conformally flat metric. This solution is localized in space and possesses a localized energy. It is shown that both the gravitational field and the nonlinearity of the scalar field are essential to the presence of such a localized solution. In recent years various classical particle models have been widely discussed which are static or time-independent solutions of nonlinear equations with localization in space and which possess a finite field energy. In particular, soliton solutions [1], solutions in the form of eddies [2], and so on have been derived and investigated. All these solutions were treated in a flat space-time. It is of interest to derive the analogous particle-like solutions with the gravitational field taken into account; in particular it is of interest to investigate the roles of the gravitational field in connection with the formation of localized objects. These problems have been discussed in [3] in the static case. We will present below a soliton-like time-dependent solution in the form of a solitary running wave as an example of the inter-action of a Born-Infeld type of nonlinear scalar field and an Einstein gravitational field in a conformally flat metric.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 5, pp. 12–17, May, 1979.     

A model of unified quantum chromodynamics and Yang-Mills gravity

Based on a generalized Yang-Mills framework, gravitational and strong interactions can be unified in analogy with the unification in the electroweak theory. By gauging T (4) × [SU (3)] color in flat space-time, we have a unified model of chromo-gravity with a new tensor gauge field, which couples universally to all gluons, quarks and anti-quarks. The space-time translational gauge symmetry assures that all wave equations of quarks and gluons reduce to a Hamilton-Jacobi equation with the same ’effective Riemann metric tensors’ in the geometric-optics (or classical) limit. The emergence of ef f ective metric tensors in the classical limit is essential for the unified model to agree with experiments. The unified model suggests that all gravitational, strong and electroweak interactions appear to be dictated by gauge symmetries in the generalized Yang-Mills framework.     

Lagrangian and Hamiltonian formalisms for relativistic dynamics of a charged particle with dipole moment

The Lagrangian and Hamiltonian formulations for the relativistic classical dynamics of a charged particle with dipole moment in the presence of an electromagnetic field are given. The differential conservation laws for the energy-momentum and angular momentum tensors of a field and particle are discussed. The Poisson brackets for basic dynamic variables, which form a closed algebra, are found. These Poisson brackets enable us to perform the canonical quantization of the Hamiltonian equations that leads to the Dirac wave equation in the case of spin 1/2. It is also shown that the classical limit of the squared Dirac equation results in equations of motion for a charged particle with dipole moment obtained from the Lagrangian formulation. The inclusion of gravitational field and non-Abelian gauge fields into the proposed formalism is discussed.Received: 4 June 2005, Published online: 27 July 2005     

Static observers in curved spaces and non-inertial frames in Minkowski spacetime

Static observers in curved spacetimes may interpret their proper acceleration as the opposite of a local gravitational field (in the Newtonian sense). Based on this interpretation and motivated by the equivalence principle, we are led to investigate congruences of timelike curves in Minkowski spacetime whose acceleration field coincides with the acceleration field of static observers of curved spaces. The congruences give rise to non-inertial frames that are examined. Specifically, we find, based on the locality principle, the embedding of simultaneity hypersurfaces adapted to the non-inertial frame in an explicit form for arbitrary acceleration fields. We also determine, from the Einstein equations, a covariant field equation that regulates the behavior of the proper acceleration of static observers in curved spacetimes. It corresponds to an exact relativistic version of the Newtonian gravitational field equation. In the specific case in which the level surfaces of the norm of the acceleration field of the static observers are maximally symmetric two-dimensional spaces, the energy?Cmomentum tensor of the source is analyzed.