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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
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,展开更多
In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric ...In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric material, is simulated numerically using the finite-difference time-domain （FDTD） algorithm. A perfectly matched layer （PML） absorbing technique is used in this simulation. The reflection and transmission coefficients of electromagnetic （EM） waves through PPCs are calculated. The characteristics of the photonic band gap （PBG） are discussed in terms of plasma density, dielectric constant ratios, number of periods, and introduced layer defect. These may provide some useful information for designing plasma photonic crystal devices.展开更多
In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration...In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration is similar to the traditional Kretschmann structure, but the metal film is pattered into a grating. As a TM-polarized laser beam impinges from the substrate, a resonant dip point in reflectance curve is produced at a certain incident angle. Our studies indicate that the sensing sensitivity and resolution are affected by the grating＇s gap and period, and after these parameters havelaeen optimized, a sensing sensitivity of 51.484°/RIU is obtained with a slightly changing resolution.展开更多
Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed...Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed to quantify the electric-field distribution within the subsurface in the presence of such a defect under normal incidence irradiation. Both impurities （inside the crack） and the chemical etching have been investigated. The results show that the maximum electric field amplitude │E│max is 9.57374 V/m when the relative dielectric constant of transparent impurity equals 8.5. And the near-field modulation will be improved if the crack is filled with the remainder polishing powders or water vapor/drops. Meanwhile, the laser-induced initial damage moves to the glass-air surface. In the etched section, the magnitude of intensification is strongly dependent on the inclination angle θ. There will be a highest modulation when θ is around π/6, and the maximum value of IEImax is 18.57314 V/m. When θ ranges from π/θ to π/4, the light intensity enhancement factor can easily be larger than 100, and the modulation follows a decreasing trend. On the other hand, the modulation curves become smooth when θ〉π/4 or θ 〈 π/8.展开更多
基金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.
基金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.
基金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.
基金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 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 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.
基金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.
基金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.
基金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,
基金supported by the Program for New Century Excellent Talents in University(No.NCET-05-0575)the Education Science Foundation of Jiangxi Province(No.Z-03510)
文摘In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric material, is simulated numerically using the finite-difference time-domain （FDTD） algorithm. A perfectly matched layer （PML） absorbing technique is used in this simulation. The reflection and transmission coefficients of electromagnetic （EM） waves through PPCs are calculated. The characteristics of the photonic band gap （PBG） are discussed in terms of plasma density, dielectric constant ratios, number of periods, and introduced layer defect. These may provide some useful information for designing plasma photonic crystal devices.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61137005 and 61178067)the Science Foundation of Shanxi Province,China(Grant No.2013021004-3/2014021021-1)+1 种基金the Pre-studied Project on Weapon Equipment,China(Grant No.201262401090404)the Specialized Research Foundation for Doctor of School,China(Grant No.20122027)
文摘In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration is similar to the traditional Kretschmann structure, but the metal film is pattered into a grating. As a TM-polarized laser beam impinges from the substrate, a resonant dip point in reflectance curve is produced at a certain incident angle. Our studies indicate that the sensing sensitivity and resolution are affected by the grating＇s gap and period, and after these parameters havelaeen optimized, a sensing sensitivity of 51.484°/RIU is obtained with a slightly changing resolution.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No.10904008)the Joint Fund of the National Natural Science Foundation of China and the China Academy of Engineering Physics (GrantNo. 11076008)the Scientific Research Foundation for the Central Universities of China (Grant Nos.ZYGX2009X007,ZYGX2010J045, and ZYGX2011J043)
文摘Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed to quantify the electric-field distribution within the subsurface in the presence of such a defect under normal incidence irradiation. Both impurities （inside the crack） and the chemical etching have been investigated. The results show that the maximum electric field amplitude │E│max is 9.57374 V/m when the relative dielectric constant of transparent impurity equals 8.5. And the near-field modulation will be improved if the crack is filled with the remainder polishing powders or water vapor/drops. Meanwhile, the laser-induced initial damage moves to the glass-air surface. In the etched section, the magnitude of intensification is strongly dependent on the inclination angle θ. There will be a highest modulation when θ is around π/6, and the maximum value of IEImax is 18.57314 V/m. When θ ranges from π/θ to π/4, the light intensity enhancement factor can easily be larger than 100, and the modulation follows a decreasing trend. On the other hand, the modulation curves become smooth when θ〉π/4 or θ 〈 π/8.
