A New Method for Calculating the Wind Speed Distribution of a Moving Tropical Cyclone

Based on gradient wind equations,including frictional force,and considering the effect of the movement of a tropical cyclone on wind speed,the Fujita Formula is improved and further simplified,and the numerical scheme for calculating the maximum wind speed radius and wind velocity distribution of a moving tropical cyclone is derived.In addition,the effect of frictional force on the internal structure of the tropical cyclone is discussed.By comparison with observational data,this numerical scheme demonstrate...     

水平井开采煤层气藏生产动态计算新数学模型

基于朗格缪尔模型、菲克第二扩散定律和达西定律,综合考虑重力和气体压缩效应的影响,推导出描述煤层气和水渗流规律的基本偏微分方程组(并给出定解条件),且利用有限差分法对该数学模型进行了数值求解。计算结果表明,该模型能够有效模拟煤层气水平井的生产动态。通过对比是否考虑非稳态扩散和气体压缩效应时不同裂隙渗透率下的产气动态可知,忽略非稳态扩散和压缩效应,将人为缩短达到气体峰值产量的时间,增大产气峰值;将降低气体膨胀能,低估气井后期产能。煤岩裂隙渗透率越高,气井产能受渗透率回归效应影响越大,峰值产量的持续时间越短且产量下降速度快。因此在煤层气藏开采过程中,应根据气井的产气动态调整井底流压,建立合理的储层压力系统,以降低渗透率回归效应带来的不利影响。     

水文重力效应改正中一维地下水模拟算法的对比

基于一维地下水渗透方程详细推导其有限差分解算过程,引入不同于显式差分的隐式差分和中心差分格式,对比分析不同差分格式对地下水模拟结果及其相应地下水重力效应的影响,并对其中的层间参数取值和非线性方程的线性化问题进行探讨。结果表明,在日本Isawa扇形地区超导台站,不同层间参数加权公式能够引起最大约0.15 μGal的重力效应差异,影响在1.9%以内;不同差分格式和线性化方法组合形式能够引起最大约0.12 μGal的重力效应差异,影响在1.5%以内。     

小振幅海洋内波的演变、破碎和所致混合

利用基于谱方法和MPI并行运算的数值模式SpectralModel,直接数值模拟了三维小振幅海洋内波的演变、破碎和所致湍流混合,指出导致其不稳定而破碎的为PSI(parametric subharmonic instability)机制;对于内波破碎所致的湍流混合过程,分析了跨等密度面扩散系数kρ、混合效率γ、浮力通量谱、动能谱以及势能谱等统计性质:内波破碎前,kρ和γ保持低值水平,浮力通量谱值为负,且集中在低波数段;内波破碎后,kρ和γ迅速增大,最大值分别约为0.9×103m2/s和0.18,浮力通量谱值在低波数段为负值,在高波数段为正值,这是因为层化湍流中势能向小尺度运动传递和动能向小尺度运动传递相比更为有效。在内波破碎、强湍流混合阶段,势能谱存在一谱段满足kz3律,P(kz)=0.2N2kz3。此外,与二维模拟结果相比较,导致内波不稳定而破碎的均为PSI机制,kρ、浮力通量谱、势能谱变化趋势大体一致;但三维数值实验中,内波破碎时间提前,湍流衰减加快;KE谱在高波数部分下降速度相对减小,更接近于kz3。     

Long‐term annual surface heat and water balances over Canada and the United States South of 60°N: Reconciliation of precipitation,run‐off and temperature fields

Abstract

The analysis compares the observed field of run‐off (assumed correct) with adjusted precipitation over North America (as amended by den Hartog and LeDrew over Canada) and derives the principal hydroclimatological ratios for each five‐degree latitude‐longitude square. The amended precipitation field yields values of the Budyko dry ness index close to values suggested by the vegetation distribution. The Priestley‐Taylor parameter, α, lies between unity (equilibrium) and potential (1.26) values over much of humid North America, but exceeds these values in the northwest Pacific squares, where advective heating may be the cause. Other regions of strong seasonal advective heating (e.g. the Great Plains) do not appear to influence the distribution strongly. A weighted convective forcing temperature is derived, varying from 298 K in the extreme south to below 285 K in the north. This function (and the Bowen ratio) achieve improbable values in northern Labrador‐ Ungava. The precipitation, run‐off and net radiation régimes appear still to be out of balance in these squares. An adjustment of either precipitation or net radiation by about a tenth corrects the imbalance, but the method is not capable of deciding which field (or both) is in error. Over the rest of the continent the adjusted precipitation field now appears to be in balance with observed run‐off and temperature distributions.     

水质大涡模拟数学模型研究

建立了基于弱可压缩流基础上的水质大涡模拟数学模型及污染物质扩散系数与平均运动的关系。数值计算采用有限体积法和预测校正法,固壁边界条件采用"部分滑移条件"。并应用该模型对重庆嘉陵江、长江两江交汇段的流场,污染物质浓度分布进行了计算,计算结果与实测值比较吻合。研究表明,横向扩散与流动状态有关,当支流流量相对较小时,横向扩散较弱,形成岸边污染带,当支流流量相对较大时,横向扩散较强,污染带向江心移动。     

McCreary和Fickian湍扩散形式差异的分析

用矩阵方法,求解McCreary和Fickian铅垂向湍扩散形式下,无界海洋条件时在风应力和深层浮力作用下的线性解。得出不同的湍扩散系数A情况下McCreary与Fickian湍扩散形式的流场结构。指出随湍扩散系数的增加,McCreary和Fickian湍扩散形式下u分量结构及其量值差异增大。在Fickian湍扩散形式下,能量穿越密度跃层的传递能力较强。当A=0.1cm~2s~(-3)时,在风应力和深层浮力作用下两种湍扩散形式造成流场的差异很小。提出在小湍扩散系数下,用McCreary湍扩散形式研究热带海洋环流,其结果和Fickian形式应相当一致。     

Ultrahigh temperature deformation microstructures in felsic granulites of the Napier Complex, Antarctica

Detailed electron microscope and microstructural analysis of two ultrahigh temperature felsic granulites from Tonagh Island, Napier Complex, Antarctica show deformation microstructures produced at  1000 °C at 8–10 kbar. High temperature orthopyroxene (Al 7 wt.% and  11 wt.%), exhibits crystallographic preferred orientation (CPO) and frequent subgrain boundaries which point to dislocation creep as the dominating deformation mechanism within opx. Two different main slip systems are observed: in opx bands with exclusively opx grains containing subgrain boundaries with traces parallel to [010] and a strong coupling of low angle misorientations (2.5°–5°) with rotation axes parallel to [010] the dominating slip system is (100)[001]. Isolated opx grains and grain clusters of 2–5 grains embedded in a qtz–fsp matrix show an additional slip system of (010)[001]. The latter slip system is harder to activate. We suggest that differences in the activation of these slip systems is a result of higher differential stresses imposed onto the isolated opx grains and grain clusters. In contrast to opx, large qtz grains (up to 200 μm) show random crystallographic orientation. This together with their elongate and cuspate shape and the lack of systematic in the rotation axes associated with the subgrain boundaries is consistent with diffusion creep as the primary deformation mechanism in quartz.Our first time detailed microstructural observations of ultrahigh temperature and medium to high pressure granulites and their interpretation in terms of active deformation mechanisms give some insight into the type of rheology that can be expect at lower crustal conditions. If qtz is the mineral phase governing the rock rheology, Newtonian flow behaviour is expected and only low differential stress can be supported. However, if the stress supporting mineral phase is opx, the flow law resulting from dislocation creep will govern the rheology of the rock unit; hence, an exponential relationship between stress and strain rate is to be expected.     

Three‐parameter discontinuous distributions for hydrological samples with zero values

A consistent approach to the frequency analysis of hydrologic data in arid and semiarid regions, i.e. the data series containing several zero values (e.g. monthly precipitation in dry seasons, annual peak flow discharges, etc.), requires using discontinuous probability distribution functions. Such an approach has received relatively limited attention. Along the lines of physically based models, the extensions of the Muskingum‐based models to three parameter forms are considered. Using 44 peak flow series from the USGS data bank, the fitting ability of four three‐parameter models was investigated: (1) the Dirac delta combined with Gamma distribution; (2) the Dirac delta combined with two‐parameter generalized Pareto distribution; (3) the Dirac delta combined with two‐parameter Weibull (DWe) distribution; (4) the kinematic diffusion with one additional parameter that controls the probability of the zero event (KD3). The goodness of fit of the models was assessed and compared both by evaluation of discrepancies between the results of both estimation methods (i.e. the method of moments (MOM) and the maximum likelihood method (MLM)) and using the log of likelihood function as a criterion. In most cases, the DWe distribution with MLM‐estimated parameters showed the best fit of all the three‐parameter models. Copyright © 2005 John Wiley & Sons, Ltd.