A Schr(o)idinger formulation research for light beam propagation

The wave equation of light beam propagation was written in the form of an axial-coordinate-dependent Schrdinger equation, and the expectation value of a dynamical variable, the trial function of variational approach and the ABCD law were discussed by use of quantum mechanics approach. In view of the evolution equations of expectation values of dynamical variables in the framework of quantum mechanics, the definition of a potential function representing the beam propagation stability and its universal formula with the quality factor, the universal formula of beam width and curvature radius for a paraxial beam and cylindrically symmetric non-paraxial beam, the general formula of second derivative of beam width with respect to the axial coordinate of beam for a paraxial beam, and the general criteria of the conservation of beam quality factor and the existence of a potential well of a potential function for a paraxial beam, were given or derived, respectively. Starting with the same trial function, the comparative research of our formulation with variational approach was done, which gave some further insight into the physical nature of a beam propagation parameters. The ABCD law of non-paraxial beam was discussed in terms of the definition of the non-paraxial expectation value of a dynamical variable for the first time. The applications to the media of constant second derivative of beam width with respect to the axial coordinate of a beam, square law media and the media of constant refractive index in the momentum representation were discussed, respectively.     

A Schr?dinger formulation research for light beam propagation

The wave equation of light beam propagation was written in the form of an axial-coordinate-dependent Schr?dinger equation, and the expectation value of a dynamical variable, the trial function of variational approach and the ABCD law were discussed by use of quantum mechanics approach. In view of the evolution equations of expectation values of dynamical variables in the framework of quantum mechanics, the definition of a potential function representing the beam propagation stability and its universal formula with the quality factor, the universal formula of beam width and curvature radius for a paraxial beam and cylindrically symmetric non-paraxial beam, the general formula of second derivative of beam width with respect to the axial coordinate of beam for a paraxial beam, and the general criteria of the conservation of beam quality factor and the existence of a potential well of a potential function for a paraxial beam, were given or derived, respectively. Starting with the same trial function, the comparative research of our formulation with variational approach was done, which gave some further insight into the physical nature of a beam propagation parameters. The ABCD law of nonparaxial beam was discussed in terms of the definition of the non-paraxial expectation value of a dynamical variable for the first time. The applications to the media of constant second derivative of beam width with respect to the axial coordinate of a beam, square law media and the media of constant refractive index in the momentum representation were discussed, respectively.     

A Schrödinger formulation research for light beam propagation through the media of complex refractive index

The Helmhotzequation oflightbeam propagating through a mediumof complex refractive indexis reduced totheaxial-coordinate-dependent Schrödingerequation ofcomplexpotential. Thenewbravector,thenewexpectationvalueof a dynamical variableandtheextended Heisenberg picturearedefinedby the inverse oftheevolutionoperator instead ofitsHermitianadjoint, andthecomplexbeam propagation parametersdefinedin terms ofthenewexpectationvalue, the complexABCD law and theABCD formulationof the Huygens’ integralarediscussedin terms of quantummechanics. Itis shown thattheevolutionequationsof the complex beampropagation parameters are the sameas those ofthebeam propagation parametersof beampropagating throughamedium ofreal refractiveindex. The researchon anopticalsystem oftheconservative complexbeam quality factorshowsthat the complexABCD lawholds, the evolution ofitscoordinate operator and the momentumoperator islinear,andtheHuygens’integral isof theABCD formulation.     

A boundary integral equation method for the transmission eigenvalue problem

We propose a new integral equation formulation to characterize and compute transmission eigenvalues for constant refractive index that play an important role in inverse scattering problems for penetrable media. As opposed to the recently developed approach by Cossonnière and Haddar [1,2] which relies on a two by two system of boundary integral equations our analysis is based on only one integral equation in terms of Dirichlet-to-Neumann or Robin-to-Dirichlet operators which results in a noticeable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further we employ the numerical algorithm for analytic non-linear eigenvalue problems that was recently proposed by Beyn [3] for the numerical computation of transmission eigenvalues via this new integral equation.     

Preconditioning techniques for the iterative solution of scattering problems

We consider a time-harmonic electromagnetic scattering problem for an inhomogeneous medium. Some symmetry hypotheses on the refractive index of the medium and on the electromagnetic fields allow to reduce this problem to a two-dimensional scattering problem. This boundary value problem is defined on an unbounded domain, so its numerical solution cannot be obtained by a straightforward application of usual methods, such as for example finite difference methods, and finite element methods. A possible way to overcome this difficulty is given by an equivalent integral formulation of this problem, where the scattered field can be computed from the solution of a Fredholm integral equation of second kind. The numerical approximation of this problem usually produces large dense linear systems. We consider usual iterative methods for the solution of such linear systems, and we study some preconditioning techniques to improve the efficiency of these methods. We show some numerical results obtained with two well known Krylov subspace methods, i.e., Bi-CGSTAB and GMRES.     

An analytic method for the initial value problem of the electric field system in vertical inhomogeneous anisotropic media

The time-dependent system of partial differential equations of the second order describing the electric wave propagation in vertically inhomogeneous electrically and magnetically biaxial anisotropic media is considered. A new analytical method for solving an initial value problem for this system is the main object of the paper. This method consists in the following: the initial value problem is written in terms of Fourier images with respect to lateral space variables, then the resulting problem is reduced to an operator integral equation. After that the operator integral equation is solved by the method of successive approximations. Finally, a solution of the original initial value problem is found by the inverse Fourier transform.     

Geometrical study of paraxial light beam transmission in free space

By introducing an imaginary space transform curvature ρ_s, a complex space called Riemannian space is constructed, in which the light propagating in free space has the trajectory of straight line while propagating. Moreover, this curvature couples with that of the wave front of the paraxial beam ρ_w, and therefore a complex curvature ρ_c is constructed, which can be employed to investigate the behavior of the light transmission and to generalize the ABCD law. Project supported by the National Hi-Tech Inertial Confinement Fusion Committee, the Guangdong Natural Science Foundation the Postdoctoral Foundation of Guangdong and National Postdoctoral Foundation of China.     

A boundary integral equation for the transmission eigenvalue problem for Maxwell equation

We propose a new integral equation formulation to characterize and compute transmission eigenvalues in electromagnetic scattering. As opposed to the approach that was recently developed by Cakoni, Haddar and Meng (2015) which relies on a two‐by‐two system of boundary integral equations, our analysis is based on only one integral equation in terms of the electric‐to‐magnetic boundary trace operator that results in a simplification of the theory and in a considerable reduction of computational costs. We establish Fredholm properties of the integral operators and their analytic dependence on the wave number. Further, we use the numerical algorithm for analytic nonlinear eigenvalue problems that was recently proposed by Beyn (2012) for the numerical computation of the transmission eigenvalues via this new integral equation.     

A Formulation for Water Waves over Topography

Many interesting free-surface flow problems involve a varying bottom. Examples of such flows include ocean waves propagating over topography, the breaking of waves on a beach, and the free surface of a uniform flow over a localized bump. We present here a formulation for such flows that is general and, from the outset, demonstrates the wave character of the free-surface evolution. The evolution of the free surface is governed by a system of equations consisting of a nonlinear wave-like partial differential equation coupled to a time-independent linear integral equation. We assume that the free-surface deformation is weakly nonlinear, but make no a priori assumption about the scale or amplitude of the topography. We also extend the formulation to include the effect of mean flows and surface tension. We show how this formulation gives some of the well-known limits for such problems once assumptions about the amplitude and scale of the topography are made.     

一类柔性梁系统的谱分析和半群生成

讨论了一类双臂三关节柔性梁系统的分析问题．首先,建立了一个与柔性梁的偏微分方程组及初值边值条件相应的希尔伯特空间中的一阶发展系统．接着讨论系统算子的谱性质和半群性质．最后借助系统算子的谱性质和半群性质提出并证明了柔性梁系统的指数稳定性．     

The operator equations of Lippmann–Schwinger type for acoustic and electromagnetic scattering problems in L 2

This article is concerned with the scattering of acoustic and electromagnetic time harmonic plane waves by an inhomogeneous medium. These problems can be translated into volume integral equations of the second kind – the most prominent example is the Lippmann–Schwinger integral equation. In this work, we study a particular class of scattering problems where the integral operator in the corresponding operator equation of Lippmann–Schwinger type fails to be compact. Such integral equations typically arise if the modelling of the inhomogeneous medium necessitates space-dependent coefficients in the highest order terms of the underlying partial differential equation. The two examples treated here are acoustic scattering from a medium with a space-dependent material density and electromagnetic medium scattering where both the electric permittivity and the magnetic permeability vary. In these cases, Riesz theory is not applicable for the solution of the arising integral equations of Lippmann–Schwinger type. Therefore, we show that positivity assumptions on the relative material parameters allow to prove positivity of the arising volume potentials in tailor-made weighted spaces of square integrable functions. This result merely holds for imaginary wavenumber and we exploit a compactness argument to conclude that the arising integral equations are of Fredholm type, even if the integral operators themselves are not compact. Finally, we explain how the solution of the integral equations in L ² affects the notion of a solution of the scattering problem and illustrate why the order of convergence of a Galerkin scheme set up in L ² does not suffer from our L ² setting, compared to schemes in higher order Sobolev spaces.     

N-group neutron transport theory: A criticality problem in slab geometry

The steady-state equation for N-group neutron transport in slab geometry is written as an integral equation. A spectral analysis is made of the integral operator and related to the criticality problem. The method depends on a representation for the resolvent kernel for a subcritical slab and on analytic continuation in a complex parameter to characterize eigenvalues in terms of singularities of the resolvent. The analytic continuation is based on a bifurcation analysis of some nonlinear matrix integral equations whose solutions provide a matrix Wiener-Hopf factorization of the Fourier transform of the kernel of the transport operator.     

Nonlinear transmission eigenvalue problem describing TE wave propagation in two-layered cylindrical dielectric waveguides

The paper is concerned with propagation of surface TE waves in a circular nonhomogeneous two-layered dielectric waveguide filled with a Kerr nonlinear medium. The problem is reduced to the analysis of a nonlinear integral equation with a kernel in the form of a Green’s function. The existence of propagating TE waves is proved using the contraction mapping method. For the numerical solution of the problem, two methods are proposed: an iterative algorithm (whose convergence is proved) and a method based on solving an auxiliary Cauchy problem (the shooting method). The existence of roots of the dispersion equation (propagation constants of the waveguide) is proved. Conditions under which k waves can propagate are obtained, and regions of localization of the corresponding propagation constants are found.     

Diffraction by a wedge at skew incidence: integral representations of Cauchy-Carleman for the electromagnetic fields

An electromagnetic diffraction problem in a wedge shaped region is reduced to a system of coupled functional difference equations by means of Sommerfeld integrals and Malyuzhinets theorem. By introducing an integral operator it is shown that the solutions of this system of functional equations can be defined in terms of integral representations whose kernels are solutions of a singular integral equation of Cauchy-Carleman type for which an explicit solution is given.     

Integral equation formulation and flutter analysis of damped non-conservative Timoshenko beams

This paper presents a mathematical modeling and a numerical methodological approach based on integral equation formulation and radial basis functions (RBF) for the dynamic behaviour damped Timoshenko beams under follower forces. Using the fundamental solution of the main operator and the RBF, the governing non-self-adjoint partial differential equation is transformed into an integral equation. Based on the harmonic assumption and internal concatenation points an eigenvalue problem is obtained and numerically solved. The flutter analysis, complex modes and the frequency load responses are investigated for beams under various subtangential follower loads, aspect ratios, internal and viscous damping.     

On the local solvability of evolution equations for the interface of two filtering fluids in the class of analytic functions

We consider a system of equations that describes the evolution of the interface of two fluids in the two-dimensional problem of their combined filtration in a porous homogeneous medium. This system contains a nonlinear integro-differential equation with a singular integral over an unknown contour together with a Fredholm integral equation of the second kind for the jump of the velocity vector potential on the movable boundary. We prove the unique solvability of such a system on a small time interval for the case in which the initial shape is defined parametrically and the functions describing the dependence of the coordinates of the point on the line on the parameter admit analytic continuations into the complex plane.     

Time harmonic acoustic scattering in anisotropic media

The scattering problem of a plane or a point source generated wave is considered for the case where both the medium of propagation and the interior of the scatterer exhibit their own anisotropies. A particular redirected gradient operator is introduced, which carries all directional characteristics of the anisotropic medium. Once the fundamental solution is obtained, integral representations for the scattered as well as for the interior and the total fields are generated. For such media even the handling of the singularities, in generating integral representations, depends on the characteristics of the particular medium. A modified, also medium dependent, radiation condition is introduced. Detailed asymptotic analysis leads to an integral representation for the scattering amplitude. The associated energy functionals are presented and the relative cross sections are also defined. Copyright © 2005 John Wiley & Sons, Ltd.     

Direct reconstruction of complex refractive index distribution from boundary measurements

We obtain the reconstruction of the refractive index distribution of body based on the intensity and normal derivative of the intensity measurements. The Helmholtz equation is inverted either directly or indirectly through repeated implementation of the forward operator and its adjoint, for recovering the complex refractive index distribution. We do not adopt the procedure of recovery of phase (normally required for complete knowledge distribution). We derive certain sensitivity relations which is used for the easy computation of the Jacobian. Our procedure successfully reconstructs the real and imaginary parts of the complex refractive index from the measurement of the two data types derived from the complex amplitude at the boundary. Our other interest is the reconstruction of the spectroscopic variations of optical absorption coefficients and visco-elastic properties of a tissue which is extremely useful in diagnostic medicines. The research is on progress and some results are available. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Symmetric Galerkin Boundary Element Method for Linear Poroelasticity

In linear poroelasticity so far only collocation boundary element methods have been available. However, in some applications, e.g., when coupling with finite elements is desired, a symmetric formulation is preferable. Choosing a Galerkin approach which involves the second boundary integral equation, such a formulation is possible. Here, a previously presented integration by part technique for the regularization of the first boundary integral equation is extended to the second boundary integral equation as well. While the weakly singular representation of the double layer operator has been presented before, the emphasis lies here on the so called hyper-singular boundary integral operator. Due to the regularization, this operator can be evaluated numerically and, hence, be used within a numerical scheme for the first time. Different numerical studies will be presented to show the behavior of the established symmetric Galerkin boundary element method, also comparing it with collocation boundary element methods. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)