地球磁层顶的漂移动力学Alfvn波不稳定性及其反常输运效应

本文用回旋动力学研究了磁层顶等离子体低频漂移动力学不稳定性.在β≥1附近发现两支不稳定的漂移动力学Alfvèn模(DKA).它们可在▽T(∥-▽n)(?)0时激发,速度剪切提供主要的自由能源.两支DKA均具有非零的(δE_y,δE_∥)和(δB_x,δB_∥).在湍动的非线性饱和状态下,δB_x的起伏可导致很强的反常输运,当|δB_x|=1nT时,D_⊥可达到10~9m~2/s的数量级.因此,DKA可能在太阳风-磁层耦合过程中起重要作用.     

地球磁层顶的漂移动力学Alfvn波不稳定性及其反常输运效应

本文用回旋动力学研究了磁层顶等离子体低频漂移动力学不稳定性.在β≥1附近发现两支不稳定的漂移动力学Alfvèn模（DKA）.它们可在▽T（∥-▽n）≠0时激发,速度剪切提供主要的自由能源.两支DKA均具有非零的（δEy,δE∥）和（δBx,δB∥）.在湍动的非线性饱和状态下,δBx的起伏可导致很强的反常输运,当|δBx|=1nT时,D⊥可达到109m2/s的数量级.因此,DKA可能在太阳风-磁层耦合过程中起重要作用.     

地球磁层顶的漂移动力学Alfvn波不稳定性及其反常输运效应

本文用回旋动力学研究了磁层顶等离子体低频漂移动力学不稳定性.在β≥1附近发现两支不稳定的漂移动力学Alfvèn模（DKA）.它们可在▽T（∥-▽n）≠0时激发,速度剪切提供主要的自由能源.两支DKA均具有非零的（δE_y,δE_∥）和（δB_x,δB_∥）.在湍动的非线性饱和状态下,δB_x的起伏可导致很强的反常输运,当|δB_x|=1nT时,D_⊥可达到10⁹m²/s的数量级.因此,DKA可能在太阳风-磁层耦合过程中起重要作用.     

磁层亚暴扩展相的数值模拟

本文用多步隐格式方法,求解包含电阻的磁流体力学方程组,对磁层亚暴扩展相作数值模拟.计算结果展现了亚暴扩展相期间磁尾变化的主要特征.它表明,球向和尾向等离子体流与准稳态重联有关.等离子体团的尾向喷发,致使中性线尾侧电流片内的密度降低,等离子体片变薄.中性线近地侧的等离子体团朝着地球运动,并合并于地球附近的重联区内.     

磁层亚暴扩展相的数值模拟

本文用多步隐格式方法,求解包含电阻的磁流体力学方程组,对磁层亚暴扩展相作数值模拟.计算结果展现了亚暴扩展相期间磁尾变化的主要特征.它表明,球向和尾向等离子体流与准稳态重联有关.等离子体团的尾向喷发,致使中性线尾侧电流片内的密度降低,等离子体片变薄.中性线近地侧的等离子体团朝着地球运动,并合并于地球附近的重联区内.     

中纬电离层大尺度扰动的数值模拟

利用计算机模拟研究了不同条件下中纬电离层F区大尺度扰动的时空演变.当存在幅度约为5％的初始扰动时,Perkins不稳定性能发展成为大尺度电离层不均匀体.大气重力波能触发Perkins不稳定性,在更短时间内增长成为大幅度电离层振荡.非均匀局部电场结构也能产生局部电离层的大幅度扰动.本义的模拟表明,中纬电离层大尺度扰动主要表现为磁通管内等离子体沿磁力线的振荡.模拟结果与观测现象符合得很好.     

中纬电离层的低频等离子体不稳定性

本文研究了中纬电离层Rayleigh-Taylor不稳定性和漂移不稳定性的基本性质.如果不存在中性风和背景电场,中纬Rayleigh-Taylor不稳定性的增长率比赤道相应条件下的小.中性风和电场对中纬Rayleigh-Taylor不稳定性有重要影响.当中性风速度达到10m/s时,其作用就超过了重力.漂移不稳定性也受到重力、中性风和电场的影响,离子-中性粒子碰撞能降低漂移不稳定性的增长率.本文的分析表明,在赤道电离层容易产生大尺度的Rayleigh-Taylor不稳定性,在中纬电离层易产生小尺度的漂移不稳定性.     

中纬电离层大尺度扰动的数值模拟

利用计算机模拟研究了不同条件下中纬电离层F区大尺度扰动的时空演变.当存在幅度约为5％的初始扰动时,Perkins不稳定性能发展成为大尺度电离层不均匀体.大气重力波能触发Perkins不稳定性,在更短时间内增长成为大幅度电离层振荡.非均匀局部电场结构也能产生局部电离层的大幅度扰动.本义的模拟表明,中纬电离层大尺度扰动主要表现为磁通管内等离子体沿磁力线的振荡.模拟结果与观测现象符合得很好.     

中纬电离层的低频等离子体不稳定性

本文研究了中纬电离层Rayleigh-Taylor不稳定性和漂移不稳定性的基本性质.如果不存在中性风和背景电场,中纬Rayleigh-Taylor不稳定性的增长率比赤道相应条件下的小.中性风和电场对中纬Rayleigh-Taylor不稳定性有重要影响.当中性风速度达到10m/s时,其作用就超过了重力.漂移不稳定性也受到重力、中性风和电场的影响,离子-中性粒子碰撞能降低漂移不稳定性的增长率.本文的分析表明,在赤道电离层容易产生大尺度的Rayleigh-Taylor不稳定性,在中纬电离层易产生小尺度的漂移不稳定性.     

中、低纬度地区电离层不同层结间的耦合对夜晚F区不规则结构的影响

电离层不规则结构的形成和演化与电离层等离子体不稳定性密切相关 .长期以来 ,在中低纬区扩展 -F现象的研究中 ,没有考虑电离层上下层结之间的相互作用 .本文从理论上全面探讨了在中纬度地区 ,E区可变的Pedersen电导率和Hall电导率 ,与F区可变的Pedersen电导率的共同作用对F区梯度漂移不稳定性的影响 ,导出了存在这种耦合时电离层梯度漂移不稳定性的统一表达式 .这一耦合理论不但解释了实际观测中发现的 ,在某些地方电离层F区顶部的不稳定性发生率要高于F区底部这一现象 ;同时还表明 ,电离层E区与F区的耦合对F区夜晚梯度漂移不稳定性的形成不仅会有阻碍作用 ,同时还使得中、低纬度地区的扰动增长具有了方向的选择性 .该理论的一个重要结论是 ,F区中某一地区能否发展出梯度漂移不稳定性 ,并不完全由当地电离层F区的状态决定 ,同一根磁力线连接的、位于不同纬度地区的E区层结对其发展和演化也会有相当大的影响 .     

木卫二尾迹的数值模拟

木卫二在木星磁层中运动,与快速共转的木星磁场相互作用形成木星极光中木卫二的尾迹,此尾迹反映了连续扰动的磁通量管对木星电离层电流注入的过程.本文应用细丝理论对此过程进行数值模拟,未扰动的上游磁通量管经过大约720 s的时间与木卫二作用形成角度约为4°且相对于木卫二静止的拉伸磁通量管,随后在磁张力的作用下加速运动.在随后大约2160 s的加速时间内,磁通量管的速度逐渐与背景磁力线接近,向木星电离层注入的电流也由9×105A逐渐衰减至0.以1/e初始电流处认为是可观测的木卫二尾迹限度,则木卫二的尾迹仅大约有一至两个经度,这与哈勃望远镜观测到的木卫二尾迹大小大致相同.木卫二尾迹与木卫一尾迹相比,输送的能量约为2.3×1013J,仅有木卫一的1/10,导致尾迹长度远小于木卫一的尾迹长度.此外,模拟过程中发现木卫二下游的压力槽对磁通量管运动的影响很小,可以忽略.     

Numerical simulations of non-ergodic solute transport in three-dimensional heterogeneous porous media

Numerical simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity, , were conducted in a three-dimensional heterogeneous and statistically isotropic aquifer. The hydraulic conductivity, K(x), is modeled as a random field which is assumed to be log-normally distributed with an exponential covariance. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume and the plume centroid variances were simulated with 1600 Monte Carlo (MC) runs for three variances of log K, _Y²=0.09, 0.23, and 0.46, and a square source normal to of three dimensionless lengths. It is showed that 1600 MC runs are needed to obtain stabilized results in mildly heterogeneous aquifers of _Y²0.5 and that large uncertainty may exist in the simulated results if less MC runs are used, especially for the transverse second spatial moments and the plume centroid variance in transverse directions. The simulated longitudinal second spatial moment and the plume centroid variance in longitudinal direction fit well to the first-order theoretical results while the simulated transverse moments are generally larger than the first-order values. The ergodic condition for the second spatial moments is far from reaching in all cases simulated and transport in transverse directions may reach ergodic condition much slower than that in longitudinal direction.     

Theory of azimuthally small-scale hydromagnetic waves in the axisymmetric magnetosphere with finite plasma pressure

The structure of monochromatic MHD-waves with large azimuthal wave number m 1 in a two-dimensional model of the magnetosphere has been investigated. A joint action of the field line curvature, finite plasma pressure, and transversal equilibrium current leads to the phenomenon that waves, standing along the field lines, are travelling across the magnetic shells. The wave propagation region, the transparency region, is bounded by the poloidal magnetic surface on one side and by the resonance surface on the other. In their meaning these surfaces correspond to the usual and singular turning points in the WKB-approximation, respectively. The wave is excited near the poloidal surface and propagates toward the resonance surface where it is totally absorbed due to the ionospheric dissipation. There are two transparency regions in a finite-beta magnetosphere, one of them corresponds to the Alfvén mode and the other to the slow magneto-sound mode.     

Comparison of laboratory-scale solute transport visualization experiments with numerical simulation using cross-bedded sandstone

Using a slab of Massillon Sandstone, laboratory-scale solute tracer experiments were carried out to test numerical simulations using the Advection–Dispersion Equation (ADE). While studies of a similar nature exist, our work differs in that we combine: (1) experimentation in naturally complex geologic media, (2) X-ray absorption imaging to visualize and quantify two-dimensional solute transport, (3) high resolution transport property characterization, with (4) numerical simulation. The simulations use permeability, porosity, and solute concentration measured to sub-centimeter resolution. While bulk breakthrough curve characteristics were adequately matched, large discrepancies exist between the experimental and simulated solute concentration fields. Investigation of potential experimental errors suggests that the failure to fit solute concentration fields may lie in loss of intricate connectivity within the cross-bedded sandstone occurring at scales finer than our property characterization measurements (i.e., sub-centimeter).     

激发态氮分子N*2在电离层中的作用的模拟研究

激发态氮分子N*2在电离层F区中起着重要的作用,它使F区占主导地位的O+离子的损失率增大,从而使该区的电子浓度减少. 本文利用理论电离层数值模型,通过考虑与不考虑N*2的作用,对包括1990年6月、1997年5月、1998年5月以及2000年4月磁暴事件在内的时间区间的电离层响应情形进行模拟研究,并与实测结果进行对比. 结果表明,N*2对电离层电子浓度的影响在太阳活动高年非常明显,在太阳活动低年虽有些影响,但效果并不明显,其程度远不如高年. 在太阳活动高年,不仅是磁暴期间,在较宁静期间也必须考虑N*2的影响. 而且,在考虑N*2的作用时,还与激发态振动温度Tν有关,在采用Tν=Tn（其中Tn为背景中性大气的温度）的简化处理时,所得结果与观测结果的符合程度不如对Tν进行精确计算时所得的结果好. 模拟结果还表明,太阳活动高年,N*2作用的结果主要是使150km高度以上的F区电离层电子浓度减少,而对150km以下高度的电离层电子浓度则影响不大. 另外,N*2基本不影响F2层峰高hmF2的值.     

均匀覆盖层中氡迁移的数值模拟

氡能从地下深部运移到地表,可带来反映地下地质状态的信息.本文以氡在均匀覆盖层中的迁移过程为对流-扩散方式作为理论基础,推导了二维氡浓度微分方程,用有限差分方法计算了无限氡源层、有限单氡源层、有限双氡源层在均匀覆盖层中氡浓度的数值解,以二维断面的形式展示了氡的分布,其结果与野外实测相符.     

Numerical simulation of dam-break flow and bed change considering the vegetation effects

Existing numerical investigations of dam-break flows rarely consider the effects of vegetation.This paper presents a depth-averaged two-dimensional model for dam-break flows over mobile and vegetated beds.In the model,both the consequences of reducing space for storing mass and momentum by the existence of vegetation and dragging the flow are considered:the former is considered by introducing a factor (1-c) to the flow depth,where c is the vegetation density;the later is considered by including an additional sink term in the momentum equations.The new governing equations are discretized by the finite volume method;and an existing second-order central-upwind scheme embedded with the hydrostatic reconstruction method for water depth,is used to estimate the fluxes;the source terms are estimated by either explicit or semi-explicit methods fulfilling the stability requirement.Laboratory experiments of dam-break flows or quasi-steady flows with/without vegetation effects/sediment transport are simulated.The good agreements between the measurements and the numerical simulations demonstrate a satisfactory performance of the model in reproducing the flow depth,velocity and bed deformation depth.Numerical case studies of six scenarios of dam-break flows over a mobile and vegetated bed are conducted.It is shown that when the area of the vegetation zone,the vegetation density,and the pattern of the vegetation distribution are varied,the resulted bed morphological change differs greatly,suggesting a great influence of vegetation on the dam-break flow evolution.Specifically,the vegetation may divert the direction of the main flow,hindering the flow and thus result in increased deposition upstream of the vegetation.     

岩石样品中超声波传播数值模拟及分析

为了深入研究超声波速度测试中超声波在岩石样品中的传播规律,对岩石超声波速度实验采用交错网格有限差分算法进行了数值模拟.选取了适合于超声波震源的差分精度、稳定性条件和吸收边界.通过模拟波形与实验记录波形的对比分析,结果表明:激发器和接收器直接位于岩心与探头接触面时,接收到的波形中干扰波较少,合成记录波形特征明显,易于分析;使用雷克子波作为超声波数值模拟的震源子波,计算得到实际首波初至时间为21.55μs,而波形图中首波初至时间为23.90μs,这是由于子波自身的长度使得波初至时间比实际波初至的时间晚;相对于岩心与探头接触面的反射波,岩心柱体与空气接触面的反射波对有效波形信号的影响更大.