油水饱和球管孔隙模型弛豫的理论计算与计算机模拟

根据Brownstein-Tarr理论,采用特征函数展开法,通过求解基于扩散效应的Bloch控制方程,给出了油水饱和球管孔隙模型弛豫的理论计算公式.理论计算结果表明,在水润湿条件下,油的弛豫过程只与含油饱和度有关,而与岩石的孔隙结构无关.根据理论计算结果对球管模型中水的弛豫进行了数值模拟,模拟结果显示,球和管的主弛豫过程是一个单指数函数,其余部分与之相比可以忽略,即球管模型中水的弛豫可近似为一个双指数衰减过程. 关键词： 核磁共振 弛豫 球管孔隙模型 Bloch控制方程     

射频场照射下多自旋体系弛豫的理论计算

根据Liouville-von Neumann方程，对射频场照射下多自旋体系的弛豫进行了理论描述，并用WBR理论推导出了体系的弛豫方程组，给出了各类弛豫速率的理论计算公式.在此基础上，编制了弛豫方程组数值解的计算程序，分别用此程序和Bloch方程计算了双自旋体系在不同情况下的稳态解，并对计算结果进行了简要的分析和讨论. 关键词： 核磁共振 弛豫 射频场 多自旋体系     

玻璃化转变的分子串模型中分子串弛豫模式的计算机模拟

基于玻璃化转变的分子串模型的哈密顿量(Hamiltonian),提出了模拟分子串第一和第二弛豫模式的蒙特卡罗(Monte Carlo)模拟方案.模拟得出的第一和第二弛豫模式的弛豫时间,与分子串模型中分子串弛豫方程所预言的结果完全一致,即理论预期和模拟结果相互印证.这不仅说明了分子串模型的分子串弛豫方程、至少是第一和第二弛豫模式的理论预言的正确性,同时也表明了本文所提出的模拟方法的有效性,并进一步明晰了分子串中分子的随机涨落和跃迁运动的图像,也为三态甚至是多态的分子串弛豫动力学研究,以及对进一步模拟分子串之 关键词： 弛豫动力学 蒙特卡罗模拟 分子串     

PEMFC气体扩散层干燥过程孔隙网络模拟

本文采用三维孔隙网络模型从孔隙尺度上对质子交换膜燃料电池气体扩散层中水分蒸发的干燥过程进行了模拟,并考虑了阴极流道的影响.计算结果表明毛细力在气体扩散层的干燥过程中起主要作用,气相从阴极流道底部开始,以沿气体扩散层厚度方向侵入为先,直到到达气体扩散层底部,随后气相才向流道肩部水平延伸侵入.蒸发速率随气体扩散层中水饱和度...     

基于ADI-FDTD模拟的岩石核磁共振横向弛豫特征分析

将复杂的骨架-孔隙系统抽象成等效双孔介质,根据Bloch方程构建数学模型,用交替隐式时域有限差分(ADI-FDTD)和联合反演迭代法(SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振T₂谱的反演,定量研究扩散系数、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响.结果表明:横向宏观磁化矢量衰减速率与扩散系数和微孔隙分量成正比,与孔隙宽度成反比,与表面弛豫速率基本无关.当扩散系数较大、孔隙宽度较小时,核磁共振T₂谱难以直观反映孔隙组分及孔隙结构.应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响.     

岩心孔隙介质中流体的核磁共振弛豫

弛豫时间是核磁共振研究中的一个重要参数,岩心孔隙介质流体的弛豫过程是自由流体弛豫机制、表面弛豫机制和流体的扩散弛豫机制共同作用的结果,它包含了丰富的孔隙和流体本身的信息. 弛豫时间和自扩散系数的测量及对弛豫时间的分析是核磁共振技术应用于岩心分析和石油勘测的重要内容.     

具有表面弛豫的液体自扩散的M0nte Carlo模拟

研究多孔介质中液体的自扩散行为能获得介质的微观结构信息,有助于了解介质中液体的传输性质.以具有不同孔隙大小的无限长圆柱体模型中的液体为对象,采用Monte Cado随机行走方法,模拟存在表面弛豫时液体的自扩散系数和核自旋磁化强度随时间的变化,导出将NMR弛豫参数和随机行走参数联系在一起的表达式.结果表明:在快扩散区,液体的核自旋磁化强度随弛豫时间呈单指数衰减,且自扩散系数在短时情况下独立于表面弛豫率;在慢扩散区,液体的核自旋磁化强度衰减和自扩散系数在短时情况下均与表面弛豫率无关.模拟结果与理论分析相吻合,可用于求解介质的表面积与体积之比及孔径大小等结构信息.     

球对称铷原子玻色-爱因斯坦凝聚中单极子模的朗道阻尼和频移

采用含时哈特里-福克-博戈留波夫近似研究球对称铷原子玻色-爱因斯坦凝聚中单极子模的朗道阻尼和频移,并用现有实验和数值模拟研究的粒子数和囚禁频率参量,解析计算了阻尼系数和频移大小及其它们的温度依赖.计算中,考虑元激发的实际弛豫及其各弛豫间的正交关系以获得阻尼和频移计算公式,把基态波函数取为高斯分布函数的一级近似以消除三模耦合矩阵元的发散.我们的计算结果与数值模拟结果和实验结果分别进行直接和间接地对比,讨论和说明了我们理论方法的合理性.     

具有表面弛豫的液体自扩散的Monte Carlo模拟

研究多孔介质中液体的自扩散行为能获得介质的微观结构信息,有助于了解介质中液体的传输性质.以具有不同孔隙大小的无限长圆柱体模型中的液体为对象,采用Monte Carlo随机行走方法,模拟存在表面弛豫时液体的自扩散系数和核自旋磁化强度随时间的变化,导出将NMR弛豫参数和随机行走参数联系在一起的表达式.结果表明:在快扩散区,液体的核自旋磁化强度随弛豫时间呈单指数衰减,且自扩散系数在短时情况下独立于表面弛豫率;在慢扩散区,液体的核自旋磁化强度衰减和自扩散系数在短时情况下均与表面弛豫率无关.模拟结果与理论分析相吻合,可用于求解介质的表面积与体积之比及孔径大小等结构信息.     

可激发气体分子声弛豫过程振动模式能量转移速率计算

研究可激发气体中振动模式能量转移速率和声弛豫过程形成的关系,将单一气体Tanczos弛豫方程理论[J.Chem.Phys.25,439(1956)]扩展应用于混合气体中振动模式的振动-振动(V-V)和振动-平动(V-T)能量转移速率的计算。在室温下CO₂,CH₄,CL₂,N₂和O₂组成的多种混合气体中,振动模式能量转移速率的计算结果表明:对于多个振动模式所形成的声复合弛豫过程,各振动模式的声激发能可由V-V能量转移相互耦合后传递给具有最快V-T转移速率的最低振动频率振动模式,再通过该振动模式的V-T转移退激发形成主弛豫过程。这种选择最快转移路径的声激发量弛豫方式,造成了大多数可激发气体中声弛豫吸收谱的实测数据只存在一个吸收峰的现象。从而提供了一个可通过计算微观振动能量转移速率分析混合气体声弛豫过程形成机理的理论模型。      

Absorption of intense regular and noise waves in relaxing media

An integro-differential equation is written down that contains terms responsible for nonlinear absorption, visco-heat-conducting dissipation, and relaxation processes in a medium. A general integral expression is obtained for calculating energy losses of the wave with arbitrary characteristics—intensity, profile (frequency spectrum), and kernel describing the internal dynamics of the medium. It is shown that for weak waves, the general integral leads to well-known results of a linear approximation. Profiles of stationary solutions are constructed both for an exponential relaxation kernel and for other types of kernels. Energy losses at the front of week shock waves are calculated. General integral formulas are obtained for energy losses of intense noise, which are determined by the form of the kernel, the structure of the noise correlation function, and the mean square of the derivative of realization of a random process.     

Hierarchical description for a non-Markoff relaxation process

A phenomenal model is proposed in this paper to describe a non-Markoff relaxation process. The main feature of the model is the fractal time concept introduced to take account of the coupling of a relaxation mode with a heat bath. The theoretical investigation indicates that a weakly non-Markoff relaxation process can be related to a set of hierarchically distributed Markoff relaxation processes and the fractal exponential law exp [?(t/τ)^1?σ] is the character of a non-Markoff relaxation process. Calculations about the temperature and time scale dependence of the exponent σ have been made, and comparison with recent experiments is also given.     

CPMG relaxation by diffusion with constant magnetic field gradient in a restricted geometry: numerical simulation and application

Carr-Purcell-Meiboom-Gill (CPMG) measurements are the primary nuclear magnetic resonance (NMR) technique used for evaluating formation properties and reservoir fluid properties in the well logging industry and laboratory sample analysis. The estimation of bulk volume irreducible (BVI), permeability, and fluid type relies on the accurate interpretation of the spin-spin relaxation time (T(2)) distribution. The interpretation is complicated when spin's self-diffusion in an inhomogeneous field and restricted geometry becomes dominant. The combined effects of field gradient, diffusion, and a restricted geometry are not easily evaluated analytically. We used a numerical method to evaluate the dependence of the free and restricted diffusion on the system parameters in the absence of surface relaxation, which usually can be neglected for the non-wetting fluids (e.g., oil or gas). The parameter space that defines the relaxation process is reduced to two dimensionless groups: D* and tau*. Three relaxation regimes: free diffusion, localization, and motionally averaging regimes are identified in the (log(10)D*, log(10)tau*) domain. The hypothesis that the normalized magnetization, M*, relaxes as a single exponential with a constant dimensionless relaxation time T*(2) is justified for most regions of the parameter space. The numerical simulation results are compared with the analytical solutions from the contour plots of T*(2). The locations of the boundaries between different relaxation regimes, derived from equalizing length scales, are challenged by observed discrepancies between numerical and analytical solutions. After adjustment of boundaries by equalizing T*(2), numerical simulation result and analytical solution match each other for every relaxation regime. The parameters, fluid diffusivity and pore length, can be estimated from analytical solutions in the free diffusion and motionally averaging regimes, respectively. Estimation of the parameters near the boundaries of the regimes may require numerical simulation.     

Effects of off-resonance irradiation, cross-relaxation, and chemical exchange on steady-state magnetization and effective spin-lattice relaxation times

In the presence of an off-resonance radiofrequency field, recovery of longitudinal magnetization to a steady state is not purely monoexponential. Under reasonable conditions with zero initial magnetization, recovery is nearly exponential and an effective relaxation rate constant R_1eff = 1/T_1eff can be obtained. Exact and approximate formulas for R_1eff and steady-state magnetization are derived from the Bloch equations for spins undergoing cross-relaxation and chemical exchange between two sites in the presence of an off-resonance radiofrequency field. The relaxation formulas require that the magnetization of one spin is constant, but not necessarily zero, while the other spin relaxes. Extension to three sites with one radiofrequency field is explained. The special cases of off-resonance effects alone and with cross-relaxation or chemical exchange, cross-relaxation alone, and chemical exchange alone are compared. The inaccuracy in saturation transfer measurements of exchange rate constants by published formulas is discussed for the creatine kinase reaction.     

An Ising model for the stress relaxation of solids

A kinetic theory of stress relaxation of solids as a cooperative process is proposed. The theory is based on a two-state model for the relaxation. It is shown that the conventional mean field approximation leads to an exponential dependence of the rate of stress relaxation on the stress while the multiplicative approximation of Vol'kenstein et al. leads to a power law. It is argued that the exponential law should be valid initially in the relaxation process while the power law is appropriate when the system is nearer equilibrium, which is in qualitative agreement with recent experimental findings.     

Relaxation processes of water confined to AlMCM-41 molecular sieves. Influence of the hydroxyl groups of the pore surface

A series of AlMCM-41 molecular sieves was prepared with constant composition (Si/Al = 14.7) and presumably same pore structure but different pore diameters (from 2.3 to 4.6 nm). The pore size distribution is narrow for each sample. The rotational fluctuations of water molecules confined inside the pores were investigated applying broadband dielectric spectroscopy (10⁻²–10⁷ Hz) over a large temperature interval (213–333K). A relaxation process, slower than that expected for bulk water, was observed which is assigned to water molecules forming a surface layer on the pore walls. The estimated relaxation time has an unusual non-monotonic temperature dependence, which is rationalized and modeled assuming two competing processes: rotational fluctuations of constrained water molecules and defect formation (Ryabov model). This paper focuses on the defects and notably the influence of the hydroxyl groups of the pore walls. The Ryabov model is fitted to the data and characteristic parameters are obtained. Their dependence on pore diameter is considered for the first time. The found results are compared with those obtained for other types of molecular sieves and related materials.     

Product operator analysis of the influence of chemical exchange on relaxation rates

Measurements of chemical-exchange processes by NMR are widely used to obtain valuable information about molecular dynamics and structure. Here, a computational method is introduced to assess the influence of chemical exchange on spin relaxation rates. The method is based on the inclusion of a random exchange process in product operator calculations on a microscopic level. This product operator approach can be applied to estimate exchange contributions when using sophisticated pulse sequences that cannot be easily described analytically. The method applies to the full range of exchange times measurable by NMR and can incorporate interference effects between exchange and other processes such as scalar coupling. To demonstrate its utility, simulated relaxation data were compared with theoretical predictions of spin-locking and Carr-Purcell spin-echo sequences with hard and adiabatic pulses, using different time scales for a two-site chemical-exchange process. Finally, simulations were used to examine a system in which a second random process is superimposed on a simple two-site exchange process. The method was found to provide a simple and robust tool to analyze pulse sequences and equations commonly used to study exchange-induced relaxation.