Regularity of solutions of the model Venttsel’ problem for quasilinear parabolic systems with nonsmooth in time principal matrices

The Venttsel’ problem in the model statement for quasilinear parabolic systems of equations with nondiagonal principal matrices is considered. It is only assumed that the principal matrices and the boundary condition are bounded with respect to the time variable. The partial smoothness of the weak solutions (Hölder continuity on a set of full measure up to the surface on which the Venttsel’ condition is defined) is proved. The proof uses the A(t)-caloric approximation method, which was also used in [1] to investigate the regularity of the solution to the corresponding linear problem.     

On the problem of determining the structure of a layered medium and the shape of an impulse source

For the equation of wave propagation in the half-space ? ₊ ² + = {(x, y) ∈ ?² | y > 0} we consider the problem of determining the speed of wave propagation that depends only on the variable y and the shape of a point impulse source on the boundary of the half-space. We show that, under some assumptions on the shape of the source and the structure of the medium, both unknown functions of one variable are uniquely determined by the displacements of boundary points of the medium. We estimate stability of a solution to the problem.     

Asymptotic Behavior of Massless Dirac Waves in Schwarzschild Geometry

In this paper, we show that massless Dirac waves in the Schwarzschild geometry decay to zero at a rate t ^?2λ , where λ = 1, 2, . . . is the angular momentum. Our technique is to use Chandrasekhar’s separation of variables whereby the Dirac equations split into two sets of wave equations. For the first set, we show that the wave decays as t ^?2λ . For the second set, in general, the solutions tend to some explicit profile at the rate t ^?2λ . The decay rate of solutions of Dirac equations is achieved by showing that the coefficient of the explicit profile is exactly zero. The key ingredients in the proof of the decay rate of solutions for the first set of wave equations are an energy estimate used to show the absence of bound states and zero energy resonance and the analysis of the spectral representation of the solutions. The proof of asymptotic behavior for the solutions of the second set of wave equations relies on careful analysis of the Green’s functions for time independent Schrödinger equations associated with these wave equations.     

On Dirichlet and Neumann problems for harmonic functions

The aim of the paper is to examine some aspects of the boundary value problems for harmonic functions in half-spaces related to approximation theory. M. V. Keldyshmentioned curious fact on richness in some sense of the solutions of Dirichlet problem in upper half-plane for a fixed continuous boundary data on the real axis. This can be considered as a model version for the Dirichlet problem with continuous boundary data, defined except a single boundary point, with no restrictions imposed on solutions near that point.Some extensions and multi-dimensional versions of Keldysh’s richness are obtained and related questions on existence, representation and richness of solutions for the Dirichlet and Neumann problems discussed.     

Polynomial decay and blow up of solutions for variable coefficients viscoelastic wave equation with acoustic boundary conditions

A variable coefficient viscoelastic wave equation with acoustic boundary conditions and nonlinear source term is considered. Under suitable conditions on the initial data and the relaxation function g, we show the polynomial decay of the energy solution and the blow up of solutions by energy methods. The estimates for the lifespan of solutions are also given.     

Boundary value problems for two types of degenerate elliptic systems in <Emphasis Type="Italic">R</Emphasis><Superscript>4</Superscript>

Firstly, the Riemann boundary value problem for a kind of degenerate elliptic system of the first order equations in R ⁴ is proposed. Then, with the help of the one-to-one correspondence between the theory of Clifford valued generalized regular functions and that of the degenerate elliptic system’s solution, the boundary value problem as stated above is transformed into a boundary value problem related to the generalized regular functions in Clifford analysis. Moreover, the solution of the Riemann boundary value problem for the degenerate elliptic system is explicitly described by using a kind of singular integral operator. Finally, the conditions for the existence of solutions of the oblique derivative problem for another kind of degenerate elliptic system of the first order equations in R ⁴ are derived.     

Existence of solutions of boundary value problems for singular fractional <Emphasis Type="Italic">q</Emphasis>-difference equations

In this paper, we study the existence and uniqueness of solutions for a class of singular three-point boundary value problems of fractional q-difference equations invovling fractional q-derivative of Riemann–Liouville type. Based on the generalization of Banach contraction principle, we obtain a sufficient condition for existence and uniqueness of solutions of the problem. By applying the Krasnoselskii’s fixed point theorem, we establish a sufficient condition for the existence of at least one solution of the problem. As applications, two examples are presented to illustrate our main results.     

On the existence of global solutions of a system of semilinear parabolic equations with nonlinear nonlocal boundary conditions

We establish conditions for the existence and nonexistence of global solutions of an initial–boundary value problem for a system of semilinear parabolic equations with nonlinear nonlocal boundary conditions. The results depend on the behavior of variable coefficients as t→∞.     

Two-Parameter Asymptotics in a Bisingular Cauchy Problem for a Parabolic Equation

The Cauchy problem for a quasilinear parabolic equation with a small parameter ε at the highest derivative is considered. The initial function, which has the form of a smoothed step, depends on a “stretched” variable x/ρ, where ρ is another small parameter. This problem statement is of interest for applications as a model of propagation of nonlinear waves in physical systems in the presence of small dissipation. In the case corresponding to a compression wave, asymptotic solutions of the problem are constructed in the parameters ε and ρ independently tending to zero. It is assumed that ε/ρ → 0. Far from the line of discontinuity of the limit solution, asymptotic solutions are constructed in the form of series in powers of ε and ρ. In a small domain of linear approximation, an asymptotic solution is constructed in the form of a series in powers of the ratio ρ/ε. The coefficients of the inner expansion are determined from a recursive chain of initial value problems. The asymptotics of these coefficients at infinity is studied. The time of reconstruction of the scale of the internal space variable is determined.     

The class of solenoidal planar-helical vector fields

The class of solenoidal vector fields whose lines lie in planes parallel to R ² is constructed by the method of mappings. This class exhausts the set of all smooth planarhelical solutions of Gromeka’s problem in some domain D ? R ³. In the case of domains D with cylindrical boundaries whose generators are orthogonal to R ², it is shown that the choice of a specific solution from the constructed class is reduced to the Dirichlet problem with respect to two functions that are harmonic conjugates in D ² = D∩R ²; i.e., Gromeka’s nonlinear problem is reduced to linear boundary value problems. As an example, a specific solution of the problem for an axisymmetric layer is presented. The solution is based on solving Dirichlet problems in the form of series uniformly convergent in \(\bar D^2\) in terms of wavelet systems that form bases of various spaces of functions harmonic in D ².     

Kinetic Monte Carlo models for the study of chemical reactions in the Earth’s upper atmosphere

A stochastic approach to study the non-equilibrium chemistry in the Earth’s upper atmosphere is presented, which has been developed over a number of years. Kinetic Monte Carlo models based on this approach are an effective tool for investigating the role of suprathermal particles both in local variations of the atmospheric chemical composition and in the formation of the hot planetary corona.     

Fractal and multifractal analysis of the rise of oxygen in Earth’s early atmosphere

The rise of oxygen in Earth’s atmosphere that occurred 2.4–2.2 billion years ago is known as the Earth’s Great Oxidation, and its impact on the development of life on Earth has been profound. Thereafter, the increase in Earth’s oxygen level persisted, though at a more gradual pace. The proposed underlying mathematical models for these processes are based on physical parameters whose values are currently not well-established owing to uncertainties in geological and biological data. In this paper, a previously developed model of Earth’s atmosphere is modified by adding different strengths of noise to account for the parameters’ uncertainties. The effects of the noise on the time variations of oxygen, carbon and methane for the early Earth are investigated by using fractal and multifractal analysis. We show that the time variations following the Great Oxidation cannot properly be described by a single fractal dimension because they exhibit multifractal characteristics. The obtained results demonstrate that the time series as obtained exhibit multifractality caused by long-range time correlations.     

Quasi-periodic water waves

We present the result and the ideas of the recent paper (Berti and Montalto, Quasi-periodic standing wave solutions of gravity-capillary water waves, http://arxiv.org/abs/1602.02411, 2016) concerning the existence of Cantor families of small-amplitude time quasi-periodic standing wave solutions (i.e. periodic and even in the space variable x) of a 2-dimensional ocean, with infinite depth, in irrotational regime, under the action of gravity and surface tension at the free boundary. These quasi-periodic solutions are linearly stable.     

Feynman and quasi-Feynman formulas for evolution equations

New methods for obtaining representations of solutions of the Cauchy problem for linear evolution equations, i.e., equations of the form u _t '(t, x) = Lu(t, x), where the operator L is linear and depends only on the spatial variable x and does not depend on time t, are proposed. A solution of the Cauchy problem, that is, the exponential of the operator tL, is found on the basis of constructions proposed by the author combined with Chernoff’s theorem on strongly continuous operator semigroups.     

To the theory of the Dirichlet and Neumann problems for strongly elliptic systems in Lipschitz domains

For strongly elliptic Systems with Douglis-Nirenberg structure, we investigate the regularity of variational solutions to the Dirichlet and Neumann problems in a bounded Lipschitz domain. The solutions of the problems with homogeneous boundary conditions are originally defined in the simplest L ₂-Sobolev spaces H ^σ . The regularity results are obtained in the potential spaces H _p ^σ and Besov spaces B _p ^σ . In the case of second-order Systems, the author’s results obtained a year ago are strengthened. The Dirichlet problem with nonhomogeneous boundary conditions is considered with the use of Whitney arrays.     

Local bifurcations of plane running waves for the generalized cubic Schrödinger equation

For the generalized cubic Schrödinger equation, we consider a periodic boundary value problem in the case of n independent space variables. For this boundary value problem, there exists a countable set of plane running waves periodic with respect to the time variable. We analyze their stability and local bifurcations under the change of stability. We show that invariant tori of dimension 2, ..., n + 1 can bifurcate from each of them. We obtain asymptotic formulas for the solutions on invariant tori and stability conditions for bifurcating tori as well as parameter ranges in which, starting from n = 3, a subcritical (stiff) bifurcation of invariant tori is possible.     

The problem of recovering the permittivity coefficient from the modulus of the scattered electromagnetic field

Under consideration is the stationary system of equations of electrodynamics relating to a nonmagnetic nonconducting medium. We study the problem of recovering the permittivity coefficient ε from given vectors of electric or magnetic intensities of the electromagnetic field. It is assumed that the field is generated by a point impulsive dipole located at some point y. It is also assumed that the permittivity differs from a given constant ε₀ only inside some compact domain Ω ? R³ with smooth boundary S. To recover ε inside Ω, we use the information on a solution to the corresponding direct problem for the system of equations of electrodynamics on the whole boundary of Ω for all frequencies from some fixed frequency ω ₀ on and for all y ∈ S. The asymptotics of a solution to the direct problem for large frequencies is studied and it is demonstrated that this information allows us to reduce the initial problem to the well-known inverse kinematic problem of recovering the refraction index inside Ω with given travel times of electromagnetic waves between two arbitrary points on the boundary of Ω. This allows us to state uniqueness theorem for solutions to the problem in question and opens up a way of its constructive solution.     

Solution of the inverse magnetotelluric sounding problem by means of a method using the synthesized magnetic field

A method of solution of the inverse magnetotelluric sounding problem is considered. The method uses the known frequency-dependent magnetic field on the Earth’s surface. The magnetic field on the Earth’s surface is found with the help of a special technique from the known impedance on the Earth’s surface. The method is applied to solve typical inverse problems and provides for efficient determination of final solutions.     

Convergence of a family of solutions to a Fujita-type equation in domains with cavities

The Dirichlet problem for a Fujita-type equation, i.e., a second-order quasilinear uniformly elliptic equation is considered in domains Ω_ε with spherical or cylindrical cavities of characteristic size ε. The form of the function in the condition on the cavities’ boundaries depends on ε. For ε tending to zero and the number of cavities increasing simultaneously, sufficient conditions are established for the convergence of the family of solutions {u _ε(x)} of this problem to the solution u(х) of a similar problem in the domain Ω with no cavities with the same boundary conditions imposed on the common part of the boundaries ?Ω and ?Ω_ε. Convergence rate estimates are given.     

Optimal homotopy analysis of a chaotic HIV-1 model incorporating AIDS-related cancer cells

The studies of nonlinear models in epidemiology have generated a deep interest in gaining insight into the mechanisms that underlie AIDS-related cancers, providing us with a better understanding of cancer immunity and viral oncogenesis. In this article, we analyze an HIV-1 model incorporating the relations between three dynamical variables: cancer cells, healthy CD4 ⁺ T lymphocytes, and infected CD4 ⁺ T lymphocytes. Recent theoretical investigations indicate that these cells interactions lead to different dynamical outcomes, for instance to periodic or chaotic behavior. Firstly, we analytically prove the boundedness of the trajectories in the system’s attractor. The complexity of the coupling between the dynamical variables is quantified using observability indices. Our calculations reveal that the highest observable variable is the population of cancer cells, thus indicating that these cells could be monitored in future experiments in order to obtain time series for attractor’s reconstruction. We identify different dynamical behaviors of the system varying two biologically meaningful parameters: r ₁, representing the uncontrolled proliferation rate of cancer cells, and k ₁, denoting the immune system’s killing rate of cancer cells. The maximum Lyapunov exponent is computed to identify the chaotic regimes. Considering very recent developments in the literature related to the homotopy analysis method (HAM), we calculate the explicit series solutions of the cancer model and focus our analysis on the dynamical variable with the highest observability index. An optimal homotopy analysis approach is used to improve the computational efficiency of HAM by means of appropriate values for the convergence control parameter, which greatly accelerate the convergence of the series solution. The approximated analytical solutions are used to compute density plots, which allow us to discuss additional dynamical features of the model.