In the light of the problem of oil pollution brought about by ships, in this paper we present the concept of backward tracing oil spills. In the course of backward calculation of the two-dimensional convection & d...In the light of the problem of oil pollution brought about by ships, in this paper we present the concept of backward tracing oil spills. In the course of backward calculation of the two-dimensional convection & diffusion equation, on the one hand, the advection term itself has the strong unilateral property, which means information in the upper reaches is transmitted downstream via the advection term; on the other hand, because of the opposite direction of calculation, it is essential for information to be conveyed upstream by means of the advection term. In addition, unlike that in the forward calculation, the diffusion term in the backward calculation is prone to accumulate errors, and thus renders the whole scheme unstable. Therefore, we adopt the central difference to deal with both the convectional term and the diffusion term. By examining two practical examples (1) under the unlimited boundary condition, and (2) under the limited boundary condition, it is proven that this method could achieve fundamentally satisfactory results not only in the open ocean but also in the closed or semi-closed bay.展开更多
A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missi...A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missing and overflowing".A sensitive heat capacity model is introduced through which the calculation errors are analyzed.Then the equation using the self-adjusted time step is presented and solved using finite differences.Through this equation,the time needed for a space cell to reach the phase change point temperature is calculated.Using this time allows the time step to be adjusted so that errors caused by "phase change missing and overflowing" are successfully eliminated.Above all,the obvious features of this method are an accelerated rate for adjusting the time step and simplifing the computations.An actual example proves that this method can accurately calculate the temperature fields during soil freezing and thawing.It is an improvement over traditional methods and can be widely used on complicated multi-dimensional phase change problems.展开更多
基金funded by NSFC 40076005 and Frontier Innovation Project L390221103 from the Chinese Academy of Sciencesthe financial support from the National Tenth Five-Year Key Project H57022113.
文摘In the light of the problem of oil pollution brought about by ships, in this paper we present the concept of backward tracing oil spills. In the course of backward calculation of the two-dimensional convection & diffusion equation, on the one hand, the advection term itself has the strong unilateral property, which means information in the upper reaches is transmitted downstream via the advection term; on the other hand, because of the opposite direction of calculation, it is essential for information to be conveyed upstream by means of the advection term. In addition, unlike that in the forward calculation, the diffusion term in the backward calculation is prone to accumulate errors, and thus renders the whole scheme unstable. Therefore, we adopt the central difference to deal with both the convectional term and the diffusion term. By examining two practical examples (1) under the unlimited boundary condition, and (2) under the limited boundary condition, it is proven that this method could achieve fundamentally satisfactory results not only in the open ocean but also in the closed or semi-closed bay.
基金Project 2006G1662-00 supported by the Key Science and Technology Project of Heilongjiang Province
文摘A new variable time step method,which is called the backwards calculating time step method,is presented in this paper.It allows numerical simulation of soil freezing and thawing while avoiding "phase change missing and overflowing".A sensitive heat capacity model is introduced through which the calculation errors are analyzed.Then the equation using the self-adjusted time step is presented and solved using finite differences.Through this equation,the time needed for a space cell to reach the phase change point temperature is calculated.Using this time allows the time step to be adjusted so that errors caused by "phase change missing and overflowing" are successfully eliminated.Above all,the obvious features of this method are an accelerated rate for adjusting the time step and simplifing the computations.An actual example proves that this method can accurately calculate the temperature fields during soil freezing and thawing.It is an improvement over traditional methods and can be widely used on complicated multi-dimensional phase change problems.