CN106202757B - A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method - Google Patents
A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method Download PDFInfo
- Publication number
- CN106202757B CN106202757B CN201610561555.2A CN201610561555A CN106202757B CN 106202757 B CN106202757 B CN 106202757B CN 201610561555 A CN201610561555 A CN 201610561555A CN 106202757 B CN106202757 B CN 106202757B
- Authority
- CN
- China
- Prior art keywords
- heat
- ladle
- temperature
- model
- numerical value
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Treatment Of Steel In Its Molten State (AREA)
Abstract
The present invention proposes that a kind of numerical value of Con casting ladle space radiant heat flux distribution determines method, belong to field of steel metallurgy, this method will rely on continuous casting technique practical, heat source temperature and heat flux distribution are obtained using research methods such as mathematical computations and numerical simulations, and based on this exploitation continuous casting steel machine process radiations heat energy model, both new way was provided for continuous casting steel machine heat recovery, and also implemented based theoretical and oil storage method method for exploitation steel and iron manufacturing whole process energy conservation and environmental protection new technology and effective " energy-saving and emission-reduction "；Enterprise's manufacturing cost will be not only greatly reduced in this, also can greatly alleviate industrial energy saving emission reduction, and then generate very big economic results in society, correlative study has important theory and realistic meaning.
Description
Technical field
The invention belongs to field of steel metallurgy, and in particular to a kind of numerical value of Con casting ladle space radiant heat flux distribution determines
Method.
Background technique
Steel and iron industry is national economy mainstay type basic industry, and the operational characteristiies such as flow manufacturing and smelting processing are determined
The characteristics of having determined its energy consumption and remaining waste heat discharge rich and influential family；Statistics shows that steel manufacturing procces energy consumption accounts for about national work
Industry total energy consumption 10%~15%, energy consumption is huge, even and if manufacturing process release remaining waste heat in ignore necessary heat loss, still have
Extremely considerable heat is wasted and fails effectively to be recycled and reused；
It is encountered with gradually attention and the steel industry sustainable development for resource, the energy and environmental problem
Ecology, economical one problem, various remaining Technogy of Waste Gas Heat have obtained quick development.Such as: sintering waste heat recycles skill
Art, pelletizing waste heat cyclic utilization technology, blast funnace hot blast stove preheating technique, electric furnace, blast furnace flue gas heat recovery dedusting skill
Art etc..But due to being limited by technology and space, effective use efficiency is 10% to 30%.Such as coke and sinter
Sensible heat could only be recycled by gas-solid heat exchange method, production steam simultaneously generates electricity, effective use efficiency only has 17%；
Recycling for dross sensible heat, although tapping temperature is up to 1500 DEG C or more, since recycling is difficult, at present in addition to blast furnace slag is adopted
It is other to be still in the experimental study stage other than water quenching recovery waste heat water.By water quenching recovery waste heat water although converted
The thermal efficiency of journey is very high but available energy is lost seriously, and effective use efficiency only has 12%.Therefore need to seek new thinking
Method or direction.
In above-mentioned techniques or methods, continuous casting steel machine waste heat is not taken seriously and conducts a research yet, continuous casting steel machine temperature span
Larger, solidifying out into the total temperature drop of continuous casting billet up to 1600 DEG C or so from high-temperature liquid state molten steel (is about 800 if continuous casting and rolling temperature drop
DEG C), in addition to portion of energy is absorbed by thermal convection by equipment and cooling water, there is quite a few energy to scatter and disappear with radiation mode
And it is largely wasted；Therefore, only clear continuous casting steel machine process radiation heat loss vector and distribution situation could be directed to its waste heat and return
It receives and carries out further investigation using technology.
Summary of the invention
In view of the deficiencies of the prior art, the present invention proposes a kind of numerical value determination side of Con casting ladle space radiant heat flux distribution
Method, the radiations heat energy of accurate description ladle different time different spaces situation of scattering and disappearing to reduce cost to reach realize section
The purpose of energy emission reduction.
A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method, comprising the following steps:
Step 1, using operating condition ladle as object, establish 1/2 three-dimensional physical model of ladle；
Step 2 carries out sliding-model control to model using finite difference calculus, using tetrahedron and hexahedron hybrid grid pair
Model is divided；
Step 3 determines cast material physical parameter, and the side of numerical simulation feedback modifiers is combined using mathematical analysis prediction
Method determines NATURAL CONVECTION COEFFICIENT OF HEAT；
Step 4 simulates the steady state heat transfer of the solid domain of model, obtains temperature field；
Step 5, on the basis of step 4, using solid domain temperature field as primary condition, to the coupling wink in the solid domain of stream of model
State is simulated, and different moments ladle surface temperature field and radiation heat loss's vector everywhere are obtained；
Step 6, according to ladle surface, temperature field and radiation heat loss's vector carry out waste heat recycling everywhere.
Cast material physical parameter described in step 3, comprising: thermal coefficient, density, specific heat, viscosity, thermal expansion coefficient and
Radiance.
Heat transfer free convection is determined using the method that mathematical analysis prediction combines numerical simulation feedback modifiers described in step 3
Coefficient, specific as follows:
Step 3-1, initial qualitative temperature is determined；
Qualitative temperature formula is as follows:
Wherein, tmIndicate qualitative temperature；twFor wall surface temperature, the value range of initial value is 200~400 DEG C；tfFor environment
Temperature；
Step 3-2, the quasi- number of Nu Saier is obtained；
Wherein, Nu indicates the quasi- number of Nu Saier；C and n is constant, is obtained according to turbulent flow property；β indicates air thermal expansion system
Number；G is acceleration of gravity；D indicates qualitative size；Δ t indicates wall surface and circumstance of temperature difference；V indicates viscosity；Pr indicates Prandtl
Number；
Step 3-3, according to the quasi- number of the Nu Saier of acquisition, NATURAL CONVECTION COEFFICIENT OF HEAT is determined:
Wherein, ααIndicate NATURAL CONVECTION COEFFICIENT OF HEAT；λ indicates thermal coefficient；
Step 3-4, obtained NATURAL CONVECTION COEFFICIENT OF HEAT is inputted in steel ladle steady state heat transfer numerical simulation, is obtained
Steady-State Thermal Field obtains ladle surface temperature everywhere；
Step 3-5, by obtained ladle, surface temperature is returned in substitution formula (1) everywhere, to NATURAL CONVECTION COEFFICIENT OF HEAT
It is modified；
Step 3-6, judge whether adjacent NATURAL CONVECTION COEFFICIENT OF HEAT relative error twice is less than setting value, if so,
Final convection transfer rate is obtained, it is no to then follow the steps 3-4.
The steady state heat transfer of the solid domain of model is solved described in step 4, the specific steps are as follows:
Step 4-1, setting numerical value is calculated as parallel computation, and grid auto-partition is arranged；
Step 4-2, selection is based on Pressure solution device, is set as steady-state process, is not provided with gravity；
Step 4-3, energy conservation equation and S2S radiation patterns are opened；
Step 4-4, boundary condition is set, specific as follows:
Ladle side wall, packet bottom, packet eaves, headroom liner and slag blanket upper surface use third boundary condition；
Environment temperature, external radiation rate and convection transfer rate are set；
At stream liquid/solid interface, the wall surface of setting solid area side is constant temperature wall, and fluid side is set as adiabatic wall；
Step 4-5, setting solver is pressure-velocity couple solution device, scheme SIMPLE, momentum, Turbulent Kinetic, rapids
It flows dissipative shock wave and energy chooses Second-order Up-wind operation, under-relaxation factor uses system default value；
Step 4-6, Initialize installation is carried out, molten steel temperature is set；
Step 4-7, it is iterated until convergence, obtains temperature field after stopping.
Described in step 5 on the basis of step 4, the coupling transient state in the solid domain of the stream of model is simulated, specific steps
It is as follows:
Step 5-1, steady-state process is revised as transient process, and opens gravity；
Step 5-2, standard k-ε turbulence model, and selective enhancement casing treatment and full buoyancy model are increased；
Step 5-3, it modifies to boundary condition, i.e., sets coupling wall surface for stream liquid/solid interface；
Step 5-4, fluid domain is arranged using patch as solid domain primary condition in the temperature field for obtaining steady state heat transfer
Primary condition, and start iteration；
Step 5-5, time step is set as adaptivity step-length, is revised as fixed step size after residual error is stablized, obtains not
Ladle surface temperature field and radiation heat loss's vector everywhere in the same time.
The invention has the advantages that
The present invention proposes that a kind of numerical value of Con casting ladle space radiant heat flux distribution determines that method, this method will rely on continuous
Technology of Steel Castings is practical, obtains heat source temperature and heat flux distribution using research methods such as mathematical computations and numerical simulations, and be based on this
Continuous casting steel machine process radiations heat energy model is developed, both provided new way for continuous casting steel machine heat recovery, is also exploitation steel
It manufactures whole process energy conservation and environmental protection new technology and effective " energy-saving and emission-reduction " implements based theoretical and oil storage method method；This is not
Enterprise's manufacturing cost will only be greatly reduced, also can greatly alleviate industrial energy saving emission reduction, and then generate very big social economy
Benefit, correlative study have important theory and realistic meaning.
Detailed description of the invention
Fig. 1 is that the numerical value that the Con casting ladle space radiant heat flux of one embodiment of the present invention is distributed determines method flow
Figure；
Fig. 2 is the 1/2 three-dimensional physical model schematic diagram of ladle of one embodiment of the present invention；
Fig. 3 is that the NATURAL CONVECTION COEFFICIENT OF HEAT of one embodiment of the present invention determines flow chart；
Fig. 4 is that certain moment of one embodiment of the present invention always radiates and wherein the heat flow density of heat loss through radiation part is illustrated
Figure；
Fig. 5 is that certain moment each position radiations heat energy of one embodiment of the present invention accounts for the ratio signal of global radiation heat
Figure.
Specific embodiment
An embodiment of the present invention is described further with reference to the accompanying drawing.
In the embodiment of the present invention, the numerical value of Con casting ladle space radiant heat flux distribution determines method, method flow diagram such as Fig. 1
It is shown, comprising the following steps:
Step 1, using operating condition ladle as object, establish 1/2 three-dimensional physical model of ladle；
In the embodiment of the present invention, 250 tons of Con casting ladles of certain factory are chosen as research object, remove the hoisting mechanism of ladle,
The influence of mainstream control mechanism and bottom gas-feeding device, ladle wall construction and thickness of slag layer use site measured data are not considered,
CAD modeling software is used to establish 1/2 three-dimensional physical model of ladle as shown in Fig. 2, 1 is molten steel；2 be working lining；3 be permanent layer；4
For heat insulation layer；5 be steel shell；6 be slag blanket；Ladle is tapered back taper bench-type in figure, to show that internal structure carries out portion to it
Divide subdivision；Ladle, which show ladle, expires bag-like state, and interior zone is molten steel, and top is slag blanket；Ladle side wall and packet bottom are all four
Layer structure, from inside to outside respectively working lining, permanent layer, heat insulation layer, steel shell；
Table 1 is ladle key dimension；
2 ladle dimensional parameters of table
Step 2 carries out sliding-model control to model using finite difference calculus, using tetrahedron and hexahedron hybrid grid pair
Model is divided；
In the embodiment of the present invention, the model of building is imported into grid dividing software and carries out grid dividing；Grid is using four sides
Body and hexahedral hybrid grid increase boundary layer, integral grid size in molten steel and working lining, slag blanket stream liquid/solid interface interface
Maximum is no more than 0.05m, grid number 2056224；
Step 3 determines cast material physical parameter, and the side of numerical simulation feedback modifiers is combined using mathematical analysis prediction
Method determines NATURAL CONVECTION COEFFICIENT OF HEAT；
In the embodiment of the present invention, the air thermal physical property parameter of corresponding temperature is checked in by " Metallurgy Transport Principle " annex two, is wrapped
It includes: thermal coefficient, density, specific heat, viscosity, thermal expansion coefficient and radiance；
In the embodiment of the present invention, ready-portioned grid is imported into limited bulk and is solved in software, and checks that mesh quality is protected
Finite elements are demonstrate,proved without negative volume, is then arranged and material property is arranged according to table 1, wherein ladle structure and slag blanket are set as solid, steel
Liquid is set as liquid, and uses Boussinesq it is assumed that initial reference temperature is 1873K；
It is 30 DEG C that ladle ambient air temperature, which is arranged, in environment temperature, and side wall and packet bottom are that steel plate setting blackness is 0.8, slag blanket
Blackness is 0.9, and headroom inner lining surface blackness is 0.9；Temperature when reaching stable state due to each surface can not determine that convection current is changed
The hot bad determination of coefficient, the method that numerical simulation feedback modifiers are combined using mathematical analysis prediction described in the embodiment of the present invention
Determine NATURAL CONVECTION COEFFICIENT OF HEAT, method flow diagram is as shown in figure 3, specific as follows:
Step 3-1, initial qualitative temperature is determined；
At 300-400 DEG C, ladle surface and air contact surface temperature when general ladle working naturally, are nature pair
Stream；For slag blanket, the free convection of the face-up horizontal wall surface of heat can be considered as to heat；In the embodiment of the present invention, if Air Temperature
Degree is 30 DEG C, and the temperature of slag blanket and packet bottom outer surface is about 300 DEG C, then qualitative temperature are as follows:
Wherein, tmIndicate qualitative temperature；twFor wall surface temperature, tfFor environment temperature；
Step 3-2, the quasi- number of Nu Saier is obtained；
Wherein, Nu indicates the quasi- number of Nu Saier；C and n is constant, as shown in table 2, is obtained according to turbulent flow property；β indicates empty
Gas thermal expansion coefficient；G is acceleration of gravity；D indicates qualitative size, in the embodiment of the present invention, d=4m；Δ t indicates wall surface and ring
The border temperature difference；V indicates viscosity, v=30.1 × 10 in the embodiment of the present invention-6m2s-1；Pr indicates Prandtl number, the embodiment of the present invention
In, Pr=0.685；
Constant C and n value in 2 formula of table (2)
Step 3-3, according to the quasi- number of the Nu Saier of acquisition, NATURAL CONVECTION COEFFICIENT OF HEAT is determined:
Wherein, ααIndicate NATURAL CONVECTION COEFFICIENT OF HEAT；λ indicates thermal coefficient, λ=3.65 × 10 in the embodiment of the present invention- 2W·m-1·℃-1；
Step 3-4, obtained NATURAL CONVECTION COEFFICIENT OF HEAT is inputted in steel ladle steady state heat transfer numerical simulation, is obtained
Steady-State Thermal Field obtains ladle surface temperature everywhere；
In the embodiment of the present invention, steel ladle steady state heat transfer numerical simulation includes: selection solver and model specification material
Expect parameter, setting boundary condition (i.e. input NATURAL CONVECTION COEFFICIENT OF HEAT) is arranged parameter and initializes, and stable state iterative calculation obtains
Obtain surface temperature；
Step 3-5, by obtained ladle, surface temperature is returned in substitution formula (1) everywhere, to NATURAL CONVECTION COEFFICIENT OF HEAT
It is modified；
Step 3-6, adjacent NATURAL CONVECTION COEFFICIENT OF HEAT relative error twice is judged whether less than 0.01, if so, obtaining
Convection transfer rate finally is obtained, it is no to then follow the steps 3-4.
In the embodiment of the present invention, gained convection transfer rate Input Software will be calculated and carry out the calculating of ladle steady state heat transfer, asked
Slag blanket surface actual temperature is obtained, is then brought into mathematical analysis process above using simulation gained temperature and is calculated again respectively
It solves, obtaining convection transfer rate is 10.24Wm-2·K-1, final result as correction after convection transfer rate be applied to
In the boundary condition setting of the present embodiment；
The accuracy of primary condition setting directly affects the convergence and accuracy of Transient Numerical Simulation；Operating condition ladle passes through
After refining, the temperature of steel ladle is non-uniform Distribution, and especially in thickness direction, there are the temperature gradients of momentary stabilization；The present invention
In embodiment, for obtain with operating condition actually closest containment wall temperature field as primary condition, first only to solid area into
The iterative calculation of row stable state, carries out instantaneous iteration for calculated result as the primary condition of fluid-solid conjugated heat transfer again；
The determination of boundary condition is also particularly important；The boundary setting of the embodiment of the present invention is divided to two processes, first to solid
Body region carry out steady state heat transfer when boundary condition setting, second is that when fluid and structural simulation boundary condition setting；
Step 4 simulates the steady state heat transfer of the solid domain of model, obtains temperature field；
In the embodiment of the present invention, for solid stable state diabatic process, the side wall of ladle, packet bottom, packet eaves, headroom liner,
Convection current and radiation effects occur for slag blanket upper surface and ambient enviroment, therefore use third boundary condition, setting environment temperature, outer
Portion's radiance and convection transfer rate.At stream liquid/solid interface, the wall surface of setting solid area side is First Boundary Condition,
Temperature is the constant temperature wall of 1873K, and fluid side is set as adiabatic wall；
Specific step is as follows:
Step 4-1, setting numerical value is calculated as parallel computation, and grid auto-partition is arranged；
In the embodiment of the present invention, since lattice number is huge, it is two using server that numerical value, which calculates, which uses parallel computation,
6 cores, 3.74GHz processor import grid and to its subregion；
Step 4-2, selection is based on Pressure solution device, is set as steady-state process, is not provided with gravity；
Step 4-3, energy conservation equation and S2S radiation patterns are opened；
Step 4-4, boundary condition is set, specific as follows:
Ladle side wall, packet bottom, packet eaves, headroom liner and slag blanket upper surface use third boundary condition；
Environment temperature, external radiation rate and convection transfer rate are set；
At stream liquid/solid interface, the wall surface of setting solid area side is constant temperature wall, and fluid side is set as adiabatic wall；
Step 4-5, solver selects pressure-velocity couple solution device, scheme SIMPLE, and pressure difference value format is selected
PRESTO！, momentum, Turbulent Kinetic, turbulence dissipation rate and energy selection Second-order Up-wind operation, under-relaxation factor use system default
Value；
Step 4-6, Initialize installation is carried out, molten steel temperature is set；In the embodiment of the present invention, setting molten steel temperature is
1873K；
Step 4-7, it is iterated until convergence, obtains temperature field after stopping.
In the embodiment of the present invention, iteration number is set as 20, remaining parameter uses default setting, starts to carry out stable state meter
It calculates, calculates convergence after 10 step of iteration and stop；Open post-processing observation temperature field effect；
Step 5, on the basis of step 4, using solid domain temperature field as primary condition, to the coupling wink in the solid domain of stream of model
State is simulated, and different moments ladle surface temperature field and radiation heat loss's vector everywhere are obtained；
After the completion of solid stable state Calculation of Heat Transfer, before fluid and structural simulation, need to reset boundary condition, i.e., it will stream
Liquid/solid interface is set as Coupled wall, so that stream be coupled admittedly；
Specific step is as follows:
Step 5-1, steady-state process is revised as transient process, and opens gravity；
Step 5-2, standard k-ε turbulence model, and selective enhancement casing treatment and full buoyancy model are increased；
Step 5-3, it modifies to boundary condition, i.e., sets Coupled wall for stream liquid/solid interface；
Step 5-4, fluid domain is arranged using patch as solid domain primary condition in the temperature field for obtaining steady state heat transfer
Primary condition, and start iteration；
In the embodiment of the present invention, solver setting is constant, selects pressure-velocity couple solution device, scheme SIMPLE, pressure
Power difference format selects PRESTO！, momentum, Turbulent Kinetic, turbulence dissipation rate and energy choose Second-order Up-wind operation, Asia relaxation because
Son uses system default value；
Step 5-5, time step is set as adaptivity step-length, is revised as fixed step size after residual error is stablized, obtains not
Ladle surface temperature field and radiation heat loss's vector everywhere in the same time；
In the embodiment of the present invention, time started step-length is set as adaptivity step-length, it is to be calculated after a certain period of time, observe residual
Step-length is changed to fixed step size, thus in the base for not influencing accuracy if residual error is stablized in a smaller range by poor curvilinear trend
Accelerate calculating speed on plinth to complete until calculating.
Step 6, according to ladle surface, temperature field and radiation heat loss's vector carry out waste heat recycling everywhere.
In the embodiment of the present invention, ladle bulk temperature field and flow field are obtained with the poster processing soft and to surface heat flow data
It is for statistical analysis as follows:
Fig. 4 is the heat flow density always to radiate with wherein heat loss through radiation part at certain moment chosen the standing phase, is seen by comparison
Out in the most areas of ladle, specific gravity of the heat loss through radiation in total heat dissipation is all very big, i.e., heat loss through radiation is ladle heat damage
The principal mode of mistake；Consider in rate of heat dispation direction, the heat flow density of slag blanket is maximum and is much larger than other positions.
Fig. 5 is the ratio that certain the moment each position radiations heat energy chosen the standing phase accounts for global radiation heat, side wall unit time
Radiations heat energy the largest loss, account for more than half of global radiation heat, slag blanket takes second place, about 31%.
As long as thus obtaining correlative study of the ladle thermal loss based on radiation heat loss, for ladle radiation heat loss
It is of great significance；The radiant heat of slag blanket and side wall by be waste heat recycling main object.
Claims (5)
1. a kind of numerical value of Con casting ladle space radiant heat flux distribution determines method, which comprises the following steps:
Step 1, using operating condition ladle as object, establish 1/2 three-dimensional physical model of ladle；
Step 2 carries out sliding-model control to model using finite difference calculus, using tetrahedron and hexahedron hybrid grid to model
It is divided；
Step 3 determines cast material physical parameter, and combines the method for numerical simulation feedback modifiers true using mathematical analysis prediction
Determine NATURAL CONVECTION COEFFICIENT OF HEAT；
Step 4 simulates the steady state heat transfer of the solid domain of model, obtains temperature field；
Step 5, on the basis of step 4, using solid domain temperature field as primary condition, to the coupling transient state in the solid domain of stream of model into
Row simulation obtains different moments ladle surface temperature field and radiation heat loss's vector everywhere；
Step 6, according to ladle surface, temperature field and radiation heat loss's vector carry out waste heat recycling everywhere.
2. the numerical value of Con casting ladle space according to claim 1 radiant heat flux distribution determines method, which is characterized in that step
Cast material physical parameter described in rapid 3, comprising: thermal coefficient, density, specific heat, viscosity, thermal expansion coefficient and radiance.
3. the numerical value of Con casting ladle space according to claim 1 radiant heat flux distribution determines method, which is characterized in that step
NATURAL CONVECTION COEFFICIENT OF HEAT is determined using the method that mathematical analysis prediction combines numerical simulation feedback modifiers described in rapid 3, specifically
It is as follows:
Step 3-1, initial qualitative temperature is determined；
Qualitative temperature formula is as follows:
Wherein, tmIndicate qualitative temperature；twFor wall surface temperature, the value range of initial value is 200~400 DEG C；tfFor environment temperature
Degree；
Step 3-2, the quasi- number of Nu Saier is obtained；
Wherein, Nu indicates the quasi- number of Nu Saier；C and n is constant, is obtained according to turbulent flow property；β indicates air thermal expansion coefficient；g
For acceleration of gravity；D indicates qualitative size；Δ t indicates wall surface and circumstance of temperature difference；V indicates viscosity；Pr indicates Prandtl number；
Step 3-3, according to the quasi- number of the Nu Saier of acquisition, NATURAL CONVECTION COEFFICIENT OF HEAT is determined:
Wherein, αaIndicate NATURAL CONVECTION COEFFICIENT OF HEAT；λ indicates thermal coefficient；
Step 3-4, obtained NATURAL CONVECTION COEFFICIENT OF HEAT is inputted in steel ladle steady state heat transfer numerical simulation, obtains stable state
Temperature field obtains ladle surface temperature everywhere；
Step 3-5, by obtained ladle, surface temperature is returned in substitution formula (1) everywhere, is carried out to NATURAL CONVECTION COEFFICIENT OF HEAT
Amendment；
Step 3-6, judge whether adjacent NATURAL CONVECTION COEFFICIENT OF HEAT relative error twice is less than setting value, if so, obtaining
Final convection transfer rate, it is no to then follow the steps 3-4.
4. the numerical value of Con casting ladle space according to claim 1 radiant heat flux distribution determines method, which is characterized in that step
The steady state heat transfer of the solid domain of model is solved described in rapid 4, the specific steps are as follows:
Step 4-1, setting numerical value is calculated as parallel computation, and grid auto-partition is arranged；
Step 4-2, selection is based on Pressure solution device, is set as steady-state process, is not provided with gravity；
Step 4-3, energy conservation equation and S2S radiation patterns are opened；
Step 4-4, boundary condition is set, specific as follows:
Ladle side wall, packet bottom, packet eaves, headroom liner and slag blanket upper surface use third boundary condition；
Environment temperature, external radiation rate and convection transfer rate are set；
At stream liquid/solid interface, the wall surface of setting solid area side is constant temperature wall, and fluid side is set as adiabatic wall；
Step 4-5, setting solver is pressure-velocity couple solution device, scheme SIMPLE, momentum, Turbulent Kinetic, turbulent flow consumption
The rate of dissipating and energy choose Second-order Up-wind operation, and under-relaxation factor uses system default value；
Step 4-6, Initialize installation is carried out, molten steel temperature is set；
Step 4-7, it is iterated until convergence, obtains temperature field after stopping.
5. the numerical value of Con casting ladle space according to claim 1 radiant heat flux distribution determines method, which is characterized in that step
Described in rapid 5 on the basis of step 4, the coupling transient state in the solid domain of the stream of model is simulated, the specific steps are as follows:
Step 5-1, steady-state process is revised as transient process, and opens gravity；
Step 5-2, standard k-ε turbulence model, and selective enhancement casing treatment and full buoyancy model are increased；
Step 5-3, it modifies to boundary condition, i.e., sets coupling wall surface for stream liquid/solid interface；
Step 5-4, the temperature field for obtaining steady state heat transfer is initial using patch setting fluid domain as solid domain primary condition
Condition, and start iteration；
Step 5-5, time step is set as adaptivity step-length, fixed step size is revised as after residual error is stablized, when obtaining different
Carve ladle surface temperature field and radiation heat loss's vector everywhere.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610561555.2A CN106202757B (en) | 2016-07-15 | 2016-07-15 | A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610561555.2A CN106202757B (en) | 2016-07-15 | 2016-07-15 | A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106202757A CN106202757A (en) | 2016-12-07 |
| CN106202757B true CN106202757B (en) | 2019-07-23 |
Family
ID=57475832
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610561555.2A Active CN106202757B (en) | 2016-07-15 | 2016-07-15 | A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106202757B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107084805B (en) * | 2017-04-06 | 2019-08-23 | 东北大学 | A kind of method of determining Con casting ladle radiant heat flux |
| CN107301276B (en) * | 2017-06-01 | 2020-08-04 | 上海理工大学 | Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner |
| CN108629145A (en) * | 2018-07-12 | 2018-10-09 | 朱林 | The determination method of damage is cut based on the pendulum-type plough plow breast chisel of heat flow piercement behavior |
| CN113128137A (en) * | 2021-04-06 | 2021-07-16 | 中国辐射防护研究院 | Design method and device of radioactive contamination ventilation protective clothing |
| CN113343428A (en) * | 2021-05-12 | 2021-09-03 | 江苏永钢集团有限公司 | Method for monitoring baking state of steel ladle in real time |
| CN113673132A (en) * | 2021-08-27 | 2021-11-19 | 东北大学 | Concave thermoelectric experimental device and convection radiation numerical simulation method thereof |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103382515A (en) * | 2013-07-19 | 2013-11-06 | 东北大学 | System and method for monitoring molten steel temperature in RH refining process in online real-time manner |
-
2016
- 2016-07-15 CN CN201610561555.2A patent/CN106202757B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103382515A (en) * | 2013-07-19 | 2013-11-06 | 东北大学 | System and method for monitoring molten steel temperature in RH refining process in online real-time manner |
Non-Patent Citations (4)
| Title |
|---|
| 100t钢包热循环传热模型的研究;李青;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20130315(第3期);第B023-97页 |
| 130T钢包包壳的温度场和应力场的有限元分析;王倩;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20150615(第6期);第B023-24页 |
| 钢包加盖过程仿真及其保温性能研究;蒋波 等;《冶金动力》;20140415(第4期);第61-65页 |
| 钢包稳态温度场的有限元模拟;杨治立 等;《特殊钢》;20070501;第28卷(第3期);第41-43页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202757A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202757B (en) | A kind of numerical value of Con casting ladle space radiant heat flux distribution determines method | |
| Feng et al. | Constructal optimization of a sinter cooling process based on exergy output maximization | |
| CN109902433B (en) | Cross-dimension coupling method for pressurized water reactor passive containment waste heat discharge system | |
| CN106957935B (en) | A kind of flexible measurement method of blast furnace inside soft heat belt shape | |
| CN112464543B (en) | Method for calculating movement of inclusions in VIM smelting process | |
| Tian et al. | Global sensitivity analysis and multi-objective optimization design of temperature field of sinter cooler based on energy value | |
| Shen et al. | Numerical study on the heat recovery and cooling effect by built-in pipes in a graphitization furnace | |
| CN106197762B (en) | A kind of method of determining Con casting ladle radiant heat flux | |
| Oumarou et al. | Transient process model of open anode baking furnace | |
| Zhang et al. | Parameter optimization of gas–solid heat transfer process in sinter packed bed based on further exergy analysis | |
| Lai et al. | Experimental study on storage performance of packed bed solar thermal energy storage system using sintered ore particles | |
| Shao et al. | Numerical simulation of heat transfer process in cement grate cooler based on dynamic mesh technique | |
| CN113111549B (en) | Erosion model modeling method and modeling system for casting repaired blast furnace hearth | |
| CN107084805A (en) | A kind of method for determining Con casting ladle radiant heat flux | |
| Oumarou et al. | A dynamic process model for predicting the performance of horizontal anode baking furnaces | |
| Shen et al. | Numerical simulation and analyses for sinter cooling process with convective and radiative heat transfer | |
| Cao et al. | Smoothed particle hydrodynamics modeling and simulation of foundry filling process | |
| Zhang et al. | The study of heat-mass transfer characteristics and optimization on electric arc furnace with the novel fence electrode | |
| Liu et al. | Asymptotic Model of Refractory and Buildup State of the Blast Furnace Hearth | |
| ABHALE et al. | COMPREHENSIVE MODEL FOR BLAST FURNACE USING OPENFOAM® | |
| Zhang et al. | Numerical investigation of gas-solid heat transfer process and parameter optimization in shaft kiln for high-purity magnesia | |
| Wang et al. | Numerical simulation of deformation during hot procedure for large hydraulic turbine runner blade | |
| Gao et al. | Estimation of Oxygen Concentration for Roller Kiln Based on Boundary Measurements | |
| Xu et al. | Structure design of aluminum alloy solution furnace based on fluid-solid coupling numerical simulation | |
| Luo et al. | Simulation analysis of temperature field and its influence factors of the new structure ladle |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |