Coupled analysis of nonlinear sloshing and ship motions

A coupled numerical model considering nonlinear sloshing flows and the linear ship motions has been developed based on a boundary element method. Hydrodynamic performances of a tank containing internal fluid under regular wave excitations in sway are investigated by the present time-domain simulation model and comparative model tests. The numerical model features well the hydrodynamic performance of a tank and its internal sloshing flows obtained from the experiments. In particular, the numerical simulations of the strong nonlinear sloshing flows at the natural frequency have been validated. The influence of the excitation wave height and wave frequency on ship motions and internal sloshing has been investigated. The magnitude of the internal sloshing increases nonlinearly as the wave excitation increases. It is observed that the asymmetry of the internal sloshing relative to still water surface becomes more pronounced at higher wave excitation. The internal sloshing-induced wave elevation is found to be amplitude-modulated. The frequency of the amplitude modulation envelope is determined by the difference between the incident wave frequency and the natural frequency of the internal sloshing. Furthermore, the coupling mechanism between ship motions and internal sloshing is discussed.     

贵州晴隆中二叠统大厂层砾岩成因研究

贵州晴隆中二叠统大厂层中—上部产出一套呈囊状、透镜状分布的砾岩层。对大厂层砾岩沉积序列、结构、构造特征和矿物组合进行研究,结果表明:沉积序列为"反粒序",砾岩成分单一,基本为玄武质,是由峨眉山玄武岩直接改造而成,未发育河道相沉积特征;大厂层砾石具特殊的"泥化边"构造,是玄武岩在海水中大量水解的直接证据,填隙物中的青磐岩化矿物组合表明砾岩遭受了低温热液的改造。填隙物具有接近凝灰岩的Zr/Hf(30.7~43.4,均值38.0),揭示玄武岩在改造形成砾石的过程中有火山碎屑(火山灰)参与成岩;砾石和填隙物相对玄武岩贫Na(Na2O)富K(K2O),三者具有相似的稀土配分模式,其中∑REE呈规律性变化(玄武岩最高,砾石次之,填隙物中最低)。沉积背景分析认为,峨眉地幔柱作用使地壳发生了差异抬升,抬升一侧暴露遭受剥蚀形成不整合面,相对沉降一侧继续接受沉积,大厂层正是峨眉山玄武岩前锋带在局限海盆边缘(沉降侧),于茅口灰岩之上继续沉积的产物。大厂层砾岩形成机制是:炽热玄武岩流在海水中急剧爆裂,同时火山碎屑参与沉积,形成的玄武岩前积层(淬碎熔岩角砾)及凝灰岩,沉积在茅口组灰岩之上;在后期的演化过程中,炽热玄武岩冷凝收缩,柱状节理发育,在潮汐流和沿岸流的相互作用下不断磨蚀形成大厂层砾岩。     

论本溪群

本溪群的原始定义属岩石地层单位,70年来的研究表明;其界线层型及岩百地层单位性质没有改变过;本溪群可以延展到东北南部及华北,并具有二分、三分的概分性。华北本溪群曾被划分为湖田组(铁铝岩组)和畔沟组,业已证实二分的可行性;张守信(1988,1992)欲以太原组(新涵义)替代本溪群(剔除下部铁铝岩组合)是不可行的。     

MM5和ETA相似理论近地层方案对农田下垫面通量模拟比较研究

近地层湍流通量计算对于中尺度数值模式有重要意义, 湍流通量的参数化是当前大气边界层研究的重要课题之一。选择青藏高原东缘大理观象台边界层通量观测系统, 离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA 相似理论迭代方案B), 并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析。在大理观象台观测场不同植被随季节交替的状况下, 根据边界层铁塔4层高度风速拟合, 发现近地层空气动力学粗糙度随季节变化特征明显。将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算。结果表明:稳定度是影响近地层参数化方案精度的重要因素, 在不稳定条件下方案B低估了动量通量, 方案A优于方案B, 而在稳定条件下方案A低估了动量通量, 方案B优于方案A, 两种方案总体来看误差不大。对于大理边界层通量观测场地农田植被交替的环境条件, 不同季节下垫面植被类型的差异, 以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显著影响。方案B考虑了空气动力学粗糙度z₀和热量粗糙度z_0h的差异, 不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A。针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象, 在使用了(Zeng, et al, 1998)提出的对于使用辐射地表温度在裸土下垫面时的订正方法后, 计算结果也有明显改善。     

Measurements of Flows Over Isolated Valleys

Wind-tunnel measurements of the flow over an isolated valley both normal and at an angle (45°) to a simulated neutrally stable atmospheric boundary layer are presented. Attention is concentrated on the nature of the flow within the valley itself. The work formed part of a wider study that included detailed field measurements around an African desert valley and some limited comparisons with that work are included. A scale of about 1:1000 was used for the laboratory work, in which an appropriate combination of hot wire and particle image velocimetry was employed. For a valley normal to the upwind flow, it is shown that the upstream influence of the valley extends to a distance of at least one half of the axial valley width upstream of the leading edge, whereas differences in mean flow and turbulence could be identified well beyond two valley widths from the downwind edge. Non-normal wind angles lead to significant along-valley flows within the valley and, even at two valley heights above the valley ridge level, there remains a significant spanwise flow component. Downwind turbulence levels are somewhat lower in this case, but are still considerably higher than in the undisturbed boundary layer. At both flow angles, there are significant recirculation regions within the valleys, starting from mean separation just beyond the leading edge, but the strong spanwise flow in the 45° case reduces the axial extent of the separated zone. The flow is shown to be in some ways analogous to flow over an isolated hill. Our results usefully enhance the field data and could be used to improve modelling of saltation processes in the field.     

The Neutral, Barotropic Planetary Boundary Layer, Capped by a Low-Level Inversion

The presence of a low-level, capping inversion layer will affect the height and structure of the planetary boundary layer (PBL). Results from models of varying levels of sophistication, including analytical, turbulent kinetic energy (TKE), second-order closure (SOC), large-eddy simulation (LES) and direct numerical simulation (DNS) models, are used to investigate this influence for the neutral, barotropic PBL. Predicted and observed profiles of stress and geostrophic departure components, and integral measures, such as the parameters of Rossby-number similarity theory, are compared for the KONTUR, Marine Stratocumulus, JASIN, Leipzig, Pre-Wangara and Upavon field experiments.Analytical models of the equilibrium value of inversion height z_i, which depend on the surface friction velocity u_*, and both the Coriolis parameter f and the free-flow Brunt-Väisälä frequency N, are found to give reasonable estimates of the PBL height. They also indicate that only the KONTUR and Marine Stratocumulus experiments were strongly influenced by N. More quantitative comparisons would require larger, more comprehensive datasets. The effects of the presence of a capping inversion on the profile structure were found to be insignificant for h_* = |f|z_i/u_* > 0.15.The simple analytical model performed quite well over all values of h_*; it predicted the profiles of the longitudinal stress component (in the direction of the surface stress) better than the lateral component. The more advanced models performed well for small values of h_* (for flow over the sea), but systematically underestimated the cross-isobaric angle for flow over land. These models predicted the profiles of the lateral stress component better than the longitudinal component. The profiles of the analytical model agreed with those of the advanced models when the constant eddy viscosity of the outer layer was increased.Agreement with DNS was achieved by increasing the eddyviscosity of the analytical model by a factor of 5.Zilitinkevich and Esau(2002, Boundary-Layer Meteorology, 104, 371–379)suggest that the neutral, barotropic values of A and B of Rossby-numbersimilarity theory are not universal constants, but depend on the ratio N/|f|. The dependence for A and B is calculated using the analytical model and TKE models. Over the sea (h_* 0.1; N/|f| 100, where we have used the Zilitinkevich-Esau relation to convert between h_* and N/|f|) there is agreement between the model predictions and observations; however over land where the equilibrium boundary-layer height is greater (h_* 0.35; N/|f| 10) the inconsistency between the advanced model predictions (TKE, SOC, LES, and DNS) and observations, as noted previously by Hess and Garratt, still exists. We attribute this disagreement to violations of the strict assumptions of steady, horizontally homogeneous, neutral, barotropic conditions implicit in the observations. At small values of z_i and a strongly stable background stratification (h_* 0.04; N/|f| 1000) both the TKE and analytical models predict that A and B depend significantly on h_*, however observations are unavailable to confirm these predictions. Zilitinkevich and Esau call this case the `long-lived near-neutral PBL', and state that it is found in cold weather at high latitudes.     

One-Equation Turbulence Modelling for Atmospheric and Engineering Applications

A new algebraic turbulent length scale model is developed, based on previous one-equation turbulence modelling experience in atmospheric flow and dispersion calculations. The model is applied to the neutral Ekman layer, as well as to fully-developed pipe and channel flows. For the pipe and channel flows examined the present model results can be considered as nearly equivalent to the results obtained using the standard k– model. For the neutral Ekman layer, the model predicts satisfactorily the near-neutral Cabauw friction velocities and a dependence of the drag coefficient versus Rossby number very close to that derived from published (G. N. Coleman) direct numerical simulations. The model underestimates the Cabauw cross-isobaric angles, but to a less degree than the cross-isobar angle versus Rossby dependence derived from the Coleman simulation. Finally, for the Cabauw data, with a geostrophic wind magnitude of 10 ms^–1, the model predicts an eddy diffusivity distribution in good agreement with semi-empirical distributions used in current operational practice.     

一次梅雨锋结构及锋区强度对锋面结构影响的数值模拟

利用中尺度模式MM5对1999年6月下旬发生在长江中下游地区的一次锋上西南低涡发展的梅雨锋暴雨过程进行了研究。通过敏感性试验, 探讨了低涡对锋强、 锋生的敏感性, 锋生、 锋消与低涡发展演变的时间空间上的联系等问题\.结果表明, 低涡对中层锋面强度的变化敏感, 中层锋强发生改变会导致整层涡度发生同位相的变化。中层系统对梅雨锋气旋的发展起主要作用, 中层锋面对低涡发展的影响比边界层锋更关键。中层锋面强度改变会激发出中低层的垂直次级环流, 低层的辐合上升运动发生改变, 最终影响到低涡强度的发展。敏感性试验证明, 锋面强度变化对低涡发展造成的影响有超前性, 中层锋面先发生变化, 6 h左右以后低涡强度发生改变。锋区强度和低涡强度发生变化的地理位置基本在同一经度上, 锋强变化的位置在涡强变化位置的北面。     

副热带北太平洋SST在ENSO冷暖位相转变中的作用

2015/2016年发生的极端El Ni?o事件,与1997/1998年El Ni?o事件具有可比拟的强度,但是2016年事件转变为弱La Ni?a,而1998年事件则为强La Ni?a。本文通过对比这两次极端El Ni?o事件,揭示其转变为不同强度La Ni?a事件的物理机制。混合层热收支分析的结果表明,在El Ni?o衰减年的4~11月,2016年平流反馈和温跃层反馈相对较弱,混合层温度衰减速率慢,其产生的主要原因是赤道中西太平洋的东风异常较弱。进一步分析表明,赤道中西太平洋的东风异常与副热带东北太平洋的海表温度异常（SSTA）有关,该地区的SST在1998年表现为冷异常,2016年为暖异常。副热带东北太平洋冷的SSTA有利于信风加强,从而加强中西太平洋的东风异常;而暖的SSTA使得赤道以北出现西南风异常,从而削弱中西太平洋的东风异常。此外,合成分析也表明,副热带东北太平洋SSTA与转变的La Ni?a的强度具有关联,El Ni?o转变为强La Ni?a的情况在位相转变期伴随着副热带北太平洋冷的SSTA,而El Ni?o转变为弱La Ni?a的情况没有明显的冷SSTA。     

边界层温度廓线遥感反演的本地化改进研究

采用MODIS资料和美国发展的MODIS大气温、湿度廓线统计反演算法,估算大气温度、湿度廓线作为初始场,应用101层快速透过率模式(PFAAST)估算了大气透过率,并采用Newton非线性迭代算法反演中国西北荒漠戈壁地区大气温度廓线。结果表明:该方法对边界层高度及以上部分的大气温度反演得比较好,误差基本都在2 K范围内,边界层范围内的温度反演误差较大,反演误差与气溶胶光学厚度增量和地表温度估算误差呈显著正相关关系,与大气水汽混合比的关系较差。文中从敏感性试验和理论分析角度阐述了地表温度和气溶胶光学厚度估算误差对大气温度反演误差的影响,发现不同光谱波段的地表温度权重均随地表温度的增加有不同程度增加,地表温度反演误差增加将增加地表温度权重,提高地表温度估算误差有助于提高地表温度权重的精度;荒漠戈壁地区大气边界层中气溶胶浓度较高,光学厚度较大,使边界层大气透过率降低,进而降低卫星红外遥感波段的地表温度权重和空气温度权重。由于该模式没有很好地考虑边界层中沙尘气溶胶的影响,使卫星反演的大气透过率偏高,以至于高估地表温度权重和大气温度权重,使得反演的表面温度和空气温度偏低。该研究结合太阳光度计获得的光学厚度资料,采用统计方法对气溶胶效应引起的大气透过率误差和表面温度估算误差进行校正,并对物理算法进行本地化改进,实现了边界层温度廓线的反演。