In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments...In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments of Nino3 SST anomalies and Tahiti-Darwin SO index. The results show that the scheme is feasible and ENSO predictable.展开更多
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.展开更多
Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fre...Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.展开更多
A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell＇s equations is introduced and analyzed. The pro- posed method combines a centered approximation for the eval...A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell＇s equations is introduced and analyzed. The pro- posed method combines a centered approximation for the evaluation of fluxes at the in- terface between neighboring elements, with a Nth-order leap-frog time scheme. More- over, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwelrs equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high- order elements show the potential of the method.展开更多
We sttidy the problem of scheduling n jobs on m parallel bounded batch machines to minimize the sum of squared machine loads. Each batch contains at most B jobs, and the processing time of a batch is equal to the long...We sttidy the problem of scheduling n jobs on m parallel bounded batch machines to minimize the sum of squared machine loads. Each batch contains at most B jobs, and the processing time of a batch is equal to the longest processing time of the jobs in this batch. We prove this problem to be NP-hard. Furthermore, we present a polynomial time approximation scheme （PTAS） and a fully polynomial time approximation scheme （FPTAS） for this problem.展开更多
During the past decade, increasing attention has been given to the development of meshless methods using radial basis functions for the numerical solution of Partial Differential Equations （PDEs）. A level set method...During the past decade, increasing attention has been given to the development of meshless methods using radial basis functions for the numerical solution of Partial Differential Equations （PDEs）. A level set method is a promising design tool for tracking, modelling and simulating the motion of free boundaries in fluid mechanics, combustion, computer animation and image processing. In the conventional level set methods, the level set equation is solved to evolve the interface using a capturing Eulerian approach. The solving procedure requires an appropriate choice of the upwind schemes, reinitialization, etc. Our goal is to include Multiquadric Radial Basis Functions （MQ RBFs） into the level set method to construct a more efficient approach and stabilize the solution process with the adaptive greedy algorithm. This paper presents an alternative approach to the conventional level set methods for solving moving-boundary problems. The solution was compared to the solution calculated by the exact explicit lime integration scheme. The examples show that MQ RBFs and adaptive greedy algorithm is a very promising calculation scheme.展开更多
In this paper a new approach for time discretization of an integro-differential equation of parabolic type is proposed. The methods are based on the backward-Euler and Crank-Nicolson Schemes but the memory and computa...In this paper a new approach for time discretization of an integro-differential equation of parabolic type is proposed. The methods are based on the backward-Euler and Crank-Nicolson Schemes but the memory and computational requirements are greatly reduced without assuming more regularities on the solution u.展开更多
One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for differ...One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth （multistep）, and second-order total variation diminishing （TVD） Runge-Kutta （multistage） time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic （rectangular, trapezoidal, triangular, and parabolic cross-sections） and non-prismatic （natural river） channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.展开更多
In fact,the popular semi-implicit time difference scheme of spectral model still includes some important linear terms using time explicit difference scheme,and the major terms are directly related to fast internal-and...In fact,the popular semi-implicit time difference scheme of spectral model still includes some important linear terms using time explicit difference scheme,and the major terms are directly related to fast internal-and external-gravity waves in the atmospheric forecasting equation. Additionally,due to using time difference on two terms at different time.the popular scheme artificially introduces unbalance between pressure gradient force and Coriolis force terms while numerically computing their small difference between large quantities.According to the computational stability analysis conducted to the linear term time difference scheme in simple harmonic motion equation,one improved semi-implicit time difference scheme is also designed in our study.By adopting a kind of revised time-explicit-difference scheme to these linear terms that still included in spectral model governing equations,the defect of spectral model which only partly using semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectral coefficients of prognostic equations,especially of Helmholtz divergence equation,can be worked out without any numerical iteration,the time-step (computation stability) can also be enlarged (enhanced) by properly introducing an adjustable coefficient.展开更多
In the field of optical interconnecting network and in super fast photonic computing system, the tree architecture and optical nonlinear materials can play a significant role. Nonlinear optical material may find impor...In the field of optical interconnecting network and in super fast photonic computing system, the tree architecture and optical nonlinear materials can play a significant role. Nonlinear optical material may find important uses in optical switching. Optical switch using nonlinear material makes it possible for one optical signal to control and switch another optical signal through nonlinear interaction in a material. In this communication such materials have been successfully exploited to design an all-optical tree-net architecture, which can be utilized for time division multiplexing scheme in all-optical domain.展开更多
We study the early work scheduling problem on identical parallel machines in order to maximize the total early work,i.e.,the parts of non-preemptive jobs that are executed before a common due date.By preprocessing and...We study the early work scheduling problem on identical parallel machines in order to maximize the total early work,i.e.,the parts of non-preemptive jobs that are executed before a common due date.By preprocessing and constructing an auxiliary instance which has several good properties,for any desired accuracy,we propose an efficient polynomial time approximation scheme with running time O(f(1/ε)n),where n is the number of jobs and f(1/ε)is exponential in 1/ε,and a fully polynomial time approximation scheme with running time O(1/ε^(2m+1)+n)when the number of machines is fixed.展开更多
Discontinuous deformation analysis(DDA) is a numerical method for analyzing the deformation of block system. It employs unified dynamic formulation for both static and dynamic analysis, in which the so-called kinetic ...Discontinuous deformation analysis(DDA) is a numerical method for analyzing the deformation of block system. It employs unified dynamic formulation for both static and dynamic analysis, in which the so-called kinetic damping is adopted for absorbing dynamic energy. The DDA dynamic equations are integrated directly by the constant acceleration algorithm of Newmark family integrators. In order to have an insight into the DDA time integration scheme, the performance of Newmark time integration scheme for dynamic equations with kinetic damping is systematically investigated, formulae of stability, bifurcation, spectral radius, critical kinetic damping and algorithmic damping are presented. Combining with numerical examples, recognition and suggestions of Newmark integration scheme application in the DDA static and dynamic analysis are proposed.展开更多
We consider several uniform parallel-machine scheduling problems in which the processing time of a job is a linear increasing function of its starting time.The objectives are to minimize the total completion time of a...We consider several uniform parallel-machine scheduling problems in which the processing time of a job is a linear increasing function of its starting time.The objectives are to minimize the total completion time of all jobs and the total load on all machines.We show that the problems are polynomially solvable when the increasing rates are identical for all jobs;we propose a fully polynomial-time approximation scheme for the standard linear deteriorating function,where the objective function is to minimize the total load on all machines.We also consider the problem in which the processing time of a job is a simple linear increasing function of its starting time and each job has a delivery time.The objective is to find a schedule which minimizes the time by which all jobs are delivered,and we propose a fully polynomial-time approximation scheme to solve this problem.展开更多
In this paper,we study high order discretization methods for solving the Maxwell equations on hybrid triangle-quad meshes.We have developed high order finite edge element methods coupled with different high order time...In this paper,we study high order discretization methods for solving the Maxwell equations on hybrid triangle-quad meshes.We have developed high order finite edge element methods coupled with different high order time schemes and we compare results and efficiency for several schemes.We introduce in particular a class of simple high order low dissipation time schemes based on a modified Taylor expansion.展开更多
In this paper,we consider the single-machine scheduling with step-deteriorating jobs and rejection.Each job is either rejected by paying a rejection penalty,or accepted and processed on the single machine,and the actu...In this paper,we consider the single-machine scheduling with step-deteriorating jobs and rejection.Each job is either rejected by paying a rejection penalty,or accepted and processed on the single machine,and the actual processing time of each accepted job is a step function of its starting time and the common deteriorating date.The objective is to minimize the makespan of the accepted jobs plus the total penalty of the rejected jobs.For the case of common deteriorating penalty,we first show that the problem is NP-hard in the ordinary sense.Then we present two pseudo-polynomial algorithms and a 2-approximation algorithm.Furthermore,we propose a fully polynomial time approximation scheme.For the case of common normal processing time,we present two pseudo-polynomial time algorithms,a 2-approximation algorithm and a fully polynomial time approximation scheme.展开更多
In this work we establish an existence theorem of regulated solutions for a class of Stieltjes equations which involve generalized fuemann kind of integrals. The general method spplied consists in considering the cont...In this work we establish an existence theorem of regulated solutions for a class of Stieltjes equations which involve generalized fuemann kind of integrals. The general method spplied consists in considering the continuous-time Stieltjes equation as limit of discrete processes. This approach will prove fruitful in the study of the controllability of Stieltjes systems, because it will be possible to get properties on the continuous time equation by transferring properties of the discrete ones.展开更多
文摘In this paper a new .mnultidimensional time series forecasting scheme based on the empirical orthogonal function (EOF) stepwise iteration process is introduced. The scheme is tested in a series of forecast experiments of Nino3 SST anomalies and Tahiti-Darwin SO index. The results show that the scheme is feasible and ENSO predictable.
基金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.
文摘Tow-phase flow mixed variational formulations of evolution filtration problems with seawater intrusion are analyzed. A dual mixed fractional flow velocity-pressure model is considered with an air-fresh water and a fresh water-seawater characterization. For analysis and computational purposes, spatial decompositions based on nonoverlapping multidomains, above and below the sea level, are variationally introduced with internal boundary fluxes dualized as weak transmission constraints. Further, parallel augmented and exactly penalized duality algorithms, and proximation semi-implicit time marching schemes, are established and analyzed.
基金supported by a grant from the French National Ministry of Education and Research(MENSR,19755-2005)
文摘A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell＇s equations is introduced and analyzed. The pro- posed method combines a centered approximation for the evaluation of fluxes at the in- terface between neighboring elements, with a Nth-order leap-frog time scheme. More- over, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwelrs equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high- order elements show the potential of the method.
文摘We sttidy the problem of scheduling n jobs on m parallel bounded batch machines to minimize the sum of squared machine loads. Each batch contains at most B jobs, and the processing time of a batch is equal to the longest processing time of the jobs in this batch. We prove this problem to be NP-hard. Furthermore, we present a polynomial time approximation scheme （PTAS） and a fully polynomial time approximation scheme （FPTAS） for this problem.
文摘During the past decade, increasing attention has been given to the development of meshless methods using radial basis functions for the numerical solution of Partial Differential Equations （PDEs）. A level set method is a promising design tool for tracking, modelling and simulating the motion of free boundaries in fluid mechanics, combustion, computer animation and image processing. In the conventional level set methods, the level set equation is solved to evolve the interface using a capturing Eulerian approach. The solving procedure requires an appropriate choice of the upwind schemes, reinitialization, etc. Our goal is to include Multiquadric Radial Basis Functions （MQ RBFs） into the level set method to construct a more efficient approach and stabilize the solution process with the adaptive greedy algorithm. This paper presents an alternative approach to the conventional level set methods for solving moving-boundary problems. The solution was compared to the solution calculated by the exact explicit lime integration scheme. The examples show that MQ RBFs and adaptive greedy algorithm is a very promising calculation scheme.
文摘In this paper a new approach for time discretization of an integro-differential equation of parabolic type is proposed. The methods are based on the backward-Euler and Crank-Nicolson Schemes but the memory and computational requirements are greatly reduced without assuming more regularities on the solution u.
文摘One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth （multistep）, and second-order total variation diminishing （TVD） Runge-Kutta （multistage） time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic （rectangular, trapezoidal, triangular, and parabolic cross-sections） and non-prismatic （natural river） channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.
基金The project is supported by the Beijing New Star Program of Science and Technology of China during 2001-2004 under Grant No.H013610330119.
文摘In fact,the popular semi-implicit time difference scheme of spectral model still includes some important linear terms using time explicit difference scheme,and the major terms are directly related to fast internal-and external-gravity waves in the atmospheric forecasting equation. Additionally,due to using time difference on two terms at different time.the popular scheme artificially introduces unbalance between pressure gradient force and Coriolis force terms while numerically computing their small difference between large quantities.According to the computational stability analysis conducted to the linear term time difference scheme in simple harmonic motion equation,one improved semi-implicit time difference scheme is also designed in our study.By adopting a kind of revised time-explicit-difference scheme to these linear terms that still included in spectral model governing equations,the defect of spectral model which only partly using semi-implicit integrating scheme can be overcome effectively.Moreover,besides all spectral coefficients of prognostic equations,especially of Helmholtz divergence equation,can be worked out without any numerical iteration,the time-step (computation stability) can also be enlarged (enhanced) by properly introducing an adjustable coefficient.
文摘In the field of optical interconnecting network and in super fast photonic computing system, the tree architecture and optical nonlinear materials can play a significant role. Nonlinear optical material may find important uses in optical switching. Optical switch using nonlinear material makes it possible for one optical signal to control and switch another optical signal through nonlinear interaction in a material. In this communication such materials have been successfully exploited to design an all-optical tree-net architecture, which can be utilized for time division multiplexing scheme in all-optical domain.
基金the National Natural Science Foundation of China(No.12071417)the Project for Innovation Team(Cultivation)of Yunnan Province.
文摘We study the early work scheduling problem on identical parallel machines in order to maximize the total early work,i.e.,the parts of non-preemptive jobs that are executed before a common due date.By preprocessing and constructing an auxiliary instance which has several good properties,for any desired accuracy,we propose an efficient polynomial time approximation scheme with running time O(f(1/ε)n),where n is the number of jobs and f(1/ε)is exponential in 1/ε,and a fully polynomial time approximation scheme with running time O(1/ε^(2m+1)+n)when the number of machines is fixed.
基金supported by the Fundamental Research Funds for Central Public Welfare Research Institutes(Grant No.CKSF2014053/CL)
文摘Discontinuous deformation analysis(DDA) is a numerical method for analyzing the deformation of block system. It employs unified dynamic formulation for both static and dynamic analysis, in which the so-called kinetic damping is adopted for absorbing dynamic energy. The DDA dynamic equations are integrated directly by the constant acceleration algorithm of Newmark family integrators. In order to have an insight into the DDA time integration scheme, the performance of Newmark time integration scheme for dynamic equations with kinetic damping is systematically investigated, formulae of stability, bifurcation, spectral radius, critical kinetic damping and algorithmic damping are presented. Combining with numerical examples, recognition and suggestions of Newmark integration scheme application in the DDA static and dynamic analysis are proposed.
基金This work was supported by the National Natural Science Foundation of China(Nos.11071142,11201259)the Natural Science Foundation of Shan Dong Province(No.ZR2010AM034)+1 种基金the Doctoral Fund of the Ministry of Education(No.20123705120001)We thank the two anonymous reviewers for their helpful and detailed comments on an earlier version of our paper.
文摘We consider several uniform parallel-machine scheduling problems in which the processing time of a job is a linear increasing function of its starting time.The objectives are to minimize the total completion time of all jobs and the total load on all machines.We show that the problems are polynomially solvable when the increasing rates are identical for all jobs;we propose a fully polynomial-time approximation scheme for the standard linear deteriorating function,where the objective function is to minimize the total load on all machines.We also consider the problem in which the processing time of a job is a simple linear increasing function of its starting time and each job has a delivery time.The objective is to find a schedule which minimizes the time by which all jobs are delivered,and we propose a fully polynomial-time approximation scheme to solve this problem.
基金This work was partially supported by the Agence Nationale de la Recherche,ANR-06-CIS6-0013.
文摘In this paper,we study high order discretization methods for solving the Maxwell equations on hybrid triangle-quad meshes.We have developed high order finite edge element methods coupled with different high order time schemes and we compare results and efficiency for several schemes.We introduce in particular a class of simple high order low dissipation time schemes based on a modified Taylor expansion.
基金supported by the National Natural Science Foundation of China(Nos.12271295 and 12001313)the Provincial Natural Science Foundation of Shandong(No.ZR2022MA019).
文摘In this paper,we consider the single-machine scheduling with step-deteriorating jobs and rejection.Each job is either rejected by paying a rejection penalty,or accepted and processed on the single machine,and the actual processing time of each accepted job is a step function of its starting time and the common deteriorating date.The objective is to minimize the makespan of the accepted jobs plus the total penalty of the rejected jobs.For the case of common deteriorating penalty,we first show that the problem is NP-hard in the ordinary sense.Then we present two pseudo-polynomial algorithms and a 2-approximation algorithm.Furthermore,we propose a fully polynomial time approximation scheme.For the case of common normal processing time,we present two pseudo-polynomial time algorithms,a 2-approximation algorithm and a fully polynomial time approximation scheme.
基金This work was developed at the Institute of Mathematics of the Czech Republic Academy of Sciences atPraha, with financial supp
文摘In this work we establish an existence theorem of regulated solutions for a class of Stieltjes equations which involve generalized fuemann kind of integrals. The general method spplied consists in considering the continuous-time Stieltjes equation as limit of discrete processes. This approach will prove fruitful in the study of the controllability of Stieltjes systems, because it will be possible to get properties on the continuous time equation by transferring properties of the discrete ones.
