Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional fin...Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.展开更多
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.展开更多
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.展开更多
The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy...The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.展开更多
The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces o...The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.展开更多
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.展开更多
In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this...In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.展开更多
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.展开更多
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.展开更多
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 Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the M...The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.展开更多
In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
Mechanical transplanting has been applied to rice cultivation to save labor costs and ease labor shortages in Asian countries, especially in China. However, little information is available related to the characteristi...Mechanical transplanting has been applied to rice cultivation to save labor costs and ease labor shortages in Asian countries, especially in China. However, little information is available related to the characteristics of agronomic performance when comparing inter-sub-specific hybrid rice(IHR) and inbred japonica rice(IJR) under mechanical transplanting method. In 2013 and 2014, field experiments were conducted using IHR(Yongyou 2640) and IJR(Wuyunjing 24) under two cultivation patterns, that is, pot seedlings mechanically transplanted(PS) and carpet seedlings mechanically transplanted(CS). Grain yield, yield components, leaf area index(LAI), leaf area duration(LAD), aboveground biomass, crop growth rate(CGR), nitrogen(N) uptake, and N accumulation were investigated. When compared with CS, PS displayed significantly increased grain yield for both varieties because the larger sink size allowed higher N accumulation from panicle initiation to maturity. Moreover, total aboveground biomass under PS increased significantly compared with that under CS; that is, higher photosynthetic productivity resulted from a greater LAI and higher LAD during the grain filling stage. Higher N absorption capacity in the middle and late growth periods resulted in significantly enhanced total N uptake under PS. When compared with IJR for both treatments, IHR generated 75.2% more grain yield. However, the characteristics creating high yield of IHR were different from those of IJR. Greater aboveground biomass production as well as higher N uptake and accumulation created higher grain yield in IHR than in IJR. These results suggest higher yield could be achieved using PS with IHR, attributing to exploit both yield superiority and productive potential.展开更多
The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in th...The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in this paper will be solved by Finite-Difference Mode, that＇ s simple and has an extensive use. The way it works is that by dividing the component into many units. Finite-difference methods （FDM） are numerical methods for anoroximating, the solutions to differential eauations usine finite difference equations to approximate derivatives.展开更多
The inter-sequence analysis method is a popular decision making method considering subjective information of the decision maker.The validity of this method is sensitive to the subjective information provided by the de...The inter-sequence analysis method is a popular decision making method considering subjective information of the decision maker.The validity of this method is sensitive to the subjective information provided by the decision maker significantly.Thus,a new weighting method based on the inter-sequence analysis method is proposed.This new method can appropriately reduce the participation of the subjective judgment of the decision maker,and then accordingly decrease the influence of the subjective judgment to the results.It only needs the decision maker to give the inter sequence of the criteria.The importance comparison of two criteria is objectively determined by considering the ratio of a criterion being superior to another one.Some numerical examples are given to illustrate the calculation process of the new method.展开更多
基金supported by National Key Research and Development Program of China(No.2022YFF0800602)Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(No.2022B1212010002)Shenzhen Science and Technology Program(No.KQTD20170810111725321).
文摘Global acoustic simulations are significant in revealing the internal and physical structure of the Earth.However,due to the limited flexibility of grids and the difficulties in handling boundaries,the traditional finite-difference method(FDM)is usually less used in global simulations.Nevertheless,these issues can be well resolved by employing a multi-block structured grid to discretize circular regions.In this paper,we propose an O-H grid approach to partition the circular region and utilize the curvilinear grid finite-difference method(CGFDM)to solve the acoustic wave equation within this circular domain.By appropriately stretching the grid,the interconnections between each grid block are sufficiently smooth for stable information exchange.To verify the efficacy of this method,we conducted three numerical experiments,by comparing results with alternative approaches.Our test results demonstrate good agreement between our findings and the reference solutions.Since the proposed algorithm can effectively solve wave propagation problems in circular regions,it can contribute to 2D global simulation,particularly in interpreting the Earth’s interior.
基金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.
基金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.
文摘The problem for calculating EM energy absorption by lossy dielectric scatterer ir-radiated by plane wave are discussed.The factors affecting the accuracy of computation arediscussed.The calculated results of EM energy absorption and its distribution in homogeneousand layered homogenous lossy dielectric spheres are presented,and a comparison of these resultswith analytical solution is given.The calculation is carried out for dielectric cylinder on conduct-ing ground as well,and the results are compared with the image theory.All the computationsshew that the finite-difference time-domain method can give satisfactory results.
基金Taif University Researchers Supporting Project Number(TURSP-2020/164),Taif University,Taif,Saudi Arabia.
文摘The present article deals with the investigation thermal stress of a magnetothermoelastic cylinder subjected to rotation,open or closed circuit,thermal and mechanical boundary conditions.The outer and inner surfaces of the cylinder are subjected to both mechanical and thermal boundary conditions.A The transient coupled thermoelasticity in an infinite cylinder with its base abruptly exposed to a heat flux of a decaying exponential function of time is devised solve by the finite-difference method.The fundamental equations’system is solved by utilizing an implicit finite-difference method.This current method is a second-order accurate in time and space;it is also unconditionally stable.To illustrate the present model’s efficiency,we consider a suitable material and acquire the numerical solution of temperature,displacement components,and the components of stresses with time t and through the radial of an infinite cylinder.The results indicate that the effect of coupled thermoelasticity,magnetic field,and rotation on the temperature,stresses,and displacement is quite pronounced.In order to illustrate and verify the analytical developments,the numerical solution of partial differential equations,stress components,displacement components and temperature is carried out and computer simulated results are presented graphically.This study is helpful in the development of piezoelectric devices.
基金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.
基金partially supported by China National Major Science and Technology Project (Subproject No:2011ZX05024-001-03)
文摘In seismic exploration, it is common practice to separate the P-wavefield from the S-wavefield by the elastic wavefield decomposition technique, for imaging purposes. However, it is sometimes difficult to achieve this, especially when the velocity field is complex. A useful approach in multi-component analysis and modeling is to directly solve the elastic wave equations for the pure P- or S-wavefields, referred as the separate elastic wave equa- tions. In this study, we compare two kinds of such wave equations： the first-order （velocity-stress） and the second- order （displacement-stress） separate elastic wave equa- tions, with the first-order （velocity-stress） and the second- order （displacement-stress） full （or mixed） elastic wave equations using a high-order staggered grid finite-differ- ence method. Comparisons are given of wavefield snap- shots, common-source gather seismic sections, and individual synthetic seismogram. The simulation tests show that equivalent results can be obtained, regardless of whether the first-order or second-order separate elastic wave equations are used for obtaining the pure P- or S-wavefield. The stacked pure P- and S-wavefields are equal to the mixed wave fields calculated using the corre- sponding first-order or second-order full elastic wave equations. These mixed equations are computationallyslightly less expensive than solving the separate equations. The attraction of the separate equations is that they achieve separated P- and S-wavefields which can be used to test the efficacy of wave decomposition procedures in multi-com- ponent processing. The second-order separate elastic wave equations are a good choice because they offer information on the pure P-wave or S-wave displacements.
基金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.
基金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.
基金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.
基金supported by the Natural Science Foundation of China(No.41374078)Geological Survey Projects of Ministry of Land and Resources of China(No.12120113086100 and 12120113101300)
文摘The Z-Axis tiPPer eiectromagnetic （ZTEM） technique is based on a frequency-domain airbome electromagnetic system that measures the natural magnetic field. A survey area was divided into several blocks by using the Maxwell＇s equations, and the magnetic components at the center of each edge of the grid cell are evaluated by applying the staggered-grid finite-difference method. The tipper and its divergence are derived to complete the 3D ZTEM forward modeling algorithm. A synthetic model is then used to compare the responses with those of 2D finite-element forward modeling to verify the accuracy of the algorithm. ZTEM offers high horizontal resolution to both simple and complex distributions of conductivity. This work is the theoretical foundation for the interpretation of ZTEM data and the study of 3D ZTEM inversion.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金the National Key Research Program of China(2016YFD0300503)the Special Fund for Agro-scientific Research in the Public Interest,China(201303102)+2 种基金the Key Research Program of Jiangsu Province,China(BE2016344)the Major Independent Innovation Project in Jiangsu Province,China(CX(15)1002)a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions,China
文摘Mechanical transplanting has been applied to rice cultivation to save labor costs and ease labor shortages in Asian countries, especially in China. However, little information is available related to the characteristics of agronomic performance when comparing inter-sub-specific hybrid rice(IHR) and inbred japonica rice(IJR) under mechanical transplanting method. In 2013 and 2014, field experiments were conducted using IHR(Yongyou 2640) and IJR(Wuyunjing 24) under two cultivation patterns, that is, pot seedlings mechanically transplanted(PS) and carpet seedlings mechanically transplanted(CS). Grain yield, yield components, leaf area index(LAI), leaf area duration(LAD), aboveground biomass, crop growth rate(CGR), nitrogen(N) uptake, and N accumulation were investigated. When compared with CS, PS displayed significantly increased grain yield for both varieties because the larger sink size allowed higher N accumulation from panicle initiation to maturity. Moreover, total aboveground biomass under PS increased significantly compared with that under CS; that is, higher photosynthetic productivity resulted from a greater LAI and higher LAD during the grain filling stage. Higher N absorption capacity in the middle and late growth periods resulted in significantly enhanced total N uptake under PS. When compared with IJR for both treatments, IHR generated 75.2% more grain yield. However, the characteristics creating high yield of IHR were different from those of IJR. Greater aboveground biomass production as well as higher N uptake and accumulation created higher grain yield in IHR than in IJR. These results suggest higher yield could be achieved using PS with IHR, attributing to exploit both yield superiority and productive potential.
文摘The use of columns on elastic foundation is very common in Civil Engineering, like bridge pier, the foundation of the buildings etc. So, it will be useful to find the critical load for the structure, the problem in this paper will be solved by Finite-Difference Mode, that＇ s simple and has an extensive use. The way it works is that by dividing the component into many units. Finite-difference methods （FDM） are numerical methods for anoroximating, the solutions to differential eauations usine finite difference equations to approximate derivatives.
基金Supported by the National Natural Science Foundation of China(71961018,71671031,71901079,71473033,41661116)the Ministry of Education Humanities and Social Sciences Research Youth Fund Project(17YJC630067)+1 种基金the Youth Foundation Project for Humanities and Social Sciences Research in Colleges and Universities of Jiangxi Province(GL18225)General Project of Social Science Planning in Jiangxi Province(19GL06).
文摘The inter-sequence analysis method is a popular decision making method considering subjective information of the decision maker.The validity of this method is sensitive to the subjective information provided by the decision maker significantly.Thus,a new weighting method based on the inter-sequence analysis method is proposed.This new method can appropriately reduce the participation of the subjective judgment of the decision maker,and then accordingly decrease the influence of the subjective judgment to the results.It only needs the decision maker to give the inter sequence of the criteria.The importance comparison of two criteria is objectively determined by considering the ratio of a criterion being superior to another one.Some numerical examples are given to illustrate the calculation process of the new method.
