A variational method for identifying a spacewise-dependent heat source

The inverse problem of determining a spacewise-dependent heat source for the parabolic heat equation using the usual conditions of the direct problem and information from one supplementary temperature measurement at a given instant of time is studied. This spacewise-dependent temperature measurement ensures that this inverse problem has a unique solution, but the solution is unstable and hence the problem is ill-posed. We propose a variational conjugate gradient-type iterative algorithm for the stable reconstruction of the heat source based on a sequence of well-posed direct problems for the parabolic heat equation which are solved at each iteration step using the boundary element method. The instability is overcome by stopping the iterative procedure at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented which have the input measured data perturbed by increasing amounts of random noise. The numerical results show that the proposed procedure yields stable and accurate numerical approximations after only a few iterations.

Determination of a spacewise dependent heat source

This paper investigates the inverse problem of determining a spacewise dependent heat source in the parabolic heat equation using the usual conditions of the direct problem and information from a supplementary temperature measurement at a given single instant of time. The spacewise dependent temperature measurement ensures that the inverse problem has a unique solution, but this solution is unstable, hence the problem is ill-posed. For this inverse problem, we propose an iterative algorithm based on a sequence of well-posed direct problems which are solved at each iteration step using the boundary element method (BEM). The instability is overcome by stopping the iterations at the first iteration for which the discrepancy principle is satisfied. Numerical results are presented for various typical benchmark test examples which have the input measured data perturbed by increasing amounts of random noise.

A probabilistic indirect adaptive control for systems with input-dependent noise

A probabilistic indirect adaptive controller is proposed for the general nonlinear multivariate class of discrete time system. The proposed probabilistic framework incorporates input–dependent noise prediction parameters in the derivation of the optimal control law. Moreover, because noise can be nonstationary in practice, the proposed adaptive control algorithm provides an elegant method for estimating and tracking the noise. For illustration purposes, the developed method is applied to the affine class of nonlinear multivariate discrete time systems and the desired result is obtained: the optimal control law is determined by solving a cubic equation and the distribution of the tracking error is shown to be Gaussian with zero mean. The efficiency of the proposed scheme is demonstrated numerically through the simulation of an affine nonlinear system.

Impulsive Partial Hyperbolic Functional Differential Equations of Fractional Order with State-Dependent Delay

MSC 2010: 26A33, 34A37, 34K37, 34K40, 35R11

Local Energy Decay in Even Dimensions for the Wave Equation with a Time-Periodic Non-Trapping Metric and Applications to Strichartz Estimates

2000 Mathematics Subject Classification: 35B40, 35L15.

Green's function for paraxial equation

The theory and experimental applications of optical Airy beams are in active development recently. The Airy beams are characterised by very special properties: they are non-diffractive and propagate along parabolic trajectories. Among the striking applications of the optical Airy beams are optical micro-manipulation implemented as the transport of small particles along the parabolic trajectory, Airy-Bessel linear light bullets, electron acceleration by the Airy beams, plasmonic energy routing. The detailed analysis of the mathematical aspects as well as physical interpretation of the electromagnetic Airy beams was done by considering the wave as a function of spatial coordinates only, related by the parabolic dependence between the transverse and the longitudinal coordinates. Their time dependence is assumed to be harmonic. Only a few papers consider a more general temporal dependence where such a relationship exists between the temporal and the spatial variables. This relationship is derived mostly by applying the Fourier transform to the expressions obtained for the harmonic time dependence or by a Fourier synthesis using the specific modulated spectrum near some central frequency. Spatial-temporal Airy pulses in the form of contour integrals is analysed near the caustic and the numerical solution of the nonlinear paraxial equation in time domain shows soliton shedding from the Airy pulse in Kerr medium. In this paper the explicitly time dependent solutions of the electromagnetic problem in the form of time-spatial pulses are derived in paraxial approximation through the Green's function for the paraxial equation. It is shown that a Gaussian and an Airy pulse can be obtained by applying the Green's function to a proper source current. We emphasize that the processes in time domain are directional, which leads to unexpected conclusions especially for the paraxial approximation.

Gain dependent pulse regimes transitions in a dissipative dispersion-managed fibre laser

For the first time, we demonstrate the possibility to switch between three distinct pulse regimes in a dissipative dispersion-managed (DM) fibre laser by solely controlling the gain saturation energy. Nonlinear Schrödinger equation based simulations show the transitions between hyper-Gaussian similaritons, parabolic similaritons, and dissipative solitons in the same laser cavity. It is also shown that such transitions exist in a wide dispersion range from all-normal to slightly net-normal dispersion. This work demonstrates that besides dispersion and filter managements gain saturation energy can be a new degree of freedom to manage pulse regimes in DM fibre lasers, which offers flexibility in designing ultrafast fibre lasers. Also, the result indicates that in contrast to conservative soliton lasers whose intensity profiles are unique, dissipative DM lasers show diversity in pulse shapes. The findings not only give a better understanding of pulse shaping mechanisms in mode-locked lasers, but also provide insight into dissipative systems.

Pulsed parabolic equation model of acoustic transmission (PPEMAT)

Underwater sound is very important in the field of oceanography where it is used for remote sensing in much the same way that radar is used in atmospheric studies. One way to mathematically model sound propagation in the ocean is by using the parabolic-equation method, a technique that allows range dependent environmental parameters. More importantly, this method can model sound transmission where the source emits either a pure tone or a short pulse of sound. Based on the parabolic approximation method and using the split-step Fourier algorithm, a computer model for underwater sound propagation was designed and implemented. This computer model differs from previous models in its use of the interactive mode, structured programming, modular design, and state-of-the-art graphics displays. In addition, the model maximizes the efficiency of computer time through synchronization of loosely coupled dual processors and the design of a restart capability. Since the model is designed for adaptability and for users with limited computer skills, it is anticipated that it will have many applications in the scientific community.

Dispersive estimates for the Schrödinger equation

In this document we explore the issue of $L^1\to L^\infty$ estimates for the solution operator of the linear Schr\"{o}dinger equation, \begin{align*} iu_t-\Delta u+Vu&=0 &u(x,0)=f(x)\in \mathcal S(\R^n). \end{align*} We focus particularly on the five and seven dimensional cases. We prove that the solution operator precomposed with projection onto the absolutely continuous spectrum of $H=-\Delta+V$ satisfies the following estimate $\|e^{itH} P_{ac}(H)\|_{L^1\to L^\infty} \lesssim |t|^{-\frac{n}{2}}$ under certain conditions on the potential $V$. Specifically, we prove the dispersive estimate is satisfied with optimal assumptions on smoothness, that is $V\in C^{\frac{n-3}{2}}(\R^n)$ for $n=5,7$ assuming that zero is regular, $|V(x)|\lesssim \langle x\rangle^{-\beta}$ and $|\nabla^j V(x)|\lesssim \langle x\rangle^{-\alpha}$, $1\leq j\leq \frac{n-3}{2}$ for some $\beta>\frac{3n+5}{2}$ and $\alpha>3,8$ in dimensions five and seven respectively. We also show that for the five dimensional result one only needs that $|V(x)|\lesssim \langle x\rangle^{-4-}$ in addition to the assumptions on the derivative and regularity of the potential. This more than cuts in half the required decay rate in the first chapter. Finally we consider a problem involving the non-linear Schr\"{o}dinger equation. In particular, we consider the following equation that arises in fiber optic communication systems, \begin{align*} iu_t+d(t) u_{xx}+|u|^2 u=0. \end{align*} We can reduce this to a non-linear, non-local eigenvalue equation that describes the so-called dispersion management solitons. We prove that the dispersion management solitons decay exponentially in $x$ and in the Fourier transform of $x$.

Centre-of-mass motion in multi-particle Schrodinger-Newton dynamics

We investigate the implication of the nonlinear and non-local multi-particle Schrodinger-Newton equation for the motion of the mass centre of an extended multi-particle object, giving self-contained and comprehensible derivations. In particular, we discuss two opposite limiting cases. In the first case, the width of the centre-of-mass wave packet is assumed much larger than the actual extent of the object, in the second case it is assumed much smaller. Both cases result in nonlinear deviations from ordinary free Schrodinger evolution for the centre of mass. On a general conceptual level we include some discussion in order to clarify the physical basis and intention for studying the Schrodinger-Newton equation.

Dynamic instabilities in scalar neural field equations with space-dependent delays

In this paper we consider a class of scalar integral equations with a form of space-dependent delay. These non-local models arise naturally when modelling neural tissue with active axons and passive dendrites. Such systems are known to support a dynamic (oscillatory) Turing instability of the homogeneous steady state. In this paper we develop a weakly nonlinear analysis of the travelling and standing waves that form beyond the point of instability. The appropriate amplitude equations are found to be the coupled mean-field Ginzburg-Landau equations describing a Turing-Hopf bifurcation with modulation group velocity of O(1). Importantly we are able to obtain the coefficients of terms in the amplitude equations in terms of integral transforms of the spatio-temporal kernels defining the neural field equation of interest. Indeed our results cover not only models with axonal or dendritic delays but those which are described by a more general distribution of delayed spatio-temporal interactions. We illustrate the predictive power of this form of analysis with comparison against direct numerical simulations, paying particular attention to the competition between standing and travelling waves and the onset of Benjamin-Feir instabilities.

Measure-valued mass evolution problems with flux boundary conditions and solution-dependent velocities

In this paper we prove well-posedness for a measure-valued continuity equation with solution-dependent velocity and flux boundary conditions, posed on a bounded one-dimensional domain. We generalize the results of an earlier paper [J. Differential Equations, 259 (2015), pp. 10681097] to settings where the dynamics are driven by interactions. In a forward-Euler-like approach, we construct a time-discretized version of the original problem and employ those results as a building block within each subinterval. A limit solution is obtained as the mesh size of the time discretization goes to zero. Moreover, the limit is independent of the specific way of partitioning the time interval [0, T]. This paper is partially based on results presented in Chapter 5 of [Evolution Equations for Systems Governed by Social Interactions, Ph.D. thesis, Eindhoven University of Technology, 2015], while a number of issues that were still open there are now resolved.

Specific shear-dependent viscoelastic third-grade fluid model

A specific modified constitutive equation for a third-grade fluid is proposed so that the model be suitable for applications where shear-thinning or shear-thickening may occur. For that, we use the Cosserat theory approach reducing the exact three-dimensional equations to a system depending only on time and on a single spatial variable. This one-dimensional system is obtained by integrating the linear momentum equation over the cross-section of the tube, taking a velocity field approximation provided by the Cosserat theory. From this reduced system, we obtain the unsteady equations for the wall shear stress and mean pressure gradient depending on the volume flow rate, Womersley number, viscoelastic coefficient and flow index over a finite section of the tube geometry with constant circular cross-section.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.