Search

In this paper, we prove a Carleman estimate for fully-discrete approximations of parabolic operators in which the discrete parameters $h$ and $\triangle t$ are connected to the large Carleman parameter. We use this estimate to obtain relaxed observability inequalities which yield, by duality, controllability results for fully-discrete linear and semilinear parabolic equations.

We develop a new model of laser-matter interaction based on Optical Bloch Equations, which includes photo-ionization, impact ionization, and various relaxation processes typical of dielectric materials. This approach is able to describe the temporal evolution of the electron dynamics in the conduction band driven by few-cycle laser pulses of any wavelength. Moreover, the nonlinear polarization response of both centrosymmetric and non-centrosymmetric materials can be described while ensuring the proper selection rules for the harmonics emission.

The purpose of this article is to introduce radial basis function, (RBFs), methods for solving null control problems for the Stokes system with few internal scalar controls and Dirichlet or Navier slip boundary conditions. To the best of our knowledge, it has not been reported in the literature any numerical approximation through RBFs to solve the direct Stokes problem with Navier slip boundary conditions. In this paper we fill this gap to show its application for solving the null control problem for the Stokes system. To achieve this goal, we introduce two radial basis function solvers, one global and the other local, to discretized the primal and adjoint systems related to the control problem. Both techniques are based on divergence free global RBFs. Stability analysis for these methods is performed in terms of the spectral radius of the corresponding Gram matrices. By using a conjugate gradient algorithm, adapted to the radial basis function setting, the control problem is solved. Several test problems in two dimensions are numerically solved by these RBFs methods to test their feasibility. The solutions to these problems are also obtained by finite elements techniques, (FE), to compare their relative performance.

PDE-constrained optimization problems have been barely solved by radial basis functions (RBFs) methods [Pearson, 2013]. It is well known that RBF methods can attain an exponential rate of convergence when $C^{\infty}$ kernels are used, also, these techniques, which are truly scattered, are known to be flexible to discretize complex boundaries in several dimensions. On the other hand, exponential convergence implies an exponential growth of the condition number of the Gram matrix associated with these meshfree methods and global collocation techniques are known to be computationally expensive. In this paper, and in the context of optimal constrained optimization problems, we aim to explore a possible answer to both problems. Specifically, we introduce two local RBF methods: LAM-DQ based in the combination of an asymmetric local method (LAM), inspired in local Hermite interpolation (LHI), with the differential quadrature method (DQ), and LAM-LAM which use two times the local asymmetric method. The efficiency of these local methods against global collocation by solving several synthetic convection-diffusion control problems is analyzed. In this article, we also propose a preconditioning technique and treat the ill-conditioning problem by using extended arithmetic precision. We think that these local methods, which are highly parallelizable, shows a possible way to solve massive optimization control problems in an efficient way.

Let $L$ be a fixed branch -- that is, an irreducible germ of curve -- on a normal surface singularity $X$. If $A,B$ are two other branches, define $u_L(A,B) := \dfrac{(L \cdot A) \: (L \cdot B)}{A \cdot B}$, where $A \cdot B$ denotes the intersection number of $A$ and $B$. Call $X$ arborescent if all the dual graphs of its resolutions are trees. In a previous paper, the first three authors extended a 1985 theorem of P{\l}oski by proving that whenever $X$ is arborescent, the function $u_L$ is an ultrametric on the set of branches on $X$ different from $L$. In the present paper we prove that, conversely, if $u_L$ is an ultrametric, then $X$ is arborescent. We also show that for any normal surface singularity, one may find arbitrarily large sets of branches on $X$, characterized uniquely in terms of the topology of the resolutions of their sum, in restriction to which $u_L$ is still an ultrametric. Moreover, we describe the associated tree in terms of the dual graphs of such resolutions. Then we extend our setting by allowing $L$ to be an arbitrary semivaluation on $X$ and by defining $u_L$ on a suitable space of semivaluations. We prove that any such function is again an ultrametric if and only if $X$ is arborescent, and without any restriction on $X$ we exhibit special subspaces of the space of semivaluations in restriction to which $u_L$ is still an ultrametric., Comment: 50 pages, 14 figures. Final version

We investigate terahertz emission from two-color fs-laser-induced microplasmas. Under strongest focusing conditions, microplasmas are shown to act as point-sources for broadband terahertz-to-far-infrared radiation, where the emission bandwidth is determined by the plasma density. Semi-analytical modeling allows us to identify scaling laws with respect to important laser parameters. In particular, we find that the optical-to-THz conversion efficiency crucially depends on the focusing conditions. We use this insight to demonstrate by means of Maxwell-consistent 3D simulations, that for only 10-$\mu$J laser energy a conversion efficiency well above $10^{-4}$ can be achieved.