A nonlinear PDE model for reconstructing a regular surface from sampled data using a level set formulation on triangular meshes

In this paper, we propose a nonlinear PDE model for reconstructing a regular surface from sampled data. At first, we show the existence and the uniqueness of a viscosity solution to this problem. Then we propose a numerical scheme for solving the nonlinear level set equation on unstructured triangulations adapted to the data sample. We show the consistency of this scheme. In addition, we show how to compute nodewise first and second order derivatives. Some application examples of curve or surface construction are provided to illustrate the potential and to demonstrate the accuracy of this method.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

A completely iterative method for the infinite domain electrostatic problem with nonlinear dielectric media

We present an iterative method for the solution of the exterior all-space electrostatic problem for nonlinear dielectric media. The electric potential is specified on interior boundaries and the electric field decays at infinity. Our approach uses a natural variational formulation based on the total energy of the nonlinear dielectric medium subject to boundary conditions. The problem is decomposed into an exterior calculation and an interior calculation with the boundary-specified electric potentials imposed as constraints between them. Together, these enable an iterative method that is based on the variational formulation. In contrast to direct solution of the electrostatic problems, we avoid the construction, storage and solution of dense and large linear systems. This provides important advantages for multiphysics problems that couple the linear electrostatic Poisson problem to nonlinear physics: the latter necessarily involves iterative approaches, and our approach replaces a large number of direct solves for the electrostatics with an iterative algorithm that can be coupled to the iterations of the nonlinear problem. We present examples applying the method to inhomogeneous, anisotropic nonlinear dielectrics. A key advantage of our variational formulation is that we require only the free-space, isotropic, homogeneous Greens function for all these settings.     

Refractive index imaging from radiative transfer equation-based reconstruction algorithm: Fundamentals

We present the first reconstruction algorithm for refractive index imaging, which is based on the radiative transfer equation (RTE). An objective function is iteratively minimized to find a solution to the problem of inversion of the refractive index field. The function describes the discrepancies of the emerging light measurements on the surface of the sample to be probed with predicted data from the corresponding numerical model. The unknown refractive index field is updated within each reconstruction iteration according to a search direction on the index distribution given by the adjoint model to the RTE. In this paper, emphasis is placed on the theoretical aspects. Preliminary tests are demonstrated on generic phantoms.     

Preconditioning techniques for iterative solvers in the Discrete Sources Method

Different preconditioning techniques for the iterative method MinRes as solver for the Discrete Sources Method (DSM) are presented. This semi-analytical method is used for light scattering computations by particles in the Mie scattering regime. Its numerical schema includes a linear least-squares problem commonly solved using the QR decomposition method. This could be the subject of numerical difficulties and instabilities for very large particles or particles with extreme geometry. In these cases, we showed that iterative methods with preconditioning techniques can provide a satisfying solution.In our previous paper, we studied four different iterative solvers (RGMRES, BiCGStab, BiCGStab(l), and MinRes) considering the performance and the accuracy of a solution. Here, we study several preconditioning techniques for the MinRes method for a variety of oblate and prolate spheroidal particles of different size and geometrical aspect ratio. Using preconditioning techniques we highly accelerated the iterative process especially for particles with a higher aspect ratio.     

Implementation and investigation of iterative solvers in the Discrete Sources Method

The implementation of iterative methods as solvers for the Discrete Sources Method (DSM) is presented. In this method, light scattering computation linear systems with dense and relative small matrices are generated. The linear systems are traditionally solved using the QR-decomposition method. For large particles or particles with extreme geometries even this commonly stable solver can fail. In these cases, we expect that iterative methods can provide a satisfying solution nevertheless.We will present our investigation in two consecutive papers. Here, we study four different iterative solvers (RGMRES, BiCGStab, BiCGStab(l), and MinRes) considering the performance and the accuracy for typical light scattering problems. Using these iterative methods we increased the quality of a solution, especially for oblate spheroids with a higher aspect ratio. Preconditioning technique is considered in the following paper.     

Inverse transient radiation analysis in one-dimensional non-homogeneous participating slabs using particle swarm optimization algorithms

Particle Swarm Optimization (PSO) algorithms, including standard PSO, Stochastic PSO, and Multi-Phase PSO, are applied to solve the time-domain inverse transient radiation problems in the present research. Time-resolved transmittance and reflectance signals of four different measuring models serve as the measurement data, which estimate absorption, scattering coefficients, and geometric position within one-dimensional non-homogeneous media by inverse simulation. To check retrieval performances and accuracies of PSO-based approaches, four different inverse transient radiation cases are investigated to deal with one homogeneous layer, two-layer, three-layer, and continuous non-homogenous media. The influences of different searching ranges, swarm sizes, and maximum fly velocities on the fitness function of PSO are discussed. Meanwhile, the effects of measurement errors on the reconstruction accuracy are also investigated. All the results confirm that radiative parameters could be estimated accurately with measurement noise using PSO-based approaches.     

Study on inverse taper based mode transformer for low loss coupling between silicon wire waveguide and lensed fiber

We have simulated the coupling loss of three types of Inverse Taper and taper-lensed fiber using three dimensional (3D) semi-vectorial beam propagation methods (BPM) respectively. Our results showed that the performances of exponential inverse taper and quadratic inverse taper were better than the commonly used linear inverse taper. Especially, for TM mode the improvement in the reduction of devices size is 53% and 136% for exponential and quadratic inverse taper compared with the linear inverse taper.     

Vector Preisach hysteresis modeling: Measurement, identification and application

The paper presents a Preisach model to simulate the vector hysteresis properties of ferromagnetic materials. The vector behavior has been studied using a single sheet tester with a disk-shaped specimen at low frequency. The locus of the magnetic flux density vector has been controlled by a digital measurement system. An inverse vector Preisach hysteresis model has been developed and identified by applying the measured data. Finally, the inverse model has been inserted into a finite element procedure through the fixed point technique and the reduced magnetic scalar potential formulation to simulate the measurement system. The applicability of the measurement system as well as the developed model has been proven by comparing measured and simulated results.     

The solution of the global relation for the derivative nonlinear Schrödinger equation on the half-line

We consider initial-boundary value problems for the derivative nonlinear Schrödinger (DNLS) equation on the half-line x>0. In a previous work, we showed that the solution q(x,t) can be expressed in terms of the solution of a Riemann-Hilbert problem with jump condition specified by the initial and boundary values of q(x,t). However, for a well-posed problem, only part of the boundary values can be prescribed; the remaining boundary data cannot be independently specified, but are determined by the so-called global relation. In general, an effective solution of the problem therefore requires solving the global relation. Here, we present the solution of the global relation in terms of the solution of a system of nonlinear integral equations. This also provides a construction of the Dirichlet-to-Neumann map for the DNLS equation on the half-line.