Finite element approximation of the Fokker-Planck equation for diffuse optical tomography

In diffuse optical tomography, light transport theory is used to describe photon propagation inside turbid medium. A commonly used simplification for the radiative transport equation is the diffusion approximation due to computational feasibility. However, it is known that the diffusion approximation is not valid close to the sources and boundary and in low-scattering regions. Fokker-Planck equation describes light propagation when scattering is forward-peaked. In this article a numerical solution of the Fokker-Planck equation using finite element method is developed. Approach is validated against Monte Carlo simulation and compared with the diffusion approximation. The results show that the Fokker-Planck equation gives equal or better results than the diffusion approximation on the boundary of a homogeneous medium and in turbid medium containing a low-scattering region when scattering is forward-peaked.     

Three-dimensional diffuse optical imaging of hand joints: System description and phantom studies

In this paper we describe a three-dimensional (3D) continuous wave (CW) diffuse optical tomography (DOT) system and present 3D volumetric reconstruction studies using this DOT system with simple phantom models that simulate hand joints. The CCD-based DOT system consists of 64×64 source/detector fiber optic channels, which are arranged in four layers, forming a cylindrical fiber optic/tissue interface. Phantom experiments are used to evaluate system performance with respective to axial spatial resolution, optical contrast and target position for detection of osteoarthritis where cartilage is the primary target region of interest. These phantom studies suggest that we are able to quantitatively resolve a 2 mm thick “cartilage” and qualitatively resolve a 1 mm thick “cartilage” using our 3D reconstruction approach. Our results also show that optical contrast of 3:1–7:1 between the “disease cartilage” and normal cartilage can be quantitatively recovered. Finally, the target position along axial direction on image reconstruction is studied. All the images are obtained using our 3D finite-element-based reconstruction algorithm.     

Measurements of optical parameters of phantom solution and bulk animal tissues in vitro at 650 nm

Optical parameters of bulk animal tissue in vitro, including absorption coefficient (μ_a), reduced scattering coefficient (μ′_s) or scattering coefficient (μ_s), total attenuation coefficient (μ_t), anisotropy factor (g) and refractive index (n) are measured at wavelength of 650 nm. Clinical Intralipid-10% is diluted in distilled water into different concentrations to use as tissue phantoms. Four types of animal tissues in vitro are studied. The relationships among the optical parameters are analyzed systemically. For animal tissues, μ_a, μ′_s or μ_s and n rely on muscle fiber orientations. μ_s and μ_t range from 10 to 20 mm⁻¹, μ_a from 10⁻² to 10⁻³ mm⁻¹ and g from 0.95 to 0.99.     

Comparative study of classifiers to mitigate intersymbol interference in diffuse indoor optical wireless communication links

The maximum data rate that can be achieved in diffuse indoor optical wireless communication (OWC) is limited due to the effect of intersymbol interference (ISI). The adverse effect of ISI on the system performance can be minimised using a channel equaliser at the receiver. In this study, digital signal detection is formulated as a classification problem and hence a classifier is adopted at the receiver. The bit error performance of classifiers with non-linear decision boundary including a multilayer perceptron (MLP), a support vector machine (SVM), the radial basis function (RBF), and the Bayesian classifier is studied along with traditional equaliser and reported here. The MLP offers the best performance; however there is trade-off between the performance and complexity especially at highly diffuse channel.     

A basic study on the sound field analysis of periodic structures by boundary integral equation method

This study deals with the development of the approximate method to analyze the sound field around equally spaced finite obstacles, using the periodic boundary condition. First, on the assumption that the equally spaced finite obstacles are the periodically arranged obstacles, the sound field is analyzed by boundary integral equation method with a Green’s function which satisfies the periodic boundary condition. Furthermore, by comparing these results and the exact solution by using the fundamental solution as Green’s function, the validity of the approximate method is also investigated. Next, in order to evaluate the applicability of the approximate method, the simple formula using some parameters, i.e., the frequency, the period, and the number of obstacles, etc., is proposed. The results of the sound field analysis applied the formula are presented.     

Image reconstruction in diffuse optical tomography using the coupled radiative transport-diffusion model

The coupled radiative transport-diffusion model can be used as light transport model in situations in which the diffusion equation is not a valid approximation everywhere in the domain. In the coupled model, light propagation is modelled with the radiative transport equation in sub-domains in which the approximations of the diffusion equation are not valid, such as within low-scattering regions, and the diffusion approximation is used elsewhere in the domain. In this paper, an image reconstruction method for diffuse optical tomography based on using the coupled radiative transport-diffusion model is developed. In the approach, absorption and scattering distributions are estimated by minimising a regularised least-squares error between the measured data and solution of the coupled model. The approach is tested with simulations. Reconstructions from different cases including domains with low-scattering regions are shown. The results show that the coupled radiative transport-diffusion model can be utilised in image reconstruction problem of diffuse optical tomography and that it produces as good quality reconstructions as the full radiative transport equation also in the presence of low-scattering regions.     

A solution of radiative transfer in isotropic plane-parallel atmosphere by integrating Milne equation

In this paper we solve the inversion problem of the radiative transfer process in the isotropic plane-parallel atmosphere by iterative integrations of the Milne integral equation. As a result, we obtain the scattering function in the form of a cubic polynomial in optical thickness. The author has already solved the same problem by iterative integrations of Chandrasekhar's integral equation. In the Milne integral equation, both the cosines of the viewing angles and the optical thickness are integral variables, while in Chandrasekhar's integral equation the cosines of the viewing angles are variables but the optical thickness is not. We derive several series of exponential-like functions as intermediate derivations. Their convergences are evaluated by the author's previous work in the solution of Chandrasekhar's integral equation. The truncated scattering function up to the third order in optical thickness thus obtained is identical to that obtained from Chandrasekhar's integral equation, though their apparent forms are different. Chandrasekhar pointed out that the solution of Chandrasekhar's integral equation does not have a uniqueness of solution. The Milne equation, in contrast, has been proven to have a unique solution. We discuss the uniqueness of the solution by these two methods.     

The novel boundary integral equation method for solving arbitrary cross-section waveguides

A novel method for solving arbitrary cross-sectionn waveguides is presented. The novel method is a modification of the eigen-weighted boundary integral equation method; the EWBIEM is modified by using the eigenfunction of a fictitious regular boundary as weighting function, whose eigenvalue may be the known value, and meanwhile using the domain-bases. To confirm the validity of the novel method, numerical analysis are presented for circular groove guide as an example.The Project Supported by National Natural Science Foundation of P.R. China     

An efficient and robust reconstruction method for optical tomography with the time-domain radiative transfer equation

An efficient and robust method based on the complex-variable-differentiation method (CVDM) is proposed to reconstruct the distribution of optical parameters in two-dimensional participating media. An upwind-difference discrete-ordinate formulation of the time-domain radiative transfer equation is well established and used as forward model. The regularization term using generalized Gaussian Markov random field model is added in the objective function to overcome the ill-posed nature of the radiative inverse problem. The multi-start conjugate gradient method was utilized to accelerate the convergence speed of the inverse procedure. To obtain an accurate result and avoid the cumbersome formula of adjoint differentiation model, the CVDM was employed to calculate the gradient of objective function with respect to the optical parameters. All the simulation results show that the CVDM is efficient and robust for the reconstruction of optical parameters.     

求解波束波导中激励-传输-辐射问题的积分方程方法

针对波束波导系统结构和馈电特性,应用体面积分方程建立了任意结构波束波导中电磁波的激励、传输和辐射问题所对应的物理模型.在激励端,采用模式展开方法描述输入口面场分布,并利用模式匹配方法精确模拟波束波导馈电端的匹配状况.通过对波束波导内部的实际结构的精确拟合,基于积分方程建立了激励-传输-辐射的电磁模型,并以高效数值解法给...     

Some new exact solutions to the Burgers--Fisher equation and generalized Burgers--Fisher equation

Some new exact solutions of the Burgers--Fisher equation and generalized Burgers--Fisher equation have been obtained by using the first integral method. These solutions include exponential function solutions, singular solitary wave solutions and some more complex solutions whose figures are given in the article. The result shows that the first integral method is one of the most effective approaches to obtain the solutions of the nonlinear partial differential equations.     

A method for multi-frequency calculation of boundary integral equation in acoustics based on series expansion

The paper presents a method to solve the problem of multi-frequency calculation of Helmholtz boundary integral equation in acoustics. Based on series expansion, system matrices are independent of wavenumber and become the matrix power series of wavenumber. As a result, all matrices in the matrix power series are only dependent on the structure geometry. In addition, an element transform method to calculate the singular integral and Cauchy singular integral is also discussed because the singular integral need to be solved using the method. The convergence of the series expansion method is also proved in this paper. The effectiveness of the method is confirmed by two numerical examples.     

Random integral equation formulation of a generalized Langevin equation

In this paper, the generalized Langevin equation introduced by Kubo and Mori is formulated as a random integral equation. We consider (1) the existence and uniqueness of the solution, (2) moments of the solution process, (3) a comparison theorem for solution processes, and (4) the Cauchy polygonal approximation to the solution.     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

半球形电极弥散放电等离子体参量的光谱研究

利用两个半球形水电极介质阻挡放电装置,在空气中实现了弥散放电,采用发射光谱法,对分子振动温度、分子转动温度及电子平均能量等随电压的变化进行了研究。实验利用氮分子第二正带系(C3Πu→B3Πg)的发射谱线得到了氮分子的振动温度;通过氮分子离子(N2+)的第一负带系(B2Σu+→X2Σg+)的发射谱线计算了氮分子的转动温度;采集了氮分子离子391.4nm和激发态的氮分子337.1nm两条发射谱线的相对强度之比,研究了电子的平均能量的变化。结果表明:随外加电压的增加,分子的转动温度呈上升趋势,分子的振动温度呈下降趋势,电子能量呈缓慢的下降趋势。     

A note on the mathematical structure of the optical cubic-quintic Schrödinger equation

The mathematical structure of the optical cubic-quintic Schrödinger equation is investigated in a special way by considering a potential depending upon the modulus of the wave-functions involved. In this context, an associated operator is defined.     

Simultaneous solution of Kompaneets equation and radiative transfer equation in the photon energy range 1-125 keV

Radiative transfer equation in plane parallel geometry and Kompaneets equation is solved simultaneously to obtain theoretical spectrum of 1-125 keV photon energy range. Diffuse radiation field are calculated using time-independent radiative transfer equation in plane parallel geometry, which is developed using discrete space theory (DST) of radiative transfer in a homogeneous medium for different optical depths. We assumed free-free emission and absorption and emission due to electron gas to be operating in the medium. The three terms n, n² and (∂n/∂x_k) where n is photon phase density and x_k=(hν/kT_e), in Kompaneets equation and those due to free-free emission are utilized to calculate the change in the photon phase density in a hot electron gas. Two types of incident radiation are considered: (1) isotropic radiation with the modified black body radiation I^MB[1] and (2) anisotropic radiation which is angle dependent. The emergent radiation at τ=0 and reflected radiation τ=τ_max are calculated by using the diffuse radiation from the medium. The emergent and reflected radiation contain the free-free emission and emission from the hot electron gas. Kompaneets equation gives the changes in photon phase densities in different types of media. Although the initial spectrum is angle dependent, the Kompaneets equation gives a spectrum which is angle independent after several Compton scattering times.     

各向异性海底地层海洋可控源电磁响应三维积分方程法数值模拟

采用积分方程法对各向异性海底地层海洋可控源电磁（MCSEM）响应进行三维数值模拟.首先利用算子理论给出各向异性海底地层中的压缩积分方程,由于满足压缩映射条件,该积分方程在任意参数下总是迭代收敛的.然后为了提高计算效率,引入分区域多重网格准线性近似技术.通过具体算例验证了所述算法在计算精度与效率方面的有效性.最后利用该算法考察并分析了海底地层的各向异性对MCSEM三维响应特征的影响. 关键词： 海洋可控源电磁法 各向异性 三维模拟 积分方程法     

Real-time optical imaging and tracking of micron-sized particles

We report real-time imaging and dynamics monitoring of micrometer predefined and random sized particles by time-space-wavelength mapping technology using a single-detector. Experimentally, we demonstrate real-time line imaging of a 5 μm polystyrene microsphere, glass powder particles and patterns such as fingerprints with up to 5 μm resolution at 1 line/50 ns capture rate. By using the same setup, real-time displacement tracking of micrometer-size glass particles with 50 ns temporal resolution and up to 5 μm spatial resolution is achieved. We also show that existing correlation spectroscopy algorithms can be adopted to extract dynamic information in a complex environment.     

Analysis of diffuse reflectance theories in the MIR region

There exist in the literature a number of theories to estimate the diffuse reflectance values of materials varying in their optical properties. In this article, an analysis of some of the theories has been made with a view to assess the dependence of diffuse reflectance of materials on their optical constants and their grain size (d).