The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic m...3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.展开更多
This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.展开更多
Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open...Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE met...To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.展开更多
A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spheric...A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.展开更多
The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order sy...The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.展开更多
A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect ...A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the...Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.展开更多
Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly mat...Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.展开更多
The accurate identification of the oil-paper insulation state of a transformer is crucial for most maintenance strategies.This paper presents a multi-feature comprehensive evaluation model based on combination weighti...The accurate identification of the oil-paper insulation state of a transformer is crucial for most maintenance strategies.This paper presents a multi-feature comprehensive evaluation model based on combination weighting and an improved technique for order of preference by similarity to ideal solution(TOPSIS)method to perform an objective and scientific evaluation of the transformer oil-paper insulation state.Firstly,multiple aging features are extracted from the recovery voltage polarization spectrum and the extended Debye equivalent circuit owing to the limitations of using a single feature for evaluation.A standard evaluation index system is then established by using the collected time-domain dielectric spectrum data.Secondly,this study implements the per-unit value concept to integrate the dimension of the index matrix and calculates the objective weight by using the random forest algorithm.Furthermore,it combines the weighting model to overcome the drawbacks of the single weighting method by using the indicators and considering the subjective experience of experts and the random forest algorithm.Lastly,the enhanced TOPSIS approach is used to determine the insulation quality of an oil-paper transformer.A verification example demonstrates that the evaluation model developed in this study can efficiently and accurately diagnose the insulation status of transformers.Essentially,this study presents a novel approach for the assessment of transformer oil-paper insulation.展开更多
In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing ...In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing the periodic metallic nanoparticle arrays. The dispersion of the metal, caused by the evanescent wave propagating along the metal-dielectric interface, is expressed by the Drude model and solved with a generalized auxiliary differential equation（ADE） technique.In the sub-gridding scheme, the ADE–FDTD is applied to the global coarse grids while the ADE–LOD–FDTD is applied to the local fine grids. The time step sizes in the fine-grid region and coarse-grid region can be synchronized, and thus obviating the temporal interpolation of the fields in the time-marching process. Numerical examples about extraordinary optical transmission through the periodic metallic nanoparticle array are provided to show the accuracy and efficiency of the proposed method.展开更多
The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） b...The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method.展开更多
In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is i...In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.展开更多
Nonlinear wave loads can induce low-frequency and high-frequency resonance motions of a moored platform in deep water. For the analysis of the nonlinear response of an offshore platform under the action of irregular w...Nonlinear wave loads can induce low-frequency and high-frequency resonance motions of a moored platform in deep water. For the analysis of the nonlinear response of an offshore platform under the action of irregular waves, the most widely used method in practice is the Cummins method, in which the second-order exciting forces in the time domain are computed by a two-term Volterra series model based on incident waves, first-order body motion response, and quadratic transfer functions(QTFs). QTFs are bichromatic waves acting on a body and are computed in the frequency domain in advance. For moving bodies, QTFs are related to the first-order body response, which is to be determined in the simulation process of body motion response but is unknown in the computation procedure of QTFs. In solving this problem, Teng and Cong(2017) proposed a method to divide the QTFs into different components,which are unrelated to the body response. With the application of the new QTF components, a modified Cummins method can be developed for the simulation of the nonlinear response of a moored floating platform. This paper presents a review of the theory.展开更多
The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulat...The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,展开更多
This paper studies the propagating characteristics of the electromagnetic waves through the coupled-resonator optical waveguides based on the two-dimensional square-lattice photonic crystals by the finite-difference t...This paper studies the propagating characteristics of the electromagnetic waves through the coupled-resonator optical waveguides based on the two-dimensional square-lattice photonic crystals by the finite-difference time-domain method. When the traditional circular rods adjacent to the centre of the cavities are replaced by the oval rods, the simulated results show that the waveguide mode region can be adjusted only by the alteration of the oval rods＇ obliquity. When the obliquity of the oval rods around one cavity is different from the obliquity of that around the adjacent cavities, the group velocities of the waveguide modes can be greatly reduced and the information of different frequencies can be shared and chosen at the same time by the waveguide branches with different structures. If the obliquities of the oval rods around two adjacent cavities are equal and they alternate between two values, the group velocities can be further reduced and a maximum value of 0.0008c （c is the light velocity in vacuum） can be acquired.展开更多
The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matchin...The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.展开更多
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
文摘The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
基金The authors thank the funds supported by the China National Nuclear Corporation under Grants Nos.WUQNYC2101 and WUHTLM2101-04National Natural Science Foundation of China(42074132,42274154).
文摘3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
文摘This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.
基金supported by the National Natural Science Foundation of China(Grant No.61231003)
文摘Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
基金Project supported by the National Natural Science Foundation of China(Grant No.51025622)
文摘To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61301063 and 41305017)
文摘A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.
基金supported by the National Natural Science Foundation of China(Grant Nos.60931002 and 61101064)the Universities Natural Science Foundation of Anhui Province,China(Grant Nos.KJ2011A002 and 1108085J01)
文摘The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.
基金National defense technical basic research project,Terahertz detection technology and application research on ceramic matrix composites(JSZL2015411C002)
文摘A method for extracting optical parameters of plastics materials based on terahertz time domain spectroscopy is presented. The transmission-type Terahertz Time-Domain Spectroscopy(THz TDS) system is adopted to detect the refractive index and extinction coefficient on different plastic materials. Then the corresponding spectral information is obtained by Fourier transform of the terahertz time domain waveform of the sampling points, including the corresponding amplitude and phase information of the waveform. The optical parameter extraction model is built. By using the simplex optimization method, the curves of the refractive index and extinction coefficient for the plastic material are obtained. The experimental samples are made of different plastic parallel plate materials. The experimental results show that the optimization of optical parameters can improve their extraction accuracy, and the error of refractive index is ±0.005. Extraction technology with the simplex optimization method of optical parameter based on THz TDS can help to extract the optical parameters of engineering plastics. It is of great significance for the research of terahertz nondestructive testing.
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金Project supported by the National Key Research and Development Program of China(Grant No.2023YFA1407000)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB0460000)+4 种基金the National Natural Science Foundation of China(Grant Nos.12322401,12127807,and 12393832)CAS Key Research Program of Frontier Sciences(Grant No.ZDBS-LY-SLH004)Beijing Nova Program(Grant No.20230484301)Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023125)CAS Project for Young Scientists in Basic Research(Grant No.YSBR-026)。
文摘Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.
基金Project supported by the National Natural Science Foundation of China (Grant No 60571058)the Specialized Research Fund for the Doctoral Program of Higher Education,China (Grant No 20070701010)
文摘Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.
基金supported by the Natural Science Foundation of the Fujian Province(2021J01109).
文摘The accurate identification of the oil-paper insulation state of a transformer is crucial for most maintenance strategies.This paper presents a multi-feature comprehensive evaluation model based on combination weighting and an improved technique for order of preference by similarity to ideal solution(TOPSIS)method to perform an objective and scientific evaluation of the transformer oil-paper insulation state.Firstly,multiple aging features are extracted from the recovery voltage polarization spectrum and the extended Debye equivalent circuit owing to the limitations of using a single feature for evaluation.A standard evaluation index system is then established by using the collected time-domain dielectric spectrum data.Secondly,this study implements the per-unit value concept to integrate the dimension of the index matrix and calculates the objective weight by using the random forest algorithm.Furthermore,it combines the weighting model to overcome the drawbacks of the single weighting method by using the indicators and considering the subjective experience of experts and the random forest algorithm.Lastly,the enhanced TOPSIS approach is used to determine the insulation quality of an oil-paper transformer.A verification example demonstrates that the evaluation model developed in this study can efficiently and accurately diagnose the insulation status of transformers.Essentially,this study presents a novel approach for the assessment of transformer oil-paper insulation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61471105 and 61331007)
文摘In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing the periodic metallic nanoparticle arrays. The dispersion of the metal, caused by the evanescent wave propagating along the metal-dielectric interface, is expressed by the Drude model and solved with a generalized auxiliary differential equation（ADE） technique.In the sub-gridding scheme, the ADE–FDTD is applied to the global coarse grids while the ADE–LOD–FDTD is applied to the local fine grids. The time step sizes in the fine-grid region and coarse-grid region can be synchronized, and thus obviating the temporal interpolation of the fields in the time-marching process. Numerical examples about extraordinary optical transmission through the periodic metallic nanoparticle array are provided to show the accuracy and efficiency of the proposed method.
基金National Natural Science Foundation of China (No. 60471002) and the Natural Science Foundation ofJiangxi Province (No. 0412014)
文摘The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education(No.20010614003)
文摘In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.
基金the National Key R&D Program of China (Grant No.2016YFE0200100)the National Natural Science Foundation of China (Grant Nos.51490672 and 51479026).
文摘Nonlinear wave loads can induce low-frequency and high-frequency resonance motions of a moored platform in deep water. For the analysis of the nonlinear response of an offshore platform under the action of irregular waves, the most widely used method in practice is the Cummins method, in which the second-order exciting forces in the time domain are computed by a two-term Volterra series model based on incident waves, first-order body motion response, and quadratic transfer functions(QTFs). QTFs are bichromatic waves acting on a body and are computed in the frequency domain in advance. For moving bodies, QTFs are related to the first-order body response, which is to be determined in the simulation process of body motion response but is unknown in the computation procedure of QTFs. In solving this problem, Teng and Cong(2017) proposed a method to divide the QTFs into different components,which are unrelated to the body response. With the application of the new QTF components, a modified Cummins method can be developed for the simulation of the nonlinear response of a moored floating platform. This paper presents a review of the theory.
文摘The characteristics of a cylindrical conformal microstrip patch antenna are analyzed by using the characteristic-based time domain （CBTD） method. A governing equation in the cylindrical coordinate system is formulated directly to facilitate the analysis of cylindrically conformal microstrip patch antennas. The algorithm has second-order accuracy both in time and space domain and has the potential to eliminate the spurious wave reflection from the numerical boundaries of the computational domain, Numerical results demonstrate the important merits and accuracy of the proposed technique in computational electromagnetics,
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 10904176 and 11004169)the Research Foundation of the State Ethnic Affairs Commission of People’s Republic of China (Grant Nos. 10ZY05 and 09ZY012)the "985 Project" and "211 Project" of the Ministry of Education of China
文摘This paper studies the propagating characteristics of the electromagnetic waves through the coupled-resonator optical waveguides based on the two-dimensional square-lattice photonic crystals by the finite-difference time-domain method. When the traditional circular rods adjacent to the centre of the cavities are replaced by the oval rods, the simulated results show that the waveguide mode region can be adjusted only by the alteration of the oval rods＇ obliquity. When the obliquity of the oval rods around one cavity is different from the obliquity of that around the adjacent cavities, the group velocities of the waveguide modes can be greatly reduced and the information of different frequencies can be shared and chosen at the same time by the waveguide branches with different structures. If the obliquities of the oval rods around two adjacent cavities are equal and they alternate between two values, the group velocities can be further reduced and a maximum value of 0.0008c （c is the light velocity in vacuum） can be acquired.
基金Project supported by the National Basic Research Program of China(Grant No.2014CB340203)the National Natural Science Foundation of China(Grant Nos.61431010 and 61501350)the Natural Science Foundation of Shaanxi Province,China(Grant Nos.2018JM6016 and 2016JM1001)
文摘The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.
