Estimates of oscillatory integrals with stationary phase and singular amplitude: Applications to propagation features for dispersive equations

In this paper, we study time‐asymptotic propagation phenomena for a class of dispersive equations on the line by exploiting precise estimates of oscillatory integrals. We propose first an extension of the van der Corput Lemma to the case of phases which may have a stationary point of real order and amplitudes allowed to have an integrable singular point. The resulting estimates provide optimal decay rates which show explicitly the influence of these two particular points. Then we apply these abstract results to solution formulas of a class of dispersive equations on the line defined by Fourier multipliers. Under the hypothesis that the Fourier transform of the initial data has a compact support or an integrable singular point, we derive uniform estimates of the solutions in space‐time cones, describing their motions when the time tends to infinity. The method permits also to show that symbols having a restricted growth at infinity may influence the dispersion of the solutions: we prove the existence of a cone, depending only on the symbol, in which the solution is time‐asymptotically localized. This corresponds to an asymptotic version of the notion of causality for initial data without compact support.     

Energy decay estimates for the dissipative wave equation with space–time dependent potential

We study the asymptotic behavior of solutions of dissipative wave equations with space–time‐dependent potential. When the potential is only time‐dependent, Fourier analysis is a useful tool to derive sharp decay estimates for solutions. When the potential is only space‐dependent, a powerful technique has been developed by Todorova and Yordanov to capture the exact decay of solutions. The presence of a space–time‐dependent potential, as in our case, requires modifications of this technique. We find the energy decay and decay of the L² norm of solutions in the case of space–time‐dependent potential. Copyright © 2010 John Wiley & Sons, Ltd.     

Cylindrical estimates for mean curvature flow of hypersurfaces in CROSSes

We consider the mean curvature flow of a closed hypersurface in the complex or quaternionic projective space. Under a suitable pinching assumption on the initial data, we prove apriori estimates on the principal curvatures which imply that the asymptotic profile near a singularity is either strictly convex or cylindrical. This result generalizes to a large class of symmetric ambient spaces the estimates obtained in the previous works on the mean curvature flow of hypersurfaces in Euclidean space and in the sphere.     

Tunneling of the hard‐core model on finite triangular lattices

We consider the hard‐core model on finite triangular lattices with Metropolis dynamics. Under suitable conditions on the triangular lattice sizes, this interacting particle system has 3 maximum‐occupancy configurations and we investigate its high‐fugacity behavior by studying tunneling times, that is, the first hitting times between these maximum‐occupancy configurations, and the mixing time. The proof method relies on the analysis of the corresponding state space using geometrical and combinatorial properties of the hard‐core configurations on finite triangular lattices, in combination with known results for first hitting times of Metropolis Markov chains in the equivalent zero‐temperature limit. In particular, we show how the order of magnitude of the expected tunneling times depends on the triangular lattice sizes in the low‐temperature regime and prove the asymptotic exponentiality of the rescaled tunneling time leveraging the intrinsic symmetry of the state space.     

The nonlinear Schrödinger equation with a potential

We consider the cubic nonlinear Schrödinger equation with a potential in one space dimension. Under the assumptions that the potential is generic, sufficiently localized, with no bound states, we obtain the long-time asymptotic behavior of small solutions. In particular, we prove that, as time goes to infinity, solutions exhibit nonlinear phase corrections that depend on the scattering matrix associated to the potential. The proof of our result is based on the use of the distorted Fourier transform – the so-called Weyl–Kodaira–Titchmarsh theory – a precise understanding of the “nonlinear spectral measure” associated to the equation, and nonlinear stationary phase arguments and multilinear estimates in this distorted setting.     

Convergence proof of the velocity field for a stokes flow immersed boundary method

The immersed boundary (IB) method is a computational framework for problems involving the interaction of a fluid and immersed elastic structures. It is characterized by the use of a uniform Cartesian mesh for the fluid, a Lagrangian curvilinear mesh on the elastic material, and discrete delta functions for communication between the two grids. We consider a simple IB problem in a two‐dimensional periodic fluid domain with a one‐dimensional force generator. We obtain error estimates for the velocity field of the IB solution for the stationary Stokes problem. We use this result to prove convergence of a simple small‐amplitude dynamic problem. We test our error estimates against computational experiments. © 2007 Wiley Periodicals, Inc.     

Asymptotic Properties of Quasi-Maximum Likelihood Estimators for ARMA Models with Time-Dependent Coefficients

For about thirty years, time series models with time-dependent coefficients have sometimes been considered as an alternative to models with constant coefficients or non-linear models. Analysis based on models with time-dependent models has long suffered from the absence of an asymptotic theory except in very special cases. The purpose of this paper is to provide such a theory without using a locally stationary spectral representation and time rescaling. We consider autoregressive-moving average (ARMA) models with time-dependent coefficients and a heteroscedastic innovation process. The coefficients and the innovation variance are deterministic functions of time which depend on a finite number of parameters. These parameters are estimated by maximising the Gaussian likelihood function. Deriving conditions for consistency and asymptotic normality and obtaining the asymptotic covariance matrix are done using some assumptions on the functions of time in order to attenuate non-stationarity, mild assumptions for the distribution of the innovations, and also a kind of mixing condition. Theorems from the theory of martingales and mixtingales are used. Some simulation results are given and both theoretical and practical examples are treated. Received 2004; Final version 23 December 2004     

A multilevel finite element method in space‐time for the Navier‐Stokes problem

A multilevel finite element method in space‐time for the two‐dimensional nonstationary Navier‐Stokes problem is considered. The method is a multi‐scale method in which the fully nonlinear Navier‐Stokes problem is only solved on a single coarsest space‐time mesh; subsequent approximations are generated on a succession of refined space‐time meshes by solving a linearized Navier‐Stokes problem about the solution on the previous level. The a priori estimates and error analysis are also presented for the J‐level finite element method. We demonstrate theoretically that for an appropriate choice of space and time mesh widths: h_j ～ h, k_j ～ k, j = 2, …, J, the J‐level finite element method in space‐time provides the same accuracy as the one‐level method in space‐time in which the fully nonlinear Navier‐Stokes problem is solved on a final finest space‐time mesh. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Two‐grid methods for time‐harmonic Maxwell equations

In this paper, we develop several two‐grid methods for the Nédélec edge finite element approximation of the time‐harmonic Maxwell equations. We first present a two‐grid method that uses a coarse space to solve the original problem and then use a fine space to solve a corresponding symmetric positive definite problem. Then, we present two types of iterative two‐grid methods, one is to add the kernel of the curl ‐operator in the fine space to a coarse mesh space to solve the original problem and the other is to use an inner iterative method for dealing with the kernel of the curl ‐operator in the fine space and the coarse space, separately. We provide the error estimates for the first two methods and present numerical experiments to show the efficiency of our methods.Copyright © 2012 John Wiley & Sons, Ltd.     

A space‐time continuous Galerkin method with mesh modification for viscoelastic wave equations

In this article, we consider the space‐time continuous Galerkin (STCG) method for the viscoelastic wave equations. It allows variable temporal step‐sizes, and the changing of the spatial grids in two adjacent time levels. The existence, uniqueness, and stability of the approximate solutions are demonstrated and the error estimates with global and local spatial mesh sizes in norm are derived without any restrictive assumptions on the space‐time meshes. If the meshes in each time level satisfy some reasonable assumptions, then we can get the optimal order error estimates both in time and space. Finally, we give a numerical example on unstructured meshes to confirm the theoretical findings. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1183–1207, 2017     

Stabilized nodal‐based finite element methods for the magnetic induction problem

We consider the time‐dependent magnetic induction model as a step towards the resistive magnetohydrodynamics model in incompressible media. Conforming nodal‐based finite element approximations of the induction model with inf‐sup stable finite elements for the magnetic field and the magnetic pseudo‐pressure are investigated. Based on a residual‐based stabilization technique proposed by Badia and Codina, SIAM J. Numer. Anal. 50 (2012), pp. 398–417, we consider a stabilized nodal‐based finite element method for the numerical solution. Error estimates are given for the semi‐discrete model in space. Finally, we present some examples, for example, for the magnetic flux expulsion problem, Shercliff's test case and singular solutions of the Maxwell problem. Copyright © 2015 John Wiley & Sons, Ltd.     

A posteriori error estimates for discontinuous Galerkin method to the elasticity problem

This work concerns with the discontinuous Galerkin (DG) method for the time‐dependent linear elasticity problem. We derive the a posteriori error bounds for semidiscrete and fully discrete problems, by making use of the stationary elasticity reconstruction technique which allows to estimate the error for time‐dependent problem through the error estimation of the associated stationary elasticity problem. For fully discrete scheme, we make use of the backward‐Euler scheme and an appropriate space‐time reconstruction. The technique here can be applicable for a variety of DG methods as well.     

Asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations

This article discusses the analyticity and the long‐time asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations in . By a Laplace transform argument, we prove that the decay rate of the solution as t→∞ is dominated by the order of the time‐fractional derivative. We consider the decay rate also in a bounded domain. Copyright © 2016 John Wiley & Sons, Ltd.     

A Cahn‐Hilliard–type equation with application to tumor growth dynamics

We consider a Cahn‐Hilliard–type equation with degenerate mobility and single‐well potential of Lennard‐Jones type. This equation models the evolution and growth of biological cells such as solid tumors. The degeneracy set of the mobility and the singularity set of the cellular potential do not coincide, and the absence of cells is an unstable equilibrium configuration of the potential. This feature introduces a nontrivial difference with respect to the Cahn‐Hilliard equation analyzed in the literature. We give existence results for different classes of weak solutions. Moreover, we formulate a continuous finite element approximation of the problem, where the positivity of the solution is enforced through a discrete variational inequality. We prove the existence and uniqueness of the discrete solution for any spatial dimension together with the convergence to the weak solution for spatial dimension d=1. We present simulation results in 1 and 2 space dimensions. We also study the dynamics of the spinodal decomposition and the growth and scaling laws of phase ordering dynamics. In this case, we find similar results to the ones obtained in standard phase ordering dynamics and we highlight the fact that the asymptotic behavior of the solution is dominated by the mechanism of growth by bulk diffusion.     

Optimal search efficiency of Barker’s algorithm with an exponential fitness function

Markov Chain Monte Carlo (MCMC) methods may be employed to search for a probability distribution over a bounded space of function arguments to estimate which argument(s) optimize(s) an objective function. This search-based optimization requires sampling the suitability, or fitness, of arguments in the search space. When the objective function or the fitness of arguments vary with time, significant exploration of the search space is required. Search efficiency then becomes a more relevant measure of the usefulness of an MCMC method than traditional measures such as convergence speed to the stationary distribution and asymptotic variance of stationary distribution estimates. Search efficiency refers to how quickly prior information about the search space is traded-off for search effort savings. Optimal search efficiency occurs when the entropy of the probability distribution over the space during search is maximized. Whereas the Metropolis case of the Hastings MCMC algorithm with fixed candidate generation is optimal with respect to asymptotic variance of stationary distribution estimates, this paper proves that Barker’s case is optimal with respect to search efficiency if the fitness of the arguments in the search space is characterized by an exponential function. The latter instance of optimality is beneficial for time-varying optimization that is also model-independent.     

Decay Estimates for Solutions of Linear Elasticity for Anisotropic Media

We analyse the time decay of solutions to the Cauchy problem for the linear hyperbolic system of elasticity for anisotropic media. As an example, we will consider media with hexagonal symmetry. First we derive decay estimates for special initial data using the method of stationary phase in several variables and degenerate phase function based on the Malgrange preparation theorem. Asymptotic expansions are given to prove the sharpness of the weaker time decay found for zinc and beryl than in the isotropic case. A method using Besov spaces leads to ℒ︁^p–ℒ︁^q-estimates.     

WEB‐spline based mesh‐free finite element analysis for the heat equation and the time‐dependent Navier–Stokes equation: A survey

In this article, we give some numerical techniques and error estimates using web‐spline based mesh‐free finite element method for the heat equation and the time‐dependent Navier–Stokes equations on bounded domains. The web‐spline method uses weighted extended B‐splines on a regular grid as basis functions and does not require any grid generation. We demonstrate the method by providing numerical results for the Poisson's and stationary Stokes equation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Advection of methane in the hydrate zone: model,analysis and examples

A two‐phase two‐component model is formulated for the advective–diffusive transport of methane in liquid phase through sediment with the accompanying formation and dissolution of methane hydrate. This free‐boundary problem has a unique generalized solution in L¹; the proof combines analysis of the stationary semilinear elliptic Dirichlet problem with the nonlinear semigroup theory in Banach space for an m‐accretive multi‐valued operator. Additional estimates of maximum principle type are obtained, and these permit appropriate maximal extensions of the phase‐change relations. An example with pure advection indicates the limitations of these estimates and of the model developed here. We also consider and analyze the coupled pressure equation that determines the advective flux in the transport model. Copyright © 2015 John Wiley & Sons, Ltd.     

A steady barycentric Lagrange interpolation method for the 2D higher‐order time‐fractional telegraph equation with nonlocal boundary condition with error analysis

In this paper, we consider the numerical solution of the time‐fractional telegraph equation with a nonlocal boundary condition. A novel barycentric Lagrange interpolation collocation method is developed to solve this equation. Two difficulties have been sorted: the singularity of the integration and the higher accuracy. At the same, we put forward a steady barycentric Lagrange interpolation technique to overcome the new “Runge” phenomenon in computation. Error estimates of the barycentric Lagrange interpolation and the time‐fractional telegraph system for the present method are presented in Sobolev spaces. High convergence rates of the proposed method are obtained and are consisted with the numerical values. Especially in the time dimension, we get the error bound, for h‐refinement and for n_t‐density in the L₂ norms. The numerical results obtained show that the proposed numerical algorithm is accurate and computationally efficient for solving time‐fractional telegraph equation. Experiments demonstrate the high convergence rates of the proposed method are consisted with the theoretical values.     

A priori error estimates of an extrapolated space‐time discontinuous galerkin method for nonlinear convection‐diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection‐diffusion equation. We employ a combination of the discontinuous Galerkin finite element (DGFE) method for the space as well as time discretization. The linear diffusive and penalty terms are treated implicitly whereas the nonlinear convective term is treated by a special higher order explicit extrapolation from the previous time step, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the L^∞(L²) –norm and the L²(H¹) –seminorm with respect to the mesh size h and time step τ. Finally, we present an efficient solution strategy and numerical examples verifying the theoretical results. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1456–1482, 2010