A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and im...A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.展开更多
The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the ...The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.展开更多
In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with...In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with adjustable time intervals, which is proved to be effective and remarkably time-saving in numerical tests and applications.展开更多
Direct simulation of 3-D MHD(magnetohydrodynamics) flows in liquid metal fusion blanket with flow channel insert(FCI) has been conducted.Two kinds of pressure equilibrium slot (PES) in FCI,which are used to balance th...Direct simulation of 3-D MHD(magnetohydrodynamics) flows in liquid metal fusion blanket with flow channel insert(FCI) has been conducted.Two kinds of pressure equilibrium slot (PES) in FCI,which are used to balance the pressure difference between the inside and outside of FCI,are considered with a slot in Hartmann wall or a slot in side wall,respectively.The velocity and pressure distribution of FCI made of SiC/SiC_f are numerically studied to illustrate the 3-D MHD flow effects,which clearly show that the flows in fusion blanket with FCI are typical three-dimensional issues and the assumption of 2-D fully developed flows is not the real physical problem of the MHD flows in dual-coolant liquid metal fusion blanket.The optimum opening location of PES has been analyzed based on the 3-D pressure and velocity distributions.展开更多
Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmos...Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmospheric and oceanic equations is analyzed and the computationally stable explicit complete square conservative difference schemes are constructed. The theoretical analysis and numerical experiment prove that the explicit complete square conservative difference schemes are computationally stable and deserve to be disseminated.展开更多
The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the ...The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.展开更多
In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical...In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical observation data.The schemes have three advantages:high-order accuracy in time,generalized square conservation,and smart use of historical observations.Numerical tests based on the one-dimensional linear advection equations suggest that reasonable consideration of accuracy,square conservation,and inclusion of historical observations is critical for good performance of a finite difference scheme.展开更多
For the conservative and non-conservative schemes of nonlinear evolution equations, by taking the two-dimensional shallow water wave equations as an example, a comparative analysis on computational stability is carrie...For the conservative and non-conservative schemes of nonlinear evolution equations, by taking the two-dimensional shallow water wave equations as an example, a comparative analysis on computational stability is carried out. The relationship between the nonlinear computational stability, the structure of the difference schemes, and the form of initial values is also discussed.展开更多
In order to meet the needs of work in numerical weather forecast and in numerical simulations for climate change and ocean current, a kind of difference scheme in high precision in the time direction developed from th...In order to meet the needs of work in numerical weather forecast and in numerical simulations for climate change and ocean current, a kind of difference scheme in high precision in the time direction developed from the completely square-conservative difference scheme in explicit way is built by means of the Taylor expansion. A numerical test with 4-wave Rossby-Haurwitz waves on them and an application of them on the monthly mean current the of South China Sea are carried out, from which, it is found that not only do the new schemes have high harmony and approximate precision but also can the time step of the schemes be lengthened and can much computational time be saved. Therefore, they are worth generalizing and applying.展开更多
In this paper, we investigate a numerical method for the generalized Novikov equation. We propose a conservative finite difference scheme and use Brouwer fixed point theorem to obtain the existence of the solution of ...In this paper, we investigate a numerical method for the generalized Novikov equation. We propose a conservative finite difference scheme and use Brouwer fixed point theorem to obtain the existence of the solution of the corresponding difference equation. We also prove the convergence and stability of the solution by using the discrete energy method. Moreover, we obtain the truncation error of the difference scheme which is .展开更多
With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor co...With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.展开更多
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.展开更多
Existing broadcasting schemes provide services for the stored videos. The basic approach in these schemes is to divide the video into segments and organize them over the channels for proper transmission. Some schemes ...Existing broadcasting schemes provide services for the stored videos. The basic approach in these schemes is to divide the video into segments and organize them over the channels for proper transmission. Some schemes use segments as a basic unit, whereas the others require segments to be further divided into subsegments. In a scheme, the number of segments/subsegments depends upon the bandwidth allocated to the video by the video server. For constructing segments, the video length should be known. If it is unknown, then the segments cannot be constructed and hence the scheme cannot be applied to provide the video services. This is an important issue especially in live broadcasting applications wherein the ending time of the video is unknown, for example, cricket match. In this paper, we propose a mechanism for the conservative staircase scheme so that it can support live video broadcasting.展开更多
A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some nume...A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.展开更多
In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated an...In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.展开更多
In this paper, two finite difference schemes are presented for initial-boundary value problems of Regularized Long-Wave(RLW) equation. They all have the advantages that there are discrete energies which are conserved....In this paper, two finite difference schemes are presented for initial-boundary value problems of Regularized Long-Wave(RLW) equation. They all have the advantages that there are discrete energies which are conserved. Convergence and stability of difference solutions with order O(h2+τ2) are proved in the energy norm. Numerical experiment results demonstrate the effectiveness of the proposed schemes.展开更多
In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws...In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.展开更多
A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourthorder central weighted essentially nonoscillatory (CW...A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourthorder central weighted essentially nonoscillatory (CWENO) reconstruction for one-dimensional cases, which is generalized to two-dimensional cases by the dimension-by-dimension approach. The large stability domain Runge-Kutta-type solver ROCK4 is used for time integration. The resulting method requires neither the use of Riemann solvers nor the computation of Jacobians and therefore it enjoys the main advantage of the relaxation schemes. The high accuracy and high-resolution properties of the present method are demonstrated in one- and two-dimensional numerical experiments.展开更多
In order to improve the accuracy of forecasts of atmospheric and oceanic phenomena which possess a wide range of space and time scales, it is crucial to design the high-order and stable schemes. On the basis of the ex...In order to improve the accuracy of forecasts of atmospheric and oceanic phenomena which possess a wide range of space and time scales, it is crucial to design the high-order and stable schemes. On the basis of the explicit square-conservative scheme, a high-order compact explicit square-conservative scheme is proposed in this paper. This scheme not only keeps the square-conservative characteristics, but also is of high accuracy. The numerical example shows that this scheme has less computing errors and better computational stability, and it could be considered to be tested and used in many atmospheric and oceanic problems.展开更多
This paper establishes the resolution finite volume scheme with slope entropy convergence of a new two-value high relaxation for conservation laws. This scheme, motivated by the general method of high resolution schem...This paper establishes the resolution finite volume scheme with slope entropy convergence of a new two-value high relaxation for conservation laws. This scheme, motivated by the general method of high resolution schemes that have high-order accuracy in smooth regions of solutions and are free of oscillations near discontinuities, unifies and evolves slopes directly with a slope relaxation equation that governs the evolution of slopes in both smooth and discontinuous regions. Proper choices of slopes are realized adaptively via a relaxation parameter. The scheme is shown to be total-variation-bounded (TVB) stable and satisfies cell-entropy inequalities.展开更多
文摘A new conservative finite difference scheme is presented based on the numerical analysis for an initialboundary value problem of a class of Schroedinger equation with the wave operator. The scheme can be linear and implicit or explicit based on the parameter choice. The initial value after discretization has second-order accuracy that is consistent with the scheme accuracy. The existence and the uniqueness of the difference solution are proved. Based on the priori estimates and an inequality about norms, the stability and the convergence of difference solutions with the second-order are proved in the energy norm. Experimental results demonstrate the efficiency of the new scheme.
文摘The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.
基金Partly supported by the State Major Key Project for Basic Researches of China
文摘In the present paper, a class of explicit forward time-difference schemes are established from a geometric view with strict analytical deductions. This class includes the schemes with a constant time interval and with adjustable time intervals, which is proved to be effective and remarkably time-saving in numerical tests and applications.
基金supported by National Natural Science Foundation of China with grant Nos.10872212,50936006National Magnetic Confinement Fusion Science Program in China with grant No.2009GB10401
文摘Direct simulation of 3-D MHD(magnetohydrodynamics) flows in liquid metal fusion blanket with flow channel insert(FCI) has been conducted.Two kinds of pressure equilibrium slot (PES) in FCI,which are used to balance the pressure difference between the inside and outside of FCI,are considered with a slot in Hartmann wall or a slot in side wall,respectively.The velocity and pressure distribution of FCI made of SiC/SiC_f are numerically studied to illustrate the 3-D MHD flow effects,which clearly show that the flows in fusion blanket with FCI are typical three-dimensional issues and the assumption of 2-D fully developed flows is not the real physical problem of the MHD flows in dual-coolant liquid metal fusion blanket.The optimum opening location of PES has been analyzed based on the 3-D pressure and velocity distributions.
基金the Outstanding State Key Laboratory Project of National Science Foundation of China (Grant No. 40023001 )the Key Innovatio
文摘Based on the forced dissipetive nonlinear evolution equations for describing the motion of atmosphere and ocean, the computational stability of the explicit difference schemes of the forced dissipotive nonlinear atmospheric and oceanic equations is analyzed and the computationally stable explicit complete square conservative difference schemes are constructed. The theoretical analysis and numerical experiment prove that the explicit complete square conservative difference schemes are computationally stable and deserve to be disseminated.
文摘The conservative form and singular perturbed ordinary differential equation with periodic boundary value problem were studied, and a conservative difference scheme was constructed. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, it is proved that the scheme converges uniformly to the solution of differential equation with order one.
基金the Ministry of Science and Technology of China for funding the National Basic Research Program of China (973 Program,Grant No.2011CB309704)
文摘In this paper,a multistep finite difference scheme has been proposed,whose coefficients are determined taking into consideration compatibility and generalized quadratic conservation,as well as incorporating historical observation data.The schemes have three advantages:high-order accuracy in time,generalized square conservation,and smart use of historical observations.Numerical tests based on the one-dimensional linear advection equations suggest that reasonable consideration of accuracy,square conservation,and inclusion of historical observations is critical for good performance of a finite difference scheme.
文摘For the conservative and non-conservative schemes of nonlinear evolution equations, by taking the two-dimensional shallow water wave equations as an example, a comparative analysis on computational stability is carried out. The relationship between the nonlinear computational stability, the structure of the difference schemes, and the form of initial values is also discussed.
文摘In order to meet the needs of work in numerical weather forecast and in numerical simulations for climate change and ocean current, a kind of difference scheme in high precision in the time direction developed from the completely square-conservative difference scheme in explicit way is built by means of the Taylor expansion. A numerical test with 4-wave Rossby-Haurwitz waves on them and an application of them on the monthly mean current the of South China Sea are carried out, from which, it is found that not only do the new schemes have high harmony and approximate precision but also can the time step of the schemes be lengthened and can much computational time be saved. Therefore, they are worth generalizing and applying.
文摘In this paper, we investigate a numerical method for the generalized Novikov equation. We propose a conservative finite difference scheme and use Brouwer fixed point theorem to obtain the existence of the solution of the corresponding difference equation. We also prove the convergence and stability of the solution by using the discrete energy method. Moreover, we obtain the truncation error of the difference scheme which is .
基金Dr.X.L.Li at the China Meteorological Administration.This study was supported by grants from the National Natural Science Foundation of China(Nos.41575103 and 91637210).
文摘With an increase in model resolution,compact high-order numerical advection scheme can improve its effectiveness and competitiveness in oceanic modeling due to its high accuracy and scalability on massive-processor computers.To provide high-quality numerical ocean simulation on overset grids,we tried a novel formulation of the fourth-order multi-moment constrained finite volume scheme to simulate continuous and discontinuous problems in the Cartesian coordinate.Utilizing some degrees of freedom over each cell and derivatives at the cell center,we obtained a two-dimensional(2D)cubic polynomial from which point values on the extended overlap can achieve fourth-order accuracy.However,this interpolation causes a lack of conservation because the flux between the regions are no longer equal;thus,a flux correction is implemented to ensure conservation.A couple of numerical experiments are presented to evaluate the numerical scheme,which confirms its approximately fourth-order accuracy in conservative transportation on overset grid.The test cases reveal that the scheme is effective to suppress numerical oscillation in discontinuous problems,which may be powerful for salinity advection computing with a sharp gradient.
基金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.
文摘Existing broadcasting schemes provide services for the stored videos. The basic approach in these schemes is to divide the video into segments and organize them over the channels for proper transmission. Some schemes use segments as a basic unit, whereas the others require segments to be further divided into subsegments. In a scheme, the number of segments/subsegments depends upon the bandwidth allocated to the video by the video server. For constructing segments, the video length should be known. If it is unknown, then the segments cannot be constructed and hence the scheme cannot be applied to provide the video services. This is an important issue especially in live broadcasting applications wherein the ending time of the video is unknown, for example, cricket match. In this paper, we propose a mechanism for the conservative staircase scheme so that it can support live video broadcasting.
文摘A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.
文摘In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.
基金The project supported by the national natural science foundation of China (10471023,10572057)
文摘In this paper, two finite difference schemes are presented for initial-boundary value problems of Regularized Long-Wave(RLW) equation. They all have the advantages that there are discrete energies which are conserved. Convergence and stability of difference solutions with order O(h2+τ2) are proved in the energy norm. Numerical experiment results demonstrate the effectiveness of the proposed schemes.
基金Project supported by the National Natural Science Foundation of China(No.11571366)the Basic Research Foundation of National Numerical Wind Tunnel Project(No.NNW2018-ZT4A08)
文摘In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.
基金the National Natural Science Foundation of China (60134010)The English text was polished by Yunming Chen.
文摘A fourth-order relaxation scheme is derived and applied to hyperbolic systems of conservation laws in one and two space dimensions. The scheme is based on a fourthorder central weighted essentially nonoscillatory (CWENO) reconstruction for one-dimensional cases, which is generalized to two-dimensional cases by the dimension-by-dimension approach. The large stability domain Runge-Kutta-type solver ROCK4 is used for time integration. The resulting method requires neither the use of Riemann solvers nor the computation of Jacobians and therefore it enjoys the main advantage of the relaxation schemes. The high accuracy and high-resolution properties of the present method are demonstrated in one- and two-dimensional numerical experiments.
文摘In order to improve the accuracy of forecasts of atmospheric and oceanic phenomena which possess a wide range of space and time scales, it is crucial to design the high-order and stable schemes. On the basis of the explicit square-conservative scheme, a high-order compact explicit square-conservative scheme is proposed in this paper. This scheme not only keeps the square-conservative characteristics, but also is of high accuracy. The numerical example shows that this scheme has less computing errors and better computational stability, and it could be considered to be tested and used in many atmospheric and oceanic problems.
基金Project supported by the National Natural Science Foundation of China(Nos.11371063,11501040,and 91530108)the Doctoral Program from the Education Ministry of China(No.20130003110004)
文摘This paper establishes the resolution finite volume scheme with slope entropy convergence of a new two-value high relaxation for conservation laws. This scheme, motivated by the general method of high resolution schemes that have high-order accuracy in smooth regions of solutions and are free of oscillations near discontinuities, unifies and evolves slopes directly with a slope relaxation equation that governs the evolution of slopes in both smooth and discontinuous regions. Proper choices of slopes are realized adaptively via a relaxation parameter. The scheme is shown to be total-variation-bounded (TVB) stable and satisfies cell-entropy inequalities.