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.展开更多
In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibi...In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.展开更多
Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investig...Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.展开更多
Detecting internal defects,particularly voids behind linings,is critical for ensuring the structural integrity of aging high-speed rail(HSR)tunnel networks.While ground-penetrating radar(GPR)is widely employed,systema...Detecting internal defects,particularly voids behind linings,is critical for ensuring the structural integrity of aging high-speed rail(HSR)tunnel networks.While ground-penetrating radar(GPR)is widely employed,systematic quantification of performance boundaries for air-coupled(A-CGPR)and ground-coupled(G-CGPR)systems within the complex electromagnetic environment of multilayer reinforced HSR tunnels remains limited.This study establishes physics-based quantitative performance limits for A-CGPR and G-CGPR through rigorously validated GPRMax finite-difference time-domain(FDTD)simulations and comprehensive field validation over a 300 m operational HSR tunnel section.Key performance metrics were quantified as functions of:(a)detection distance(A-CGPR:2.0–4.5 m;G-CGPR:≤0.1 m),(b)antenna frequency(A-CGPR:300 MHz;G-CGPR:400/900 MHz),(c)reinforcement configuration(unreinforced,single-layer,multilayer rebar),and(d)void geometry(axial length:0.1–1.0 m;radial depth:0.1–0.5 m).Key findings demonstrate:a.A-CGPR(300 MHz):Reliably detects axial voids≥0.3 m at distances≤3 m in minimally reinforced(single-layer rebar)linings(field R2=0.89).Performance degrades significantly at distances>3 m(>60%signal attenuation at 4.5 m)or under multilayer rebar interference,causing 25%–40%accuracy loss for voids<0.3 m.Optimal distance:2.0–2.5 m.b.G-CGPR(900 MHz):Achieves<5%size measurement error for axial voids≥0.1 m and radial voids≥0.2 m in unreinforced linings.Resolution degrades under multilayer reinforcement due to severe signal attenuation,increasing axial void detection error to 10%–20%for voids≥0.3 m and constraining radial size measurement.c.Synergistic Framework:A hybrid inspection protocol is proposed,integrating A-CGPR(20 km/h)for rapid large-area screening and targeted G-CGPR(3 km/h)for high-resolution verification of identified anomalies.This framework enhances NDT efficiency while reducing estimated lifecycle inspection costs by 34%compared to G-CGPR alone.This research provides the first physics-derived quantitative detection thresholds for A-CGPR and G-CGPR in multi-rebar HSR tunnels,validated through field-correlated simulations.Future work will focus on multi-frequency antenna arrays and deep learning algorithms to mitigate reinforcement interference.The established performance boundaries and hybrid framework offer critical guidance for optimizing tunnel lining inspection strategies in extensive HSR networks.展开更多
In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a mess...In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.展开更多
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.展开更多
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.展开更多
To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However...To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.展开更多
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.展开更多
Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circ...Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circular nanotube arrays, two photonic band gaps are emerged in the transmission spectra offing-shaped nanotube arrays, the two band gaps and transmission spectra are adjusted by the length, inner radius, intertube spacing and the dielectric constants of the core and embedding medium, and magnitude modification, redshift and blueshift of the resonance modes are observed. A metallic ring-shaped nanotube arrays for subwavelength band-stop filter in the range of visible light can be achieved. To understand its physical origin, field-interference mechanism was suggested by the field distributions. The proposed nanostructures and results may have great potential applications in subwavelength near-field optics.展开更多
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 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.展开更多
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.展开更多
The voltage stability is substantially a dynamic stability, but the primary method which is more mature and engineering practical to analyze the stability of voltage is still static analysis. The time-domain simulatio...The voltage stability is substantially a dynamic stability, but the primary method which is more mature and engineering practical to analyze the stability of voltage is still static analysis. The time-domain simulation is an important measure in research of complex power grid. With the development of full dynamic simulation technology, the research of dynamic voltage stability by using full dynamic simulation program which is based on time-domain simulation can be carried out. This paper uses full dynamic simulation program in dynamic voltage stability research, lays special stress on research in how generator over-excitation limiter functioned and influence in dynamic voltage stability research, and raise 2 methods and steps to figure out dynamic stable voltage in both over-excitation counted and not counted. The simulation results of examples indicate the correctness and effectiveness of these methods, and also fully verify the function and influence of generator over-excitation limiter in full dynamic voltage stability research.展开更多
In this study, the terahertz time-domain spectroscopy (THz-TDS) of crystalline methedrine, which is one of the illegal drugs, is performed using molecular dynamics simulation by the Fourier transform of time derivativ...In this study, the terahertz time-domain spectroscopy (THz-TDS) of crystalline methedrine, which is one of the illegal drugs, is performed using molecular dynamics simulation by the Fourier transform of time derivative auto-correlation functions of the dipole moment. In order to accurately detect the drugs from samples, it is necessary to build a complete database for terahertz spectra under different external conditions from theoretical calculation, which are hardly obtained from the experiments directly. Our results show remarkable consistency with the available experimental data in the frequency range of 10 - 100 cm-1 indicating that the presented method has significant capability to simulate terahertz spectra at various conditions. We investigated the effects of temperature and pressure on THz-TDS by simulating the system at temperature range between 78.4 K and 400 K at pressures up to 100 atm. Results show the spectral features of THz-TDS both in intensity and profile are highly sensitive to the variation of temperature and with a lower magnitude to the variation of pressure. The vanishing, rebuilding and shifting of spectral peaks are due to the complex mechanisms such as the anharmonicity, shifting in the vibration energy levels, formation and destruction of hydrogen-binding and the deformation of the potential energy surface during the environment changing. This improved our understanding for complicated THz-TDS of crystalline methedrine and would be useful for assignment of the practical measurements.展开更多
An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness wav...An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.展开更多
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.展开更多
With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving larg...With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.展开更多
The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) appr...The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.展开更多
基金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 National Natural Science Foundation of China(No.40774056)Program of Excellent Team in Harbin Institute of Technology
文摘In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.
基金financially supported by Lloyd’s Register Foundation(LRF),a UK-registered charity and sole shareholder of Lloyd’s Register Group Ltd.the Youth Innovation Fund of State Key Laboratory of Ocean Engineering(Grant No.GKZD010059-21)
文摘Simulating the coupled motions of multiple bodies in the time domain is a complex problem because of the strong hydrodynamic interactions and coupled effect of various mechanical connectors. In this study, we investigate the hydrodynamic responses of three barges moored side-by-side in a floatover operation in the frequency and time domains. In the frequency domain, the damping lid method is adopted to improve the overestimated hydrodynamic coefficients calculated from conventional potential flow theory. A time-domain computing program based on potential flow theory and impulse theory is compiled for analyses that consider multibody hydrodynamic interactions and mechanical effects from lines and fenders. Correspondingly, an experiment is carried out for comparison with the numerical results. All statistics, time series, and power density spectra from decay and irregular wave tests are in a fairly good agreement.
基金funded by the Key Project of Science&Technology Research ofChina Academy of Railway Sciences,grant number 2023YJ022.
文摘Detecting internal defects,particularly voids behind linings,is critical for ensuring the structural integrity of aging high-speed rail(HSR)tunnel networks.While ground-penetrating radar(GPR)is widely employed,systematic quantification of performance boundaries for air-coupled(A-CGPR)and ground-coupled(G-CGPR)systems within the complex electromagnetic environment of multilayer reinforced HSR tunnels remains limited.This study establishes physics-based quantitative performance limits for A-CGPR and G-CGPR through rigorously validated GPRMax finite-difference time-domain(FDTD)simulations and comprehensive field validation over a 300 m operational HSR tunnel section.Key performance metrics were quantified as functions of:(a)detection distance(A-CGPR:2.0–4.5 m;G-CGPR:≤0.1 m),(b)antenna frequency(A-CGPR:300 MHz;G-CGPR:400/900 MHz),(c)reinforcement configuration(unreinforced,single-layer,multilayer rebar),and(d)void geometry(axial length:0.1–1.0 m;radial depth:0.1–0.5 m).Key findings demonstrate:a.A-CGPR(300 MHz):Reliably detects axial voids≥0.3 m at distances≤3 m in minimally reinforced(single-layer rebar)linings(field R2=0.89).Performance degrades significantly at distances>3 m(>60%signal attenuation at 4.5 m)or under multilayer rebar interference,causing 25%–40%accuracy loss for voids<0.3 m.Optimal distance:2.0–2.5 m.b.G-CGPR(900 MHz):Achieves<5%size measurement error for axial voids≥0.1 m and radial voids≥0.2 m in unreinforced linings.Resolution degrades under multilayer reinforcement due to severe signal attenuation,increasing axial void detection error to 10%–20%for voids≥0.3 m and constraining radial size measurement.c.Synergistic Framework:A hybrid inspection protocol is proposed,integrating A-CGPR(20 km/h)for rapid large-area screening and targeted G-CGPR(3 km/h)for high-resolution verification of identified anomalies.This framework enhances NDT efficiency while reducing estimated lifecycle inspection costs by 34%compared to G-CGPR alone.This research provides the first physics-derived quantitative detection thresholds for A-CGPR and G-CGPR in multi-rebar HSR tunnels,validated through field-correlated simulations.Future work will focus on multi-frequency antenna arrays and deep learning algorithms to mitigate reinforcement interference.The established performance boundaries and hybrid framework offer critical guidance for optimizing tunnel lining inspection strategies in extensive HSR networks.
基金Project supported by the National Defense Laboratory Foundation (Grant No.51444020103QT0601)the Shanghai Leading Academic Discipline Project (Grant No.T0102)
文摘In this work, we treat scattering objects, water, surface and bottom in a truly unified manner in a parallel finitedifference time-domain （FDTD） scheme, which is suitable for distributed parallel computing in a message passing interface （MPI） programming environment. The algorithm is implemented on a cluster-based high performance computer system. Parallel computation is performed with different division methods in 2D and 3D situations. Based on analysis of main factors affecting the speedup rate and parallel efficiency, data communication is reduced by selecting a suitable scheme of task division. A desirable scheme is recommended, giving a higher speedup rate and better efficiency. The results indicate that the unified parallel FDTD algorithm provides a solution to the numerical computation of acoustic scattering.
基金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.
基金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 Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+1 种基金the Funds of Creative Research Groups of China(No.41821002)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003).
文摘To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.
基金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.
基金Projects(11164007,61275174)supported by the National Natural Science Foundation of ChinaProject(20100162110068)supported by the Doctoral Program of Higher Education of China+1 种基金Project(20132BAB212007)supported by the Natural Science Foundation of Jiangxi Province,ChinaProject(GJJ11107)supported by Scientific Foundation of Jiangxi Education Department,China
文摘Metallic ring-shaped nanotube arrays are proposed and its optical transmission properties are studied by using finite-difference time-domain （FDTD） method. Compared with the transmission spectra of conventional circular nanotube arrays, two photonic band gaps are emerged in the transmission spectra offing-shaped nanotube arrays, the two band gaps and transmission spectra are adjusted by the length, inner radius, intertube spacing and the dielectric constants of the core and embedding medium, and magnitude modification, redshift and blueshift of the resonance modes are observed. A metallic ring-shaped nanotube arrays for subwavelength band-stop filter in the range of visible light can be achieved. To understand its physical origin, field-interference mechanism was suggested by the field distributions. The proposed nanostructures and results may have great potential applications in subwavelength near-field optics.
基金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.
基金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.
基金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 voltage stability is substantially a dynamic stability, but the primary method which is more mature and engineering practical to analyze the stability of voltage is still static analysis. The time-domain simulation is an important measure in research of complex power grid. With the development of full dynamic simulation technology, the research of dynamic voltage stability by using full dynamic simulation program which is based on time-domain simulation can be carried out. This paper uses full dynamic simulation program in dynamic voltage stability research, lays special stress on research in how generator over-excitation limiter functioned and influence in dynamic voltage stability research, and raise 2 methods and steps to figure out dynamic stable voltage in both over-excitation counted and not counted. The simulation results of examples indicate the correctness and effectiveness of these methods, and also fully verify the function and influence of generator over-excitation limiter in full dynamic voltage stability research.
文摘In this study, the terahertz time-domain spectroscopy (THz-TDS) of crystalline methedrine, which is one of the illegal drugs, is performed using molecular dynamics simulation by the Fourier transform of time derivative auto-correlation functions of the dipole moment. In order to accurately detect the drugs from samples, it is necessary to build a complete database for terahertz spectra under different external conditions from theoretical calculation, which are hardly obtained from the experiments directly. Our results show remarkable consistency with the available experimental data in the frequency range of 10 - 100 cm-1 indicating that the presented method has significant capability to simulate terahertz spectra at various conditions. We investigated the effects of temperature and pressure on THz-TDS by simulating the system at temperature range between 78.4 K and 400 K at pressures up to 100 atm. Results show the spectral features of THz-TDS both in intensity and profile are highly sensitive to the variation of temperature and with a lower magnitude to the variation of pressure. The vanishing, rebuilding and shifting of spectral peaks are due to the complex mechanisms such as the anharmonicity, shifting in the vibration energy levels, formation and destruction of hydrogen-binding and the deformation of the potential energy surface during the environment changing. This improved our understanding for complicated THz-TDS of crystalline methedrine and would be useful for assignment of the practical measurements.
基金Supported by the National Natural Science Fbundation of China(No.69931030)
文摘An Improved Locally Conformal Finite-Difference Time-Domain (ILC-FDTD) method is presented in this paper, which is used to analyze the edge inclined slots penetrating adjacent broadwalls of a finite wall thickness waveguide. ILC-FDTD not only removes tile instability of the original locally conformal FDTD algorithm, but also improves the computational accuracy by locally modifying magnetic field update equations and the virtual iterative electric fields accordlng to the complexity of tile slot fringe fields. The mutual coupling between two edge inclined slots can also be analyzed by ILC-FDTD effectively.
基金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(2025ZD1009704)supported by the National Science and Technology Major Project of ChinaProjects(2023JJ30659,2022JJ30706)supported by Hunan Provincial Natural Science Foundation,China。
文摘With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.
基金Project(41174061) supported by the National Natural Science Foundation of ChinaProject(2011QNZT011) supported by the Free Exploration Program of Central South University,China
文摘The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.
