The precursors of dipole blocking are obtained by a numerical approach based upon a quasi-geostrophic barotropic planetary- to synoptic-scale interaction model without topography and with a localized synoptic-scale wa...The precursors of dipole blocking are obtained by a numerical approach based upon a quasi-geostrophic barotropic planetary- to synoptic-scale interaction model without topography and with a localized synoptic-scale wave-maker. The optimization problem related to the precursors of blocking is formulated and the nonlinear optimization method is used to examine the optimal synoptic-scale initial field successfully. The results show that the prominent characteristics of the optimal synoptic-scale initial field are that the synoptic-scale wave train structures exist upstream of the incipient blocking. In addition, the large-scale low/high eddy-forcing pattern upstream of the incipient blocking is an essential precondition for the onset of dipole blocking.展开更多
The computation of the radiative transfer equation is expensive mainly due to two stiff terms:the transport term and the collision operator.The stiffness in the former comes from the fact that particles(such as photon...The computation of the radiative transfer equation is expensive mainly due to two stiff terms:the transport term and the collision operator.The stiffness in the former comes from the fact that particles(such as photons)travel at the speed of light,while that in the latter is due to the strong scattering in the optically thick region.We study the fully implicit scheme for this equation to account for the stiffness.The main challenge in the implicit treatment is the coupling between the spacial and angular coordinates that requires the large size of the to-be-inverted matrix,which is also ill-conditioned and not necessarily symmetric.Our main idea is to utilize the spectral structure of the ill-conditioned matrix to construct a pre-conditioner,which,along with an exquisite split of the spatial and angular dependence,significantly improve the condition number and allows a matrix-free treatment.We also design a fast solver to compute this pre-conditioner explicitly in advance.Our method is shown to be efficient in both diffusive and free streaming limit,and the computational cost is comparable to the state-of-the-art method.Various examples including anisotropic scattering and two-dimensional problems are provided to validate the effectiveness of our method.展开更多
基金This paper was supported by the Outstanding Youth Natural Science Foundation of China(Grant No.40325016)the Teaching and Research Award Program for Outstanding Young Teachers in Higher Education Institutions of MOE,PRC(TRAPOYT)the National Natural Science Foundation of China(Grant No.40175011).
文摘The precursors of dipole blocking are obtained by a numerical approach based upon a quasi-geostrophic barotropic planetary- to synoptic-scale interaction model without topography and with a localized synoptic-scale wave-maker. The optimization problem related to the precursors of blocking is formulated and the nonlinear optimization method is used to examine the optimal synoptic-scale initial field successfully. The results show that the prominent characteristics of the optimal synoptic-scale initial field are that the synoptic-scale wave train structures exist upstream of the incipient blocking. In addition, the large-scale low/high eddy-forcing pattern upstream of the incipient blocking is an essential precondition for the onset of dipole blocking.
基金The work of Q.Li is supported in part by a start-up fund from UW-Madison and National Science Foundation under the grant DMS-1619778The work of L.Wang is supported in part by the National Science Foundation under the grant DMS-1620135Both authors would like to express gratitude to the support from the NSF research network grant RNMS11-07444(KI-Net).We also thank Professors Shi Jin,Jim Morel and Cory Hauck for fruitful discussions.
文摘The computation of the radiative transfer equation is expensive mainly due to two stiff terms:the transport term and the collision operator.The stiffness in the former comes from the fact that particles(such as photons)travel at the speed of light,while that in the latter is due to the strong scattering in the optically thick region.We study the fully implicit scheme for this equation to account for the stiffness.The main challenge in the implicit treatment is the coupling between the spacial and angular coordinates that requires the large size of the to-be-inverted matrix,which is also ill-conditioned and not necessarily symmetric.Our main idea is to utilize the spectral structure of the ill-conditioned matrix to construct a pre-conditioner,which,along with an exquisite split of the spatial and angular dependence,significantly improve the condition number and allows a matrix-free treatment.We also design a fast solver to compute this pre-conditioner explicitly in advance.Our method is shown to be efficient in both diffusive and free streaming limit,and the computational cost is comparable to the state-of-the-art method.Various examples including anisotropic scattering and two-dimensional problems are provided to validate the effectiveness of our method.
