A Godunov-Type Scheme for Atmospheric Flows on Unstructured Grids: Euler and Navier-Stokes Equations

In recent years there has been a growing interest in using Godunov-type methods for atmospheric flow problems. Godunov's unique approach to numerical modeling of fluid flow is characterized by introducing physical reasoning in the development of the numerical scheme (van Leer, 1999). The construction of the scheme itself is based upon the physical phenomenon described by the equation sets. These finite volume discretizations are conservative and have the ability to resolve regions of steep gradients accurately, thus avoiding dispersion errors in the solution. Positivity of scalars (an important factor when considering the transport of microphysical quantities) is also guaranteed by applying the total variation diminishing condition appropriately. This paper describes the implementation of a Godunov-type finite volume scheme based on unstructured adaptive grids for simulating flows on the meso-, micro- and urban-scales. The Harten-Lax-van Leer-Contact (HLLC) approximate Riemann solver used to calculate the Godunov fluxes is described in detail. The higher-order spatial accuracy is achieved via gradient reconstruction techniques after van Leer and the total variation diminishing condition is enforced with the aid of slope-limiters. A multi-stage explicit Runge-Kutta time marching scheme is used for maintaining higher-order accuracy in time. The scheme is conservative and exhibits minimal numerical dispersion and diffusion. The subgrid scale diffusion in the model is parameterized via the Smagorinsky-Lilly turbulence closure. The scheme uses a non-staggered mesh arrangement of variables (all quantities are cell-centered) and requires no explicit filtering for stability. A comparison with exact solutions shows that the scheme can resolve the different types of wave structures admitted by the atmospheric flow equation set. A qualitative evaluation for an idealized test case of convection in a neutral atmosphere is also presented. The scheme was able to simulate the onset of Kelvin-Helmholtz type instability and shows promise in simulating atmospheric flows characterized by sharp gradients without using explicit filtering for numerical stability.     

自适应网格技术在数值模式中的应用研究 II.二维问题

用一个具有解析解的二维动力锋生成过程算例,对比自适应网格方案和固定网格方案的优劣,结果表明:同等误差要求下,固定网格方案的网格点数为自适应网格方案的3倍。从时间演化上看,自适应网格对温带气旋的预报,在同等误差要求下,可比固定网格的预报延长10小时以上。文章对加权系数选取及对解的影响进行了分析,从几个切面的计算结果展示了自适应网格对网格的优良安排并能抓住锋面特征结构。文章分析了光滑性、正交性对结果的影响。结果表明:网格的光滑性影响有一个优化值;在网格适当安排情况下,要适度考虑正交性。     

利用泰森多边形和格网平滑的人口密度空间化研究——以徐州市为例

在运用泰森多边形作边界替代得到徐州市区人口密度的分布图后,采用GCAWI法实现了徐州市区人口数据的空间化,提出了采用大小两种格网错置的平滑方法,并运用此方法得到徐州市2000年和2005年的500 m×500 m的格网人口密度分布图。结果表明:①通过对比两个年份的人口密度分布图可以发现,徐州市区人口密度在2000年~2005年期间具有南移东扩的趋势;②利用大小两种格网错置平滑的方法实现了对人口密度空间化与平滑两个功能,平滑后的结果表明,徐州市区人口密度的空间分布具有块状集聚、轴状延伸的格局,整体来看呈现出"大"字形的态势;③通过与直接运用核函数内插得到的人口密度分布图的对比可以发现,前者对反映人口密度分布的细部差异比较敏感,而后者仅能反映出人口密度空间分布的整体格局。     

P–T modelling of the andalusite–kyanite–andalusite sequence and related assemblages in high-Al graphitic pelites. Prograde and retrograde paths in a late kyanite belt in the Variscan Iberia

The exceptional andalusite–kyanite–andalusite sequence occurs in Al‐rich graphitic slates in a narrow pelite belt on the hangingwall of a ductile normal fault in NW Variscan Iberia. Early chiastolite is replaced by Ky–Ms–Pg aggregates, which are overgrown by pleochroic andalusite near granites intruded along the fault. Slates plot in AKFM above the chloritoid‐chlorite tie‐line. Their P–T grids are modelled with Thermocalc v2.7 and the 1998 databases in the NaKFMASH and KFMASH systems. The univariant reaction Ctd + And/Ky = St + Chl + Qtz + H₂O ends at progressively lower pressure as F/FM increases and A/AFM decreases, shrinking the assemblage Cld–Ky–Chl, and opening a chlorite‐free Cld–Ky trivariant field on the low temperature reaction side. This modelling matches the observed absence of chlorite in high F/FM rocks, which is restricted to low pressure in the andalusite stability field. The P–T path deduced from modelling shows a first prograde event in the andalusite field followed by retrogression into the kyanite field, most likely coupled with a slight pressure increase. The final prograde evolution into the andalusite field can be explained by two different prograde paths. Granite intrusion caused the first prograde part of the path with andalusite growth. The subsequent thermal relaxation, together with a_H2O decrease, generated the retrograde andalusite–kyanite transformation, plus chlorite consumption and chloritoid growth. This transformation could have been related to folding in the beginning, and aided later by downthrowing due to normal faulting. Heat supplied by syntectonic granite intrusion explains the isobaric part of the path in the late stages of evolution, causing the prograde andalusite growth after the assemblage St–Ky–Chl. Near postectonic granites, a prograde path with pressure decrease originated the assemblage St–And–Chl.     

地学现象的分维与标度关系——以中国大陆山系为例

分形理论产生以后,已经在许多领域里取得了较为广泛的应用,但是在分形维数的计算方面仍存在一定的问题,最根本的问题之一是大量相关研究往往忽略了标度问题的影响,而将计算出的分维结果随意比较。本文则是以中国大陆山系为例,对地学现象的分维与标度的关系进行了研究,得出:(1)中国大陆山系的分形性质是客观存在的,其分形性质是不以测算标度而出现变化的,是始终存在的;(2)虽然在不同的测算标度下,中国大陆山系的分形性质始终是客观存在的,但是随着测算标度范围的变化,中国大陆山系的分维也会相应地出现一定的变化,即中国大陆山系客观存在的分形性质是绝对的,而分维则是一个相对的值,是随着测算标度而出现变化的;(3)分维不仅随着测算标度的变化而出现一定的变化,而且,分维是测算平均标度的函数,即D=f(r)(式中:r为测算平均标度,D为分维)关系成立;(4)因为分维值会随着不同标度范围的变化而出现相应的变化,所以,如果要进行分维之间的比较,则应选择相同的标度范围和标度值点。     

Global sustainability,innovation and governance dynamics of national smart electricity meter transitions

Smart electricity meters are a central feature of any future smart grid, and therefore represent a rapid and significant household energy transition, growing by our calculations from less than 23.5 million smart meters in 2010 to an estimated 729.1 million in 2019, a decadal growth rate of 3013%. What are the varying economic, governance, and energy and climate sustainability aspects associated with the diffusion of smart meters for electricity? What lessons can be learned from the ongoing rollouts of smart meters around the world? Based on an original dataset twice as comprehensive as the current state of the art, this study examines smart meter deployment across 41 national programs and 61 subnational programs that collectively target 1.49 billion installations involving 47 countries. In addition to rates of adoption and the relative influence of factors such as technology costs, we examine adoption requirements, modes of information provision, patterns of incumbency and management, behavioral changes and energy savings, emissions reductions, policies, and links to other low-carbon transitions such as energy efficiency or renewable energy. We identify numerous weak spots in the literature, notably the lack of harmonized datasets as well as inconsistent scope and quality within national cost-benefit analyses of smart meter programs. Most smart meters have a lifetime of only 20 years, leading to future challenges concerning repair, care, and waste. National-scale programs (notably China) account for a far larger number of installations than subnational ones, and national scale programs also install smart meters more affordably, i.e. with lower general costs. Finally, the transformative effect of smart meters may be oversold, and we find that smart electricity meters are a technology that is complementary, rather than disruptive or transformative, one that largely does not challenge the dominant practices and roles of electricity suppliers, firms, or network operators.     

机载LiDAR旁向重叠区域冗余数据的消除算法

本文基于点云格网化的思想,提出了一种去除机载LiDAR航带间旁向重叠区域冗余数据的算法。该算法使用规则方格表示点云覆盖的区域,通过给网格赋值的方式提取重叠区域的冗余点云,最后删除靠近条带边缘的冗余点云。经过实验测试,该方法能够快速有效地去除重叠区的冗余数据;无需依赖航迹信息,因此也适用于地面激光点云的消冗。     

Detecting and Enumerating New Building Structures Utilizing Very‐High Resolution Imaged Data and Image Processing

Abstract

Information on the number and type of new building structures is required by urban and transportation planners and the real estate industry. The goal of this paper is to explore the potential of high resolution imagery for meeting public and private sector demands for information on new buildings. The value of 1 m, 5 m, and 10 m panchromatic and 1 m color scanned aerial photography images acquired in 1997 and 1998 for a study area within the City of San Diego, California is assessed for general change detection and building enumeration. Both semi‐automated and interactive change‐detection approaches are evaluated. We demonstrate that interactive, visual‐based approaches appear to be the most accurate (within 1% of actual count) and efficient approach for generating information on the number of new buildings associated with single family residential land use. More automated approaches to detecting and enumerating image microfeatures may be useful as enhancements for visual‐based assessments and may be practical in areas composed mostly of large buildings associated with commercial and industrial land use. The highest accuracy for automated approaches was an undercounting of 11% for residential buildings and overcounting of 20% for those associated with commercial and industrial land use.