Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics ...Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics needs such numerical methods. These methods must be precise and useful in terms of memory capacity and time duration. In this paper, we study basic analytical and finite difference methods to determine the effective refractive index of AIGaAs-GaAs slab waveguide. Also, appropriate effective refractive index value is obtained with respect to number of grid points and number of matrix sizes. Finally, the validity of the obtained values by both methods is compared to using waveguide type.展开更多
In this paper, we propose a nearly analytic exponential time difference （NETD） method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approxima...In this paper, we propose a nearly analytic exponential time difference （NETD） method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approximate the high-order spatial differential operators and transform the seismic wave equations into semi-discrete ordinary differential equations （ODEs）. Then, the converted ODE system is solved by the exponential time difference （ETD） method. We investigate the properties of NETD in detail, including the stability condition for 1-D and 2-D cases, the theoretical and relative errors, the numerical dispersion relation for the 2-D acoustic case, and the computational efficiency. In order to further validate the method, we apply it to simulating acoustic/elastic wave propagation in mul- tilayer models which have strong contrasts and complex heterogeneous media, e.g., the SEG model and the Mar- mousi model. From our theoretical analyses and numerical results, the NETD can suppress numerical dispersion effectively by using the displacement and gradient to approximate the high-order spatial derivatives. In addition, because NETD is based on the structure of the Lie group method which preserves the quantitative properties of differential equations, it can achieve more accurate results than the classical methods.展开更多
With the development of molecular imaging,Cherenkov optical imaging technology has been widely concerned.Most studies regard the partial boundary flux as a stochastic variable and reconstruct images based on the stead...With the development of molecular imaging,Cherenkov optical imaging technology has been widely concerned.Most studies regard the partial boundary flux as a stochastic variable and reconstruct images based on the steadystate diffusion equation.In this paper,time-variable will be considered and the Cherenkov radiation emission process will be regarded as a stochastic process.Based on the original steady-state diffusion equation,we first propose a stochastic partial differential equationmodel.The numerical solution to the stochastic partial differential model is carried out by using the finite element method.When the time resolution is high enough,the numerical solution of the stochastic diffusion equation is better than the numerical solution of the steady-state diffusion equation,which may provide a new way to alleviate the problem of Cherenkov luminescent imaging quality.In addition,the process of generating Cerenkov and penetrating in vitro imaging of 18 F radionuclide inmuscle tissue are also first proposed by GEANT4Monte Carlomethod.The result of the GEANT4 simulation is compared with the numerical solution of the corresponding stochastic partial differential equations,which shows that the stochastic partial differential equation can simulate the corresponding process.展开更多
文摘Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics needs such numerical methods. These methods must be precise and useful in terms of memory capacity and time duration. In this paper, we study basic analytical and finite difference methods to determine the effective refractive index of AIGaAs-GaAs slab waveguide. Also, appropriate effective refractive index value is obtained with respect to number of grid points and number of matrix sizes. Finally, the validity of the obtained values by both methods is compared to using waveguide type.
文摘In this paper, we propose a nearly analytic exponential time difference （NETD） method for solving the 2D acoustic and elastic wave equations. In this method, we use the nearly analytic discrete operator to approximate the high-order spatial differential operators and transform the seismic wave equations into semi-discrete ordinary differential equations （ODEs）. Then, the converted ODE system is solved by the exponential time difference （ETD） method. We investigate the properties of NETD in detail, including the stability condition for 1-D and 2-D cases, the theoretical and relative errors, the numerical dispersion relation for the 2-D acoustic case, and the computational efficiency. In order to further validate the method, we apply it to simulating acoustic/elastic wave propagation in mul- tilayer models which have strong contrasts and complex heterogeneous media, e.g., the SEG model and the Mar- mousi model. From our theoretical analyses and numerical results, the NETD can suppress numerical dispersion effectively by using the displacement and gradient to approximate the high-order spatial derivatives. In addition, because NETD is based on the structure of the Lie group method which preserves the quantitative properties of differential equations, it can achieve more accurate results than the classical methods.
基金National Science Foundation of China(NSFC)(61671009,12171178).
文摘With the development of molecular imaging,Cherenkov optical imaging technology has been widely concerned.Most studies regard the partial boundary flux as a stochastic variable and reconstruct images based on the steadystate diffusion equation.In this paper,time-variable will be considered and the Cherenkov radiation emission process will be regarded as a stochastic process.Based on the original steady-state diffusion equation,we first propose a stochastic partial differential equationmodel.The numerical solution to the stochastic partial differential model is carried out by using the finite element method.When the time resolution is high enough,the numerical solution of the stochastic diffusion equation is better than the numerical solution of the steady-state diffusion equation,which may provide a new way to alleviate the problem of Cherenkov luminescent imaging quality.In addition,the process of generating Cerenkov and penetrating in vitro imaging of 18 F radionuclide inmuscle tissue are also first proposed by GEANT4Monte Carlomethod.The result of the GEANT4 simulation is compared with the numerical solution of the corresponding stochastic partial differential equations,which shows that the stochastic partial differential equation can simulate the corresponding process.
