Global existence for the nonlinear equations of crystal optics
In this paper we give a proof of the existence and uniqueness of smooth solutions for the nonlinear semiconductor Boltzmann equation. The method used allows to obtain global existence in time and uniqueness for dimensions 1 and 2. For dimension 3 we can only assure local existence in time and uniqueness. First, we define a sequence of solutions for a linearized equation and then, we prove the strong convergence of the sequence in a suitable space. The metod relies on the use of interpolation estimates...
We describe both the classical lagrangian and the Eulerian methods for first order Hamilton–Jacobi equations of geometric optic type. We then explain the basic structure of the software and how new solvers/models can be added to it. A selection of numerical examples are presented.
We describe both the classical Lagrangian and the Eulerian methods for first order Hamilton–Jacobi equations of geometric optic type. We then explain the basic structure of the software and how new solvers/models can be added to it. A selection of numerical examples are presented.
Poisson and Poisson-Boltzmann equations (PE and PBE) are widely used in molecular modeling to estimate the electrostatic contribution to the free energy of a system. In such applications, PE often needs to be solved multiple times for a large number of system configurations. This can rapidly become a highly demanding computational task. To accelerate such calculations we implemented a graphical processing unit (GPU) PE solver described in this work. The GPU solver performance is compared to that...
We discuss the propagation of electromagnetic waves of a special form through an inhomogeneous isotropic medium which has a cylindrical symmetry and a nonlinear dielectric response. For the case where this response is of self-focusing type the problem is treated in [1]. Here we continue this study by dealing with a defocusing dielectric response. This tends to inhibit the guidance properties of the medium and so guidance can only be expected provided that the cylindrical stratification is such that...
A simple and effective method based on Haar wavelets is proposed for the solution of Pocklington’s integral equation. The properties of Haar wavelets are first given. These wavelets are utilized to reduce the solution of Pocklington’s integral equation to the solution of algebraic equations. In order to save memory and computation time, we apply a threshold procedure to obtain sparse algebraic equations. Through numerical examples, performance of the present method is investigated concerning the...
We derive the high frequency limit of the Helmholtz equations in terms of quadratic observables. We prove that it can be written as a stationary Liouville equation with source terms. Our method is based on the Wigner Transform, which is a classical tool for evolution dispersive equations. We extend its use to the stationary case after an appropriate scaling of the Helmholtz equation. Several specific difficulties arise here; first, the identification of the source term ( which does not share the...
Low order edge elements are widely used for electromagnetic field problems. Higher order edge approximations are receiving increasing interest but their definition become rather complex. In this paper we propose a simple definition for Whitney edge elements of polynomial degree higher than one. We give a geometrical localization of all degrees of freedom over particular edges and provide a basis for these elements on simplicial meshes. As for Whitney edge elements of degree one, the basis is...
We propose transmission conditions of order 1, 2 and 3 approximating the shielding behaviour of thin conducting curved sheets for the magneto-quasistatic eddy current model in 2D. This model reduction applies to sheets whose thicknesses ε are at the order of the skin depth or essentially smaller. The sheet has itself not to be resolved, only its midline is represented by an interface. The computation is directly in one step with almost no additional cost. We prove the well-posedness w.r.t. to...
We propose transmission conditions of order 1, 2 and 3 approximating the shielding behaviour of thin conducting curved sheets for the magneto-quasistatic eddy current model in 2D. This model reduction applies to sheets whose thicknesses ε are at the order of the skin depth or essentially smaller. The sheet has itself not to be resolved, only its midline is represented by an interface. The computation is directly in one step with almost no additional cost. We prove the well-posedness w.r.t. to...
Cell motility is an integral part of a diverse set of biological processes. The quest for mathematical models of cell motility has prompted the development of automated approaches for gathering quantitative data on cell morphology, and the distribution of molecular players involved in cell motility. Here we review recent approaches for quantifying cell motility, including automated cell segmentation and tracking. Secondly, we present our own novel...
We study the Maxwell-Landau-Lifshitz system for highly oscillating initial data, with characteristic frequencies O(1 / ε) and amplitude O(1), over long time intervals O(1 / ε), in the limit ε → 0. We show that a nonlinear Schrödinger equation gives a good approximation for the envelope of the solution in the time interval under consideration. This extends previous results of Colin and Lannes [1]. This text is a short version of the article [5].