CN110472343A - A kind of simulation method of roadbed ice-melt characteristic - Google Patents
A kind of simulation method of roadbed ice-melt characteristic Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling 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
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.
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)
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 |
-
2019
- 2019-08-19 CN CN201910763904.2A patent/CN110472343B/en active Active
Patent Citations (3)
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)
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 |