Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system

Simulations of flow for a discrete-fracture model in fractured porous rocks have gradually become more practical, as a consequence of increased computer power and improved simulation and characterization techniques. Discrete-fracture models can be formulated in a lower-dimensional framework, where the fractures are modeled in a lower dimension than the matrix, or in an equi-dimensional form, where the fractures and the matrix have the same dimension.     

Design and commissioning of a laminar soil container for use on small shaking tables

This paper describes the design, fabrication and commissioning of a single axis laminar shear box for use in seismic soil–structure interaction studies. A laminar shear box is a flexible container that can be placed on a shaking table to simulate vertical shear-wave propagation during earthquakes through a soil layer of finite thickness. The laminar shear box described in this paper was designed to overcome the base shear limitations of a small shaking table at The University of Western Ontario. The design details of the box are provided in addition to results of dynamic tests performed to commission the box. A synthetic clay comprising sodium bentonite mixed with diluted glycerin was used as the model soil and 1-G similitude theory was employed to maintain model to prototype similarity. The model soil was compacted into the container in lifts to achieve soil stiffness that increased with depth. A series of shaking table tests and numerical analyses that were performed to study the performance of the laminar box and non-linear seismic behavior of the model clay are described. The results of this study show that the laminar box does not impose significant boundary effects and is able to maintain 1-D soil column behavior. In addition, the dynamic behavior of the model clay during scaled model tests was found to be consistent with the behavior measured during cyclic laboratory tests.     

Wave-induced steady streaming,mass transport and net sediment transport in rough turbulent ocean bottom boundary layers

The interaction between two important mechanisms which causes streaming has been investigated by numerical simulations of the seabed boundary layer beneath both sinusoidal waves and Stokes second order waves, as well as horizontally uniform bottom boundary layers with asymmetric forcing. These two mechanisms are streaming caused by turbulence asymmetry in successive wave half-cycles (beneath asymmetric forcing), and streaming caused by the presence of a vertical wave velocity within the seabed boundary layer as earlier explained by Longuet-Higgins. The effect of wave asymmetry, wave length to water depth ratio, and bottom roughness have been investigated for realistic physical situations. The streaming induced sediment dynamics near the ocean bottom has been investigated; both the resulting suspended load and bedload are presented. Finally, the mass transport (wave-averaged Lagrangian velocity) has been studied for a range of wave conditions. The streaming velocities beneath sinusoidal waves (Longuet-Higgins streaming) is always in the direction of wave propagation, while the streaming velocities in horizontally uniform boundary layers with asymmetric forcing are always negative. Thus the effect of asymmetry in second order Stokes waves is either to reduce the streaming velocity in the direction of wave propagation, or, for long waves relative to the water depth, to induce a streaming velocity against the direction of wave propagation. It appears that the Longuet-Higgins streaming decreases as the wave length increases for a given water depth, and the effect of wave asymmetry can dominate, leading to a steady streaming against the wave propagation. Furthermore, the asymmetry of second order Stokes waves reduces the mass transport (wave-averaged Lagrangian velocity) as compared with sinusoidal waves. The boundary layer streaming leads to a wave-averaged transport of suspended sediments and bedload in the direction of wave propagation.     

First-principles derivation of reactive transport modeling parameters for particle tracking and PDE approaches

Both Eulerian and Lagrangian reactive transport simulations in natural media require selection of a parameter that controls the “promiscuity” of the reacting particles. In Eulerian models, measurement of this parameter may be difficult because its value will generally differ between natural (diffusion-limited) systems and batch experiments, even though both are modeled by reaction terms of the same form. And in Lagrangian models, there previously has been no a priori way to compute this parameter. In both cases, then, selection is typically done by calibration, or ad hoc. This paper addresses the parameter selection problem for Fickian transport by deriving, from first principles and D (the diffusion constant) the reaction-rate-controlling parameters for particle tracking (PT) codes and for the diffusion–reaction equation (DRE). Using continuous time random walk analysis, exact reaction probabilities are derived for pairs of potentially reactive particles based on D and their probability of reaction provided that they collocate. Simultaneously, a second PT scheme directly employing collocation probabilities is derived. One-to-one correspondence between each of D, the reaction radius specified for a PT scheme, and the DRE decay constant are then developed. These results serve to ground reactive transport simulations in their underlying thermodynamics, and are confirmed by simulations.     

地表沉陷空间预报及其有理分式解法的应用

通过阐述地面沉陷机理,分析地下岩体内应力状态和地下开挖空间的形成过程.引入数学积分模型建立地面沉陷理论模型,重点说明单元下沉盆地和特征点下沉,推导出地面沉陷预报公式.并利用有理分式解法求解地面沉陷预报数值.此法的计算精度高,适用面广,有效地解决了常规解法工作量大、计算精度低的缺点.为一种有效、实用的计算地面沉陷预报值的方法.     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.     

比例边界有限元子结构法研究

比例边界有限元法最初应用于土-结构的相互作用分析，经过近几年的完善和发展，如今已经能够应用到其他很多领域。但是因为比例边界有限元理论是基于相似性要求的，使得其在处理几何形状复杂的结构时，会有很大的局限性，从而在某些领域的应用仍旧受到限制。同时由于其全时空耦合，导致大量计算量和工作量，也是其应用受限的一个原因。采用子结构法，打破这些局限性，并且分别针对有限域、无限域的问题，对比例边界有限元子结构法进行了研究，得出了有利于比例边界有限元法在工程实践中应用的结论，为其在实际工程应用中提供了可靠的依据和规律。     

Theorem of turbulent intensity and macroscopic mechanism of the turbulence development

Turbulence is one of the most common nature phenomena in everyday experience, but that is not adequately understood yet. This article reviews the history and present state of development of the turbulence theory and indicates the necessity to probe into the turbulent features and mechanism with the different methods at different levels. Therefore this article proves a theorem of turbulent transpor- tation and a theorem of turbulent intensity by using the theory of the nonequilibrium thermodynamics, and that the Reynolds turbulence and the Rayleigh-Bénard turbulence are united in the theorems of the turbulent intensity and the turbulent transportation. The macroscopic cause of the development of fluid turbulence is a result from shearing effect of the velocity together with the temperature, which is also the macroscopic cause of the stretch and fold of trajectory in the phase space of turbulent field. And it is proved by the observed data of atmosphere that the phenomenological coefficient of turbulent in- tensity is not only a function of the velocity shear but also a function of temperature shear, viz the sta- bility of temperature stratification, in the atmosphere. Accordingly, authenticity of the theorem, which is proved by the theory of nonequilibrium thermodynamics, of turbulent intensity is testified by the facts of observational experiment.     

北京延庆钠层上边界的平均特性研究

本文对北京延庆(116.0°E,40.5°N)钠荧光激光雷达2018年11月—2019年12月夜间的数据以1 h和1 km的时空分辨率进行了分析.获得了当地钠层上边界可达到的高度范围、其相应的季节变化、夜间演化等：延庆顶部钠层可达到110 km、120 km和130 km的概率分别为99.94%、84.46%和40.34%,由此得出,顶部钠层在大多数情况下可以达到120 km,突破了顶部钠层在105 km(或者110 km)的传统观点;上边界高度范围的季节变化规律为5—6月份最高,2—3月份最低;后半夜顶部钠层密度会有所增加.我们比较了偶发钠层、温度和流星注入与上边界高度范围的相关性,推测出偶发钠层对钠层上边界拓展有很大贡献,流星注入次之,温度影响最低.后半夜流星注入量增多跟此时顶部钠层密度增加有一定的相关性.对顶部钠层能够提高到120 km的这项研究,可以为拓展风温的探测高度提供理论依据.

     

On determining the noon polar cap boundary from SuperDARN HF radar backscatter characteristics

Previous work has shown that ionospheric HF radar backscatter in the noon sector can be used to locate the footprint of the magnetospheric cusp particle precipitation. This has enabled the radar data to be used as a proxy for the location of the polar cap boundary, and hence measure the flow of plasma across it to derive the reconnection electric field in the ionosphere. This work used only single radar data sets with a field of view limited to 2 h of local time. In this case study using four of the SuperDARN radars, we examine the boundary determined over 6 h of magnetic local time around the noon sector and its relationship to the convection pattern. The variation with longitude of the latitude of the radar scatter with cusp characteristics shows a bay-like feature. It is shown that this feature is shaped by the variation with longitude of the poleward flow component of the ionospheric plasma and may be understood in terms of cusp ion time-of-flight effects. Using this interpretation, we derive the time-of-flight of the cusp ions and find that it is consistent with approximately 1 keV ions injected from a subsolar reconnection site. A method for deriving a more accurate estimate of the location of the open-closed field line boundary from HF radar data is described.