CN110472343A - A kind of simulation method of roadbed ice-melt characteristic - Google Patents

A kind of simulation method of roadbed ice-melt characteristic Download PDF

Info

Publication number
CN110472343A
CN110472343A CN201910763904.2A CN201910763904A CN110472343A CN 110472343 A CN110472343 A CN 110472343A CN 201910763904 A CN201910763904 A CN 201910763904A CN 110472343 A CN110472343 A CN 110472343A
Authority
CN
China
Prior art keywords
ice
roadbed
temperature
indicate
melt
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910763904.2A
Other languages
Chinese (zh)
Other versions
CN110472343B (en
Inventor
苏文涛
赵文可
杨皓
张亚宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Liaoning Shihua University
Original Assignee
Shenyang Institute of Engineering
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenyang Institute of Engineering filed Critical Shenyang Institute of Engineering
Priority to CN201910763904.2A priority Critical patent/CN110472343B/en
Publication of CN110472343A publication Critical patent/CN110472343A/en
Application granted granted Critical
Publication of CN110472343B publication Critical patent/CN110472343B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Road Paving Structures (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The present invention proposes a kind of simulation method of roadbed ice-melt characteristic, Step 1: having ice-covered situation for road bed, establishes the numerical value computational geometry model of roadbed ice-melt, and obtain geometrical model parameter;Step 2: carrying out grid dividing to the geometrical model in step 1, and carry out the verifying of grid independence;Step 3: establishing the energy conservation equation of roadbed deicing processes;Step 4: setting boundary condition;Step 5: carrying out the unstable state two-dimensional surface numerical simulation calculation based on pressure using process of the SIMPLE algorithm to roadbed ice-melt, wherein energy conservation equation uses Second-order Up-wind discrete scheme, and the condition of convergence is defined as energy conservation equation convergence residual error.The simulation method calculation amount is small, and method is simple.

Description

A kind of simulation method of roadbed ice-melt characteristic
Technical field
The invention belongs to roadbed hot fluid ice-melt characteristic technical fields, more particularly to a kind of simulation of roadbed ice-melt characteristic Calculation method.
Background technique
Research of the hot fluid for roadbed snow melt originates in the 1950s, successively experienced using underground heat, solar energy The change procedure of hold over system.Japan is more early in the research starting of solar heat-preservation system this respect, in the side such as principle and application There is biggish progress in face.1997, Jilin University researcher took the lead in proposing that the snow melt of road Solar use and accumulation of energy exists Application in north China is imagined, and the work that conducts a research in terms of earth energy and deicing or snow melting always.When identical, University Of Tianjin is also opened The related work of road snow melt experiment and model calculating is opened up.Also the expansion in succession such as Harbin Institute of Technology, North China University of Tech Correlative study.But current some Model Calculating Methods are relative complex, computationally intensive, computational accuracy is low.
Summary of the invention
The invention aims to solve the problems of the prior art, the simulation for proposing a kind of roadbed ice-melt characteristic is calculated Method.
The present invention is achieved by the following technical solutions, and the present invention proposes a kind of simulation calculating side of roadbed ice-melt characteristic Method the described method comprises the following steps:
Step 1: having ice-covered situation for road bed, the numerical value computational geometry model of roadbed ice-melt is established, and obtain Take geometrical model parameter;
Step 2: carrying out grid dividing to the geometrical model in step 1, and carry out the verifying of grid independence;
Step 3: establishing energy conservation equation:
Wherein, E indicates hot fluid gross energy, hjIndicate the enthalpy of component j, keffIndicate effective thermal conductivity, JjIndicate component j Diffusion flux, ShIndicating that chemical reaction heat and volumetric sources item, ρ indicate density, t indicates the time, and ▽ is Hamilton operator, For indicating gradient and divergence, μ indicates dynamic viscosity, and p indicates pressure, and T indicates temperature, τeffIndicate virtual viscosity stress;
Step 4: setting boundary condition;
Step 5: carrying out the unstable state two-dimensional surface numerical value based on pressure using process of the SIMPLE algorithm to roadbed ice-melt Simulation calculates, and wherein energy conservation equation uses Second-order Up-wind discrete scheme, and the condition of convergence is defined as energy conservation equation convergence Residual error.
Further, the geometrical model parameter includes width, height, pipe laying depth and pipe laying radius.
Further, for roadbed ice-melt, boundary condition is specifically configured to: roadbed upper surface is contacted with ice cube, ice cube Upper surface is directly contacted with external environment, so geometrical model upper surface uses third boundary condition;Convection transfer rate h takes 18;The environment temperature that free convection temperature takes UDF to import;Concrete sample and ice cube two sides in roadbed use symmetrical border Condition;Geometrical model bottom surface uses isothermy, and temperature uses original ambient temperature;Geometrical model heating wall uses isothermal item Part, temperature is using experiment supply water temperature;Concrete sample and ice cube coupling face thickness in roadbed take 0.
Further, the convection transfer rate h is calculated with following formula:
Wherein a, b, n indicate that design factor, v indicate wind speed.
Detailed description of the invention
Fig. 1 is roadbed ice-melt geometrical model schematic diagram;
Fig. 2 is that the grid independence of ice-melt simulation verifies schematic diagram;
Fig. 3 is roadbed ice-melt pavement temperature experiment value and analogue value comparison diagram;
Fig. 4 is the influence schematic diagram of supply water temperature road pavement temperature in roadbed ice-melt;
Fig. 5 is temperature change cloud atlas in roadbed ice-melt;
Fig. 6 is the influence schematic diagram of pipe laying spacing road pavement temperature in roadbed ice-melt;
Fig. 7 is the influence schematic diagram of pipe laying depth road pavement temperature in roadbed ice-melt;
Fig. 8 is the influence schematic diagram of roadbed ice-melt medium floe quality road pavement temperature;
Fig. 9 is the influence schematic diagram of environment temperature road pavement temperature in roadbed ice-melt.
Specific embodiment
Technical solution in the embodiment of the present invention that following will be combined with the drawings in the embodiments of the present invention carries out clear, complete Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on this Embodiment in invention, every other reality obtained by those of ordinary skill in the art without making creative efforts Example is applied, shall fall within the protection scope of the present invention.
The present invention proposes a kind of simulation method of roadbed ice-melt characteristic, the described method comprises the following steps:
Step 1: having ice-covered situation for road bed, the numerical value computational geometry model of roadbed ice-melt is established;And it obtains Take geometrical model parameter;The geometrical model parameter includes width, height, pipe laying depth and pipe laying radius.It arranges and advises because of pipe laying Rule, by the heat flow density at two pipe laying centers be regarded as it is identical cancel out each other, can be as shown in Figure 1 by geometric model simplification.
Each roadbed model is the ideal heating model of each concrete sample, has different width, height, pipe laying depth, Specific geometrical model parameter is as shown in table 1.
1 roadbed model parameter list of table
Step 2: carrying out grid dividing to the geometrical model in step 1, and carry out the verifying of grid independence;
As shown in Fig. 2, being the pavement temperature of concrete sample in roadbed 4 with the increased change curve of grid number, ice cube thickness Maximum gauge 13mm when degree takes experiment.When grid number is 370,1782,3357,5016,7980,14573,33684, road surface temperature Degree first increases to be basically unchanged afterwards, and after grid number increases to 7980, analog result is held essentially constant, so coagulation in roadbed 4 Native test specimen selects the grid of grid number 7980 as size of mesh opening is calculated, at this time grid node interval 0.75mm.
By taking ice thickness 8mm as an example, grid node interval takes 0.75mm, in this case the number of nodes of each roadbed and Unit number is as shown in table 2.
The number of nodes and unit number of each roadbed of table 2
Step 3: establishing energy conservation equation:
Wherein, E indicates hot fluid gross energy, hjIndicate the enthalpy of component j, keffIndicate effective thermal conductivity, JjIndicate component j Diffusion flux, ShIndicating that chemical reaction heat and volumetric sources item, ρ indicate density, t indicates the time, and ▽ is Hamilton operator, For indicating gradient and divergence, μ indicates dynamic viscosity, and p indicates pressure, and T indicates temperature, τeffIndicate virtual viscosity stress;
Step 4: setting boundary condition;For roadbed ice-melt, boundary condition is specifically configured to: roadbed upper surface and ice Block contact, ice cube upper surface are directly contacted with external environment, so geometrical model upper surface uses third boundary condition;Convection current Coefficient of heat transfer h takes 18;The environment temperature that free convection temperature takes UDF to import;Concrete sample and ice cube two sides in roadbed Using symmetrical boundary condition;Geometrical model bottom surface uses isothermy, and temperature uses original ambient temperature;Geometrical model heated wall Face uses isothermy, and temperature is using experiment supply water temperature;Concrete sample and ice cube coupling face thickness in roadbed take 0.
The convection transfer rate h is calculated with following formula:
Wherein a, b, n indicate that design factor, v indicate wind speed.
Step 5: carrying out the unstable state two-dimensional surface numerical value based on pressure using process of the SIMPLE algorithm to roadbed ice-melt Simulation calculates, and wherein energy conservation equation uses Second-order Up-wind discrete scheme, and the condition of convergence is defined as energy conservation equation convergence Residual error 10-6
Simulation operating condition is consistent with experiment condition, and supply water temperature 323K, total working medium flow 18L/min, environment temperature is at any time Variation.It simulates under operating condition, the pavement temperature of the concrete sample in roadbed 1 is as shown in Figure 3 with the change curve of runing time.For Guarantee the reliability that numerical value calculates, the mathematical model of roadbed ice-melt will be verified.From the figure 3, it may be seen that under simulation operating condition Pavement temperature is identical as the variation tendency that experiment condition measures, numerical value is close.Period error is larger in simulations, and max value of error is 1.8%, after appearing in experiment operation 275 minutes, i.e. 12:43, solar irradiation is obvious at this time, therefore pavement temperature increases.It is overall From the point of view of, the error of experiment value and the analogue value is within zone of reasonableness, it was demonstrated that the reliability of mathematical model ensure that next step Numerical simulation.
The influence of supply water temperature road pavement temperature
Pavement temperature is as shown in Figure 4 with the change curve of runing time under different supply water temperatures.Operating condition is simulated between pipe laying Away from 120mm, pipe laying depth 80mm, environment temperature 263K, ice cube quality 1.5kg, supply water temperature is respectively 303,313,323K, fortune The row time 570 minutes.The temperature cloud picture of different time nodes roadbed and ice cube such as Fig. 5 institute in supply water temperature 313K simulation process Show, roadbed with the direct contact area of pipe laying to borderline region from being gradually warmed up.As shown in Figure 4, under identical supply water temperature, road surface temperature Degree is gradually increased with the increase of runing time, and temperature rise rate is gradually reduced with the increase of runing time.This is because with operation The temperature of the increase of time, concrete sample is gradually increasing, and the temperature difference of heating agent and refrigerant is gradually reduced, and heat-transfer capability weakens, institute It is reduced with temperature rise rate.Under different supply water temperatures, supply water temperature is higher, and the temperature rise rate of same time point is bigger.First 180 minutes The stage supply water temperature 303K heating rate that is rapidly heated is 0.050K/min, and supply water temperature is promoted to after 313K and 323K, temperature Raising speed rate has respectively reached 0.063K/min and 0.077K/min, and 27% and 21% has been respectively increased in temperature rise rate.This is because Supply water temperature is higher, and the temperature difference of heating agent and refrigerant is bigger, and heat-transfer capability is stronger, so temperature rise rate is higher.In addition, with water supply Temperature rises in gradient, and the pavement temperature difference of same time point also rises in gradient, in a linear relationship with supply water temperature difference, Temperature gap at 570 minutes has reached 2.9K.
The influence of pipe laying spacing road pavement temperature
Pavement temperature is as shown in Figure 6 with the change curve of runing time under different pipe laying spacing.Simulating operating condition is for water temperature 313K is spent, pipe laying depth 80mm, ice cube quality 1.5kg, environment temperature 263K, pipe laying spacing is respectively 120,140,160mm, fortune The row time 570 minutes.It will be appreciated from fig. 6 that pavement temperature is gradually increased with the increase of runing time, temperature rise under identical pipe laying spacing Rate is gradually reduced with runing time increase.It is rapidly heated the stage within first 180 minutes, the heating rate of pipe laying spacing 160mm is 0.049K/min, when pipe laying spacing is reduced to 140mm and 120mm, temperature rise rate then respectively reached 0.055K/min and 0.063K/min, temperature rise rate have been respectively increased 12.7% and 15.2%, have not decayed.Different pipe laying spacing, pipe laying spacing Smaller, the temperature rise rate of same time point is first big after small.Because pipe laying spacing is smaller, single pipe laying heating region is smaller, heating effect Fruit is better, and the variation of pavement temperature is faster, and the time for reaching dynamic equilibrium is shorter.In addition, same time point, pipe laying spacing The pavement temperature difference of 120mm and 140mm be greater than pipe laying spacing 140mm and 160mm pavement temperature difference, but with operation when Between increase tend towards stability.
The influence of pipe laying depth road pavement temperature
Pavement temperature is as shown in Figure 7 with the change curve of runing time under different pipe laying depths.Simulating operating condition is for water temperature 313K is spent, pipe laying spacing 120mm, ice cube quality 1.5kg, environment temperature 263K, pipe laying depth is respectively 60,80,100mm, fortune The row time 570 minutes.As shown in Figure 7, under identical pipe laying depth, pavement temperature is gradually increased with the increase of runing time, temperature rise Rate is gradually reduced with the increase of runing time.It is rapidly heated the stage within first 180 minutes, the temperature rise rate of pipe laying depth 100mm is 0.043K/min, and as pipe laying depth is by gradient reduction 20mm, temperature rise rate is respectively 0.063K/min and 0.087K/min, 47% and 37% has been respectively increased in temperature rise rate.Different pipe laying depths, pipe laying depth is smaller, and the variation of pavement temperature is faster, arrives Time up to dynamic equilibrium is shorter.In addition, the pavement temperature difference of pipe laying spacing 60mm and 80mm be greater than pipe laying spacing 80mm and The pavement temperature difference of 100mm, mid-term are reduced close to difference, and the later period continues to expand again.
The influence of ice cube quality road pavement temperature
Pavement temperature is as shown in Figure 8 with the change curve of runing time under different ice cube quality.Simulating operating condition is for water temperature 313K, pipe laying spacing 120mm, pipe laying depth 80mm, environment temperature 263K are spent, ice cube quality is respectively 1.5,2.0,2.5kg, fortune The row time 570 minutes.As shown in Figure 8, under identical ice cube quality, pavement temperature is gradually increased with the increase of runing time, temperature rise Rate is gradually reduced with the increase of runing time.It is rapidly heated the stage within first 180 minutes, the temperature rise rate of ice cube quality 1.5kg is 0.0632K/min, when ice cube quality increases to 2.0kg and 2.5kg, temperature rise rate is then respectively 0.0633K/min and 0.063K/ Min, temperature rise rate have been respectively increased 0.1% and have reduced 0.5%.Insulation effect is not played in ice cube covering.Different ice cube matter Amount, ice cube quality is smaller, and the variation of pavement temperature is faster, and the time for reaching dynamic equilibrium is shorter.
The influence of environment temperature road pavement temperature
Pavement temperature is as shown in Figure 9 with the change curve of runing time at a temperature of varying environment.Simulating operating condition is for water temperature 313K, pipe laying spacing 120mm, pipe laying depth 80mm, ice cube quality 1.5kg are spent, environment temperature is respectively 268,263,258K, fortune The row time 570 minutes.As shown in Figure 9, under identical environment temperature, pavement temperature is gradually increased with the increase of runing time, temperature rise Rate is gradually reduced with the increase of runing time.It is rapidly heated the stage within first 180 minutes, the temperature rise rate of environment temperature 268K is 0.058K/min, when environment temperature is reduced to 263K and 258K, temperature rise rate is respectively 0.063K/min and 0.069K/ 9.4% and 8.6% has been respectively increased in min, temperature rise rate.At a temperature of varying environment, environment temperature is lower, the temperature of same time point Raising speed rate is bigger.In addition, the pavement temperature difference of same time point then rises in gradient as environment temperature declines in gradient, It is in a linear relationship with environment temperature difference.
Above to a kind of simulation method of roadbed ice-melt characteristic proposed by the invention, it is described in detail, this Apply that a specific example illustrates the principle and implementation of the invention in text, the explanation of above example is only intended to It facilitates the understanding of the method and its core concept of the invention;At the same time, for those skilled in the art, think of according to the present invention Think, there will be changes in the specific implementation manner and application range, in conclusion the content of the present specification should not be construed as pair Limitation of the invention.

Claims (4)

1. a kind of simulation method of roadbed ice-melt characteristic, it is characterised in that: the described method comprises the following steps:
Step 1: having ice-covered situation for road bed, the numerical value computational geometry model of roadbed ice-melt is established, and obtain several What model parameter;
Step 2: carrying out grid dividing to the geometrical model in step 1, and carry out the verifying of grid independence;
Step 3: establishing energy conservation equation:
Wherein, E indicates hot fluid gross energy, hjIndicate the enthalpy of component j, keffIndicate effective thermal conductivity, JjIndicate the expansion of component j Dissipate flux, ShIndicate that chemical reaction heat and volumetric sources item, ρ indicate density, t indicates the time, and ▽ is Hamilton operator, is used to Indicate gradient and divergence, μ indicates dynamic viscosity, and p indicates pressure, and T indicates temperature, τeffIndicate virtual viscosity stress;
Step 4: setting boundary condition;
Step 5: carrying out the unstable state two-dimensional surface numerical simulation based on pressure using process of the SIMPLE algorithm to roadbed ice-melt It calculates, wherein energy conservation equation uses Second-order Up-wind discrete scheme, and the condition of convergence is defined as energy conservation equation convergence residual error.
2. according to the method described in claim 1, it is characterized by: the geometrical model parameter includes width, height, pipe laying depth Degree and pipe laying radius.
3. according to the method described in claim 1, it is characterized by: boundary condition is specifically configured to: road for roadbed ice-melt Base upper surface is contacted with ice cube, and ice cube upper surface is directly contacted with external environment, so geometrical model upper surface uses third class Boundary condition;Convection transfer rate h takes 18;The environment temperature that free convection temperature takes UDF to import;Concrete sample in roadbed Symmetrical boundary condition is used with ice cube two sides;Geometrical model bottom surface uses isothermy, and temperature uses original ambient temperature;It is several What model heating wall uses isothermy, and temperature is using experiment supply water temperature;Concrete sample and ice cube coupling in roadbed Face thickness takes 0.
4. according to the method described in claim 3, it is characterized by: the convection transfer rate h is calculated with following formula:
Wherein a, b, n indicate that design factor, v indicate wind speed.
CN201910763904.2A 2019-08-19 2019-08-19 Simulation calculation method for ice melting characteristics of roadbed Active CN110472343B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910763904.2A CN110472343B (en) 2019-08-19 2019-08-19 Simulation calculation method for ice melting characteristics of roadbed

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910763904.2A CN110472343B (en) 2019-08-19 2019-08-19 Simulation calculation method for ice melting characteristics of roadbed

Publications (2)

Publication Number Publication Date
CN110472343A true CN110472343A (en) 2019-11-19
CN110472343B CN110472343B (en) 2023-09-08

Family

ID=68511125

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910763904.2A Active CN110472343B (en) 2019-08-19 2019-08-19 Simulation calculation method for ice melting characteristics of roadbed

Country Status (1)

Country Link
CN (1) CN110472343B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007179501A (en) * 2005-12-28 2007-07-12 Sharp Corp Thermal fluid simulation device, thermal fluid simulation method, thermal fluid simulation program, and storage medium recording the program
JP2015215772A (en) * 2014-05-12 2015-12-03 株式会社東芝 Heat transfer simulation device and heat transfer simulation method
CN107784161A (en) * 2017-09-27 2018-03-09 北京理工大学 A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007179501A (en) * 2005-12-28 2007-07-12 Sharp Corp Thermal fluid simulation device, thermal fluid simulation method, thermal fluid simulation program, and storage medium recording the program
JP2015215772A (en) * 2014-05-12 2015-12-03 株式会社東芝 Heat transfer simulation device and heat transfer simulation method
CN107784161A (en) * 2017-09-27 2018-03-09 北京理工大学 A kind of analysis method of the compressible supercavity flow dynamic characteristic of high speed

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
姜宝成等: "内融冰式蓄冰管融冰过程中自然对流的影响", 《哈尔滨工业大学学报》 *
徐慧宁等: "流体加热道路融雪系统温-湿耦合融雪模型", 《中国公路学报》 *

Also Published As

Publication number Publication date
CN110472343B (en) 2023-09-08

Similar Documents

Publication Publication Date Title
Bauer et al. Transient 3D analysis of borehole heat exchanger modeling
Jiang et al. A novel three-dimensional transient model for subsurface heat exchange in enhanced geothermal systems
Al‐Khoury et al. Efficient finite element formulation for geothermal heating systems. Part II: transient
Kandaswamy et al. Magnetoconvection in an enclosure with partially active vertical walls
Cheng Soret and Dufour effects on natural convection heat and mass transfer from a vertical cone in a porous medium
CN104899351B (en) Bitumen mixture specimen temperature field analysis method based on discrete element simulation
CN109992846B (en) Simulation method for solar cross-season buried pipe heat storage
Ghasemi et al. Natural convection between a circular enclosure and an elliptic cylinder using control volume based finite element method
Lhendup et al. In-situ measurement of borehole thermal properties in Melbourne
Wołoszyn et al. Modelling of a borehole heat exchanger using a finite element with multiple degrees of freedom
Bartrons et al. A finite volume method to solve the frost growth using dynamic meshes
CN107391789A (en) Silicon melt reconstructing temperature field method based on free surface measured temperature and characteristic function interpolation
Zhou et al. Determination and analysis of parameters for an in-situ thermal response test
CN110472343A (en) A kind of simulation method of roadbed ice-melt characteristic
CN104594879A (en) Temperature control device and method for oil and gas field physical simulation experiments
Schiavi 3D simulation of the thermal response test in a U-tube borehole heat exchanger
Terzi et al. Experimental investigation on the evaporation of a wet porous layer inside a vertical channel with resolution of the heat equation by inverse method
Jahangir Heat and mass transfer analysis in unsaturated ground soils around the buried heating pipe
Liu et al. Analysis and experimental study on thermal dispersion effect of small scale saturated porous aquifer
CN113468743B (en) Medium-deep buried pipe fluid temperature field analysis method considering groundwater seepage
Van de Wal The importance of thermodynamics for modeling the volume of the Greenland ice sheet
Vasilyev et al. Estimation of the thermal effect of ground moisture condensation on heat transfer outside a geothermal borehole
CN111400833A (en) Verilog-A macro model modeling method for thermal wind speed sensor
Maestre et al. Experimental validation of a hybrid model for vertical ground heat exchangers under on-off operation conditions
Tabarsa et al. A numerical investigation of the effect of the temperature on the seepage calculation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230323

Address after: Wanghua District West Dandong Road 113001 Liaoning city of Fushun province No. 1

Applicant after: Liaoming Petrochemical University

Address before: 110136 Pu Chang Road, Shenbei New Area, Shenyang, Liaoning Province, No. 18

Applicant before: SHENYANG INSTITUTE OF ENGINEERING

GR01 Patent grant
GR01 Patent grant