In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald...In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.展开更多
In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formul...In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.展开更多
The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the ...The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.展开更多
Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploi...Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.展开更多
The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important iss...The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important issues being concerned in past researches.An efficient method for determining the interaction of piles is introduced in present research.First the wave force is calculated by the exact method using the diffraction theory,then in the finite difference numerical method the force is calculated by adding the velocity potentials of each pile and integration of pressure on their surface.The results showed that the ratio of the wave force on each of the double piles to a single pile has a damped oscillation around unity in which the amplitude of oscillation decreases with the increase in the spacing parameter.Also different wave incident directions and diffraction parameters were used and the results showed that the numerical solution has acceptable accuracy when the diffraction parameter is larger than unity.展开更多
An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solutio...An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.展开更多
The Zeldovich-von Neumann-Doring model and the Chapman-Jouguet model for a simplified combustion model-Majda＇s model is studied.The author proves a uniform maximum norm estimate,then proves that as the rate of chemic...The Zeldovich-von Neumann-Doring model and the Chapman-Jouguet model for a simplified combustion model-Majda＇s model is studied.The author proves a uniform maximum norm estimate,then proves that as the rate of chemical reaction tends to infinity the solutions to the Zeldovich-von Neumann-Doring model tend to that of the Chapman-Jouguet model.The type of combustion waves is studied.This result is compared with the result of the projection and finite difference method for the same model.展开更多
Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface...Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.展开更多
An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient....An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient. The cross-correlation algorithm with sub-pixel accuracy is introduced and used to calculate the background-element displacement of a disturbed image and a fourth-order difference scheme is also developed to solve the Poisson equation. An experiment for a disturbed flow field caused by a burning candle is performed to validate the built BOS system and the results indicate that density and temperature distribution of the disturbed flow field can be reconstructed accurately. A notable conclusion is that in order to make the reconstructed results have a satisfactory accuracy, the inquiry step length should be less than the size of the interrogation window.展开更多
In this paper, a computer visualization approach is proposed for electromagnetic wave interaction with structures by mains of finite difference-time doain method (F-D) and computer graphics. By visualization of FDTD, ...In this paper, a computer visualization approach is proposed for electromagnetic wave interaction with structures by mains of finite difference-time doain method (F-D) and computer graphics. By visualization of FDTD, Phenomena such as wave propagation, penetration through structures, renection and absorption by structures are observed. Visualization of electromagnetic wave interactions with two wing-shaped structures is demonstrated. These examples indicate that the approach describe in the paper offers an effective way for investigating electromagnetic wave phenomena and is helpful to the engineers in controlling radar signature of the targets.展开更多
This article describes a local error estimator for Glimm＇s scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm＇s scheme by an optimal choice. As a by-product of...This article describes a local error estimator for Glimm＇s scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm＇s scheme by an optimal choice. As a by-product of the local error estimator, the procedure provides a global error estimator that is shown numerically to be a very accurate estimate of the error in L1 （R） for all times. Although there is partial mathematical evidence for the error estimator proposed, at this stage the error estimator must be considered ad- hoc. Nonetheless, the error estimator is simple to compute, relatively inexpensive, without adjustable parameters and at least as accurate as other existing error estimators. Numerical experiments in 1-D for Burgers＇ equation and for Euler＇s system are performed to measure the asymptotic accuracy of the resulting scheme and of the error estimator.展开更多
Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the h...Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the help of intrinsically safe borehole video instrument.And according to the geology of this working face,a discrete element model was calculated by UDEC.Combining in situ experimental data with numerical results,the relationship between the fractal dimension of boreholes'wall and the distribution of advanced abutment pressure was studied under the condition of mining advance.The results show that the variation tendency of fractal dimension and the abutment pressure has the same characteristic value.The distance between working face and the peak value of the abutment pressure has a slight increasing trend with the advancing of mining-face.When the working face is set as the original point,the trend of fractal dimension from the far place to the origin can be divided into three phases:constant,steady increasing and constant.And the turning points of these phases are the max-influencing distance(50 m)and peak value(15 m)of abutment pressure.展开更多
The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used ...The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used to evaluate the azimuthal character of borehole formation. In this paper, a 3D finite difference method was used to simulate the acoustic fields in a fluid-filled borehole generated by a traditional monopole source and a phased arc array. Acoustic waveforms were presented for both cases. The analysis of the simulated waveforms showed that different from the monopole source, the acoustic energy generated by the phased arc array transmitter mainly radiated to the borehole in a narrow azimuthal range, which was the key technique to implement azimuthal acoustic well logging. Similar to the monopole source, the waveforms generated by the phased arc array in the fluid-filled borehole also contain compressional （P） waves and shear （S） waves refracted along the borehole, which is the theoretical foundation of phased arc array acoustic well logging.展开更多
In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numer...In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.展开更多
The cluster-shaped plasmonic nanostructures are used to manage the incident light inside an ultra-thin silicon solar cell.Here we simulate spherical,conical,pyramidal,and cylindrical nanoparticles in a form of a clust...The cluster-shaped plasmonic nanostructures are used to manage the incident light inside an ultra-thin silicon solar cell.Here we simulate spherical,conical,pyramidal,and cylindrical nanoparticles in a form of a cluster at the rear side of a thin silicon cell,using the finite difference time domain(FDTD)method.By calculating the optical absorption and hence the photocurrent,it is shown that the clustering of nanoparticles significantly improves them.The photocurrent enhancement is the result of the plasmonic effects of clustering the nanoparticles.For comparison,first a cell with a single nanoparticle at the rear side is evaluated.Then four smaller nanoparticles are put around it to make a cluster.The photocurrents of 20.478 mA/cm2,23.186 mA/cm2,21.427 mA/cm2,and 21.243 mA/cm2 are obtained for the cells using clustering conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.These values are 13.987 mA/cm2,16.901 mA/cm2,16.507 mA/cm2,17.926 mA/cm2 for the cell with one conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.Therefore,clustering can significantly improve the photocurrents.Finally,the distribution of the electric field and the generation rate for the proposed structures are calculated.展开更多
A general numerical tool, based on thermal diffusion equation and full-vectorial eigen-mode equation, has been presented for the systematic analysis of graded index channel waveguide fabricated by ion exchange on Er^3...A general numerical tool, based on thermal diffusion equation and full-vectorial eigen-mode equation, has been presented for the systematic analysis of graded index channel waveguide fabricated by ion exchange on Er^3＋ doped glass. Finite difference method with full-vectorial formulation （FV-FDM） is applied to solving the full-vectorial modes of graded index channel waveguide for the first time. The coupled difference equations based on magnetic fields in FV-FDM are derived from the Taylor series expansion and accurate formulation of boundary conditions. Hybrid nature of vectorial guided modes for both pump （980 nm） and signal light （1550 nm） are demonstrated by the simulation. Results show that the fabrication parameters of ion exchange, such as channel opening width and time ratio of second step to first step in ion exchange, have large influence on the properties of waveguide. By optimizing the fabrication parameters, maintenance of monomode for signal light and improvement of the gain dynamics can be achieved in Er^3＋ doped waveguide amplifier （EDWA） fabricated by ion exchange technique. This theoretical model is significant for the design and fabrication of EDWA with ion exchange technique. Furthermore, a single polarization EDWA, which operates at wavelength from 1528 nm to 1541 nm for HE polarization, is numerically designed.展开更多
In this paper,we consider a delayed diffusive SVEIR model with general incidence.We first establish the threshold dynamics of this model.Using a Nonstandard Finite Difference(NSFD) scheme,we then give the discretizati...In this paper,we consider a delayed diffusive SVEIR model with general incidence.We first establish the threshold dynamics of this model.Using a Nonstandard Finite Difference(NSFD) scheme,we then give the discretization of the continuous model.Applying Lyapunov functions,global stability of the equilibria are established.Numerical simulations are presented to validate the obtained results.The prolonged time delay can lead to the elimination of the infectiousness.展开更多
Under higher temperatures the charge within a warhead will be subjected to physical and chemical changes, which will influence the security of the warhead launching process. In this paper the problem is studied adopti...Under higher temperatures the charge within a warhead will be subjected to physical and chemical changes, which will influence the security of the warhead launching process. In this paper the problem is studied adopting finite difference method for the case of a rocket powered missile. Temperature distribution tables are given through quadrature experiments, and the results fit the actual measured values very well. The results also show that to ensure the launching security of rocket powered missiles having warhead charge critical temperature close to about 60?℃, the designer can select appropriate heat insulation structural parameters.展开更多
The factors affecting one dimensional （1D） and two dimensional （2D） photonic crystals （PhCs） are systemically analyzed in this paper by numerical simulation. Transfer matrix method （TMM） is employed for 1D PCs...The factors affecting one dimensional （1D） and two dimensional （2D） photonic crystals （PhCs） are systemically analyzed in this paper by numerical simulation. Transfer matrix method （TMM） is employed for 1D PCs, both finite difference time domain method （FDTD） and plane wave expansion method （PWE） are employed for 2D PCs. The result shows that the photonic bandgaps （PBG） are directly affected by crystal type, crystal lattice constant, modulation of refractive index and periodicity, and it is should be useful for design of different type photonic crystals with the required PBG and functional devices. Finally, as an example, a near-IR 1D PCs narrow filter was designed.展开更多
In this study,we investigate the influence of doping on the charge transfer and device characteristics parameters in the bulk heterojunction solar cells based on poly（3-hexylthiophene）（P3HT） and a methanofuUerene ...In this study,we investigate the influence of doping on the charge transfer and device characteristics parameters in the bulk heterojunction solar cells based on poly（3-hexylthiophene）（P3HT） and a methanofuUerene derivative（PCBM）.Organic semiconductors are also known to be not pure and they have defects and impurities,some of them are being charged and act as p-type or n-type dopants.Calculations of the solar cell characteristics parameters versus the p-doping level have been done at three different n-dopings（N_d） that consist of 5 × 10^（17） cm^（-3）,10^（18） cm^（-3）,and 5 × 10^（18） cm^（-3）.We perform the analysis of the doping concentration through the drift-diffusion model,and calculate the current and voltage doping dependency.We find that at three different n-dopant levels,optimum p-type doping is about N_p = 6 × 10^（18） cm^（-3）.Simulation results have shown that by increasing doping level,V_（oc） monotonically increases by doping.Cell efficiency reaches its maximum at somewhat higher doping as FF has its peak at N_p = 3 × 10^（18） cm^（-3）.Moreover,this paper demonstrates that the optimum value for the p-doping is about N_p = 6 × 10^（18） cm^（-3） and optimum value for n-dopant is N_d = 10^（18） cm^（-3）,respectively.The simulated results confirm that doping considerably affects the performance of organic solar cells.展开更多
文摘In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.
基金partly supported by the National Natural Science Foundation of China(No.11875070)the Natural Science Foundation of Anhui Province(No.1908085MA16)
文摘In this study,single-particle energy was examined using the finite difference method by taking 208Pb as an example.If the first derivative term in the spherical Dirac equation is discretized using a three-point formula,a one-to-one correspondence occurs between the physical and spurious states.Although these energies are exactly the same,the wave functions of the spurious states exhibit a much faster staggering than those of the physical states.Such spurious states can be eliminated when applying the finite difference method by introducing an extra Wilson term into the Hamiltonian.Furthermore,it was also found that the number of spurious states can be reduced if we improve the accuracy of the numerical differential formula.The Dirac equation is then solved in a momentum space in which there is no differential operator,and we found that the spurious states can be completely avoided in the momentum space,even without an extra Wilson term.
基金Project supported partly by the Open Research Program in State Key Laboratory of Millimeter Waves of China(Grant No.K200802)partly by the National Natural Science Foundation of China(Grant No.60971122)
文摘The transmission coefficients of electromagnetic （EM） waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain （SO-FDTD） method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.
基金supported by the Major State BasicResearch Program of China(19990328)the National Tackling Key Problem Programs(20050200069)+4 种基金the National Natural Science Foundation of China(1077112410372052)the Doctorate Foundation of the Ministryof Education of China(20030422047)Shandong Provance Natural Science Foundation(2R2009AQ12)the Independent Innovation Foundation of Shandong University(2010TS031)
文摘Coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.The upwind finite difference schemes applicable to parallel arithmetic are put forward and two-dimensional and three-dimensional schemes are used to form a complete set.Some techniques,such as change of variables,calculus of variations, multiplicative commutation rule of difference operators,decomposition of high order difference operators and prior estimates,are adopted.The estimates in l~2 norm are derived to determine the error in the approximate solution.This method was already applied to the numerical simulation of migration-accumulation of oil resources.
文摘The wave force exerted on vertical piles of offshore structures is the main criterion in designing them.In structures with more than one large pile,the influence of piles on each other is one of the most important issues being concerned in past researches.An efficient method for determining the interaction of piles is introduced in present research.First the wave force is calculated by the exact method using the diffraction theory,then in the finite difference numerical method the force is calculated by adding the velocity potentials of each pile and integration of pressure on their surface.The results showed that the ratio of the wave force on each of the double piles to a single pile has a damped oscillation around unity in which the amplitude of oscillation decreases with the increase in the spacing parameter.Also different wave incident directions and diffraction parameters were used and the results showed that the numerical solution has acceptable accuracy when the diffraction parameter is larger than unity.
文摘An explicit,time-dependent variable grid finite difference method is introduced and analyzed for approximating the solution of a scalar conservation law in two dimension. The scheme is stable,and the numerical solution is proved to converges to the relevant physical solution.
基金the China State Major Key Project for Basic Researchesthe Science Fund of the Ministry of Education of China
文摘The Zeldovich-von Neumann-Doring model and the Chapman-Jouguet model for a simplified combustion model-Majda＇s model is studied.The author proves a uniform maximum norm estimate,then proves that as the rate of chemical reaction tends to infinity the solutions to the Zeldovich-von Neumann-Doring model tend to that of the Chapman-Jouguet model.The type of combustion waves is studied.This result is compared with the result of the projection and finite difference method for the same model.
基金supported in part by R & D Department of China National Petroleum Corporation (2022DQ0604-01)National Natural Science Foundation of China (42204132)+3 种基金the China Postdoctoral Science Foundations (2020M680667, 2021T140661)Harvard-CUP Joint Laboratory on Petroleum Science“111” project (B13010)the financial support from the CAS Special Research Assistant Project。
文摘Seismic wave propagation in fluid-solid coupled media is currently a popular topic. However, traditional wave equation-based simulation methods have to consider complex boundary conditions at the fluid-solid interface. To address this challenge, we propose a novel numerical scheme that integrates the lattice Boltzmann method(LBM) and lattice spring model(LSM). In this scheme, LBM simulates viscoacoustic wave propagation in the fluid area and LSM simulates elastic wave propagation in the solid area. We also introduce three different LBM-LSM coupling strategies, a standard bounce back scheme, a specular reflection scheme, and a hybrid scheme, to describe wave propagation across fluid-solid boundaries. To demonstrate the accuracy of these LBM-LSM coupling schemes, we simulate wave propagation in a two-layer model containing a fluid-solid interface. We place excitation sources in the fluid layer and the solid layer respectively, to observe the wave phenomena when seismic waves propagate to interface from different sides. The simulated results by LBM-LSM are compared with the reference wavefields obtained by the finite difference method(FDM) and the analytical solution(ANA).Our LBM-LSM coupling scheme was verified effective, as the relative errors between the LBM-LSM solutions and reference solutions were within an acceptable range, sometimes around 1.00%. The coupled LBM-LSM scheme is further used to model seismic wavefields across a more realistic rugged seabed,which reveals the potential applications of the coupled LBM-LSM scheme in marine seismic imaging techniques, such as reverse-time migration and full-waveform inversion. The method also has potential applications in simulating wave propagation in complex two-and multi-phase media.
基金supported by the Key Program of the National Natural Science Foundation of China(Grant No.NSFC 91441205)
文摘An experimental system based on the background-oriented schlieren（BOS） technique is built to reconstruct the density and temperature distribution of a flame-induced distorted flow field which has a density gradient. The cross-correlation algorithm with sub-pixel accuracy is introduced and used to calculate the background-element displacement of a disturbed image and a fourth-order difference scheme is also developed to solve the Poisson equation. An experiment for a disturbed flow field caused by a burning candle is performed to validate the built BOS system and the results indicate that density and temperature distribution of the disturbed flow field can be reconstructed accurately. A notable conclusion is that in order to make the reconstructed results have a satisfactory accuracy, the inquiry step length should be less than the size of the interrogation window.
文摘In this paper, a computer visualization approach is proposed for electromagnetic wave interaction with structures by mains of finite difference-time doain method (F-D) and computer graphics. By visualization of FDTD, Phenomena such as wave propagation, penetration through structures, renection and absorption by structures are observed. Visualization of electromagnetic wave interactions with two wing-shaped structures is demonstrated. These examples indicate that the approach describe in the paper offers an effective way for investigating electromagnetic wave phenomena and is helpful to the engineers in controlling radar signature of the targets.
基金supported by a Korea Research Foundation Grant from the Korean Government(MOEHRD)(KRF-2007-331-C00053)supported by the National Science and Engineering Council of Canada and the Canadian Foundation for Innovation
文摘This article describes a local error estimator for Glimm＇s scheme for hyperbolic systems of conservation laws and uses it to replace the usual random choice in Glimm＇s scheme by an optimal choice. As a by-product of the local error estimator, the procedure provides a global error estimator that is shown numerically to be a very accurate estimate of the error in L1 （R） for all times. Although there is partial mathematical evidence for the error estimator proposed, at this stage the error estimator must be considered ad- hoc. Nonetheless, the error estimator is simple to compute, relatively inexpensive, without adjustable parameters and at least as accurate as other existing error estimators. Numerical experiments in 1-D for Burgers＇ equation and for Euler＇s system are performed to measure the asymptotic accuracy of the resulting scheme and of the error estimator.
基金financial support from the State Key Basic Research Program of China(Nos.2011CB201201and 2010CB226802)the National Natural Science Foundation of China(No.51204112)
文摘Based on the geological conditions of coal mining face No.15-14120 at No.8 mine of Pingdingshan coal mining group,the real-time evolution of coal-roof crack network with working face advancing was collected with the help of intrinsically safe borehole video instrument.And according to the geology of this working face,a discrete element model was calculated by UDEC.Combining in situ experimental data with numerical results,the relationship between the fractal dimension of boreholes'wall and the distribution of advanced abutment pressure was studied under the condition of mining advance.The results show that the variation tendency of fractal dimension and the abutment pressure has the same characteristic value.The distance between working face and the peak value of the abutment pressure has a slight increasing trend with the advancing of mining-face.When the working face is set as the original point,the trend of fractal dimension from the far place to the origin can be divided into three phases:constant,steady increasing and constant.And the turning points of these phases are the max-influencing distance(50 m)and peak value(15 m)of abutment pressure.
基金supported by the National Natural ScienceFoundation of China(Grant Nos.10534040,40574049 and 40874097)the Research Fund for the Doctoral Programof Higher Education(Grant No.20070425028)the Foundation of State Key Laboratory of Petroleum Resourceand Prospecting,China University of Petroleum(Grant No.PRPDX2008-08).
文摘The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used to evaluate the azimuthal character of borehole formation. In this paper, a 3D finite difference method was used to simulate the acoustic fields in a fluid-filled borehole generated by a traditional monopole source and a phased arc array. Acoustic waveforms were presented for both cases. The analysis of the simulated waveforms showed that different from the monopole source, the acoustic energy generated by the phased arc array transmitter mainly radiated to the borehole in a narrow azimuthal range, which was the key technique to implement azimuthal acoustic well logging. Similar to the monopole source, the waveforms generated by the phased arc array in the fluid-filled borehole also contain compressional （P） waves and shear （S） waves refracted along the borehole, which is the theoretical foundation of phased arc array acoustic well logging.
文摘In this paper,various aspects of the 2D and 3D nonlinear liquid sloshing problems in vertically excited containers have been studied numerically along with the help of a modified-transformation.Based on this new numerical algorithm,a numerical study on a regularly and randomly excited container in vertical direction was conducted utilizing four different cases： The first case was performed utilizing a 2D container with regular excitations.The next case examined a regularly excited 3D container with two different initial conditions for the liquid free surface,and finally,3D container with random excitation in the vertical direction.A grid independence study was performed along with a series of validation tests.An iteration error estimation method was used to stop the iterative solver（used for solving the discretized governing equations in the computational domain） upon reaching steady state of results at each time step.In the present case,this method was found to produce quite accurate results and to be more time efficient as compared to other conventional stopping procedures for iterative solvers.The results were validated with benchmark results.The wave elevation time history,phase plane diagram and surface plots represent the wave nonlinearity during its motion.
文摘The cluster-shaped plasmonic nanostructures are used to manage the incident light inside an ultra-thin silicon solar cell.Here we simulate spherical,conical,pyramidal,and cylindrical nanoparticles in a form of a cluster at the rear side of a thin silicon cell,using the finite difference time domain(FDTD)method.By calculating the optical absorption and hence the photocurrent,it is shown that the clustering of nanoparticles significantly improves them.The photocurrent enhancement is the result of the plasmonic effects of clustering the nanoparticles.For comparison,first a cell with a single nanoparticle at the rear side is evaluated.Then four smaller nanoparticles are put around it to make a cluster.The photocurrents of 20.478 mA/cm2,23.186 mA/cm2,21.427 mA/cm2,and 21.243 mA/cm2 are obtained for the cells using clustering conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.These values are 13.987 mA/cm2,16.901 mA/cm2,16.507 mA/cm2,17.926 mA/cm2 for the cell with one conical,spherical,pyramidal,cylindrical NPs at the backside,respectively.Therefore,clustering can significantly improve the photocurrents.Finally,the distribution of the electric field and the generation rate for the proposed structures are calculated.
基金supported by the Foundation for Development of Science and Technology of Shanghai (Grant No 022261002)
文摘A general numerical tool, based on thermal diffusion equation and full-vectorial eigen-mode equation, has been presented for the systematic analysis of graded index channel waveguide fabricated by ion exchange on Er^3＋ doped glass. Finite difference method with full-vectorial formulation （FV-FDM） is applied to solving the full-vectorial modes of graded index channel waveguide for the first time. The coupled difference equations based on magnetic fields in FV-FDM are derived from the Taylor series expansion and accurate formulation of boundary conditions. Hybrid nature of vectorial guided modes for both pump （980 nm） and signal light （1550 nm） are demonstrated by the simulation. Results show that the fabrication parameters of ion exchange, such as channel opening width and time ratio of second step to first step in ion exchange, have large influence on the properties of waveguide. By optimizing the fabrication parameters, maintenance of monomode for signal light and improvement of the gain dynamics can be achieved in Er^3＋ doped waveguide amplifier （EDWA） fabricated by ion exchange technique. This theoretical model is significant for the design and fabrication of EDWA with ion exchange technique. Furthermore, a single polarization EDWA, which operates at wavelength from 1528 nm to 1541 nm for HE polarization, is numerically designed.
文摘In this paper,we consider a delayed diffusive SVEIR model with general incidence.We first establish the threshold dynamics of this model.Using a Nonstandard Finite Difference(NSFD) scheme,we then give the discretization of the continuous model.Applying Lyapunov functions,global stability of the equilibria are established.Numerical simulations are presented to validate the obtained results.The prolonged time delay can lead to the elimination of the infectiousness.
文摘Under higher temperatures the charge within a warhead will be subjected to physical and chemical changes, which will influence the security of the warhead launching process. In this paper the problem is studied adopting finite difference method for the case of a rocket powered missile. Temperature distribution tables are given through quadrature experiments, and the results fit the actual measured values very well. The results also show that to ensure the launching security of rocket powered missiles having warhead charge critical temperature close to about 60?℃, the designer can select appropriate heat insulation structural parameters.
基金Sponsored by the Ministerial Level Advanced Research Foundation (1006513002060)
文摘The factors affecting one dimensional （1D） and two dimensional （2D） photonic crystals （PhCs） are systemically analyzed in this paper by numerical simulation. Transfer matrix method （TMM） is employed for 1D PCs, both finite difference time domain method （FDTD） and plane wave expansion method （PWE） are employed for 2D PCs. The result shows that the photonic bandgaps （PBG） are directly affected by crystal type, crystal lattice constant, modulation of refractive index and periodicity, and it is should be useful for design of different type photonic crystals with the required PBG and functional devices. Finally, as an example, a near-IR 1D PCs narrow filter was designed.
文摘In this study,we investigate the influence of doping on the charge transfer and device characteristics parameters in the bulk heterojunction solar cells based on poly（3-hexylthiophene）（P3HT） and a methanofuUerene derivative（PCBM）.Organic semiconductors are also known to be not pure and they have defects and impurities,some of them are being charged and act as p-type or n-type dopants.Calculations of the solar cell characteristics parameters versus the p-doping level have been done at three different n-dopings（N_d） that consist of 5 × 10^（17） cm^（-3）,10^（18） cm^（-3）,and 5 × 10^（18） cm^（-3）.We perform the analysis of the doping concentration through the drift-diffusion model,and calculate the current and voltage doping dependency.We find that at three different n-dopant levels,optimum p-type doping is about N_p = 6 × 10^（18） cm^（-3）.Simulation results have shown that by increasing doping level,V_（oc） monotonically increases by doping.Cell efficiency reaches its maximum at somewhat higher doping as FF has its peak at N_p = 3 × 10^（18） cm^（-3）.Moreover,this paper demonstrates that the optimum value for the p-doping is about N_p = 6 × 10^（18） cm^（-3） and optimum value for n-dopant is N_d = 10^（18） cm^（-3）,respectively.The simulated results confirm that doping considerably affects the performance of organic solar cells.
