A multicomponent model for computing the thermal structure of collapsing protostellar clouds

A model for simulating the thermal and dynamical evolution of protostellar clouds is presented. In the model, the dust and gas temperatures are treated separately, making it possible to more precisely describe the initial stages of the cloud’s gravitational contraction and collapse. The model is fast enough to be applied in hydrodynamical computations, and has a high enough accuracy for the results to be used to compute emission spectra and comparing them with observational data. Two problems are considered as test examples and simple applications: calculation of the structure of clouds in thermal and hydrostatic equilibrium, and modeling the evolution of a protostellar cloud in a spherically symmetric approximation, including the formation of the first hydrostatic core.     

A three-dimensional heat transfer and thermal cracking model considering the effect of cracks on heat transfer

A simple three-dimensional heat transfer model is developed to consider the hindering effect of cracks on heat transfer. The 3D heat transfer model can also be applied to numerical methods such as the combined finite-discrete element method (FDEM), discrete element method (DEM), discontinuous deformation analysis (DDA), the numerical manifold method (NMM), and the finite element method (FEM) to construct thermo-mechanical coupling models that allow these methods to solve thermal cracking problems and dynamically consider the hindering effect of cracks on heat transfer. In the 3D heat transfer model, the continuous-discontinuous medium is discretized into independent tetrahedral elements, and joint elements are inserted between adjacent tetrahedral elements. Heat transfer calculations for continuous-discontinuous media are converted to heat conduction in tetrahedral elements and the heat exchange between the adjacent tetrahedral elements through the joint element. If the joint element between adjacent tetrahedral elements breaks (ie, a crack generates), the heat exchange coefficient of the joint element is reduced to account for the hindering effect of cracks on heat conduction. Then the model and the FDEM are combined to build a thermo-mechanical coupling model to simulate thermal cracking. The thermally induced deformation, stress, and cracking are investigated by the thermo-mechanical coupling model, and the numerical results are compared with analytical solutions or experimental results. The 3D heat transfer model and thermo-mechanical model can provide a powerful tool for simulating heat transfer and thermal cracking in a continuous-discontinuous medium.     

An Improved Method for Calculating Paleoheat Flow from Vitrinite Reflectance Profiles

lNTRODUCTIONIt is wel1 established that the vitrinite reflectance (R<, ) de-pends on the temperature history that a basin has undergone.In general, Ro is a function of maximum paleotemperature andeffective heating time. Thus, the paleotemperature of sedi-ments and the value of vitrinite ref1ectance are controlled by theflow of heat Q(t) into the basin and by the thermal parametersof the sediments. Therefore, in principle, it is possible to cal-culate the thermal vitrinite reflectance of …     

工作荷载下温度循环对桩基变形与应力的影响分析

目前针对工作荷载下温度循环对桩基承载力特性的研究相对较少,基于室内模型试验方法,对饱和砂土中工作荷载下桩体在多次冷热循环作用时的承载特性和传热特性进行研究,测得温度循环作用下桩体和桩周土体温度、桩周水平土压力、桩体应变以及桩顶位移随时间的变化规律。试验结果表明,施加温度循环作用后桩体及桩周土体温度变化不大,桩周水平土压力也能基本恢复到初值,但在桩体内部则会产生较少残余应变,桩顶下沉并随着循环次数的增加不断累积,从而影响结构的整体承载能力。     

含水层热量运移中自然热对流作用的数值模拟

较严密地推导了能够描述含水层热量运移中自然热对流作用的水流方程和热量运移方程,并在此基础上建立了一个三维非线性模型.根据模型的特点,提出综合采用Douglas Brian ADI差分法、特征线法返回技术及改进Picard迭代法求解,并给出了计算步骤.模型及数值方法用于模拟上海第二承压含水层中的一个群井储能试验.试验中的自然热对流作用得到了较好的体现,反映出模型合理可靠,方法有效,结果令人满意.     

Modeling of protostellar clouds and their observational properties

A physical model and two-dimensional numerical method for computing the evolution and spectra of protostellar clouds are described. The physical model is based on a system of magneto-gas-dynamical equations, including ohmic and ambipolar diffusion, and a scheme for calculating the thermal and ionization structure of a cloud. The dust and gas temperatures are determined when calculating the thermal structure of the cloud. The results of computing the dynamical and thermal structure of the cloud are used to model the transfer of continuum and molecular-line radiation in the cloud. Results are presented for clouds in hydrostatic and thermal equilibrium. The evolution of a rotating magnetic protostellar cloud that is compressed starting from a quasi-static equilibrium state is also considered. Spectral maps for optically thick lines of linear molecules are analyzed. The influence of the magnetic field and rotation can lead to a redistribution of angular momentum in the cloud and the formation of a characteristic rotational-velocity structure. As a result, the distribution of the velocity centroid of the molecular lines can acquire an hourglass shape. It is planned in future to use the developed program package and a model for the chemical evolution of clouds to interpret and model in detail observed starless and protostellar cores.     

Three-dimensional random network model for thermal conductivity in particulate materials

This paper describes a three-dimensional random network model to evaluate the thermal conductivity of particulate materials. The model is applied to numerical assemblies of poly-dispersed spheres generated using the discrete element method (DEM). The grain size distribution of Ottawa 20–30 sand is modeled using a logistic function in the DEM assemblies to closely reproduce the gradation of physical specimens. The packing density and inter-particle contact areas controlled by confining stress are explored as variables to underscore the effects of micro- and macro-scales on the effective thermal conductivity in particulate materials. It is assumed that skeletal structure of 3D granular system consists of the web of particle bodies interconnected by thermal resistor at contacts. The inter-particle contact condition (e.g., the degree of particle separation or overlap) and the particle radii determine the thermal conductance between adjacent particles. The Gauss–Seidel method allows evaluation of the evolution of temperature variation in the linear system. Laboratory measurements of thermal conductivity of Ottawa 20–30 sand corroborate the calculated results using the proposed network model. The model is extended to explore the evolution of thermal conduction depending on the nucleation habits of secondary solid phase as an anomalous material in the pore space. The proposed network model highlights that the coordination number, packing density and the inter-particle contact condition are integrated together to dominate the heat transfer characteristics in particulate materials, and allows fundamental understanding of particle-scale mechanism in macro-scale manifestation.     

一种约束盆地低温热历史的裂变径迹技术

裂变径迹(FT)技术是根据矿物中铀裂变产生的辐射损伤特性进行分析的低温热年代技术。随着对FT退火行为和实验退火模型研究的深入,使得这一技术成为约束沉积盆地低温热历史的重要手段。Laslett等、Crowley等和Ketcham等先后提出3个重要的磷灰石FT实验退火模型,其中以Ketcham等的退火模型研究最为深入,它分析了磷灰石类型、时间、温度和化学成分对其径迹退火的影响,使用c轴投射径迹长度和Dpar等参数,形成了描述磷灰石FT多元动力学退火的数学模型。锆石FT与U-Th/He技术、Ro值、地表温度和地层年龄等,均是约束磷灰石FT热历史重建的重要约束条件,HeFTy(2009)是进行低温热历史模拟的主要软件之一。     

面接触条件下胶结颗粒物导热系数模型优化

为了更好地探究岩土类胶结颗粒材料的导热系数变化规律,减少造成的热害问题,根据胶结颗粒物的结构特征,考虑颗粒物和胶结基质为混合物时的状态,并结合有效介质导热系数模型和面接触导热系数模型,构建了适用于岩土类胶结颗粒材料的有效导热系数计算模型。采用瞬态热线法测试了实验试块的导热系数,将预测值与实测结果和Maxwell-Eucken模型的计算结果对比后发现,两模型的导热系数计算值都随着孔隙率的增加而逐渐减小。4组测试数据结果显示,Maxwell-Eucken模型的导热系数计算结果平均值为0.171 72 W/（m·K）,优化的导热系数计算模型结果平均值为0.171 39 W/（m·K）,说明优化的导热系数计算模型得出的结论与实际测试结果平均值（0.169 22 W/（m·K））更接近,能够更好预测胶结颗粒物质的导热系数。      

深层地热能开发及其对地热水流场的影响

为了研究深层地热能开发对热储层流场造成的影响,选取兰考县深层地热能开发为研究对象,分析研究区新生代以来的地质结构和区域构造特征,调查地区地热开采情况及存在的问题。根据地热能类型分析了第四系Q、新近系明化镇组N₂、新近系馆陶组N₁等3个热储层特征,利用可控源音频大地电磁测深查明区内的断裂分布、热储层发育、隔水边界等条件,通过现场试验与测试手段查明不同部位热储层的参数指标,分区计算各层的地热资源量。根据开采方式、边界条件、水文地质参数建立数学模型,通过数值模拟预测出2023年Q热储层流场将恢复到原始流场模式,N₂热储层流场漏斗效应持续减弱,N₁热储层流场漏斗效应明显且分布不规则,为控制地热资源开发对环境造成的影响,提出封闭循环利用地热、同层加压回灌的开采方式及合理布局开采井、综合利用多层地热等措施,为本地及周边地区地热资源开发及环境保护提供参考。      

西宁盆地热储构造概念模型的建立

基于西宁盆地地热地质条件的研究,在Donaldson管状模型的基础上,提出的西宁盆地热储构造概念模型是大地热流为热源一低热导率岩层聚热-深循环逐渐加热受迫对流为机制-构造控水控热。进而揭示出,尽管西宁盆地地热异常分布范围较广,但也并非“遍地有热”。在生产实践中,关键是要较为准确地确定地热井的最佳构造部位和找到高渗透率的热储。     

Modelling of the thermal dependence of structural and elastic properties of calcite, CaCO3

A computational method, based on the quasiharmonic approximation, has been computer-coded to calculate the temperature dependence of elastic constants and structural features of crystals. The model is applied to calcite, CaCO₃; an interatomic potential based on a C-O Morse function and Ca-O and O-O Borntype interactions, including a shell model for O, has been used. Equilibrations in the range 300–800 K reproduce the experimental unit-cell edges and bond lengths within 1%. The simulated thermal expansion coefficients are 22.3 (//c) and 2.6 (⊥ c), against 25.5 and-3.7×10^?6K^?1 experimental values, respectively. The thermal coefficients of elastic constants tend to be underestimated; for the bulk modulus, -2.3 against-3.7×10^?4K^?1 is obtained.     

基于微观结构重塑的非饱和冻土导热系数预测

为预测非饱和冻土的导热性能,基于土体微观结构,提出了非饱和冻土特征结构识别算法和多元素生成算法,并将该算法与传统有限单元法组合,建立非饱和冻土导热系数蒙特卡洛预测模型。通过土体SEM电镜图像,采用逆向四参数增长识别法识别土体中各组分含量、大小以及各方向分布概率;改进传统的四参数随机增长法,提出了考虑土、水、冰和气的多元素生成算法;基于生成的非饱和冻土模型,通过蒙特卡洛方法获得非饱和冻土导热系数,并与规范中冻土导热系数进行对比,验证了蒙特卡洛法预测模型的合理性（平均误差<4%）;通过多因素分析研究孔隙率、颗粒大小、土体导热性、饱和度以及结冰率对非饱和冻土导热性影响,各因素与导热系数的相关系数依次为：-0.352、-0.098、0.641、0.520和0.060,影响大小为：土颗粒导热性>饱和度>孔隙率>土颗粒大小>结冰率。各影响因素对非饱和冻土导热系数影响可以归纳为对热通量形成“热链”密度、宽度、连通性、热流承载力以及对“热桥”通量的影响。     

A 2D coupled hydro-thermal model for the combined finite-discrete element method

Based on the combined finite-discrete element method (FDEM), a two-dimensional coupled hydro-thermal model is proposed. This model can simulate fluid flow and heat transfer in rock masses with arbitrary complex fracture networks. The model consists of three parts: a heat conduction model of the rock matrix, a heat-transfer model of the fluid in the fracture (including the heat conduction and convection of fluid), and a heat exchange model between the fluid and rock at the fracture surface. Three examples with analytical solutions are given to verify the correctness of the coupled model. Finally, the coupled model is applied to hydro-thermal coupling simulations of a rock mass with a fracture network. The temperature field evolution, the effect of thermal conductivity of the rock matrix thermal conductivity and the fracture aperture on the outlet temperature are studied. The coupled model presented in this paper will enable the application of FDEM to study rock rupture driven by the effect of hydro-thermo-mechanical coupling in geomaterials such as in geothermal systems, petroleum engineering, environmental engineering and nuclear waste geological storage.

     

碳质物拉曼光谱变质温度计及其在造山带热结构重建与演化中的应用

变质沉积岩普遍含有碳质物,其源自沉积母岩中有机质。在变质过程中这些有机质逐渐转变为碳质物或石墨,且碳质物结晶程度可以作为变质等级的可靠指示标志。拉曼光谱是表征碳质物结晶度的有效工具,Beyssac et al.(2002a)基于碳质物拉曼光谱参数(R1=D1/G,D1和G为碳质物拉曼光谱峰强;R2=D1/(G+D1+D2),G、D1和D2为碳质物拉曼光谱峰面积)与寄主岩变质温度之间的线性关系构建了碳质物拉曼光谱温度计:T(℃)=-445(R2)+641,其简单且实用,并被应用到阿尔卑斯和喜马拉雅造山带热结构与折返机制研究中;随后,Rahl et al.(2005)对该变质温度计进行修订,修订后温度估算表达式为:T(℃)=737.3+320.9R1-1067R2-80.638R12,并将变质温度估算范围扩展为100~700℃。本文对碳质物拉曼光谱变质温度计的基本原理、方法、应用条件及其在造山带热结构重建与演化方面的研究进展进行了综述,并对碳质物拉曼光谱温度计与传统温度计估算的变质温度进行了系统的对比分析,结果表明碳质物拉曼光谱温度计代表峰期变质温度,不会受后期退变质影响,当传统温度计结果代表峰期变质温度时,二者计算结果一致。碳质物拉曼光谱温度计已被用于造山带热结构重建、折返过程的热演化,以及高应变带、流体相关热异常等研究。尽管变质过程的压力、变质持续时间、碳质物前体类型等因素对于碳质物拉曼光谱温度计的影响尚待研究,但与传统矿物组合温压计相结合,该方法可以有效评价峰期变质条件和造山多期热演化。     

Influence of luminosity bursts on properties of protostellar disks

A (2+ 1)-dimensional numerical model for the formation and evolution of young stellar objects with sub-solar masses is presented. The numerical hydrodynamicall code describing the formation and evolution of a pølar disk in a two-dimensional approximation is supplemented by one-dimensional code for the evolution of the star and an algorithm for establishing the vertical structure of the disk. This code is used to investigate the influence of luminosity bursts with intensities similar to those observed in FU Orionis objects (FUors) on the properties and thermal balance of pølar disks. A model with gravitational instability and fragmentation of the disk, with subsequent migration of the fragments onto the protostar, is used as a basic model for FUors. Typical FUor bursts (L ～ 100 L ^⊙) can appreciably influence the thermal balance of their disks and parent envelopes, leading to an increase in the disk temperature by more than a factor of two. On the other hand, massive fragments in the disk are only weakly perturbed by such bursts, partially due to screening by the disk and partially due to their high temperature brought about by adiabatic heating. Apart from massive fragments, the characteristic thermal time scales are appreciably shorter than the dynamical time scales throughout the radial extent of the disk and envelope; this enables the use of a stationary radiative-transfer equation when determining the vertical structure of the disk.     

The Thermal Structure of the Upper Mantle in Eastern China——Inferred from the Petrological Model

Cenozoic basalt in eastern China contains abundant ultramafic xenoliths which are specimens of pyrolitesreleased during basaltic magma eruption. A total of 405 P-T data of pyroxene in the ultramafic rocks have beencollected, which present a more precise pyroxene geotherm. The average geothermal gradient in the upper man-tle represented by the pyroxene geotherm is about 3.3℃ / km, which is much less than that derived from theconductive thermal model (≈14℃ / km), implying the great significance of convective heat transfer. The calcu-lation shows that the contributions of convective and conductive heat transfers are 79% and 21%, respectively.The perturbation in the thermal structure of the upper mantle is an important manifestation of thetectonothermal event of Cenozoic continental rifting and intense basaltic volcanism in eastern China. Based onthe pyroxene geotherm and its comparison with the current geothermal field derived from the measurements ofthe surface heat flows, it is suggested that the Moho may be a secondary thermal boundary. The currentgeothermal field and the thermal structure of the lithosphere in eastern China may mainly reflect the result ofthe tectonothermal disturbance in the Neogene-Quaternary, in other words, the lithosphere has just begun toCool.     

一个网格型松散结构分布式水文模型的构建

根据流域降雨径流的基本过程,以蓄满产流理论为基础,建立了一个网格型松散结构的分布式流域水文模型。模型将流域离散为包含河道与不包含河道两种类型的单元格,以协克里金方法插值得到空间离散的降雨输入,考虑的产汇流物理过程包括降雨、植被截留及蒸散发、单元格产流、单元格汇流及河网汇流。模型结构简单,参数较少,在充分利用植被覆盖类型图及土地利用类型图的基础上,能够获得大部分参数的选用值。通过在长江三峡区间沿渡河流域的实际应用,模型计算成果令人满意。     

On the thermal regime of some tectonic units in a continental collision environment in Romania

An analysis of the evolution of the main tectonic units in Romania shows that the thermal regime of the lithosphere should be derived according, on the one hand, to the particular tectonic interactions the various tectonic units have been involved in and, on the other hand, to the investigated depth interval. A steady-state conduction model of the crustal temperature field based on the heat flow distribution and information on the structure is presented for the Romanian territory. It shows large lateral thermal variations between tectonic units, as a result of different geological and thermal histories. Within the frame of a complex modelling of the thermal evolution of the lithosphere in the East Carpathians, as proposed in this study, the zero-order thermal effects of the pre-Miocene oceanic subduction of the Eurasian plate are evaluated. The deep thermal structure of the subducted slab is derived and shown to be compatible with the velocity structure of the lithosphere and the intermediate-depth seismicity of the Vrancea area.