CN109530648A - A kind of method of crystallizer slag interface fluctuation in prediction continuous casting - Google Patents
A kind of method of crystallizer slag interface fluctuation in prediction continuous casting Download PDFInfo
- Publication number
- CN109530648A CN109530648A CN201910078981.4A CN201910078981A CN109530648A CN 109530648 A CN109530648 A CN 109530648A CN 201910078981 A CN201910078981 A CN 201910078981A CN 109530648 A CN109530648 A CN 109530648A
- Authority
- CN
- China
- Prior art keywords
- crystallizer
- model
- boundary condition
- slag interface
- continuous casting
- 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.)
- Pending
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B22—CASTING; POWDER METALLURGY
- B22D—CASTING OF METALS; CASTING OF OTHER SUBSTANCES BY THE SAME PROCESSES OR DEVICES
- B22D11/00—Continuous casting of metals, i.e. casting in indefinite lengths
- B22D11/16—Controlling or regulating processes or operations
- B22D11/18—Controlling or regulating processes or operations for pouring
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Continuous Casting (AREA)
Abstract
The present invention provides a kind of method for predicting crystallizer slag interface fluctuation in continuous casting, is related to metal casting technique field.This method is first according to steel mill's continuous casting installation for casting process, crystallizer 3-D geometric model is established with modeling software, gridding is carried out to model, formulates model boundary condition, export .msh file, and it is conducted into fluid calculation software and carries out simulation calculating, extract the data in crystallizer slag interface position, shape and neighbouring flow field, and in the poster processing soft, the liquid fluctuating height rectangular coordinate system scattergram of slab median plane radial distance corresponding thereto is established, judges that slag interface fluctuates situation accordingly；Finally adjustment slab mouth of a river insertion depth and mouth of a river shape again predict slag interface fluctuation situation.The method that crystallizer slag interface fluctuates in prediction continuous casting provided by the invention, it can be determined that slab Flow Field Distribution and slag interface in continuous casting process fluctuate behavior, to improve slab quality.
Description
Technical field
The present invention relates to crystallizer slag interface waves in metal casting technique field more particularly to a kind of prediction continuous casting
Dynamic method.
Background technique
It is higher and higher to the performance requirement of steel products with the increasingly fierceness of rapid development of economy and market competition,
Therefore, to continuous casting production process, more stringent requirements are proposed.The factors such as cleanliness, surface defect and the internal flaw of slab are to production
The quality problems of product serve vital, and crystallizer operation exactly controls the last ring that these factors reach required standard
Section.Therefore, while guaranteeing conticaster production efficiency, the slab of high-quality is obtained, has become the pass for restricting high-quality steel
Key step.
In continuous casting production process, there is the quality of continuous casting billet highly important steel slag liquid fluctuating situation in crystallizer
It influences.The flow behavior of Mold not only influences the heat transfer of Mold and the floating of field trash, but also with
The surfaces such as slab crackle, segregation and internal soundness have close relationship, and reasonable crystallizer molten steel flow field is rolled up to molten steel is reduced
The defects of slag, bubble, crackle, improves slab quality, there is good facilitation.Therefore study and control the steel in crystallizer
Raising steel product quality is significant in turn for improving cc billet surface quality for slag interface fluctuation.
If slag interface fluctuation is excessive in crystallizer, covering slag is easy to cause to be involved in molten steel, causes dew steel phenomenon, and
And casting billet surface can be brought to be mingled with and subsurface inclusion, seriously affect slab quality.If slag interface fluctuation is excessively in crystallizer
Calmness then will lead to molten steel at meniscus and update slowly, the phenomenon that being unfavorable for the fusing of covering slag, lead to insufficient lubrication.
Summary of the invention
The technical problem to be solved by the present invention is in view of the above shortcomings of the prior art, provide in a kind of prediction continuous casting
The method of crystallizer slag interface fluctuation improves Flow Field Distribution and slag interface fluctuation in continuous casting process, to improve slab matter
Amount.
In order to solve the above technical problems, the technical solution used in the present invention is: a kind of crystallizer in prediction continuous casting
The method of slag interface fluctuation, comprising the following steps:
Step 1, according to steel mill's continuous casting installation for casting process, establish crystallizer three-dimensional geometry mould with Gambit modeling software
Type, and gridding is carried out to model, formulates model boundary condition, then export .msh file；The boundary condition includes entrance
Boundary condition, export boundary condition, model top surface boundary condition and crystallizer wall surface and mouth of a river wall boundary condition；
Step 2 will carry out simulation calculating in the importing fluid calculation software of .msh file derived from step 1, method particularly includes:
Step 2.1, grid import and processing: the .msh file saved in read step 1 checks grid, guarantees grid most
Net region size is arranged less than 0 in small size；
Step 2.2, setting physical model and material: negative for Y-axis according to the type of actual setting solver, acceleration of gravity
Direction, select turbulence model and multiphase flow model for computation model, setting material be fluid, be arranged slag interface surface tension；
Step 2.3, setting operating environment and boundary condition；
The setting boundary condition method particularly includes:
The entrance boundary condition includes inlet velocity vinlet, tubulence energy k and tubulence energy dissipative shock wave ε, respectively according to formula
(1)-(3) it obtains:
Wherein, vcastFor drawing speed, SoutFor mold exit area of section, SinFor submersed nozzle inlet section face
Product；
Wherein, RnozFor nozzle inlet radius；
The model top surface is set as free surface, and the velocity component and every other each variable perpendicular to liquid level are along level surface method
The gradient in line direction is set as zero, meets formula (4):
Wherein, v is the velocity magnitude of model Free Surface x and the direction z, h=CpTinFor molten steel sensible heat, CpFor specific heat, Tin
For inlet temperature, u is throwing direction speed；
Using the interface between VOF model following phase for multiphase flow model is the volume by solving a phase or multiphase
The continuity equation of ratio is completed；To q phase, the multiphase flow conservation of mass is controlled by following formula:
Wherein, SaqFor quality source item, αqFor q phase volume fraction, ρqFor q phase density, t is the time；In default scenario
Lower formula (5) right end source item is zero, but when to each phase specified constant or user-defined Mass Sources, then right end is not zero；
The export boundary condition needs the condition met are as follows: sufficiently development, i.e., each physical quantity edge are flowed at mold exit
The normal derivative in the section is zero；
For the crystallizer wall boundary condition at crystallizer wall surface, the velocity component perpendicular to wall surface is zero, and parallel
In wall surface component using no slip boundary condition, meet formula (6):
Voutlet=Vcast, Vwall=kwall=0 (6)
Wherein, Voutlet is muzzle velocity, and Vcast is pulling rate, and Vwall is crystallizer wall face velocity, and kwall is crystallization
Wall face tubulence energy；
Step 2.4, setting method for solving and control parameter；
Step 2.5, setting monitor window and initialization；When initialization, set in range at top liquid level 15mm
For steel slag, other regions are defaulted as molten steel；
Step 2.6, setting iterative steps and time step, carry out parametric solution；
Shown in the restrictive condition of the time step such as formula (7):
Courant=u Δ t/ Δ x (7)
Wherein, u is fluid velocity, and Δ t is time step, and Δ x is size of mesh opening, and Courant is Krona spy number, is indicated
A fluid particle can pass through how many a grids in a time step；As Krona spy number < 1, calculates and stablize, but when
Between step-length it is small, calculate overlong time；As 1 < Kronas of spy number < 5, computational stability still seldom occurs calculating diverging very well；
When Krona spy number > 10, calculating is interrupted it is easy to appear diverging；Therefore, suitable time step is selected in Krona spy number assistance；
The parameter for finally solving acquisition includes the entrance velocity vinlet, tubulence energy k, tubulence energy dissipative shock wave ε and time step
Long Δ t；
Step 3, simulation extract data after calculating；In fluid calculation software extract crystallizer slag interface position,
Shape and the data in neighbouring flow field, and in the poster processing soft Tecplot, establish slab median plane radial distance corresponding thereto
The liquid fluctuating height rectangular coordinate system scattergram answered；
Step 4 analyzes post-processing result: the rectangular coordinate system scattergram obtained according to step 3 judges steel
Slag interface fluctuation situation；
Step 5, adjustment slab mouth of a river insertion depth and mouth of a river shape re-execute the steps 1-4 and fluctuate feelings to slag interface
Condition is predicted, and is compared and analyzed to the slag interface fluctuation behavior of adjustment front and back.
The beneficial effects of adopting the technical scheme are that being tied in a kind of prediction continuous casting provided by the invention
The method of brilliant device slag interface fluctuation, combined by modeling software Gambit and fluid calculation software, is being built to molten steel
Flowing is simulated in 3-D geometric model, extracted after having been calculated crystallizer slag interface position, shape and neighbouring flow field into
Row analysis can predict that crystallizer slag interface fluctuates situation, judge whether phenomena such as will appear dew coil of strip slag.Through the invention
The slag interface that prediction technique obtains fluctuates situation, it can be determined that slab Flow Field Distribution and slag interface fluctuation in continuous casting process
Behavior, to improve slab quality.Meanwhile being conducive in continuous casting process apart from each spot speed ladder in meniscus sustained height face center
Degree is in a certain range of permission；Velocity gradient is not much different available more stable molten steel face, can preferably prevent
The generation of slag inhibits the generation of the subcutaneous negative segregation of slab, with the subcutaneous negative segregation band of coordinated regulation slab and center segregation, realizes
The raising of slab homogenieity；Be conducive to slag interface undulated control in the reasonable scope, be more advantageous to slag blanket temperature and uniformly rise
Height melts slag blanket.
Detailed description of the invention
Fig. 1 is the stream of the method for crystallizer slag interface fluctuation in a kind of prediction continuous casting provided in an embodiment of the present invention
Cheng Tu；
Fig. 2 is the initial 3-D geometric model schematic diagram of crystallizer provided in an embodiment of the present invention, wherein (a) is whole mould
Type (b) is initial mouth of a river partial enlarged view；
Fig. 3 is the schematic diagram of the slag interface distribution and Flow Field Distribution situation of prediction provided in an embodiment of the present invention, wherein
(a) it is the slag interface distribution and Flow Field Distribution situation predicted using initial 3-D geometric model, (b) is to be inserted using the change mouth of a river
The slag interface distribution and Flow Field Distribution situation of 3-D geometric model prediction after entering depth；
Fig. 4 is the schematic diagram of the crystallizer 3-D geometric model after change mouth of a river shape provided in an embodiment of the present invention；
Fig. 5 is the signal of slag interface distribution and Flow Field Distribution situation after change mouth of a river shape provided in an embodiment of the present invention
Figure；
Fig. 6 is the liquid of slab median plane radial distance corresponding thereto after change mouth of a river shape provided in an embodiment of the present invention
Surface wave dynamic height rectangular coordinate system figure.
In figure, 1, slab eddy flow downspout；2, slag interface in crystallizer；3, slab eddy flow downspout exports；4, slab；5, electric
Magnetic mixing plant；6, continuous cast mold；
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
The present embodiment is by taking the original round billet design drawing of certain steel mill as an example, using crystallizer steel in prediction continuous casting of the invention
The method of slag interface fluctuation predicts the fluctuation situation of crystallizer slag interface in continuous casting.
The corresponding radius of the original round billet design drawing of the steel mill and height are respectively 250mm, 2000mm, and wherein crystallizer section has
Height 780mm is imitated, is secondary cooling zone remaining as 1220mm.The vertical disc in the mouth of a river, immersion depth 120mm, mouth of a river internal diameter 42mm, outside
Three outlets of diameter 100mm, the mouth of a river are evenly distributed on away from mouth of a river bottom 10mm, 15 ° of inclination that Open Side Down.Its production is radius
250mm round billet, nozzle inlet internal diameter are 42mm, and crystallizer effective height is 0.78m, and the mouth of a river immerses mold liquid level depth and is
0.12m, molten steel entrance velocity is 0.32m/min, if the following table 1 is the physical parameter that practical institute produces steel grade Q235 steel.
1 steel grade parameter of table
A kind of method of crystallizer slag interface fluctuation in prediction continuous casting, as shown in Figure 1, comprising the following steps:
Step 1, according to steel mill's continuous casting installation for casting process, establish crystallizer three-dimensional geometry mould with Gambit modeling software
Type, as shown in Fig. 2, and to model carry out gridding, definition formulate model boundary condition, then export .msh file；The side
Boundary's condition includes entrance boundary condition, export boundary condition, model top surface boundary condition and crystallizer wall surface and mouth of a river wall surface side
Boundary's condition；
In the present embodiment, the specific size parameter of crystallizer 3-D geometric model is as shown in table 2:
The parameter value of 2 3-D geometric model of table
| Parameter | Numerical value |
| Casting blank cross-section size, mm | R250mm |
| Model computational length, mm | 2000mm |
| Mouth of a river immersion depth, mm | 120mm |
| Crystallizer effective height, mm | 780mm |
| Distance of the water-gap outlets away from mouth of a river bottom, mm | 10mm |
| The mouth of a river is angled downward, ° | 15° |
Specifically, Gambit modeling software is opened, creating geometrical model according to above-mentioned parameter, (Y negative direction is gravity side
To), Mesh button is clicked in control panel, grid dividing is carried out to 3-D geometric model, then carries out the definition of boundary condition,
The boundary condition of definition has nozzle inlet boundary condition inlet, free surface boundary condition free surface at the top of model, mould
Type wall boundary condition moldwall and export boundary condition outlet；Modeling finishes, and three-dimensional mathematical model is saved as
Mold.msh file；
.msh file derived from step 1 is imported in fluid analysis simulation softward Ansys Fluent and is simulated by step 2
It calculates, method particularly includes:
Step 2.1, grid import and processing: the .msh file saved in read step 1 checks grid, guarantees grid most
Small size minimum volume less than 0 (General → Check), be arranged net region size (General → Scale →
Mesh Length Unit In→mm)；
Step 2.2, setting physical model and material: according to type Pressure Based of actual setting solver, again
Power acceleration is Y-axis negative direction, size 9.81m/s2, select turbulence model and multiphase flow model for computation model (Models
→Viscous Laminar→K-e/Energy；Model → Multiphase → Volume of Fluid), setting material be
Fluid (Materials → steel；Materials → Slag), the surface tension that slag interface is arranged is 1.4N/m；
Step 2.3, setting operating environment (Define → operating conditions keeps default) and boundary condition
(Boundary conditions→inlet/outlet/free surface/mold wall)；
In boundary condition, the entrance boundary condition includes inlet velocity vinlet, tubulence energy k, tubulence energy dissipative shock wave ε,
It is obtained respectively according to formula (1)-(3):
Wherein, vcastFor drawing speed, unit m/min；SoutFor mold exit area of section, unit m2；SinFor leaching
Enter area of section at formula nozzle inlet, unit m2, the vinletUnit is m/min；
Wherein, RnozFor nozzle inlet radius, unit m；
The model top surface is set as free surface, and the velocity component and every other each variable perpendicular to liquid level are along level surface method
The gradient in line direction is set as zero, meets formula (4):
Wherein, v is the velocity magnitude of model Free Surface x and the direction z, unit m/min；H is molten steel sensible heat, and unit is
J/kg, h=CpTin, CpFor specific heat, TinFor inlet temperature, unit J/kg/k；U is throwing direction speed, unit m/min；
Using the interface between VOF model following phase for multiphase flow model is the volume by solving a phase or multiphase
The continuity equation of ratio is completed；To q phase, the multiphase flow conservation of mass is controlled by following formula:
Wherein, SaqFor quality source item, αqFor q phase volume fraction, ρqFor q phase density, t is the time；In default scenario
Lower formula (5) right end source item is zero, but when to each phase specified constant or user-defined Mass Sources, then right end is not zero；
The export boundary condition needs the condition met are as follows: sufficiently development, i.e., each physical quantity edge are flowed at mold exit
The normal derivative in the section is zero；
For the crystallizer wall boundary condition at crystallizer wall surface, the velocity component perpendicular to wall surface is zero, and parallel
In wall surface component using no slip boundary condition, meet formula (6):
Voutlet=Vcast, Vwall=kwall=0 (6)
Wherein, Voutlet is muzzle velocity, and Vcast is pulling rate, and Vwall is crystallizer wall face velocity, and kwall is crystallization
Wall face tubulence energy；
In the present embodiment, inlet velocity vinletIt can be calculated V according to formula (1)inlet=0.756m/s；Tubulence energy k can be by
Formula (2) obtains, and k=5.71536e-03m2/s；Tubulence energy dissipative shock wave ε=0.020575296m is obtained by formula (3)2/s3；Mould
Type top surface is set as free surface free surface, and the velocity component and every other each variable perpendicular to liquid level are along liquid level normal
The gradient in direction is set as zero, meets following condition:Export boundary condition outlet needs full
The condition of foot are as follows: sufficiently development is flowed at mold exit, i.e., normal derivative of each physical quantity along the section is zero；Crystallizer wall surface
Boundary condition moldwall is at crystallizer wall surface, and the velocity component perpendicular to wall surface is zero, and the component for being parallel to wall surface is adopted
With no slip boundary condition, meet following condition: Vout=Vcast=0.00533m/s, Vwall=kwall=0.
Step 2.4, setting method for solving and control parameter: setting method for solving opens Solution Methods panel,
Pressure PISO algorithm corresponding with speed coupled modes, gradient select Least Squares Cell Based, and pressure uses Body
Force Weighted format, the equation of momentum select Second OrderUpwind format；Setting solves control parameter, clicks
Solution controls panel keeps default；
Step 2.5, setting monitor window and initialization；When initialization, set in range at top liquid level 15mm
Path is steel slag, other regions are defaulted as molten steel；
Step 2.6, setting iterative steps and time step progress row parametric solution；
Shown in the restrictive condition of the time step such as formula (7):
Courant=u Δ t/ Δ x (7)
Wherein, u is fluid velocity, unit m/s；Δ t is time step, unit s；Δ x is size of mesh opening, and unit is
m；Courant is Krona spy number, indicates a fluid particle can pass through how many a grids in a time step；When gram
It when Lang Te number < 1, calculates and stablizes, but time step is small, calculate overlong time；As 1 < Kronas of spy number < 5, computational stability
Still seldom occur calculating diverging very well；When Krona spy number > 10, calculating is interrupted it is easy to appear diverging；Therefore, Krona is special
Number assists to select suitable time step；
In the present embodiment, Run Calculation panel is opened, inputs iteration step at Number of Iterations
Number, the input time step-length at Time Step Size, time step can obtain Δ t=0.04s by formula (7), be drawn by step 1
Subnetting lattice acquire Courant=1.2096 in the reasonable scope it is found that Δ x=0.025m, u=0.756m/s；It clicks
Calculate button starts to calculate；
The parameter for finally solving acquisition includes the entrance velocity vinnlet, tubulence energy k, tubulence energy dissipative shock wave ε and time step
Long Δ t, as shown in table 3:
The data of the acquisition parameter of table 3
Step 3, simulation extract data after calculating；In Ansys Fluent software click File → Export →
Solution extracts the data in crystallizer slag interface position, shape and neighbouring flow field, and in the poster processing soft Tecplot
In, establish the schematic diagram of the slag interface distribution and Flow Field Distribution situation of prediction as shown in Figure 3；
Step 4 analyzes post-processing result: obtaining image according to step 3, judges that slag interface fluctuates situation；
In the present embodiment, the slag interface that analysis finds the prediction of the straight hole mouth of a river is carried out to post-processing result as shown in Figure 3
It is very steady.
Step 5, adjustment slab mouth of a river insertion depth and mouth of a river shape re-execute the steps 1-4 and fluctuate feelings to slag interface
Condition is predicted, and is compared and analyzed to the slag interface fluctuation behavior of adjustment front and back.
In the present embodiment, mouth of a river insertion depth is increased into 20mm or the straight hole mouth of a river is changed to the three hole mouths of a river, with this size
It re-execute the steps 1-4 and carries out 3 d geometric modeling, calculating, post-processing.
After the slab mouth of a river is changed to three hole eddy flow downspouts, as shown in Figure 4, wherein 1 is slab eddy flow downspout, and internal diameter is
42mm, outer diameter 100mm；2 be slag interface in crystallizer；3 export for slab eddy flow downspout, and internal diameter 25mm, drift angle is downward
15 °, be 10mm apart from mouth of a river bottom；4 be R be 250mm slab；5 be electromagnetic stirring equipment；6 be continuous cast mold；It obtains at this time
Schematic diagram and as shown in FIG. 6 change of the slag interface distribution with Flow Field Distribution situation after changing mouth of a river shape as shown in Figure 5
The rectangular coordinate system figure of the liquid fluctuating height of slab median plane radial distance corresponding thereto after the shape of the mouth of a river.
As seen from Figure 6, when changing mouth of a river shape, liquid fluctuating height increases to 4.5mm from 0mm, illustrates mouth of a river shape
When shape changes into three holes by straight hole, steel slag liquid fluctuating increases.Crystallizer steel in continuous casting can be predicted by this method
Slag interface fluctuation.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal
Replacement；And these are modified or replaceed, model defined by the claims in the present invention that it does not separate the essence of the corresponding technical solution
It encloses.
Claims (3)
1. a kind of method of crystallizer slag interface fluctuation in prediction continuous casting, it is characterised in that: the following steps are included:
Step 1, according to steel mill's continuous casting installation for casting process, establish crystallizer 3-D geometric model with Gambit modeling software,
And gridding is carried out to model, formulates model boundary condition, then export .msh file；The boundary condition includes entrance boundary
Condition, export boundary condition, model top surface boundary condition and crystallizer wall surface and mouth of a river wall boundary condition；
Step 2 will carry out simulation calculating in the importing fluid calculation software of .msh file derived from step 1, method particularly includes:
Step 2.1, grid import and processing: the .msh file saved in read step 1 checks grid, guarantees grid most corpusculum
Net region size is arranged less than 0 in product；
Step 2.2, setting physical model and material: being Y-axis losing side according to the type of actual setting solver, acceleration of gravity
To, select turbulence model and multiphase flow model for computation model, setting material be fluid, be arranged slag interface surface tension；
Step 2.3, setting operating environment and boundary condition；
Step 2.4, setting method for solving and control parameter；
Step 2.5, setting monitor window and initialization；
When initialization, it is being set as steel slag in range at top liquid level 15mm, other regions are defaulted as molten steel；
Step 2.6, setting iterative steps and time step, carry out parametric solution；
The parameter for finally solving acquisition includes the entrance velocity vinlet, tubulence energy k, tubulence energy dissipative shock wave ε and time step Δ
t；
Step 3, simulation extract data after calculating；Crystallizer slag interface position, shape are extracted in fluid calculation software
And the neighbouring data in flow field establish slab median plane radial distance corresponding thereto and in the poster processing soft Tecplot
Liquid fluctuating height rectangular coordinate system scattergram；
Step 4 analyzes post-processing result: the rectangular coordinate system scattergram obtained according to step 3 judges steel slag circle
Surface wave emotionally condition；
Step 5, adjustment slab mouth of a river insertion depth and mouth of a river shape, re-execute the steps 1-4 to slag interface fluctuation situation into
Row prediction, and the slag interface fluctuation behavior of adjustment front and back is compared and analyzed.
2. the method for crystallizer slag interface fluctuation, feature exist in a kind of prediction continuous casting according to claim 1
In: setting boundary condition described in step 2.3 method particularly includes:
The entrance boundary condition includes inlet velocity vinlet, tubulence energy k and tubulence energy dissipative shock wave ε, respectively according to formula (1)-
(3) it obtains:
Wherein, vcastFor drawing speed, SoutFor mold exit area of section, SinFor submersed nozzle inlet area of section；
Wherein, RnozFor nozzle inlet radius；
The model top surface is set as free surface, and the velocity component and every other each variable perpendicular to liquid level are along liquid level normal side
To gradient be set as zero, meet formula (4):
Wherein, v is the velocity magnitude of model Free Surface x and the direction z, h=CpTinFor molten steel sensible heat, CpFor specific heat, TinTo enter
Mouth temperature, u are throwing direction speed；
Using the interface between VOF model following phase for multiphase flow model is the volume ratio by solving a phase or multiphase
Continuity equation complete；To q phase, the multiphase flow conservation of mass is controlled by following formula:
Wherein, SaqFor quality source item, αqFor q phase volume fraction, ρqFor q phase density, t is the time；The formula under default scenario
(5) right end source item is zero, but when to each phase specified constant or user-defined Mass Sources, then right end is not zero；
The export boundary condition needs the condition met are as follows: sufficiently development is flowed at mold exit, i.e., each physical quantity is along this section
The normal derivative in face is zero；
For the crystallizer wall boundary condition at crystallizer wall surface, the velocity component perpendicular to wall surface is zero, and is parallel to wall
The component in face meets formula (6) using no slip boundary condition:
Voutlet=Vcast, Vwall=kwall=0 (6)
Wherein, Voutlet is muzzle velocity, and Vcast is pulling rate, and Vwall is crystallizer wall face velocity, and kwall is crystallizer wall
Face tubulence energy.
3. the method for crystallizer slag interface fluctuation, feature exist in a kind of prediction continuous casting according to claim 1
In: shown in the restrictive condition of time step described in step 2.6 such as formula (7):
Courant=u Δ t/ Δ x (7)
Wherein, u is fluid velocity, and Δ t is time step, and Δ x is size of mesh opening, and Courant is Krona spy number, is indicated one
A fluid particle can pass through how many a grids in a time step；As Krona spy number < 1, calculates and stablize, but time step
Length is small, calculates overlong time；As 1 < Kronas of spy number < 5, computational stability still seldom occurs calculating diverging very well；When gram
When Lang Te number > 10, calculating is interrupted it is easy to appear diverging；Therefore, suitable time step is selected in Krona spy number assistance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910078981.4A CN109530648A (en) | 2019-01-28 | 2019-01-28 | A kind of method of crystallizer slag interface fluctuation in prediction continuous casting |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910078981.4A CN109530648A (en) | 2019-01-28 | 2019-01-28 | A kind of method of crystallizer slag interface fluctuation in prediction continuous casting |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109530648A true CN109530648A (en) | 2019-03-29 |
Family
ID=65838734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910078981.4A Pending CN109530648A (en) | 2019-01-28 | 2019-01-28 | A kind of method of crystallizer slag interface fluctuation in prediction continuous casting |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109530648A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111507035A (en) * | 2020-04-17 | 2020-08-07 | 东北大学 | Design method for traveling wave electromagnetic stirring at solidification tail end of continuous casting billet |
| CN112115634A (en) * | 2020-09-21 | 2020-12-22 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
| CN112182852A (en) * | 2020-09-08 | 2021-01-05 | 山东莱钢永锋钢铁有限公司 | Method and device for predicting service life of crystallizer copper pipe |
| CN113198994A (en) * | 2021-04-25 | 2021-08-03 | 东北大学 | Method for improving subcutaneous negative segregation in large round billet crystallizer in continuous casting process |
| CN113305277A (en) * | 2021-03-16 | 2021-08-27 | 北部湾大学 | Method for judging slag coiling in plate blank crystallizer casting powder shearing |
| CN113828746A (en) * | 2021-09-22 | 2021-12-24 | 日照钢铁控股集团有限公司 | Method for evaluating crystallizer flow field by utilizing vibration trace distribution of casting blank |
| CN113926865A (en) * | 2020-06-29 | 2022-01-14 | 宝山钢铁股份有限公司 | Casting blank slag inclusion forecasting method, machine cleaning control method, computing device and storage medium |
| CN114528787A (en) * | 2022-04-24 | 2022-05-24 | 北京科技大学 | Method and device for evaluating flow uniformity of outlet of long nozzle and impact area of tundish |
| CN114918394A (en) * | 2022-04-22 | 2022-08-19 | 首钢集团有限公司 | Method and device for controlling flow field bias of crystallizer |
| CN115106499A (en) * | 2022-06-30 | 2022-09-27 | 北京科技大学 | Crystallizer liquid level abnormal fluctuation distinguishing method and system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS55149761A (en) * | 1979-05-10 | 1980-11-21 | Mishima Kosan Co Ltd | Continuous casting control device |
| JP2003145257A (en) * | 2001-11-09 | 2003-05-20 | Shinkawa Denki Kk | Method for judging completion in pouring of molten steel and judging instrument therefor |
| CN108345737A (en) * | 2018-02-05 | 2018-07-31 | 东北大学 | A kind of design method of bloom continuous casting eddy flow downspout |
| CN108920821A (en) * | 2018-06-30 | 2018-11-30 | 马鞍山钢铁股份有限公司 | A kind of steel-making continuous casting intrusive mood flowing water mouth and its method for numerical simulation |
| CN108959749A (en) * | 2018-06-22 | 2018-12-07 | 东北大学 | A kind of design method of round billet continuous casting eddy flow downspout |
-
2019
- 2019-01-28 CN CN201910078981.4A patent/CN109530648A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS55149761A (en) * | 1979-05-10 | 1980-11-21 | Mishima Kosan Co Ltd | Continuous casting control device |
| JP2003145257A (en) * | 2001-11-09 | 2003-05-20 | Shinkawa Denki Kk | Method for judging completion in pouring of molten steel and judging instrument therefor |
| CN108345737A (en) * | 2018-02-05 | 2018-07-31 | 东北大学 | A kind of design method of bloom continuous casting eddy flow downspout |
| CN108959749A (en) * | 2018-06-22 | 2018-12-07 | 东北大学 | A kind of design method of round billet continuous casting eddy flow downspout |
| CN108920821A (en) * | 2018-06-30 | 2018-11-30 | 马鞍山钢铁股份有限公司 | A kind of steel-making continuous casting intrusive mood flowing water mouth and its method for numerical simulation |
Non-Patent Citations (1)
| Title |
|---|
| 庄迎等: "不锈钢方坯连铸结晶器内钢渣界面行为的数值模拟研究", 《连铸》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111507035A (en) * | 2020-04-17 | 2020-08-07 | 东北大学 | Design method for traveling wave electromagnetic stirring at solidification tail end of continuous casting billet |
| CN113926865B (en) * | 2020-06-29 | 2024-03-08 | 宝山钢铁股份有限公司 | Casting blank slag inclusion prediction method, machine cleaning control method, computing equipment and storage medium |
| CN113926865A (en) * | 2020-06-29 | 2022-01-14 | 宝山钢铁股份有限公司 | Casting blank slag inclusion forecasting method, machine cleaning control method, computing device and storage medium |
| CN112182852A (en) * | 2020-09-08 | 2021-01-05 | 山东莱钢永锋钢铁有限公司 | Method and device for predicting service life of crystallizer copper pipe |
| CN112115634A (en) * | 2020-09-21 | 2020-12-22 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
| CN112115634B (en) * | 2020-09-21 | 2021-04-06 | 哈尔滨理工大学 | Three-dimensional numerical prediction method for grain structure in unidirectional solidification process of molten metal |
| CN113305277A (en) * | 2021-03-16 | 2021-08-27 | 北部湾大学 | Method for judging slag coiling in plate blank crystallizer casting powder shearing |
| CN113198994A (en) * | 2021-04-25 | 2021-08-03 | 东北大学 | Method for improving subcutaneous negative segregation in large round billet crystallizer in continuous casting process |
| CN113828746A (en) * | 2021-09-22 | 2021-12-24 | 日照钢铁控股集团有限公司 | Method for evaluating crystallizer flow field by utilizing vibration trace distribution of casting blank |
| CN114918394A (en) * | 2022-04-22 | 2022-08-19 | 首钢集团有限公司 | Method and device for controlling flow field bias of crystallizer |
| CN114528787A (en) * | 2022-04-24 | 2022-05-24 | 北京科技大学 | Method and device for evaluating flow uniformity of outlet of long nozzle and impact area of tundish |
| CN115106499A (en) * | 2022-06-30 | 2022-09-27 | 北京科技大学 | Crystallizer liquid level abnormal fluctuation distinguishing method and system |
| CN115106499B (en) * | 2022-06-30 | 2024-02-20 | 北京科技大学 | Method and system for judging abnormal fluctuation of liquid level of crystallizer |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109530648A (en) | A kind of method of crystallizer slag interface fluctuation in prediction continuous casting | |
| Jin et al. | Modeling and measurements of multiphase flow and bubble entrapment in steel continuous casting | |
| Liu et al. | Transient asymmetric flow and bubble transport inside a slab continuous-casting mold | |
| Chen et al. | Numerical simulation of steel and argon gas two-phase flow in continuous casting using LES+ VOF+ DPM model | |
| Liu et al. | Modeling of quasi-four-phase flow in continuous casting mold using hybrid Eulerian and Lagrangian approach | |
| Liu et al. | Euler-Euler-Lagrangian modeling for two-phase flow and particle transport in continuous casting mold | |
| Li et al. | Investigation of bubble-slag layer behaviors with hybrid Eulerian–Lagrangian modeling and large eddy simulation | |
| Javurek et al. | Simulation of nonmetallic inclusions in a continuous casting strand | |
| Wang et al. | Transient fluid flow phenomena during continuous casting: part I—cast start | |
| Chaudhary et al. | Effect of stopper-rod misalignment on fluid flow in continuous casting of steel | |
| Cwudziński | Numerical simulation of liquid steel flow and behaviour of non-metallic inclusions in one-strand slab tundish with subflux turbulence controller and gas permeable barrier | |
| CN108959749A (en) | A kind of design method of round billet continuous casting eddy flow downspout | |
| Yang et al. | Transient thermo-fluid and solidification behaviors in continuous casting mold: evolution phenomena | |
| Liu et al. | Effect of vertical length on asymmetric flow and inclusion transport in vertical-bending continuous caster | |
| CN106227954B (en) | A kind of Aluminum alloy gravity gravity die casting process optimization method | |
| Chattopadhyay et al. | Physical and mathematical modeling of inert gas-shrouded ladle nozzles and their role on slag behavior and fluid flow patterns in a delta-shaped, four-strand tundish | |
| CN108345737B (en) | Design method of bloom continuous casting rotational flow water gap | |
| CN108920821A (en) | A kind of steel-making continuous casting intrusive mood flowing water mouth and its method for numerical simulation | |
| Zhu et al. | Influence of large-scale vortex movement in lower recirculation zone on instable flow field in the mold | |
| Liu et al. | Transient turbulent flow simulation with water model validation and application to slide gate dithering | |
| Li et al. | In-situ analysis and numerical study of inclusion distribution in a vertical-bending caster | |
| Chaudhary et al. | Differences between physical water models and steel continuous casters: A theoretical evaluation | |
| Egole et al. | Micro-macro model for the transient heat and fluid transport in solidification structure evolution during static casting processes | |
| Zhang et al. | Mathematical and physical modeling of metal delivery system during top side-pouring twin-roll casting | |
| Song et al. | Numerical modeling of grade mixing and inclusion entrapment in eight strand billet tundish |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190329 |