CN116187107B - Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium - Google Patents
Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium Download PDFInfo
- Publication number
- CN116187107B CN116187107B CN202310469070.0A CN202310469070A CN116187107B CN 116187107 B CN116187107 B CN 116187107B CN 202310469070 A CN202310469070 A CN 202310469070A CN 116187107 B CN116187107 B CN 116187107B
- Authority
- CN
- China
- Prior art keywords
- field
- temperature
- abnormal
- space
- domain
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Abstract
The application discloses a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which are used for carrying out space discretization and time discretization on a target area, and iteratively calculating the total field temperature of a space domain by combining a background field control equation and an abnormal field control equation in space; the space-wave number domain abnormal field temperature of the next time node can be obtained by combining the ground temperature field recursion formula of the explicit differential format according to the space-wave number domain abnormal field temperature of the current time node, and then the space domain total field temperature of the next time node is calculated in an iterative mode; and so on to get the total field temperature of the spatial domain for all time nodes. The three-dimensional ground temperature field dynamic numerical simulation is realized, the dynamic change of the ground temperature field can be truly reflected, the background field control equation, the abnormal field control equation and the ground temperature field recurrence technology are combined, the calculation and storage cost is reduced while the fine simulation of the temperature field dynamic numerical value is ensured, and the calculation efficiency is improved.
Description
Technical Field
The application relates to the technical field of ground temperature field numerical simulation, in particular to a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium.
Background
Geothermal is a renewable energy source with abundant reserves, wide distribution and green and clean, and plays an important role in the sustainable development of society. The research of a fine and efficient three-dimensional ground temperature field numerical simulation method is of great significance in grasping the temperature distribution in the earth, evaluating the geothermal resource quantity and developing geothermal energy.
At present, most of numerical simulation methods for three-dimensional ground temperature fields are steady-state numerical simulation, for example, chinese patent application CN202111423020.6 discloses a three-dimensional ground temperature field numerical simulation method based on a heat conductivity anisotropic medium, and Chinese patent application CN202111423009.X discloses a forward modeling method, device, equipment and medium of a three-dimensional steady-state heat conduction ground temperature field, which are steady-state numerical simulation methods. However, the ground temperature field is dynamically changing in nature, i.e., transient, and dynamic simulation of the ground temperature field is of great practical significance to geothermal research. At present, although researches on transient ground temperature fields are carried out, the problems of large calculated amount and high storage requirement in the dynamic simulation of the three-dimensional ground temperature field numerical values under large-scale complex geological conditions generally exist in the technology, and the cost is high. Aiming at the problem, the application provides a fine and efficient three-dimensional ground temperature field dynamic numerical simulation technology, which provides important technical support for fine inversion of large-scale ground temperature data.
Disclosure of Invention
Aiming at the defects in the prior art, the application aims to provide a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which can realize fine and efficient simulation of the three-dimensional ground temperature field dynamic numerical value.
In a first aspect, a method for dynamically and numerically simulating a three-dimensional geothermal field is provided, including:
s1: selecting a target area containing an abnormal body, and constructing a target model;
s2: performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes;
s3: giving a thermophysical parameter value on each sampling node, wherein the thermophysical parameter value comprises a background thermophysical parameter value at a sampling node without an abnormal body and an abnormal background thermophysical parameter value at a sampling node with an abnormal body;
s4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value;
s5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation;
s6: according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value, calculating to obtain the abnormal field temperature of the space-wave number domain by combining a space-wave number domain abnormal field control equation;
s7: performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain;
s8: obtaining a new spatial domain total field temperature based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
s9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, taking the new space domain total field temperature as the current round space domain total field temperature in the next round of iteration, and returning to the step S6;
s10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
According to the first aspect, in a possible implementation manner, in the step S2, when the spatial discretization is performed, the target area is spatially subjected to three-dimensional mesh subdivision, and the x direction and the y direction are uniformly subdivided, and the z direction is uniformly or non-uniformly subdivided.
In a possible implementation manner according to the first aspect, the thermophysical parameter values include thermal conductivity, thermal capacity, heat generation rate.
In one possible implementation form according to the first aspect, the background field control equation is expressed as follows:
in the method, in the process of the application,for background field temperature, +.>For background thermal conductivity, +.>For background heat generation rate, < >>Is background heat capacity, t is time.
In one possible implementation manner according to the first aspect, the step S5 includes:
the abnormal field control equation is expressed as follows:
in the method, in the process of the application,for abnormal field temperature, +.>For background thermal conductivity, +.>For background heat capacity, +.>For an abnormal thermal conductivity,for abnormal heat generation rate, < >>Is the abnormal heat capacity, T is the total field temperature, < >>Is the heat capacity of the fluid, v is the flow rate of the fluid, and t is the time;
performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation, wherein the space-wave number domain abnormal field control equation is expressed as follows:
in the method, in the process of the application,abnormal field temperature for the space-wavenumber domain, +.>Abnormal heat generation rate for space-wave number domain, < >>Respectively->Wave number in direction, +.>Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
In one possible implementation manner, according to the first aspect, the iterative convergence condition is that an error is smaller than a preset value; the expression of the error is as follows:
wherein e is an error,the number of sampling nodes in the x, y and z directions is +.>、/>The spatial domain total field temperature for the nth and n+1th iterations, respectively.
According to the first aspect, in a possible implementation manner, in the step S10, a ground temperature field recurrence formula is expressed as follows:
in the method, in the process of the application,for the last time node->Is a spatially-wavenumber domain abnormal field temperature, +.>For the next time nodeK is a coefficient matrix related to the abnormal field of the space-wave number domain, G is a coefficient matrix related to the time derivative term of the abnormal field of the space-wave number domain, and P is a source term.
According to the first aspect, in a possible implementation manner, in the step S10, a catch-up method is used to solve the abnormal field temperature in the space-wave number domain of the next time node.
In a second aspect, there is provided an electronic device comprising:
a memory storing a computer program;
and the processor is used for loading and executing the computer program to realize the steps of the three-dimensional ground temperature field dynamic numerical simulation method.
In a third aspect, a computer readable storage medium is provided, storing a computer program which, when executed by a processor, implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method as described above.
The application provides a three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium, which are used for carrying out space discretization and time discretization on a target area, and iteratively calculating the total field temperature of a space domain by combining a background field control equation and an abnormal field control equation in space; the space-wave number domain abnormal field temperature of the next time node can be obtained by combining the ground temperature field recursion formula of the explicit differential format according to the space-wave number domain abnormal field temperature of the current time node, and then the space domain total field temperature of the next time node is calculated in an iterative mode; in this way, the total field temperature of the spatial domain of all time nodes can be obtained. The application realizes three-dimensional ground temperature field dynamic numerical simulation, is ground temperature field transient simulation, can truly reflect the dynamic change of a ground temperature field, combines a background field control equation, an abnormal field control equation and a ground temperature field recurrence technology, reduces calculation and storage costs while guaranteeing the fine simulation of the temperature field dynamic numerical value, improves calculation efficiency, and provides important technical support for fine inversion of large-scale ground temperature data.
Drawings
In order to more clearly illustrate the embodiments of the application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for simulating dynamic numerical values of a three-dimensional ground temperature field provided by an embodiment of the application;
FIG. 2 is a cross-sectional view of a target model according to an embodiment of the present application, wherein (a) is a cross-sectional view of the target model XOY, and (b) is a cross-sectional view of the target model XOZ;
FIG. 3 is a schematic diagram of the calculation results of the method and COMSOL Multiphysics software at a sampling node according to the embodiment of the present application;
FIG. 4 is a graph showing the relative error between the calculation results of the method and COMSOL Multiphysics software according to the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the technical solutions of the present application will be described in detail below. It will be apparent that the described embodiments are only some, but not all, embodiments of the application. All other embodiments, based on the examples herein, which are within the scope of the application as defined by the claims, will be within the scope of the application as defined by the claims.
The application simulates the dynamic value of the ground temperature field of each time node based on the dynamic simulation equation of the ground temperature field, and the dynamic simulation equation of the ground temperature field is expressed as follows:
(1)
in the method, in the process of the application,for heat conductivity, < >>Is the total field temperature, +.>Is the heat capacity of the fluid, v is the flow rate of the fluid, < >>The heat-generating rate of the heat-generating material,is the heat capacity of the medium, t is the time.
Based on the superposition principle, the formula (1) is split into a background field control equation and an abnormal field control equation, wherein the background field is set as a uniform layered medium, and the background field control equation is expressed as follows:
(2)
in the method, in the process of the application,for background field temperature, +.>For background thermal conductivity, +.>For background heat generation rate, < >>Is background heat capacity.
The abnormal field control equation is expressed as follows:
(3)
in the method, in the process of the application,for abnormal field temperature, +.>For abnormal thermal conductivity +.>For abnormal heat generation rate, < >>Is an abnormal heat capacity.
Based on the background field control equation and the abnormal field control equation, the embodiment of the application provides a three-dimensional ground temperature field dynamic numerical simulation method, as shown in fig. 1, comprising the following steps:
s1: and selecting a target area containing an abnormal body, and constructing a target model.
S2: and performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes.
In the embodiment, when the space discretization is performed, the target area is spatially subjected to three-dimensional grid subdivision, the x and y directions are uniform subdivision, and the z direction is uniform or non-uniform subdivision; the number of sampling nodes in the x direction isThe number of sampling nodes in the y direction is +.>The number of sampling nodes in the z direction is +.>. When time discretization is carried out, a series of time nodes are selected in time, and the number of the time nodes is +.>。
S3: the thermophysical parameter values on each sampling node are given, including a background thermophysical parameter value at a non-outlier sampling node and an outlier background thermophysical parameter value at an outlier sampling node.
In this embodiment, the thermophysical parameter values include thermal conductivityHeat capacity->And a heat generation rate A, wherein the abnormal heat physical parameter value is different from the surrounding heat physical parameter value. The heat conductivity, the heat capacity and the heat generation rate at the sampling node without abnormal body are background heat conductivity +.>Background Heat capacity->Background Heat production Rate->The background thermophysical parameter value is unchanged in the horizontal direction and only in the vertical direction. The thermal conductivity, the heat capacity and the heat generation rate at the abnormal body sampling node are abnormal thermal conductivity +.>Abnormal heat capacity->Abnormal heat generation rate->。
S4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value.
In particular, the boundary condition refers to that the temperature value or heat flux density value or heat exchange coefficient at the boundary is knownIt is known that. The initial condition is given initial timeThe initial temperature of each sampling node in the target area can be given as original data, and can be calculated based on boundary conditions and model parameters. On the premise of loading initial conditions and boundary conditions, solving based on the background field control equation of the formula (2) and the background thermophysical parameter values of all sampling nodes to obtain the background field temperature of the space domain.
S5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation (namely a one-dimensional space domain abnormal field control equation under a two-dimensional wave number domain), wherein the method is expressed as follows:
in the method, in the process of the application,abnormal field temperature for the space-wavenumber domain, +.>Abnormal heat generation rate for space-wave number domain, < >>Respectively->Wave number in direction, +.>Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
S6: and solving by adopting a one-dimensional finite element method based on quadratic interpolation according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value and by combining a space-wave number domain abnormal field control equation to obtain the space-wave number domain abnormal field temperature.
S7: and performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain.
S8: and adding the space domain background field temperature and the space domain abnormal field temperature to obtain a new space domain total field temperature.
S9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, the new spatial domain total field temperature is taken as the current round of spatial domain total field temperature in the next round of iteration, and the step S6 is returned.
In this embodiment, the iteration convergence condition is that an error is smaller than a preset value; the expression of the error is as follows:
wherein e is an error,the number of sampling nodes in the x, y and z directions is +.>、/>The spatial domain total field temperature for the nth and n+1th iterations, respectively. If the error e between the space domain total field temperature obtained by the current iteration round and the space domain total field temperature obtained by the previous iteration round is smaller than a set value, meeting a convergence condition, ending the iteration and outputting the space domain total field temperature; otherwise, if the convergence condition is not satisfied, repeating the steps S6-S9 to continue the iterative computation until convergence.
S10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
Specifically, after the final spatial domain total field temperature of the current time node is calculated by utilizing the steps S6-S9, extracting the spatial-wave number domain abnormal field temperature in the last round of iterative calculation of the spatial domain total field temperature of the current time node, then solving the spatial-wave number domain abnormal field temperature of the next time node by combining a ground temperature field recurrence formula of an explicit differential format, performing two-dimensional Fourier inverse transformation on the spatial-wave number domain abnormal field temperature of the next time node to obtain the spatial domain abnormal field temperature of the next time node, adding the spatial domain abnormal field temperature of the next time node and the spatial domain background field temperature to obtain the spatial domain total field temperature of the current round of spatial domain total field temperature of the next time node, returning to the step S6, and executing the iterative calculation process of S6-S9 to obtain the final spatial domain total field temperature (namely the ground temperature field temperature) of the next time node. Step S10 is then performed and so on until the ground temperature field temperature of all time nodes is obtained.
Wherein, the ground temperature field recurrence formula is expressed as follows:
(5)
in the method, in the process of the application,for the last time node->Is a spatially-wavenumber domain abnormal field temperature, +.>For the next time nodeIs the temperature of the abnormal field in the space-wave number domain, K is the coefficient matrix related to the abnormal field in the space-wave number domain, G is the temperature of the abnormal field in the space-wave number domainCoefficient matrix related to abnormal field time derivative term in number domain, P is source term; in the process of solving the abnormal field control equation of the space-wave number domain by adopting the finite element method to obtain the abnormal field temperature of the space-wave number domain, the value of K, G, P can be obtained by carrying out unit analysis and overall synthesis. Equation (5) can be reduced to a diagonal system of equations +.>Wherein Q is a five-diagonal matrix, X is the unknown quantity to be solved, and B is the right-hand term. The diagonal equation set is quickly solved by adopting the catch-up method, and the ordinary differential equation is solved by adopting the catch-up method, so that the method has high parallelism, and the calculation time and the memory occupation are greatly reduced.
The embodiment of the application also provides electronic equipment, which comprises:
a memory storing a computer program;
and the processor is used for loading and executing the computer program to realize the steps of the three-dimensional ground temperature field dynamic numerical simulation method according to the embodiment.
The embodiment of the application also provides a computer readable storage medium storing a computer program which, when executed by a processor, implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method as described in the above embodiment.
The effect of the solution provided by the application is examined in connection with a specific example.
The computer tested was configured as an Intel (R) Core (TM) i7-11800H, main frequency 2.30GHz, memory 16GB,64 bit win11 system.
The projection of the target model XOY plane and XOZ cross section is shown in fig. 2, the background is a full space medium, the background thermal conductivity is 2W/(m.deg.c), and the thermal capacity is 5×10 5 J/(m 3 DEG C), a fluid abnormal body is arranged, the abnormal heat conductivity is 1.0W/(m DEG C), the fluid flow direction is vertical upwards, and the flow velocity is 1 multiplied by 10 -9 m/s, abnormal heat capacity of 1X 10 6 J/(m 3 at..degree.C.) no internal heat was generated in the background, and the abnormal heat generation rate was 4X 10 -6 W/m 3 . Calculating the range of x-5000 to 5000m, y-5000 to 5000m and z-0 to 8000m, the fluid flow range is from x direction to 1000m, the y direction is from 1000m to 1000m, and the z direction is from 3000 m to 700 m. The upper boundary adopts a first type boundary condition (the temperature value is known), the initial temperature value of the boundary is 10 ℃, the lower boundary adopts a second type boundary condition (the heat flux density value is known), and the heat flux density value is 41.86 mW/m 2 . The spatial grid subdivision is 51 x 51, the number of time nodes is 11. In the single-threaded case, a standard Fourier transform is adopted, and the iterative convergence accuracy (preset value) is 10 -4 . Calculation by the method and COMSOL Multiphysics software of the application respectively, and FIG. 3 isThe calculation results of the method and COMSOL Multiphysics software of the application at the measuring points are shown in fig. 4, which is a graph of the relative errors of the two algorithms, and the relative errors of the two calculation results are less than 0.05%, thus illustrating the correctness of the method of the application. Under the same grid condition, the method takes 23 seconds and occupies 0.14GB of memory; COMSOL Multiphysics takes 426 seconds and occupies 10.48GB of memory, which illustrates the high efficiency of the method of the application.
It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While embodiments of the present application have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the application, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the application.
Claims (10)
1. The method for simulating the dynamic numerical value of the three-dimensional ground temperature field is characterized by comprising the following steps of:
s1: selecting a target area containing an abnormal body, and constructing a target model;
s2: performing spatial discretization and time discretization on the target model to obtain a series of sampling nodes and a series of time nodes;
s3: giving a thermophysical parameter value on each sampling node, wherein the thermophysical parameter value comprises a background thermophysical parameter value at a sampling node without an abnormal body and an abnormal background thermophysical parameter value at a sampling node with an abnormal body;
s4: loading initial conditions and boundary conditions, and taking the background field temperature of the space domain as the total field temperature of the initial space domain based on a background field control equation and a background thermophysical parameter value;
s5: performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation;
s6: according to the total field temperature of the current round space domain and the abnormal background thermophysical parameter value, calculating to obtain the abnormal field temperature of the space-wave number domain by combining a space-wave number domain abnormal field control equation;
s7: performing two-dimensional Fourier inverse transformation on the abnormal field temperature in the space-wave number domain to obtain the abnormal field temperature in the space domain;
s8: obtaining a new spatial domain total field temperature based on the spatial domain background field temperature and the spatial domain abnormal field temperature;
s9: judging whether the new space domain total field temperature meets the iteration convergence condition, if so, outputting the new space domain total field temperature as the ground temperature of the current time node; otherwise, taking the new space domain total field temperature as the current round space domain total field temperature in the next round of iteration, and returning to the step S6;
s10: combining a ground temperature field recursion formula in an explicit differential format, obtaining space-wave number domain abnormal field temperature of the next time node from space-wave number domain abnormal field temperature in the last iteration of the current time node, performing two-dimensional Fourier inverse transformation on the space-wave number domain abnormal field temperature, combining space domain background field temperature to obtain space domain total field temperature as the current round space domain total field temperature of the next time node, and returning to the step S6; and the like until the ground temperature field temperature of all time nodes is obtained.
2. The method according to claim 1, wherein in the step S2, the target region is spatially three-dimensionally mesh-divided, the x-direction and the y-direction are uniformly divided, and the z-direction is uniformly or non-uniformly divided.
3. The method of claim 1, wherein the thermophysical parameter values include thermal conductivity, thermal capacity, and heat generation rate.
4. The method of claim 1, wherein the background field control equation is expressed as follows:
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->For background field temperature, +.>For background thermal conductivity, +.>For background heat generation rate, < >>Is background heat capacity, t is time.
5. The method of three-dimensional geothermal field dynamic numerical simulation according to claim 1, wherein the step S5 comprises:
the abnormal field control equation is expressed as follows:
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->For abnormal field temperature, +.>For background thermal conductivity, +.>For background heat capacity, +.>For abnormal thermal conductivity +.>For abnormal heat generation rate, < >>Is the abnormal heat capacity, T is the total field temperature, < >>Is the heat capacity of the fluid, v is the flow rate of the fluid, and t is the time;
performing horizontal two-dimensional Fourier transform on the abnormal field control equation to obtain a space-wave number domain abnormal field control equation, wherein the space-wave number domain abnormal field control equation is expressed as follows:
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->Abnormal field temperature for the space-wavenumber domain, +.>Abnormal heat generation rate for space-wave number domain, < >>Respectively->Wave number in direction, +.>Is a two-dimensional Fourier transform symbol, i is an imaginary unit, v x 、v y 、v z The flow rates of the fluid in the x, y, and z directions, respectively.
6. The method for simulating the dynamic numerical value of the three-dimensional ground temperature field according to claim 1, wherein the iterative convergence condition is that the error is smaller than a preset value; the expression of the error is as follows:
the method comprises the steps of carrying out a first treatment on the surface of the Wherein e is error, ">The number of sampling nodes in the x, y and z directions is +.>、/>The spatial domain total field temperature for the nth and n+1th iterations, respectively.
7. The method according to any one of claims 1 to 6, wherein in the step S10, the ground temperature field recurrence formula is as follows:
the method comprises the steps of carrying out a first treatment on the surface of the In (1) the->For the last time node->Is a spatially-wavenumber domain abnormal field temperature, +.>For the next timeInter-node->K is a coefficient matrix related to the abnormal field of the space-wave number domain, G is a coefficient matrix related to the time derivative term of the abnormal field of the space-wave number domain, and P is a source term.
8. The method according to claim 7, wherein in the step S10, the space-wave number domain abnormal field temperature of the next time node is solved by a catch-up method.
9. An electronic device, comprising:
a memory storing a computer program;
a processor for loading and executing the computer program to implement the steps of the three-dimensional earth temperature field dynamic numerical simulation method as claimed in any one of claims 1 to 8.
10. A computer readable storage medium storing a computer program, characterized in that the computer program when executed by a processor implements the steps of the three-dimensional earth temperature field dynamic numerical simulation method according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310469070.0A CN116187107B (en) | 2023-04-27 | 2023-04-27 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310469070.0A CN116187107B (en) | 2023-04-27 | 2023-04-27 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116187107A CN116187107A (en) | 2023-05-30 |
| CN116187107B true CN116187107B (en) | 2023-08-11 |
Family
ID=86452646
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310469070.0A Active CN116187107B (en) | 2023-04-27 | 2023-04-27 | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116187107B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012027848A1 (en) * | 2010-09-03 | 2012-03-08 | Geosoft Inc. | Method and system for modeling anomalous density zones in geophysical exploration |
| WO2015110599A1 (en) * | 2014-01-24 | 2015-07-30 | Ledaflow Technologies Da | Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines |
| CN108595762A (en) * | 2018-03-26 | 2018-09-28 | 中国科学院计算技术研究所 | A kind of three-dimensional dividing method and system of general circulation model dynamical frame |
| CN113792445A (en) * | 2021-11-15 | 2021-12-14 | 中南大学 | Three-dimensional magnetotelluric numerical simulation method based on integral equation method |
| WO2022000132A1 (en) * | 2020-06-28 | 2022-01-06 | 大连理工大学 | Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface |
| CN114036805A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field |
| CN114036806A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium |
| CN114707385A (en) * | 2022-04-12 | 2022-07-05 | 中国地质大学（武汉） | Deep stratum heat conductivity coefficient three-dimensional prediction method and device based on Krylov subspace |
| CN115659579A (en) * | 2022-09-05 | 2023-01-31 | 中南大学 | Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20220261517A1 (en) * | 2020-10-29 | 2022-08-18 | Jilin University | Simulation method for electromagnetic multi-scale diffusion |
-
2023
- 2023-04-27 CN CN202310469070.0A patent/CN116187107B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012027848A1 (en) * | 2010-09-03 | 2012-03-08 | Geosoft Inc. | Method and system for modeling anomalous density zones in geophysical exploration |
| WO2015110599A1 (en) * | 2014-01-24 | 2015-07-30 | Ledaflow Technologies Da | Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines |
| CN108595762A (en) * | 2018-03-26 | 2018-09-28 | 中国科学院计算技术研究所 | A kind of three-dimensional dividing method and system of general circulation model dynamical frame |
| WO2022000132A1 (en) * | 2020-06-28 | 2022-01-06 | 大连理工大学 | Method for designing and optimizing three-dimensional porous heat dissipation structure on basis of three-cycle minimal curved surface |
| CN113792445A (en) * | 2021-11-15 | 2021-12-14 | 中南大学 | Three-dimensional magnetotelluric numerical simulation method based on integral equation method |
| CN114036805A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field |
| CN114036806A (en) * | 2021-11-26 | 2022-02-11 | 中南大学 | Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium |
| CN114707385A (en) * | 2022-04-12 | 2022-07-05 | 中国地质大学（武汉） | Deep stratum heat conductivity coefficient three-dimensional prediction method and device based on Krylov subspace |
| CN115659579A (en) * | 2022-09-05 | 2023-01-31 | 中南大学 | Three-dimensional gravity field numerical simulation method and system based on 3D AS-FT |
Non-Patent Citations (1)
| Title |
|---|
| 基于有限差分法的大线能量焊接热模拟试件的温度分布;李矗东;魏强;李玉中;黄宏虎;;工业加热(01);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116187107A (en) | 2023-05-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hickel et al. | An adaptive local deconvolution method for implicit LES | |
| Selimefendigil et al. | Numerical study and pod-based prediction of natural convection in a ferrofluids–filled triangular cavity with generalized neural networks | |
| Guo et al. | Three-dimensional transient heat conduction analysis by Laplace transformation and multiple reciprocity boundary face method | |
| Dehghan et al. | Numerical study of three-dimensional Turing patterns using a meshless method based on moving Kriging element free Galerkin (EFG) approach | |
| Dalal et al. | Natural convection in a rectangular cavity heated from below and uniformly cooled from the top and both sides | |
| Chang et al. | Review of computational schemes in inverse heat conduction problems | |
| WO2021068901A1 (en) | Material microstructure evolution simulation solved on the basis of exponential time difference format | |
| Mirzaei et al. | MLPG method for transient heat conduction problem with MLS as trial approximation in both time and space domains | |
| CN113158538A (en) | Soft measurement method for heat flux density of complex-structure boiling surface | |
| CN115758841A (en) | Finite element reduced modeling method for digital twinning application-oriented transformer temperature field | |
| Li et al. | Meshless local Petrov Galerkin method for 2D/3D nonlinear convection–diffusion equations based on LS-RBF-PUM | |
| CN105760586A (en) | Fuzzy temperature response subordinating degree function solving method based on point collocation theory | |
| Zang et al. | Isogeometric boundary element for analyzing steady-state heat conduction problems under spatially varying conductivity and internal heat source | |
| CN116187107B (en) | Three-dimensional ground temperature field dynamic numerical simulation method, equipment and medium | |
| CN105677995A (en) | Method for numerical solution of fuzzy steady state heat conduction problem based on full grid point collocation theory | |
| Danko | Functional or operator representation of numerical heat and mass transport models | |
| Kadantsev et al. | The formulation and implementation of analytic energy gradients for periodic density functional calculations with STO/NAO Bloch basis set | |
| Rahideh et al. | Two-dimensional inverse transient heat conduction analysis of laminated functionally graded circular plates | |
| Wang et al. | Solving of two-dimensional unsteady inverse heat conduction problems based on boundary element method and sequential function specification method | |
| CN114036805A (en) | Forward modeling method, device, equipment and medium for three-dimensional steady-state heat conduction geothermal field | |
| CN105808508A (en) | Random orthogonal expansion method for solving uncertain heat conduction problem | |
| CN114036806A (en) | Three-dimensional geothermal field numerical simulation method based on thermal conductivity anisotropic medium | |
| CN116341279A (en) | Anisotropic medium three-dimensional transient earth temperature field forward modeling method, equipment and medium | |
| CN112084723B (en) | Optical fiber preform one-time stretching process simulation method and device | |
| Ala et al. | Smoothed particle electromagnetics modelling on HPC-GRID environment |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |