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.     

Lie-Mei symmetry and conserved quantities of the Rosenberg problem

The Rosenberg problem is a typical but not too complicated problem of nonholonomic mechanical systems.The Lie-Mei symmetry and the conserved quantities of the Rosenberg problem are studied.For the Rosenberg problem,the Lie and the Mei symmetries for the equation are obtained,the conserved quantities are deduced from them and then the definition and the criterion for the Lie-Mei symmetry of the Rosenberg problem are derived.Finally,the Hojman conserved quantity and the Mei conserved quantity are deduced from the Lie-Mei symmetry.     

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.     

Positive solution of singularly perturbed Dirichlet problem with singularities

A class of singularly perturbed boundary value problem with singularities is considered. Introducing the stretched variables, the boundary layer corrective terms near x = 0 and x = 1 are constructed. Under suitable conditions, by using the theory of differential inequalities the existence and asymptotic behavior of solution for boundary value problem are proved, uniformly valid asymptotic expansion of solution with boundary layers are obtained,     

A Remark on the Cauchy Problem of 1D Compressible Navier-Stokes Equations with Density-dependent Viscosity Coefcients

This paper is concerned with the Cauchy problems of one-dimensional compressible Navier-Stokes equations with density-dependent viscosity coefcients.By assumingρ0∈L1(R),we will prove the existence of weak solutions to the Cauchy problems forθ0.This will improve results in Jiu and Xin’s paper(Kinet.Relat.Models,1(2):313–330(2008))in whichθ12is required.In addition,We will study the large time asymptotic behavior of such weak solutions.