A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avo...A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.展开更多
Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order tim...Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.展开更多
In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entrop...In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.展开更多
A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of ...A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.展开更多
The test scheme of static postponing time given in our actual national military test standard on shrapnels used to control riot is a nine-point test scheme on the combined action of three kinds of temperatures and thr...The test scheme of static postponing time given in our actual national military test standard on shrapnels used to control riot is a nine-point test scheme on the combined action of three kinds of temperatures and three kinds of pressures, the consumed ammunitions are more excessive. Statistic analysis and tentative checkout about a lot of test data are done, feasibility gists are put forward for optimizing of the test design scheme. The optimizing design and data analysis of test scheme of the item are done by managing uniformity design theory, two scientific and reasonable six-point test schemes are confirmed. The feasibility and reliability of the optimizing design schemes put forward above are proved ulteriorly by test validating. The gained schemes not only have good design uniformity and little ammunition wastage and meet the test demand, but also have better forecast ability for the result data of other points using the mathematic models from the actual test points.展开更多
A model for fast electron-driven high-density plasma is proposed to describe the effect of injected fast electrons on the temperature and inner pressure of the plasma in the fast heating process of the double-cone ign...A model for fast electron-driven high-density plasma is proposed to describe the effect of injected fast electrons on the temperature and inner pressure of the plasma in the fast heating process of the double-cone ignition(DCI)scheme.Due to the collision of the two low-density plasmas,the density and volume of the high-density plasma vary.Therefore,the ignition temperature and energy requirement of the high-density plasma vary at different moments,and the required energy for hot electrons to heat the plasma also changes.In practical experiments,the energy input of hot electrons needs to be considered.To reduce the energy input of hot electrons,the optimal moment and the shortest time for injecting hot electrons with minimum energy are analyzed.In this paper,it is proposed to inject hot electrons for a short time to heat the high-density plasma to a relatively high temperature.Then,the alpha particles with the high heating rate and PdV work heat the plasma to the ignition temperature,further reducing the energy required to inject hot electrons.The study of the injection time of fast electrons can reduce the energy requirement of fast electrons for the high-density plasma and increase the probability of successful ignition of the high-density plasma.展开更多
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 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.展开更多
We derive in this paper a time stable seventh-order dissipative compact finite difference scheme with simultaneous approximation terms(SATs) for solving two-dimensional Euler equations. To stabilize the scheme, the ch...We derive in this paper a time stable seventh-order dissipative compact finite difference scheme with simultaneous approximation terms(SATs) for solving two-dimensional Euler equations. To stabilize the scheme, the choice of penalty coefficients for SATs is studied in detail. It is demonstrated that the derived scheme is quite suitable for multi-block problems with different spacial steps. The implementation of the scheme for the case with curvilinear grids is also discussed.Numerical experiments show that the proposed scheme is stable and achieves the design seventh-order convergence rate.展开更多
In this work, a new numerical technique is proposed for the resolution of a fluid model based on three Boltzmann moments. The main purpose of this technique is to calculate electric and physical properties in the non-...In this work, a new numerical technique is proposed for the resolution of a fluid model based on three Boltzmann moments. The main purpose of this technique is to calculate electric and physical properties in the non-equilibrium electric discharge at low pressure. The transport and Poisson＇s equations form a self-consistent model. This equation system is written in cylindrical coordinates following the geometric shape of a plasma reactor. Our transport equation system is discretized using the finite volume approach and resolved by the N-BEE explicit scheme coupled to the time splitting method. This programming structure reduces computation time considerably. The 2D code is carried out and tested by comparing our results with those found in literature.展开更多
By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improv...By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.展开更多
基金Supported by the National Basic Research Program of China(2012CB825601)the National Natural Science Foundationof China(41031066,41231068,41274192,41074121,41204127)+1 种基金the Knowledge Innovation Program of the ChineseAcademy of Sciences(KZZD-EW-01-4)the Specialized Research Fund for State Key Laboratories
文摘A new hybrid numerical scheme of combining an E-CUSP(Energy-Convective Upwind and Split Pressure) method for the fluid part and the Constrained Transport(CT) for the magnetic induction part is proposed.In order to avoid the occurrence of negative pressure in the reconstructed profiles and its updated value,a positivity preserving method is provided.Furthermore,the MHD equations are solved at each physical time step by advancing in pseudo time.The use of dual time stepping is beneficial in the computation since the use of dual time stepping allows the physical time step not to be limited by the corresponding values in the smallest cell and to be selected based on the numerical accuracy criterion.This newly established hybrid scheme combined with positivity preserving method and dual time technique has demonstrated the accurateness and robustness through numerical experiments of benchmark problems such as the 2D Orszag-Tang vortex problem and the3 D shock-cloud interaction problem.
基金Supported by the Discipline Construction and Teaching Research Fund of LUTcte(20140089)
文摘Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.
基金Supported in part by the National Natural Science of China, NSF Grant No. DMS-8657319.
文摘In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.
基金The Project Supported by National Natural Science Foundation of China.
文摘A natural generalization of random choice finite difference scheme of Harten and Lax for Courant number larger than 1 is obtained. We handle interactions between neighboring Riemann solvers by linear superposition of their conserved quantities. We show consistency of the scheme for arbitrarily large Courant numbers. For scalar problems the scheme is total variation diminishing.A brief discussion is given for entropy condition.
文摘The test scheme of static postponing time given in our actual national military test standard on shrapnels used to control riot is a nine-point test scheme on the combined action of three kinds of temperatures and three kinds of pressures, the consumed ammunitions are more excessive. Statistic analysis and tentative checkout about a lot of test data are done, feasibility gists are put forward for optimizing of the test design scheme. The optimizing design and data analysis of test scheme of the item are done by managing uniformity design theory, two scientific and reasonable six-point test schemes are confirmed. The feasibility and reliability of the optimizing design schemes put forward above are proved ulteriorly by test validating. The gained schemes not only have good design uniformity and little ammunition wastage and meet the test demand, but also have better forecast ability for the result data of other points using the mathematic models from the actual test points.
基金Project supported by the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA_(2)5051000)the National Key R&D Program of China(Grant No.2023YFA1608400)+1 种基金the National Natural Science Foundation of China(Grant No.12005008)the Natural Science Foundation of Top Talent of SZTU(Grant No.GDRC202209).
文摘A model for fast electron-driven high-density plasma is proposed to describe the effect of injected fast electrons on the temperature and inner pressure of the plasma in the fast heating process of the double-cone ignition(DCI)scheme.Due to the collision of the two low-density plasmas,the density and volume of the high-density plasma vary.Therefore,the ignition temperature and energy requirement of the high-density plasma vary at different moments,and the required energy for hot electrons to heat the plasma also changes.In practical experiments,the energy input of hot electrons needs to be considered.To reduce the energy input of hot electrons,the optimal moment and the shortest time for injecting hot electrons with minimum energy are analyzed.In this paper,it is proposed to inject hot electrons for a short time to heat the high-density plasma to a relatively high temperature.Then,the alpha particles with the high heating rate and PdV work heat the plasma to the ignition temperature,further reducing the energy required to inject hot electrons.The study of the injection time of fast electrons can reduce the energy requirement of fast electrons for the high-density plasma and increase the probability of successful ignition of the high-density plasma.
基金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 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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11601517)the Basic Research Foundation of National University of Defense Technology(Grant No.ZDYYJ-CYJ20140101)
文摘We derive in this paper a time stable seventh-order dissipative compact finite difference scheme with simultaneous approximation terms(SATs) for solving two-dimensional Euler equations. To stabilize the scheme, the choice of penalty coefficients for SATs is studied in detail. It is demonstrated that the derived scheme is quite suitable for multi-block problems with different spacial steps. The implementation of the scheme for the case with curvilinear grids is also discussed.Numerical experiments show that the proposed scheme is stable and achieves the design seventh-order convergence rate.
文摘In this work, a new numerical technique is proposed for the resolution of a fluid model based on three Boltzmann moments. The main purpose of this technique is to calculate electric and physical properties in the non-equilibrium electric discharge at low pressure. The transport and Poisson＇s equations form a self-consistent model. This equation system is written in cylindrical coordinates following the geometric shape of a plasma reactor. Our transport equation system is discretized using the finite volume approach and resolved by the N-BEE explicit scheme coupled to the time splitting method. This programming structure reduces computation time considerably. The 2D code is carried out and tested by comparing our results with those found in literature.
基金Project supported by the National Natural Science Foundation of China （Grant No 70271069）.
文摘By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number （Re） flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50 000 in the literature; however, this paper can readily extend the maximum to 1000 000 with the above combination.
