This paper presents an efficient algorithm for reducing RLC power/ground network complexities by exploitation of the regularities in the power/ground networks. The new method first builds the equivalent models for man...This paper presents an efficient algorithm for reducing RLC power/ground network complexities by exploitation of the regularities in the power/ground networks. The new method first builds the equivalent models for many series RLC-current chains based on their Norton's form companion models in the original networks,and then the precondition conjugate gradient based iterative method is used to solve the reduced networks,which are symmetric positive definite. The solutions of the original networks are then back solved from those of the reduced networks.Experimental results show that the complexities of reduced networks are typically significantly smaller than those of the original circuits, which makes the new algorithm extremely fast. For instance, power/ground networks with more than one million branches can be solved in a few minutes on modern Sun workstations.展开更多
To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Pad6 approximation with Baker's algorithm and the program a...To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Pad6 approximation with Baker's algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000, the intensity spectrum obtained by the Padé approximation from a 2^8-item sequence output is more exact than that obtained by fast Fourier transformation from a 2^20-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method, and then the band diagrams are obatined. In addition, mode frequencies and Q-factors are calculated for photonic crystal microcavity.展开更多
We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular...We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.展开更多
The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split...The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split mode wavelengths. We find that the coupled fundamental and first order WG modes can have anticrossing mode coupling as their wavelengths approach the same value in some case,which prevents the cross of the coupled mode wavelengths. The anticrossing mode coupling greatly reduces the coupled mode Q-factor,because the coupled mode field distribution transfers between the fundamental and the first order WG modes.展开更多
The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduct...The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.展开更多
A hybrid method combining finite difference time domain(FDTD)with topology network was presented to treat with electromagnetic couplings and transmissions in large spaces A generalized matrix euqation expressing th...A hybrid method combining finite difference time domain(FDTD)with topology network was presented to treat with electromagnetic couplings and transmissions in large spaces A generalized matrix euqation expressing the relations among wave vectors at every port of the network nodes was give Scattering characteristics and electromagnetic distributions of every node was calculated independently using FDTD A structure of irises in a waveguide was taken as numerical examples This hybrid method has more advantages than the traditional FDTD method which includes saving calculation time,saving memory spaces and being flexible in setting up FDTD grids展开更多
Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and ...Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components (PCs) to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network (ANN) to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody's methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.展开更多
In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid colu...In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.展开更多
With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) meth...With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.展开更多
The mode frequencies and the quality factors for the equilateral triangular resonator (ETR), the square resonator (SR) and the rhombus resonator (RR) are numerically calculated by the finite difference time domai...The mode frequencies and the quality factors for the equilateral triangular resonator (ETR), the square resonator (SR) and the rhombus resonator (RR) are numerically calculated by the finite difference time domain technique and the Padé approximation. The numerical results show that the resonant modes confined in an equilateral triangular cavity have much higher quality factors than those in the square or the rhombus cavities. The modes in the ETR are totally confined in transverse direction while those in the SR and RR are only partly confined. For the ETR with the side length of 4μm and the refractive index of 3 2, the mode quality factor of about 5 5×10 3 at the wavelength of 1 55μm has been obtained.展开更多
A rigorous analysis of a broad wall slot coupler composed of two crossed rectangular waveguides by using FDTD method is presented. Two types of slots, a longitudinal/transverse slot and a centered-inclined slot, are a...A rigorous analysis of a broad wall slot coupler composed of two crossed rectangular waveguides by using FDTD method is presented. Two types of slots, a longitudinal/transverse slot and a centered-inclined slot, are analyzed. Coupling slot characteristics are obtained, including scattering parameters and the resonant length for different waveguide dimensions and frequencies. The numerical results are all in good agreement with those by the method of moments. A coupler system with a short circuit port is also analyzed to help designing the slot array.展开更多
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
The scalar two-dimensional finite difference time domain (FDTD) method is applied to simulate the mode field distribution of TE 0 of the waveguide grating coupler. Computer simulation shows that the same stable mode f...The scalar two-dimensional finite difference time domain (FDTD) method is applied to simulate the mode field distribution of TE 0 of the waveguide grating coupler. Computer simulation shows that the same stable mode field distribution pattern is obtained through the different kinds of driving sources. It is found that the optical field mode is determined by waveguide structure and optical wavelength other than the driving source.According to the mode field distribution, the optimum coupling efficiency can be predicted. Compared with another numerical methods,the CPU-time and memory elements of computer used by FDTD are much less.展开更多
This letter investigates an efficient design procedure integrating the Genetic Algorithm (GA) with the Finite Difference Time Domain (FDTD) for the fast optimal design of Smart Antenna Arrays (SAA). The FDTD is used t...This letter investigates an efficient design procedure integrating the Genetic Algorithm (GA) with the Finite Difference Time Domain (FDTD) for the fast optimal design of Smart Antenna Arrays (SAA). The FDTD is used to analyze SAA with mutual coupling. Then,on the basis of the Maximal Signal to Noise Ratio (MSNR) criteria, the GA is applied to the optimization of weighting elements and structure of SAA. Finally, the effectiveness of the analysis is evaluated by experimental antenna arrays.展开更多
The transport processes of solutes in two soil columns filled with undisturbed soil material collected from an unsaturated sandy aquifer formation in Belgium subjected to a variable upper boundary condition were ident...The transport processes of solutes in two soil columns filled with undisturbed soil material collected from an unsaturated sandy aquifer formation in Belgium subjected to a variable upper boundary condition were identified from breakthrough curves measured by means of time domain reflectometry (TDR). Solute breakthrough was measured with 3 TDR probes inserted into each soil column at three different depths at a 10 minutes time interval. In addition, soil water content and pressure head were measured at 3 different depths. Analytical solute transport models were used to estimate the solute dispersion coefficient and average pore-water velocity from the observed breakthrough curves. The results showed that the analytical solutions were suitable in fitting the observed solute transport. The dispersion coefficient was found to be a function of the soil depth and average pore-water velocity, imposed by the soil water flux. The mobile moisture content on the other hand was not correlated with the average pore-water velocity and the dispersion coefficient.展开更多
Based on couple wave equation and finite-difference time-domain (FDTD) algorithm,the strong couple characteristic of 2-D photonic crystal couplers is calculated.Theoretical analysis and numerical simulated results ind...Based on couple wave equation and finite-difference time-domain (FDTD) algorithm,the strong couple characteristic of 2-D photonic crystal couplers is calculated.Theoretical analysis and numerical simulated results indicate that the energy in a 2-D photonic crystal coupler can not be totally transferred between two wave-guides.Compared with the result of weak coupling theory,our result is more accurate.展开更多
The band structure of 2D photonic crystals (PCs) and localized states resulting from defects are analyzed by finite-difference time-domain (FDTD) technique and Padé approximation.The effect of dielectric constant...The band structure of 2D photonic crystals (PCs) and localized states resulting from defects are analyzed by finite-difference time-domain (FDTD) technique and Padé approximation.The effect of dielectric constant contrast and filling factor on photonic bandgap (PBG) for perfect PCs and localized states in PCs with point defects are investigated.The resonant frequencies and quality factors are calculated for PCs with different defects.The numerical results show that it is possible to modulate the location,width and number of PBGs and frequencies of the localized states only by changing the dielectric constant contrast and filling factor.展开更多
A new singularity extraction technique is presented to calculate accurately the singular integrals in Time Domain Electric Field Integral Equation (TDEFIE).In singularity extraction pro- cedure,through the aid of the ...A new singularity extraction technique is presented to calculate accurately the singular integrals in Time Domain Electric Field Integral Equation (TDEFIE).In singularity extraction pro- cedure,through the aid of the first order Taylor series of time base function including time-retardation,the singularity of the integrand can be removed.The surface current density and backscattered far-field response of a conducting cube illuminated by a Gaussian plane wave is com- puted using the presented technique.Comparisons are made with the results obtained by the Inverse Discrete Fourier Transform (IDFT) of the frequency domain and the results obtained by using Ve- chinski's time averaging technique,which demonstrate that the presented method with this new time domain singularity extraction technique to solve TDEFIE is very accurate and stable.展开更多
文摘This paper presents an efficient algorithm for reducing RLC power/ground network complexities by exploitation of the regularities in the power/ground networks. The new method first builds the equivalent models for many series RLC-current chains based on their Norton's form companion models in the original networks,and then the precondition conjugate gradient based iterative method is used to solve the reduced networks,which are symmetric positive definite. The solutions of the original networks are then back solved from those of the reduced networks.Experimental results show that the complexities of reduced networks are typically significantly smaller than those of the original circuits, which makes the new algorithm extremely fast. For instance, power/ground networks with more than one million branches can be solved in a few minutes on modern Sun workstations.
文摘To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Pad6 approximation with Baker's algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000, the intensity spectrum obtained by the Padé approximation from a 2^8-item sequence output is more exact than that obtained by fast Fourier transformation from a 2^20-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method, and then the band diagrams are obatined. In addition, mode frequencies and Q-factors are calculated for photonic crystal microcavity.
基金supported by the Young Scientists Fund of the National Natural Science Foundation of China(No.41304082)the China Postdoctoral Science Foundation(No.2016M590731)+2 种基金the Young Scientists Fund of the Natural Science Foundation of Hebei Province(No.D2014403011)the Program for Young Excellent Talents of Higher Education Institutions of Hebei Province(No.BJ2016046)the Geological survey project of China Geological Survey(No.1212011121197)
文摘We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.
文摘The mode characteristics for twin coupled microdisks are investigated by finite-difference time-domain technique. In the coupled microdisks,the same order whispering-gallery(WG) modes can form coupled modes with split mode wavelengths. We find that the coupled fundamental and first order WG modes can have anticrossing mode coupling as their wavelengths approach the same value in some case,which prevents the cross of the coupled mode wavelengths. The anticrossing mode coupling greatly reduces the coupled mode Q-factor,because the coupled mode field distribution transfers between the fundamental and the first order WG modes.
文摘The transmission and dispersive characteristics of slotline are calculated in this paper. The tail of Gaussion pulse is improved because a modified dispersive boundary condition (DBC) is adopted. It leads to a reduction in computer memory requirements and computational time. The computational domain is greatly reduced to enable performance in personal computer. At the same time because edges of a boundary and summits are treated well, the computational results is more accurate and more collector.
文摘A hybrid method combining finite difference time domain(FDTD)with topology network was presented to treat with electromagnetic couplings and transmissions in large spaces A generalized matrix euqation expressing the relations among wave vectors at every port of the network nodes was give Scattering characteristics and electromagnetic distributions of every node was calculated independently using FDTD A structure of irises in a waveguide was taken as numerical examples This hybrid method has more advantages than the traditional FDTD method which includes saving calculation time,saving memory spaces and being flexible in setting up FDTD grids
基金supported by the National Natural Science Foundation of China (Grant No. 40974039)High-Tech Research and Development Program of China (Grant No.2006AA06205)Leading Strategic Project of Science and Technology, Chinese Academy of Sciences (XDA08020500)
文摘Traditionally, airborne time-domain electromagnetic (ATEM) data are inverted to derive the earth model by iteration. However, the data are often highly correlated among channels and consequently cause ill-posed and over-determined problems in the inversion. The correlation complicates the mapping relation between the ATEM data and the earth parameters and thus increases the inversion complexity. To obviate this, we adopt principal component analysis to transform ATEM data into orthogonal principal components (PCs) to reduce the correlations and the data dimensionality and simultaneously suppress the unrelated noise. In this paper, we use an artificial neural network (ANN) to approach the PCs mapping relation with the earth model parameters, avoiding the calculation of Jacobian derivatives. The PC-based ANN algorithm is applied to synthetic data for layered models compared with data-based ANN for airborne time-domain electromagnetic inversion. The results demonstrate the PC-based ANN advantages of simpler network structure, less training steps, and better inversion results over data-based ANN, especially for contaminated data. Furthermore, the PC-based ANN algorithm effectiveness is examined by the inversion of the pseudo 2D model and comparison with data-based ANN and Zhody's methods. The results indicate that PC-based ANN inversion can achieve a better agreement with the true model and also proved that PC-based ANN is feasible to invert large ATEM datasets.
基金supported by NSFC(No.41174118)one of the major state S&T special projects(No.2008ZX05020-004)+1 种基金a Postdoctoral Fellowship of China(No.2013M530106)China Scholarship Council(No.2010644006)
文摘In acoustic logging-while-drilling (ALWD) finite difference in time domain (FDTD) simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns (radius -27 mm). Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole (the difference is 〉30 times), the stability and efficiency of the perfectly matched layer (PML) scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain (FDTD) in the ALWD simulation, including field-splitting PML (SPML), multiaxial PML(M- PML), non-splitting PML (NPML), and complex frequency-shifted PML (CFS-PML). The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor (a0) and the scaling factor (β0) depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.
基金The National Natural Science Foundation of China(No.60702027)the Free Research Fund of the National Mobile Communications Research Laboratory of Southeast University (No.2008B07)the National Basic Research Program of China(973 Program)(No.2007CB310603)
文摘With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.
文摘The mode frequencies and the quality factors for the equilateral triangular resonator (ETR), the square resonator (SR) and the rhombus resonator (RR) are numerically calculated by the finite difference time domain technique and the Padé approximation. The numerical results show that the resonant modes confined in an equilateral triangular cavity have much higher quality factors than those in the square or the rhombus cavities. The modes in the ETR are totally confined in transverse direction while those in the SR and RR are only partly confined. For the ETR with the side length of 4μm and the refractive index of 3 2, the mode quality factor of about 5 5×10 3 at the wavelength of 1 55μm has been obtained.
文摘A rigorous analysis of a broad wall slot coupler composed of two crossed rectangular waveguides by using FDTD method is presented. Two types of slots, a longitudinal/transverse slot and a centered-inclined slot, are analyzed. Coupling slot characteristics are obtained, including scattering parameters and the resonant length for different waveguide dimensions and frequencies. The numerical results are all in good agreement with those by the method of moments. A coupler system with a short circuit port is also analyzed to help designing the slot array.
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
文摘The scalar two-dimensional finite difference time domain (FDTD) method is applied to simulate the mode field distribution of TE 0 of the waveguide grating coupler. Computer simulation shows that the same stable mode field distribution pattern is obtained through the different kinds of driving sources. It is found that the optical field mode is determined by waveguide structure and optical wavelength other than the driving source.According to the mode field distribution, the optimum coupling efficiency can be predicted. Compared with another numerical methods,the CPU-time and memory elements of computer used by FDTD are much less.
文摘This letter investigates an efficient design procedure integrating the Genetic Algorithm (GA) with the Finite Difference Time Domain (FDTD) for the fast optimal design of Smart Antenna Arrays (SAA). The FDTD is used to analyze SAA with mutual coupling. Then,on the basis of the Maximal Signal to Noise Ratio (MSNR) criteria, the GA is applied to the optimization of weighting elements and structure of SAA. Finally, the effectiveness of the analysis is evaluated by experimental antenna arrays.
基金Project supported by the European Economic Community Research Program STEP.
文摘The transport processes of solutes in two soil columns filled with undisturbed soil material collected from an unsaturated sandy aquifer formation in Belgium subjected to a variable upper boundary condition were identified from breakthrough curves measured by means of time domain reflectometry (TDR). Solute breakthrough was measured with 3 TDR probes inserted into each soil column at three different depths at a 10 minutes time interval. In addition, soil water content and pressure head were measured at 3 different depths. Analytical solute transport models were used to estimate the solute dispersion coefficient and average pore-water velocity from the observed breakthrough curves. The results showed that the analytical solutions were suitable in fitting the observed solute transport. The dispersion coefficient was found to be a function of the soil depth and average pore-water velocity, imposed by the soil water flux. The mobile moisture content on the other hand was not correlated with the average pore-water velocity and the dispersion coefficient.
基金Supported by the Youth Science and Technology InnovationItem Foundation of Fujian Province(2001J008)
文摘Based on couple wave equation and finite-difference time-domain (FDTD) algorithm,the strong couple characteristic of 2-D photonic crystal couplers is calculated.Theoretical analysis and numerical simulated results indicate that the energy in a 2-D photonic crystal coupler can not be totally transferred between two wave-guides.Compared with the result of weak coupling theory,our result is more accurate.
文摘The band structure of 2D photonic crystals (PCs) and localized states resulting from defects are analyzed by finite-difference time-domain (FDTD) technique and Padé approximation.The effect of dielectric constant contrast and filling factor on photonic bandgap (PBG) for perfect PCs and localized states in PCs with point defects are investigated.The resonant frequencies and quality factors are calculated for PCs with different defects.The numerical results show that it is possible to modulate the location,width and number of PBGs and frequencies of the localized states only by changing the dielectric constant contrast and filling factor.
文摘A new singularity extraction technique is presented to calculate accurately the singular integrals in Time Domain Electric Field Integral Equation (TDEFIE).In singularity extraction pro- cedure,through the aid of the first order Taylor series of time base function including time-retardation,the singularity of the integrand can be removed.The surface current density and backscattered far-field response of a conducting cube illuminated by a Gaussian plane wave is com- puted using the presented technique.Comparisons are made with the results obtained by the Inverse Discrete Fourier Transform (IDFT) of the frequency domain and the results obtained by using Ve- chinski's time averaging technique,which demonstrate that the presented method with this new time domain singularity extraction technique to solve TDEFIE is very accurate and stable.