Modeling of Hydrogen Blending on the Leakage and Diffusion of Urban Buried Hydrogen-Enriched Natural Gas Pipeline

With the introduction of various carbon reduction policies around the world,hydrogen energy,as a kind of clean energy with zero carbon emission,has attracted much attention.The safe and economical transportation of hydrogen is of great significance to the development of hydrogen energy industries.Utilizing natural gas pipelines to transport hydrogen is considered to be an efficient and economical way.However,hydrogen has a higher risk of leakage due to its strong diffusion capacity and lower explosive limit than conventional natural gas.Therefore,it is of great significance to study the leakage and diffusion law of hydrogen-enriched natural gas(HENG)pipelines for the safe transportation of hydrogen energy.In this study,the leakage and diffusion characteristics of urban buried HENG pipelines are investigated numerically,and the dangerous degree of leakage is analyzed based on the time and area when the gas concentration reaches the lower explosive limit.The influences of hydrogen blending ratio(HBR),operating pressure,leakage hole size and direction,as well as soil type on the leakage and diffusion law of HENG are analyzed.Results show that the hydrogen mixing is not the key factor in increasing the degree of risk after gas leakage for urban buried HENG pipelines.When the HBR is 5%,10%,15% and 20%,the corresponding first dangerous time is 1053,1041,1019 and 998 s,respectively.Thiswork is expected to provide a valuable reference for the safe operation and risk prevention of HENG pipelines in the future.

Coupled heat-fluid-solid numerical study on heat extraction potential of hot dry rocks based on discrete fracture network model

Fracture networks within hot dry rock(HDR)geothermal reservoirs are complex,and heat extraction via water injection is thus a coupled process of heat-fluid-solid multifield.In this paper,utilizing the theory of normally distributed random functions,we develop a corresponding pre-processing subprogram to establish a discrete network model of complex fracture distribution in HDR reservoirs;then construct a heat-fluid-solid finite element model for heat extraction via water injection and compare the numerical solution with the analytical solution of the one-dimensional non-isothermal consolidation problem for verification.The numerical simulation results show that the main factors affecting the heat extraction efficiency of HDR reservoirs include fracture width,fracture density,fracture permeability,and matrix permeability.When a HDR reservoir is injected with water for heat extraction,there is a certain threshold value of these influential parameters,beyond which the outlet temperature drops significantly,resulting in an obvious thermal breakthrough.When injecting water for heat extraction,the values of these parameters should be controlled and kept at a reasonable level,otherwise,the HDR reservoir may enter a thermal breakthrough stage in advance,which is not conducive for long-period heat extraction.Influenced by the random distribution of complex fractures,the leading edge of the cold front may present an irregular distribution.During the process of heat extraction,the stress gradually changes from a compressional state to a tensile state,which induces further damage to the HDR reservoir.

Failure Patterns and Mechanisms of Hydraulic Fracture Propagation Behavior in the Presence of Naturally Cemented Fractures

In this study,we use the extended finite element method(XFEM)with a consideration of junction enrichment functions to investigate the mechanics of hydraulic fractures related to naturally cemented fractures.In the proposed numerical model,the lubrication equation is adopted to describe the fluid flow within fractures.The fluid-solid coupling systems of the hydraulic fracturing problem are solved using the Newton-Raphson method.The energy release rate criterion is used to determine the cross/arrest behavior between a hydraulic fracture(HF)and a cemented natural fracture(NF).The failure patterns and mechanisms of crack propagation at the intersection of natural fractures are discussed.Simulation results show that after crossing an NF,the failure mode along the cemented NF path may change from the tensile regime to the shear or mixed-mode regime.When an advancing HF kinks back toward the matrix,the failure mode may gradually switch back to the tensile-dominated regime.Key factors,including the length of the upper/lower portion of the cemented NF,horizontal stress anisotropy,and the intersection angle of the crack propagation are investigated in detail.An uncemented or partially cemented NF will form a more complex fracture network than a cemented NF.This study provides insight into the formation mechanism of fracture networks in formations that contain cemented NF.

A new analytical-numerical method for calculating interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses

Based on the elementary solutions and new integral equations,a new analytical-numerical method is proposed to calculate the interacting stresses of multiple circular holes in an infinite elastic plate under both remote stresses and arbitrarily distributed stresses applied to the circular boundaries.The validity of this new analytical-numerical method is verified by the analytical solution of the bi-harmonic stress function method,the numerical solution of the finite element method,and the analytical-numerical solutions of the series expansion and Laurent series methods.Some numerical examples are presented to investigate the effects of the hole geometry parameters(radii and relative positions)and loading conditions(remote stresses and surface stresses)on the interacting tangential stresses and interacting stress concentration factors(SCFs).The results show that whether the interference effect is shielding(k<1)or amplifying(k>1)depends on the relative orientation of holes(α)and remote stresses(σ^∞x,σ^∞y).When the maximum principal stress is aligned with the connecting line of two-hole centers andσ^∞y<0.5σ^∞x,the plate containing two circular holes has greater stability than that containing one circular hole,and the smaller circular hole has greater stability than the bigger one.This new method not only has a simple formulation and high accuracy,but also has an advantage of wide applications over common analytical methods and analytical-numerical methods in calculating the interacting stresses of a multi-hole problem under both remote and arbitrary surface stresses.

Study on a Dual Embedded Discrete Fracture Model for Fluid Flow in Fractured Porous Media

Simulation of fluid flow in the fractured porous media is very important and challenging.Researchers have developed some models for fractured porous media.With the development of related research in recent years,the prospect of embedded discrete fracture model(EDFM)is more and more bright.However,since the size of the fractures in the actual reservoir varies greatly,a very fine grid should be used which leads to a huge burden to the computing resources.To address this challenge,in the present paper,an upscaling based model is proposed.In this model,the flow in large-scale fractures is directly described by the EDFM while that in the small-scale fractures is upscaled through local simulation by EDFM.The EDFM is used to simulate the large-and small-scale fractures independently two times,so the new model is called dual embedded discrete fracture model(D-EDFM).In this paper,the detailed implementation process of D-EDFM is introduced and,through test cases,it is found the proposed model is a feasible method to simulate the flow in fractured porous media.

LES Investigation of Drag-Reducing Mechanism of Turbulent Channel Flow with Surfactant Additives

In this work,the drag-reducing mechanism of high-Reynoldsnumber turbulent channel flow with surfactant additives is investigated by using large eddy simulation(LES)method.An N-parallel finitely extensible nonlinear elastic model with Peterlin’s approximation(FENE-P)is used to describe the rheological behaviors of non-Newtonian fluid with surfactant.To close the filtered LES equations,a hybrid subgrid scale(SGS)model coupling the spatial filter and temporal filter is applied to compute the subgrid stress and other subfilter terms.The finite difference method and projection algorithm are adopted to solve the LES governing equations.To validate the correctness of our LES method and in-house code,the particle image velocimetry(PIV)experiment is carried out and representative measured results are compared with LES results in detail.Then the flow characteristics and drag-reducing mechanism of turbulent channel flow with surfactant are investigated from the perspective of drag reduction rate,mean velocity,fluctuation of deformation rate,shear stress,transport and dissipation of turbulent kinetic energy,and turbulent coherent structures.This research can shed a light on the application of turbulent drag reduction technique in district heating,petroleum transport,etc.

智能化背景下中华优秀传统文化传播路径研究

在这个信息大爆炸的年代,信息科技已不只是少数人的专属,它已成为推动社会发展、文化繁荣的一股不可忽视的力量。在人工智能和大数据等前沿技术飞速发展的今天,中华传统文化的发展迎来了空前的机遇和挑战。文章从“智慧”这个崭新的角度,对“智慧”时代下中华优秀传统文化的传承与发扬进行深入研究。本文就现代信息技术对我国传统文化的冲击进行了剖析,并就其存在的问题进行了探讨,并提出了相应的对策与建议。力求在新形势下,把中华优秀传统文化创造性地发挥出来,再创辉煌。

基于非结构化VOSET方法的沸腾传热

将非结构化VOSET方法推广至求解气液界面存在相变的流动传热问题.为准确计算界面两侧因相变导致的能量跳跃,给出了含界面非结构网格控制单元温度的计算方法,并对非界面网格温度场采用隐式求解以提高计算精度.为验证所构建相变模型的正确性,编写程序分别模拟了恒壁面过热度和恒热流密度两种边界条件下的水平表面膜态沸腾,计算结果与Klimenko经验关联式吻合良好.通过圆弧表面膜态沸腾问题,验证了本文基于带相变非结构VOSET方法对不规则区域沸腾相变问题的适用性.通过模拟近临界压力水的膜态沸腾问题,并与Berenson和Klimenko经验关联式对比,验证了本文方法对实际沸腾问题的适用性.

Condensation heat transfer enhancement mechanism for vertical upflows by the phase separation concept at small gravity

In the field of aerospace, minimum and seal of equipments cause the increase in the thermal loading sharply. Due to the lack of driving force, the performance of conventional condenser deteriorates greatly under the small gravity environment, which leads to reduction in the service life of equipments. In this study, a passive condenser, developed on basis of the phase separation concept,is utilized to improve the performance of the condensation heat transfer under the small gravity environment. As a result of the limitation of experiments, the mechanisms of heat transfer enhancement of the phase separation condenser tube are revealed through numerical simulation based on the volume-of-fluid(VOF) method. The following conclusions could be obtained:(1) A novel phase distribution of ‘‘gas near the tube wall and liquid in the tube core'' is formed. The thin liquid film is indeed created after the flow pattern modulation by inserting mesh cylinder.(2)The condensation quantity for single bubble in the annular region increases about 16 times greater than that in the bare tube region in the case of Jl= 0.0574 m/s and Jg= 0.0229 m/s.(3) Gas volume fraction affects the parameters of liquid film thickness, bubble length and liquid bridge length. The increase in the gas volume fraction results in the decrease in the evaluation index from21.56 to 12.82. The evaluation index is defined as the ratio of the condensation quantities per unit tube length of the annular region and the bare tube region.

The bubble leakage mechanism for vertical upflows by the phase separation concept

Recently,our research group proposed the phase separation condenser tube,in which a mesh cylinder was inserted to form the flow structure of‘‘gas near the tube wall and liquid in the tube core’’,significantly enhance the condensation heat transfer.But the bubble leakage towards the core region may worsen the heat transfer enhancement.In order to prevent the bubble leakage,the critical criterion was proposed based on the Young–Laplace equation,considering the inertia force,viscous force and pulsating flow.It was found that the critical criterion depends on the dimensionless parameter G*,the We number and a coefficient C.The numerical model was developed in terms of the volume of fluid method to predict the two-phase laminar flow in the phase separation condenser tube.The results show that the bubble leakage takes place at the bubble tip,which is agreed with the experimental observations.The critical curve distinguishing the non-bubble-breaking and bubble-breaking was obtained by comparing the bubble dynamics at different G*and We.The coefficient C was determined.The critical criterion for the bubble leakage is given as G We0:22We0:99349 G4:7 103 We0:00651t196:39We0:006514cosa DHW,providing the design and operation guidance for the phase separation condenser tube.