不同方案求解非线性化学动力方程组的比较

在一个简化的气相色化学模式的基础上，比较了ＨｙｂｒｉｄＱＳＳＡ，Ｓｋｌａｒｅｗ三种计算方案求解非线性化学动力学方程组的差异，并引入误差分配和线性组合两种质量守恒技术，分析它们对模拟结果的影响，研究结果表明，不同方案对有要烃类浓度的计算几乎没有影响，但对其它气态物浓度有一定的影响，用ＱＳＳＡ和Ｓｌａｒｅｗ方案预测的结果更为相近，而采用Ｈｙｂｒｉｄ方案得到的气态物浓度的平衡时间略为提前，从计算效率来看     

非线性大气化学动力学方程组数值解法的比较

大气化学模式中描述化学反应动力学的是一组高度耦合、刚性、非线性的常微分方程组,高效精确的计算方法在空气质量模拟预测、气候化学相互作用研究中尤为重要。本文介绍了求解此类非线性常微分方程组(ODEs)的七种常见的数值计算方法,包括QSSA、MQSSA、HYBRID、MHYBRID、LSODE、GongCho以及TWOSTEP方法,并基于同一反应机理进行24 h的数值模拟,通过对模拟结果的讨论分析,比较了各种方法的优劣。研究结果表明,就计算效率而言,QSSA方法最为省时。就计算精度而言,LSODE方法最优。若综合考虑计算精度和效率,HYBRID、MHY-BRID是比较适中的数值解法。     

一种非线性大气化学动力学方程组的新算法

提出了一种求解非线性大气化学动力学方程组的新算法,即预测／校正拟稳定近似法（ＰＣＱＳＳＡ）该方案集中了传统的ＱＳＳＡ方案和Ｈｙｂｒｉｄ方案的优点,具有较高的计算效率和精度,对于复杂的大气化学模式和区域污染模式较适用。     

强迫耗散非线性临界层动力学的初步研究

研究了强迫耗散非线性临界层动力学，得到了强迫耗散非线性临界层的解析解。解析解表明，在含有强迫耗散的条件下，非线性临界仍有“Ｋｅｌｖｉｎ猫眼”，当时间很大时，“猫眼”趋于定常，强迫耗散因子对非线性临界层激发的“Ｋｅｌｖｉｎ猫眼”的涡度场有重要的调制作用。     

非线性大气动力学的进展

总结了自新中国成立以来中国科学院大气物理研究所在非线性大气动力学领域的进展,其中主要介绍了在非线性适应过程、非线性稳定与不稳定、全局分析理论、可预报性、低阶谱方法与多平衡态动力学、非线性波动、波流相互作用、阻塞高压非线性动力学、中小尺度非线性动力学等几个方面的研究成果及最新进展.     

大气动力学方程组的定性理论及其应用

基于完整的湿大气动力学方程组,利用无穷维动力系统的新理论和新方法,系统讨论了强迫耗散的非线性大气系统的定性理论及其应用。将完整的强迫耗散非线性湿大气动力学方程组化为Hilbert空间中一个等价的算子方程,研究了算子的性质及其物理意义,在此基础上得到湿大气系统全局吸引子的存在性定理,揭示出系统向外源的非线性适应特征,并把结果推广到有地形动力作用和非定常外源强迫的情形。同时探讨了大气方程组惯性流形的存在,大气多平衡态产生的根源以及强迫、耗散和非线性对系统解的渐近行为的影响。在理论结果的基础上,提出强迫耗散的非线性动力系统中存在三类时间边界层、方程组简化准则、分解算法的算子约束原则以及支撑吸引子基底的少数自由度的构造方法,探讨了理论在非线性发展方程差分格式的设计和计算稳定性分析、多平衡态的数值分析、数值模式延伸预报的改进、短期气候预测以及一类中尺度系统分析与预测中的应用,指出描述长期过程动力学模式的必备条件,给出初值与模式相协调的合理解释。最后,对今后的研究方向作了展望。     

中国年降水气候振动的非线性动力学特征

从观测资料出发,利用反演方法,获得了描述我国17个站年降水气候振动的二阶非线性动力学方程组中的系数,并在无外源强迫的情况下对该非线性动力学方程组作了定性的分析,对17个站的年降水气候振动进行气候模拟试验。结果表明:我国大部分地区的年降水气候振动存在弱的非线性振动特征,17个站的年降水气候振动都可以看作是一个正阻尼的动力系统,在运动变化过程中,自身能量的变化dEdt<0,17个站年降水气候振动的固有频率的变化范围为0.6≤ω0≤1.417,其固有周期的变化范围为10.4≥T≥4.4年。从空间分布特征上可以看到,其固有周期的长短分布有明显的东西向特征。在17个站中,兰州、成都、重庆、南宁、福州、上海、温州地区的年降水气候振动是一个类似于具有渐软弹簧特性的非线性动力系统,而昆明、长春、哈尔滨、北京、广州、西安、厦门、长沙、南京、杭州地区的年降水气候振动是一个类似于具有渐硬弹簧特性的非线性动力系统,从其空间分布特征上看也有东西向特征。全国17个站年降水气候振动的动力系统都存在一个稳定的焦点,在无外源强迫作用时,各站的年降水气候振动将逐渐趋向于各自的稳定的焦点,但都不存在极限环。在仅考虑降水自身振动的情况下,利用我们反演得到的非线性动力学方程对年降水气候振动做模拟,位于内陆地区的模拟情况可以通过显著性水平检验,而靠近海洋的地区模拟效果没有内陆好。这一结果说明在降水气候振动中,靠近海洋地区的外源强迫影响作用比位于内陆地区大。     

Three-Step Difference Scheme for Solving Nonlinear Time-Evolution Partial Differential Equations

In this paper, a special three-step difference scheme is applied to the solution of nonlinear time-evolution equations, whose coefficients are determined according to accuracy constraints, necessary conditions of square conservation, and historical observation information under the linear supposition. As in the linear case, the schemes also have obvious superiority in overall performance in the nonlinear case compared with traditional finite difference schemes, e.g., the leapfrog（LF） scheme and the complete square conservation difference（CSCD） scheme that do not use historical observations in determining their coefficients, and the retrospective time integration（RTI） scheme that does not consider compatibility and square conservation. Ideal numerical experiments using the one-dimensional nonlinear advection equation with an exact solution show that this three-step scheme minimizes its root mean square error（RMSE） during the first 2500 integration steps when no shock waves occur in the exact solution, while the RTI scheme outperforms the LF scheme and CSCD scheme only in the first 1000 steps and then becomes the worst in terms of RMSE up to the 2500th step. It is concluded that reasonable consideration of accuracy, square conservation, and historical observations is also critical for good performance of a finite difference scheme for solving nonlinear equations.     

两种风场分解方法计算动能收支对比研究

使用两种风场分解方法，计算了旋转风和辐散风功能收支方程各项时空平均值，结果表明两者差异不大，具有一定程度的可比性。     

非线性发展方程的非守恒格式的计算稳定性问题

针对非线性发展方程的非守恒格式,以一维浅水波方程为例,对非守恒格式的计算稳定性进行了研究分析,探讨了非线性发展方程的非守恒格式与初值的关系.理论分析和数值试验表明,在格式结构已经确定的情况下,非守恒格式的计算稳定性主要由初值的形式所决定.     

Kinetic Energy Budget Equations of Rotational and Divergent Flow in Terrain-following Coordinates

In this study, kinetic energy budget equations of rotational and divergent flow in pressure coordinates are derived on terrain-following coordinates. The new formulation explicitly shows the terrain effects and can be applied directly to model-simulated dynamic and thermodynamic fields on the model’s original vertical grid. Such application eliminates interpolation error and avoids errors in virtual weather systems in mountainous areas. These advantages and their significance are demonstrated by a numerical study in terrain-following coordinates of a developing vortex after it moves over the Tibetan Plateau in China.     

多普勒雷达资料格点化方案的比较研究

利用高阶多项式拟合实际雷达速度观测的方法生成模拟资料。据此比较了3种雷达资料格点化的插值方案。它们是先锥面后垂直内插方案(CVI)。3维Barnes插值方案(3D—Branes)和变分分析方案(VAR)。结果表明。3D-Barnes方案明显优于其它两种,特别是在中层;CVI方案的误差略高于前者,但是光滑性较差。在高层和低层空缺区较大;VAR方案的误差最大。     

Kinetic equations and stationary energy spectra of weakly nonlinear internal gravity waves

An ensemble of random-phase internal gravity waves is considered in the dynamical framework of the Euler–Boussinesq equations. For flows with zero mean potential vorticity, a kinetic equation for the mean spectral energy density of the waves is obtained under hypothesis of Gaussian statistics with zero correlation length. Stationary scaling solutions of this equation are found for almost vertically propagating waves. The resulting spectra are anisotropic in vertical and horizontal wave numbers. For flows with small but non-zero mean potential vorticity, under the same statistical hypothesis applied to the wave part of the flow, it is shown that the vortex part and the wave part decouple. The vortex part obeys a limiting slow dynamics equation exhibiting vertical collapse and layering which may contaminate the wave-part spectra. Relation of these results to the in situ atmospheric measurements and previous work on oceanic gravity waves is discussed.     

An Analysis and Implications of Alternative Methods of Deriving the Density (WPL) Terms for Eddy Covariance Flux Measurements

We explore some of the underlying assumptions used to derive the density or WPL terms (Webb et al. (1980) Quart J Roy Meteorol Soc 106:85–100) required for estimating the surface exchange fluxes by eddy covariance. As part of this effort we recast the origin of the density terms as an assumption regarding the density fluctuations rather than as a (dry air) flux assumption. This new approach, which is similar to the expansion/compression approach of Liu (Boundary-Layer Meteorol 115:151–168, 2005), eliminates the dry-air mean advective vertical velocity from the development of the WPL terms and allows us to directly compare Liu’s assumptions for deriving the WPL terms with the analogous assumptions appropriate to the original expression of the WPL terms. We suggest, (i) that the main difference between these two approaches lies in the interpretation of the turbulent exchange flux, and (ii) that the original WPL formulation is the more appropriate approach. Given the importance of the WPL terms to accurate and reliable measurements of surface exchange fluxes, a careful analysis of their origins and their proper mathematical expression and interpretation is warranted.     

Nonlinear pulsations of horizontal jets

A shallow-water model with horizontally nonuniform density is used to study the dynamics of jet flows that arise under the influence of buoyancy and the Coriolis force. Within this approach, the jet is described by a self-similar compactly-localized solution and interpreted as a band of shear flow having a temperature contrast with the ambient fluid. In addition to stationary states, the dynamics of such jets admit cyclonic rotation with a constant angular velocity and transverse nonlinear pulsations. The phase portrait corresponding to this model shows that regimes with pulsating jets develop along closed trajectories bounded by the separatrix loop. The theory predicts that the period for warm jet pulsations is longer than the inertial oscillation period caused by the Earth’s rotation, while for cold jet pulsations, it is shorter. Thus, only warm jets can have a noticeable effect on the atmospheric dynamics in the synoptic range. In particular, they may well be responsible for additional spectral peaks that appear in this range of wind speed fluctuations.     

地形区两种典型气压梯度力计算方法的比较

试验比较了两种较为典型的地形区气压梯度力计算方法：一是先在准水平的等压面上计算,再垂直插值到地形坐标面上的方法,二是差-微-差方法。初始化分析和预报结果检验都表明,在减轻地形影响上,前者的优势在大气上层,后者则在中下层。