Tag: Physical Modeling

Abstract

This paper presents the results of tests on the suitability of designed heads (impellers) for aluminum refining. The research was carried out on a physical model of the URO-200, followed by numerical simulations in the FLOW 3D program. Four design variants of impellers were used in the study. The degree of dispersion of the gas phase in the model liquid was used as a criterion for evaluating the performance of each solution using different process parameters, i.e., gas flow rate and impeller speed. Afterward, numerical simulations in Flow 3D software were conducted for the best solution. These simulations confirmed the results obtained with the water model and verified them.

Keywords: aluminum, impeller construction, degassing process, numerical modeling, physical modeling

1. Introduction

Constantly increasing requirements concerning metallurgical purity in terms of hydrogen content and nonmetallic inclusions make casting manufacturers use effective refining techniques. The answer to this demand is the implementation of the aluminum refining technique making use of a rotor with an original design guaranteeing efficient refining [1,2,3,4]. The main task of the impeller (rotor) is to reduce the contamination of liquid metal (primary and recycled aluminum) with hydrogen and nonmetallic inclusions. An inert gas, mainly argon or a mixture of gases, is introduced through the rotor into the liquid metal to bring both hydrogen and nonmetallic inclusions to the metal surface through the flotation process. Appropriately and uniformly distributed gas bubbles in the liquid metal guarantee achieving the assumed level of contaminant removal economically. A very important factor in deciding about the obtained degassing effect is the optimal rotor design [5,6,7,8]. Thanks to the appropriate geometry of the rotor, gas bubbles introduced into the liquid metal are split into smaller ones, and the spinning movement of the rotor distributes them throughout the volume of the liquid metal bath. In this solution impurities in the liquid metal are removed both in the volume and from the upper surface of the metal. With a well-designed impeller, the costs of refining aluminum and its alloys can be lowered thanks to the reduced inert gas and energy consumption (optimal selection of rotor rotational speed). Shorter processing time and a high degree of dehydrogenation decrease the formation of dross on the metal surface (waste). A bigger produced dross leads to bigger process losses. Consequently, this means that the choice of rotor geometry has an indirect impact on the degree to which the generated waste is reduced [9,10].

Another equally important factor is the selection of process parameters such as gas flow rate and rotor speed [11,12]. A well-designed gas injection system for liquid metal meets two key requirements; it causes rapid mixing of the liquid metal to maintain a uniform temperature throughout the volume and during the entire process, to produce a chemically homogeneous metal composition. This solution ensures effective degassing of the metal bath. Therefore, the shape of the rotor, the arrangement of the nozzles, and their number are significant design parameters that guarantee the optimum course of the refining process. It is equally important to complete the mixing of the metal bath in a relatively short time, as this considerably shortens the refining process and, consequently, reduces the process costs. Another important criterion conditioning the implementation of the developed rotor is the generation of fine diffused gas bubbles which are distributed throughout the metal volume, and whose residence time will be sufficient for the bubbles to collide and adsorb the contaminants. The process of bubble formation by the spinning rotors differs from that in the nozzles or porous molders. In the case of a spinning rotor, the shear force generated by the rotor motion splits the bubbles into smaller ones. Here, the rotational speed, mixing force, surface tension, and fluid density have a key effect on the bubble size. The velocity of the bubbles, which depends mainly on their size and shape, determines their residence time in the reactor and is, therefore, very important for the refining process, especially since gas bubbles in liquid aluminum may remain steady only below a certain size [13,14,15].

The impeller designs presented in the article were developed to improve the efficiency of the process and reduce its costs. The impellers used so far have a complicated structure and are very pricey. The success of the conducted research will allow small companies to become independent of external supplies through the possibility of making simple and effective impellers on their own. The developed structures were tested on the water model. The results of this study can be considered as pilot.

2. Materials and Methods

Rotors were realized with the SolidWorks computer design technique and a 3D printer. The developed designs were tested on a water model. Afterward, the solution with the most advantageous refining parameters was selected and subjected to calculations with the Flow3D package. As a result, an impeller was designed for aluminum refining. Its principal lies in an even distribution of gas bubbles in the entire volume of liquid metal, with the largest possible participation of the bubble surface, without disturbing the metal surface. This procedure guarantees the removal of gaseous, as well as metallic and nonmetallic, impurities.

2.1. Rotor Designs

The developed impeller constructions, shown in Figure 1, Figure 2, Figure 3 and Figure 4, were printed on a 3D printer using the PLA (polylactide) material. The impeller design models differ in their shape and the number of holes through which the inert gas flows. Figure 1, Figure 2 and Figure 3 show the same impeller model but with a different number of gas outlets. The arrangement of four, eight, and 12 outlet holes was adopted in the developed design. A triangle-shaped structure equipped with three gas outlet holes is presented in Figure 4.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g001.jpg

Figure 1

A 3D model—impeller with four holes—variant B4.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g002.jpg

Figure 2

A 3D model—impeller with eight holes—variant B8.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g003.jpg

Figure 3

A 3D model—impeller with twelve holes—variant B12.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g004.jpg

Figure 4

A 3D model—‘red triangle’ impeller with three holes—variant RT3.

2.2. Physical Models

Investigations were carried out on a water model of the URO 200 reactor of the barbotage refining process (see Figure 5).

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g005.jpg

Figure 5

A schematic of the water model of reactor URO 200.

The URO 200 reactor can be classified as a cyclic reactor. The main element of the device is a rotor, which ends the impeller. The whole system is attached to a shaft via which the refining gas is supplied. Then, the shaft with the rotor is immersed in the liquid metal in the melting pot or the furnace chamber. In URO 200 reactors, the refining process lasts 600 s (10 min), the gas flow rate that can be obtained ranges from 5 to 20 dm³·min⁻¹, and the speed at which the rotor can move is 0 to 400 rpm. The permissible quantity of liquid metal for barbotage refining is 300 kg or 700 kg [8,16,17]. The URO 200 has several design solutions which improve operation and can be adapted to the existing equipment in the foundry. These solutions include the following [8,16]:

URO-200XR—used for small crucible furnaces, the capacity of which does not exceed 250 kg, with no control system and no control of the refining process.
URO-200SA—used to service several crucible furnaces of capacity from 250 kg to 700 kg, fully automated and equipped with a mechanical rotor lift.
URO-200KA—used for refining processes in crucible furnaces and allows refining in a ladle. The process is fully automated, with a hydraulic rotor lift.
URO-200KX—a combination of the XR and KA models, designed for the ladle refining process. Additionally, refining in heated crucibles is possible. The unit is equipped with a manual hydraulic rotor lift.
URO-200PA—designed to cooperate with induction or crucible furnaces or intermediate chambers, the capacity of which does not exceed one ton. This unit is an integral part of the furnace. The rotor lift is equipped with a screw drive.

Studies making use of a physical model can be associated with the observation of the flow and circulation of gas bubbles. They require meeting several criteria regarding the similarity of the process and the object characteristics. The similarity conditions mainly include geometric, mechanical, chemical, thermal, and kinetic parameters. During simulation of aluminum refining with inert gas, it is necessary to maintain the geometric similarity between the model and the real object, as well as the similarity related to the flow of liquid metal and gas (hydrodynamic similarity). These quantities are characterized by the Reynolds, Weber, and Froude numbers. The Froude number is the most important parameter characterizing the process, its magnitude is the same for the physical model and the real object. Water was used as the medium in the physical modeling. The factors influencing the choice of water are its availability, relatively low cost, and kinematic viscosity at room temperature, which is very close to that of liquid aluminum.

The physical model studies focused on the flow of inert gas in the form of gas bubbles with varying degrees of dispersion, particularly with respect to some flow patterns such as flow in columns and geysers, as well as disturbance of the metal surface. The most important refining parameters are gas flow rate and rotor speed. The barbotage refining studies for the developed impeller (variants B4, B8, B12, and RT3) designs were conducted for the following process parameters:

Rotor speed: 200, 300, 400, and 500 rpm,
Ideal gas flow: 10, 20, and 30 dm³·min⁻¹,
Temperature: 293 K (20 °C).

These studies were aimed at determining the most favorable variants of impellers, which were then verified using the numerical modeling methods in the Flow-3D program.

2.3. Numerical Simulations with Flow-3D Program

Testing different rotor impellers using a physical model allows for observing the phenomena taking place while refining. This is a very important step when testing new design solutions without using expensive industrial trials. Another solution is modeling by means of commercial simulation programs such as ANSYS Fluent or Flow-3D [18,19]. Unlike studies on a physical model, in a computer program, the parameters of the refining process and the object itself, including the impeller design, can be easily modified. The simulations were performed with the Flow-3D program version 12.03.02. A three-dimensional system with the same dimensions as in the physical modeling was used in the calculations. The isothermal flow of liquid–gas bubbles was analyzed. As in the physical model, three speeds were adopted in the numerical tests: 200, 300, and 500 rpm. During the initial phase of the simulations, the velocity field around the rotor generated an appropriate direction of motion for the newly produced bubbles. When the required speed was reached, the generation of randomly distributed bubbles around the rotor was started at a rate of 2000 per second. Table 1 lists the most important simulation parameters.

Table 1

Values of parameters used in the calculations.

Parameter	Value	Unit
Maximum number of gas particles	1,000,000	–
Rate of particle generation	2000	1·s⁻¹
Specific gas constant	287.058	J·kg⁻¹·K⁻¹
Atmospheric pressure	1.013 × 10⁵	Pa
Water density	1000	kg·m⁻³
Water viscosity	0.001	kg·m⁻¹·s⁻¹
Boundary condition on the walls	No-slip	–
Size of computational cell	0.0034	m

In the case of the CFD analysis, the numerical solutions require great care when generating the computational mesh. Therefore, computational mesh tests were performed prior to the CFD calculations. The effect of mesh density was evaluated by taking into account the velocity of water in the tested object on the measurement line A (height of 0.065 m from the bottom) in a characteristic cross-section passing through the object axis (see Figure 6). The mesh contained 3,207,600, 6,311,981, 7,889,512, 11,569,230, and 14,115,049 cells.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g006.jpg

Figure 6

The velocity of the water depending on the size of the computational grid.

The quality of the generated computational meshes was checked using the criterion skewness angle Q_EAS [18]. This criterion is described by the following relationship:

QEAS=max{βmax−βeq180−βeq,βeq−βminβeq},

(1)

where β_max, β_min are the maximal and minimal angles (in degrees) between the edges of the cell, and β_eq is the angle corresponding to an ideal cell, which for cubic cells is 90°.

Normalized in the interval [0;1], the value of Q_EAS should not exceed 0.75, which identifies the permissible skewness angle of the generated mesh. For the computed meshes, this value was equal to 0.55–0.65.

Moreover, when generating the computational grids in the studied facility, they were compacted in the areas of the highest gradients of the calculated values, where higher turbulence is to be expected (near the impeller). The obtained results of water velocity in the studied object at constant gas flow rate are shown in Figure 6.

The analysis of the obtained water velocity distributions (see Figure 6) along the line inside the object revealed that, with the density of the grid of nodal points, the velocity changed and its changes for the test cases of 7,889,512, 11,569,230, and 14,115,049 were insignificant. Therefore, it was assumed that a grid containing not less than 7,900,000 (7,889,512) cells would not affect the result of CFD calculations.

A single-block mesh of regular cells with a size of 0.0034 m was used in the numerical calculations. The total number of cells was approximately 7,900,000 (7,889,512). This grid resolution (see Figure 7) allowed the geometry of the system to be properly represented, maintaining acceptable computation time (about 3 days on a workstation with 2× CPU and 12 computing cores).

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g007.jpg

Figure 7

Structured equidistant mesh used in numerical calculations: (a) mesh with smoothed, surface cells (the so-called FAVOR method) used in Flow-3D; (b) visualization of the applied mesh resolution.

The calculations were conducted with an explicit scheme. The timestep was selected by the program automatically and controlled by stability and convergence. From the moment of the initial velocity field generation (start of particle generation), it was 0.0001 s.

When modeling the degassing process, three fluids are present in the system: water, gas supplied through the rotor head (impeller), and the surrounding air. Modeling such a multiphase flow is a numerically very complex issue. The necessity to overcome the liquid backpressure by the gas flowing out from the impeller leads to the formation of numerical instabilities in the volume of fluid (VOF)-based approach used by Flow-3D software. Therefore, a mixed description of the analyzed flow was used here. In this case, water was treated as a continuous medium, while, in the case of gas bubbles, the discrete phase model (DPM) model was applied. The way in which the air surrounding the system was taken into account is later described in detail.

The following additional assumptions were made in the modeling:

—The liquid phase was considered as an incompressible Newtonian fluid.
—The effect of chemical reactions during the refining process was neglected.
—The composition of each phase (gas and liquid) was considered homogeneous; therefore, the viscosity and surface tension were set as constants.
—Only full turbulence existed in the liquid, and the effect of molecular viscosity was neglected.
—The gas bubbles were shaped as perfect spheres.
—The mutual interaction between gas bubbles (particles) was neglected.

2.3.1. Modeling of Liquid Flow

The motion of the real fluid (continuous medium) is described by the Navier–Stokes Equation [20].

dudt=−1ρ∇p+ν∇2u+13ν∇(∇⋅ u)+F,

(2)

where du/dt is the time derivative, u is the velocity vector, t is the time, and F is the term accounting for external forces including gravity (unit components denoted by X, Y, Z).

In the simulations, the fluid flow was assumed to be incompressible, in which case the following equation is applicable:

∂u∂t+(u⋅∇)u=−1ρ∇p+ν∇2u+F.

(3)

Due to the large range of liquid velocities during flows, the turbulence formation process was included in the modeling. For this purpose, the k–ε model turbulence kinetic energy k and turbulence dissipation ε were the target parameters, as expressed by the following equations [21]:

∂(ρk)∂t+∂(ρkvi)∂xi=∂∂xj[(μ+μtσk)⋅∂k∂xi]+Gk+Gb−ρε−Ym+Sk,

(4)

∂(ρε)∂t+∂(ρεui)∂xi=∂∂xj[(μ+μtσε)⋅∂k∂xi]+C1εεk(Gk+G3εGb)+C2ερε2k+Sε,

(5)

where ρ is the gas density, σ_κ and σ_ε are the Prandtl turbulence numbers, k and ε are constants of 1.0 and 1.3, and G_k and G_b are the kinetic energy of turbulence generated by the average velocity and buoyancy, respectively.

As mentioned earlier, there are two gas phases in the considered problem. In addition to the gas bubbles, which are treated here as particles, there is also air, which surrounds the system. The boundary of phase separation is in this case the free surface of the water. The shape of the free surface can change as a result of the forming velocity field in the liquid. Therefore, it is necessary to use an appropriate approach to free surface tracking. The most commonly used concept in liquid–gas flow modeling is the volume of fluid (VOF) method [22,23], and Flow-3D uses a modified version of this method called TrueVOF. It introduces the concept of the volume fraction of the liquid phase f_l. This parameter can be used for classifying the cells of a discrete grid into areas filled with liquid phase (f_l = 1), gaseous phase, or empty cells (f_l = 0) and those through which the phase separation boundary (f_l ∈ (0, 1)) passes (free surface). To determine the local variations of the liquid phase fraction, it is necessary to solve the following continuity equation:

dfldt=0.

(6)

Then, the fluid parameters in the region of coexistence of the two phases (the so-called interface) depend on the volume fraction of each phase.

ρ=flρl+(1−fl)ρg,

(7)

ν=flνl+(1−fl)νg,

(8)

where indices l and g refer to the liquid and gaseous phases, respectively.

The parameter of fluid velocity in cells containing both phases is also determined in the same way.

u=flul+(1−fl)ug.

(9)

Since the processes taking place in the surrounding air can be omitted, to speed up the calculations, a single-phase, free-surface model was used. This means that no calculations were performed in the gas cells (they were treated as empty cells). The liquid could fill them freely, and the air surrounding the system was considered by the atmospheric pressure exerted on the free surface. This approach is often used in modeling foundry and metallurgical processes [24].

2.3.2. Modeling of Gas Bubble Flow

As stated, a particle model was used to model bubble flow. Spherical particles (gas bubbles) of a given size were randomly generated in the area marked with green in Figure 7b. In the simulations, the gas bubbles were assumed to have diameters of 0.016 and 0.02 m corresponding to the gas flow rates of 10 and 30 dm³·min⁻¹, respectively.

Experimental studies have shown that, as a result of turbulent fluid motion, some of the bubbles may burst, leading to the formation of smaller bubbles, although merging of bubbles into larger groupings may also occur. Therefore, to be able to observe the behavior of bubbles of different sizes (diameter), the calculations generated two additional particle types with diameters twice smaller and twice larger, respectively. The proportion of each species in the system was set to 33.33% (Table 2).

Table 2

Data assumed for calculations.

No	Rotor Speed (Rotational Speed) rpm	Bubbles Diameter m	Corresponding Gas Flow Rate dm³·min⁻¹	No	Rotor Speed (Rotational Speed) rpm	Bubbles Diameter m	Corresponding Gas Flow Rate dm³·min⁻¹
A	200	0.016	10	D	200	0.02	30
0.008	0.01
0.032	0.04
B	300	0.016	10	E	300	0.02	30
0.008	0.01
0.032	0.04
C	500	0.016	10	F	500	0.02	30
0.008	0.01
0.032	0.04

The velocity of the particle results from the generated velocity field (calculated from Equation (3) in the liquid ul around it and its velocity resulting from the buoyancy force ub. The effect of particle radius r on the terminal velocity associated with buoyancy force can be determined according to Stokes’ law.

ub=29 (ρg−ρl)μlgr2,

(10)

where g is the acceleration (9.81).

The DPM model was used for modeling the two-phase (water–air) flow. In this model, the fluid (water) is treated as a continuous phase and described by the Navier–Stokes equation, while gas bubbles are particles flowing in the model fluid (discrete phase). The trajectories of each bubble in the DPM system are calculated at each timestep taking into account the mass forces acting on it. Table 3 characterizes the DPM model used in our own research [18].

Table 3

Characteristic of the DPM model.

Method	Equations
Euler–Lagrange	Balance equation: dugdt=FD(u−ug)+g(ϱg−ϱ)ϱg+F. F_D (u − u_p) denotes the drag forces per mass unit of a bubble, and the expression for the drag coefficient F_D is of the form FD=18μCDReϱ⋅gd2g24. The relative Reynolds number has the form Re≡ρdg\|ug−u\|μ. On the other hand, the force resulting from the additional acceleration of the model fluid has the form F=12dρdtρg(u−ug), where u_g is the gas bubble velocity, u is the liquid velocity, d_g is the bubble diameter, and C_D is the drag coefficient.

Method

Equations

Euler–Lagrange

Balance equation:
dugdt=FD(u−ug)+g(ϱg−ϱ)ϱg+F.
F_D (u − u_p) denotes the drag forces per mass unit of a bubble, and the expression for the drag coefficient F_D is of the form
FD=18μCDReϱ⋅gd2g24.
The relative Reynolds number has the form
Re≡ρdg|ug−u|μ.
On the other hand, the force resulting from the additional acceleration of the model fluid has the form
F=12dρdtρg(u−ug),
where u_g is the gas bubble velocity, u is the liquid velocity, d_g is the bubble diameter, and C_D is the drag coefficient.

3. Results and Discussion

3.1. Calculations of Power and Mixing Time by the Flowing Gas Bubbles

One of the most important parameters of refining with a rotor is the mixing power induced by the spinning rotor and the outflowing gas bubbles (via impeller). The mixing power of liquid metal in a ladle of height (h) by gas injection can be determined from the following relation [15]:

pgVm=ρ⋅g⋅uB,

(11)

where p_g is the mixing power, V_m is the volume of liquid metal in the reactor, ρ is the density of liquid aluminum, and u_B is the average speed of bubbles, given below.

uB=n⋅R⋅TAc⋅Pm⋅t,

(12)

where n is the number of gas moles, R is the gas constant (8.314), A_c is the cross-sectional area of the reactor vessel, T is the temperature of liquid aluminum in the reactor, and P_m is the pressure at the middle tank level. The pressure at the middle level of the tank is calculated by a function of the mean logarithmic difference.

Pm=(Pa+ρ⋅g⋅h)−Paln(Pa+ρ⋅g⋅h)Pa,

(13)

where P_a is the atmospheric pressure, and h is the the height of metal in the reactor.

Themelis and Goyal [25] developed a model for calculating mixing power delivered by gas injection.

pg=2Q⋅R⋅T⋅ln(1+m⋅ρ⋅g⋅hP),

(14)

where Q is the gas flow, and m is the mass of liquid metal.

Zhang [26] proposed a model taking into account the temperature difference between gas and alloy (metal).

pg=QRTgVm[ln(1+ρ⋅g⋅hPa)+(1−TTg)],

(15)

where T_g is the gas temperature at the entry point.

Data for calculating the mixing power resulting from inert gas injection into liquid aluminum are given below in Table 4. The design parameters were adopted for the model, the parameters of which are shown in Figure 5.

Table 4

Data for calculating mixing power introduced by an inert gas.

Parameter	Value	Unit
Height of metal column	0.7	m
Density of aluminum	2375	kg·m⁻³
Process duration	20	s
Gas temperature at the injection site	940	K
Cross-sectional area of ladle	0.448	m²
Mass of liquid aluminum	546.25	kg
Volume of ladle	0.23	M³
Temperature of liquid aluminum	941.15	K

Table 5 presents the results of mixing power calculations according to the models of Themelis and Goyal and of Zhang for inert gas flows of 10, 20, and 30 dm³·min⁻¹. The obtained calculation results significantly differed from each other. The difference was an order of magnitude, which indicates that the model is highly inaccurate without considering the temperature of the injected gas. Moreover, the calculations apply to the case when the mixing was performed only by the flowing gas bubbles, without using a rotor, which is a great simplification of the phenomenon.

Table 5

Mixing power calculated from mathematical models.

Mathematical Model	Mixing Power (W·t⁻¹) for a Given Inert Gas Flow (dm³·min⁻¹)
	10	20	30
Themelis and Goyal	11.49	23.33	35.03
Zhang	0.82	1.66	2.49

The mixing time is defined as the time required to achieve 95% complete mixing of liquid metal in the ladle [27,28,29,30]. Table 6 groups together equations for the mixing time according to the models.

Table 6

Models for calculating mixing time.

Authors	Model	Remarks
Szekely [31]	τ=800ε−0.4	ε—W·t⁻¹
Chiti and Paglianti [27]	τ=CVQl	V—volume of reactor, m³ Ql—flow intensity, m³·s⁻¹
Iguchi and Nakamura [32]	τ=1200⋅Q−0.4D1.97h−1.0υ0.47	υ—kinematic viscosity, m²·s⁻¹ D—diameter of ladle, m h—height of metal column, m Q—liquid flow intensity, m³·s⁻¹

Figure 8 and Figure 9 show the mixing time as a function of gas flow rate for various heights of the liquid column in the ladle and mixing power values.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g008.jpg

Figure 8

Mixing time as a function of gas flow rate for various heights of the metal column (Iguchi and Nakamura model).

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g009.jpg

Figure 9

Mixing time as a function of mixing power (Szekly model).

3.2. Determining the Bubble Size

The mechanisms controlling bubble size and mass transfer in an alloy undergoing refining are complex. Strong mixing conditions in the reactor promote impurity mass transfer. In the case of a spinning rotor, the shear force generated by the rotor motion separates the bubbles into smaller bubbles. Rotational speed, mixing force, surface tension, and liquid density have a strong influence on the bubble size. To characterize the kinetic state of the refining process, parameters k and A were introduced. Parameters k, A, and u_B can be calculated using the below equations [33].

k=2D⋅uBdB⋅π−−−−−−√,

(16)

A=6Q⋅hdB⋅uB,

(17)

uB=1.02g⋅dB,−−−−−√

(18)

where D is the diffusion coefficient, and d_B is the bubble diameter.

After substituting appropriate values, we get

dB=3.03×104(πD)−2/5g−1/5h4/5Q0.344N−1.48.

(19)

According to the last equation, the size of the gas bubble decreases with the increasing rotational speed (see Figure 10).

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g010.jpg

Figure 10

Effect of rotational speed on the bubble diameter.

In a flow of given turbulence intensity, the diameter of the bubble does not exceed the maximum size d_max, which is inversely proportional to the rate of kinetic energy dissipation in a viscous flow ε. The size of the gas bubble diameter as a function of the mixing energy, also considering the Weber number and the mixing energy in the negative power, can be determined from the following equations [31,34]:

—Sevik and Park:

dBmax=We0.6kr⋅(σ⋅103ρ⋅10−3)0.6⋅(10⋅ε)−0.4⋅10−2.

(20)

—Evans:

dBmax=⎡⎣Wekr⋅σ⋅1032⋅(ρ⋅10−3)13⎤⎦35 ⋅(10⋅ε)−25⋅10−2.

(21)

The results of calculating the maximum diameter of the bubble d_Bmax determined from Equation (21) are given in Table 7.

Table 7

The results of calculating the maximum diameter of the bubble using Equation (21).

Model	Mixing Energy ĺ (m²·s⁻³)	Weber Number (We_kr)
0.59	1.0	1.2
Zhang and Taniguchi d_max	0.1	0.0167	0.0230	0.026
0.5	0.0088	0.0121	0.013
1.0	0.0067	0.0091	0.010
1.5	0.0057	0.0078	0.009
Sevik and Park d_B_max	0.1	0.265	0.36	0.41
0.5	0.139	0.19	0.21
1.0	0.106	0.14	0.16
1.5	0.090	0.12	0.14
Evans d_B_max	0.1	0.247	0.340	0.38
0.5	0.130	0.178	0.20
1.0	0.098	0.135	0.15
1.5	0.084	0.115	0.13

3.3. Physical Modeling

The first stage of experiments (using the URO-200 water model) included conducting experiments with impellers equipped with four, eight, and 12 gas outlets (variants B4, B8, B12). The tests were carried out for different process parameters. Selected results for these experiments are presented in Figure 11, Figure 12, Figure 13 and Figure 14.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g011.jpg

Figure 11

Impeller variant B4—gas bubbles dispersion registered for a gas flow rate of 10 dm³·min⁻¹ and rotor speed of (a) 200, (b) 300, (c) 400, and (d) 500 rpm.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g012.jpg

Figure 12

Impeller variant B8—gas bubbles dispersion registered for a gas flow rate of 10 dm³·min⁻¹ and rotor speed of (a) 200, (b) 300, (c) 400, and (d) 500 rpm.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g013.jpg

Figure 13

Gas bubble dispersion registered for different processing parameters (impeller variant B12).

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g014.jpg

Figure 14

Gas bubble dispersion registered for different processing parameters (impeller variant RT3).

The analysis of the refining variants presented in Figure 11, Figure 12, Figure 13 and Figure 14 reveals that the proposed impellers design model is not useful for the aluminum refining process. The number of gas outlet orifices, rotational speed, and flow did not affect the refining efficiency. In all the variants shown in the figures, very poor dispersion of gas bubbles was observed in the object. The gas bubble flow had a columnar character, and so-called dead zones, i.e., areas where no inert gas bubbles are present, were visible in the analyzed object. Such dead zones were located in the bottom and side zones of the ladle, while the flow of bubbles occurred near the turning rotor. Another negative phenomenon observed was a significant agitation of the water surface due to excessive (rotational) rotor speed and gas flow (see Figure 13, cases 20; 400, 30; 300, 30; 400, and 30; 500).

Research results for a ‘red triangle’ impeller equipped with three gas supply orifices (variant RT3) are presented in Figure 14.

In this impeller design, a uniform degree of bubble dispersion in the entire volume of the modeling fluid was achieved for most cases presented (see Figure 14). In all tested variants, single bubbles were observed in the area of the water surface in the vessel. For variants 20; 200, 30; 200, and 20; 300 shown in Figure 14, the bubble dispersion results were the worst as the so-called dead zones were identified in the area near the bottom and sidewalls of the vessel, which disqualifies these work parameters for further applications. Interestingly, areas where swirls and gas bubble chains formed were identified only for the inert gas flows of 20 and 30 dm³·min⁻¹ and 200 rpm in the analyzed model. This means that the presented model had the best performance in terms of dispersion of gas bubbles in the model liquid. Its design with sharp edges also differed from previously analyzed models, which is beneficial for gas bubble dispersion, but may interfere with its suitability in industrial conditions due to possible premature wear.

3.4. Qualitative Comparison of Research Results (CFD and Physical Model)

The analysis (physical modeling) revealed that the best mixing efficiency results were obtained with the RT3 impeller variant. Therefore, numerical calculations were carried out for the impeller model with three outlet orifices (variant RT3). The CFD results are presented in Figure 15 and Figure 16.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g015.jpg

Figure 15

Simulation results of the impeller RT3, for given flows and rotational speeds after a time of 1 s: simulation variants (a) A, (b) B, (c) C, (d) D, (e) E, and (f) F.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g016.jpg

Figure 16

Simulation results of the impeller RT3, for given flows and rotational speeds after a time of 5.4 s.: simulation variants (a) A, (b) B, (c) C, (d) D, (e) E, and (f) F.

CFD results are presented for all analyzed variants (impeller RT3) at two selected calculation timesteps of 1 and 5.40 s. They show the velocity field of the medium (water) and the dispersion of gas bubbles.

Figure 15 shows the initial refining phase after 1 s of the process. In this case, the gas bubble formation and flow were observed in an area close to contact with the rotor. Figure 16 shows the phase when the dispersion and flow of gas bubbles were advanced in the reactor area of the URO-200 model.

The quantitative evaluation of the obtained results of physical and numerical model tests was based on the comparison of the degree of gas dispersion in the model liquid. The degree of gas bubble dispersion in the volume of the model liquid and the areas of strong turbulent zones formation were evaluated during the analysis of the results of visualization and numerical simulations. These two effects sufficiently characterize the required course of the process from the physical point of view. The known scheme of the below description was adopted as a basic criterion for the evaluation of the degree of dispersion of gas bubbles in the model liquid.

Minimal dispersion—single bubbles ascending in the region of their formation along the ladle axis; lack of mixing in the whole bath volume.
Accurate dispersion—single and well-mixed bubbles ascending toward the bath mirror in the region of the ladle axis; no dispersion near the walls and in the lower part of the ladle.
Uniform dispersion—most desirable; very good mixing of fine bubbles with model liquid.
Excessive dispersion—bubbles join together to form chains; large turbulence zones; uneven flow of gas.

The numerical simulation results give a good agreement with the experiments performed with the physical model. For all studied variants (used process parameters), the single bubbles were observed in the area of water surface in the vessel. For variants presented in Figure 13 (200 rpm, gas flow 20 and dm³·min⁻¹) and relevant examples in numerical simulation Figure 16, the worst bubble dispersion results were obtained because the dead zones were identified in the area near the bottom and sidewalls of the vessel, which disqualifies these work parameters for further use. The areas where swirls and gas bubble chains formed were identified only for the inert gas flows of 20 and 30 dm³·min⁻¹ and 200 rpm in the analyzed model (physical model). This means that the presented impeller model had the best performance in terms of dispersion of gas bubbles in the model liquid. The worst bubble dispersion results were obtained because the dead zones were identified in the area near the bottom and side walls of the vessel, which disqualifies these work parameters for further use.

Figure 17 presents exemplary results of model tests (CFD and physical model) with marked gas bubble dispersion zones. All variants of tests were analogously compared, and this comparison allowed validating the numerical model.

An external file that holds a picture, illustration, etc.
Object name is materials-15-05273-g017.jpg

Figure 17

Compilations of model research results (CFD and physical): A—single gas bubbles formed on the surface of the modeling liquid, B—excessive formation of gas chains and swirls, C—uniform distribution of gas bubbles in the entire volume of the tank, and D—dead zones without gas bubbles, no dispersion. (a) Variant B; (b) variant F.

It should be mentioned here that, in numerical simulations, it is necessary to make certain assumptions and simplifications. The calculations assumed three particle size classes (Table 2), which represent the different gas bubbles that form due to different gas flow rates. The maximum number of particles/bubbles (Table 1) generated was assumed in advance and related to the computational capabilities of the computer. Too many particles can also make it difficult to visualize and analyze the results. The size of the particles, of course, affects their behavior during simulation, while, in the figures provided in the article, the bubbles are represented by spheres (visualization of the results) of the same size. Please note that, due to the adopted Lagrangian–Eulerian approach, the simulation did not take into account phenomena such as bubble collapse or fusion. However, the obtained results allow a comprehensive analysis of the behavior of gas bubbles in the system under consideration.

The comparative analysis of the visualization (quantitative) results obtained with the water model and CFD simulations (see Figure 17) generated a sufficient agreement from the point of view of the trends. A precise quantitative evaluation is difficult to perform because of the lack of a refraction compensating system in the water model. Furthermore, in numerical simulations, it is not possible to determine the geometry of the forming gas bubbles and their interaction with each other as opposed to the visualization in the water model. The use of both research methods is complementary. Thus, a direct comparison of images obtained by the two methods requires appropriate interpretation. However, such an assessment gives the possibility to qualitatively determine the types of the present gas bubble dispersion, thus ultimately validating the CFD results with the water model.

A summary of the visualization results for impellers RT3, i.e., analysis of the occurring gas bubble dispersion types, is presented in Table 8.

Table 8

Summary of visualization results (impeller RT3)—different types of gas bubble dispersion.

No Exp.	A	B	C	D	E	F
Gas flow rate, dm³·min⁻¹	10	30
Impeller speed, rpm	200	300	500	200	300	500
Type of dispersion	Accurate	Uniform	Uniform/excessive	Minimal	Excessive	Excessive

Tests carried out for impeller RT3 confirmed the high efficiency of gas bubble distribution in the volume of the tested object at a low inert gas flow rate of 10 dm³·min⁻¹. The most optimal variant was variant B (300 rpm, 10 dm³·min⁻¹). However, the other variants A and C (gas flow rate 10 dm³·min⁻¹) seemed to be favorable for this type of impeller and are recommended for further testing. The above process parameters will be analyzed in detail in a quantitative analysis to be performed on the basis of the obtained efficiency curves of the degassing process (oxygen removal). This analysis will give an unambiguous answer as to which process parameters are the most optimal for this type of impeller; the results are planned for publication in the next article.

It should also be noted here that the high agreement between the results of numerical calculations and physical modelling prompts a conclusion that the proposed approach to the simulation of a degassing process which consists of a single-phase flow model with a free surface and a particle flow model is appropriate. The simulation results enable us to understand how the velocity field in the fluid is formed and to analyze the distribution of gas bubbles in the system. The simulations in Flow-3D software can, therefore, be useful for both the design of the impeller geometry and the selection of process parameters.

4. Conclusions

The results of experiments carried out on the physical model of the device for the simulation of barbotage refining of aluminum revealed that the worst results in terms of distribution and dispersion of gas bubbles in the studied object were obtained for the black impellers variants B4, B8, and B12 (multi-orifice impellers—four, eight, and 12 outlet holes, respectively).

In this case, the control of flow, speed, and number of gas exit orifices did not improve the process efficiency, and the developed design did not meet the criteria for industrial tests. In the case of the ‘red triangle’ impeller (variant RT3), uniform gas bubble dispersion was achieved throughout the volume of the modeling fluid for most of the tested variants. The worst bubble dispersion results due to the occurrence of the so-called dead zones in the area near the bottom and sidewalls of the vessel were obtained for the flow variants of 20 dm³·min⁻¹ and 200 rpm and 30 dm³·min⁻¹ and 200 rpm. For the analyzed model, areas where swirls and gas bubble chains were formed were found only for the inert gas flow of 20 and 30 dm³·min⁻¹ and 200 rpm. The model impeller (variant RT3) had the best performance compared to the previously presented impellers in terms of dispersion of gas bubbles in the model liquid. Moreover, its design differed from previously presented models because of its sharp edges. This can be advantageous for gas bubble dispersion, but may negatively affect its suitability in industrial conditions due to premature wearing.

The CFD simulation results confirmed the results obtained from the experiments performed on the physical model. The numerical simulation of the operation of the ‘red triangle’ impeller model (using Flow-3D software) gave good agreement with the experiments performed on the physical model. This means that the presented model impeller, as compared to other (analyzed) designs, had the best performance in terms of gas bubble dispersion in the model liquid.

In further work, the developed numerical model is planned to be used for CFD simulations of the gas bubble distribution process taking into account physicochemical parameters of liquid aluminum based on industrial tests. Consequently, the obtained results may be implemented in production practice.

Funding Statement

This paper was created with the financial support grants from the AGH-UST, Faculty of Foundry Engineering, Poland (16.16.170.654 and 11/990/BK_22/0083) for the Faculty of Materials Engineering, Silesian University of Technology, Poland.

Author Contributions

Conceptualization, K.K. and D.K.; methodology, J.P. and T.M.; validation, M.S. and S.G.; formal analysis, D.K. and T.M.; investigation, J.P., K.K. and S.G.; resources, M.S., J.P. and K.K.; writing—original draft preparation, D.K. and T.M.; writing—review and editing, D.K. and T.M.; visualization, J.P., K.K. and S.G.; supervision, D.K.; funding acquisition, D.K. and T.M. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Zhang L., Xuewei L., Torgerson A.T., Long M. Removal of Impurity Elements from Molten Aluminium: A Review. Miner. Process. Extr. Metall. Rev. 2011;32:150–228. doi: 10.1080/08827508.2010.483396. [CrossRef] [Google Scholar]

2. Saternus M. Impurities of liquid aluminium-methods on their estimation and removal. Met. Form. 2015;23:115–132. [Google Scholar]

3. Żak P.L., Kalisz D., Lelito J., Gracz B., Szucki M., Suchy J.S. Modelling of non-metallic particle motion process in foundry alloys. Metalurgija. 2015;54:357–360. [Google Scholar]

4. Kalisz D., Kuglin K. Efficiency of aluminum oxide inclusions rmoval from liquid steel as a result of collisions and agglomeration on ceramic filters. Arch. Foundry Eng. 2020;20:43–48. [Google Scholar]

5. Kuglin K., Kalisz D. Evaluation of the usefulness of rotors for aluminium refining. IOP Conf. Ser. Mater. Sci. Eng. 2021;1178:012036. doi: 10.1088/1757-899X/1178/1/012036. [CrossRef] [Google Scholar]

6. Saternus M., Merder T. Physical modeling of the impeller construction impact o the aluminium refining process. Materials. 2022;15:575. doi: 10.3390/ma15020575. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

7. Saternus M., Merder T. Physical modelling of aluminum refining process conducted in batch reactor with rotary impeller. Metals. 2018;8:726. doi: 10.3390/met8090726. [CrossRef] [Google Scholar]

8. Saternus M., Merder T., Pieprzyca J. The influence of impeller geometry on the gas bubbles dispersion in uro-200 reactor—RTD curves. Arch. Metall. Mater. 2015;60:2887–2893. doi: 10.1515/amm-2015-0461. [CrossRef] [Google Scholar]

9. Hernández-Hernández M., Camacho-Martínez J., González-Rivera C., Ramírez-Argáez M.A. Impeller design assisted by physical modeling and pilot plant trials. J. Mater. Process. Technol. 2016;236:1–8. doi: 10.1016/j.jmatprotec.2016.04.031. [CrossRef] [Google Scholar]

10. Mancilla E., Cruz-Méndez W., Garduño I.E., González-Rivera C., Ramírez-Argáez M.A., Ascanio G. Comparison of the hydrodynamic performance of rotor-injector devices in a water physical model of an aluminum degassing ladle. Chem. Eng. Res. Des. 2017;118:158–169. doi: 10.1016/j.cherd.2016.11.031. [CrossRef] [Google Scholar]

11. Michalek K., Socha L., Gryc K., Tkadleckova M., Saternus M., Pieprzyca J., Merder T. Modelling of technological parameters of aluminium melt refining in the ladle by blowing of inert gas through the rotating impeller. Arch. Metall. Mater. 2018;63:987–992. [Google Scholar]

12. Walek J., Michalek K., Tkadlecková M., Saternus M. Modelling of Technological Parameters of Aluminium Melt Refining in the Ladle by Blowing of Inert Gas through the Rotating Impeller. Metals. 2021;11:284. doi: 10.3390/met11020284. [CrossRef] [Google Scholar]

13. Michalek K., Gryc K., Moravka J. Physical modelling of bath homogenization in argon stirred ladle. Metalurgija. 2009;48:215–218. [Google Scholar]

14. Michalek K. The Use of Physical Modeling and Numerical Optimization for Metallurgical Processes. VSB; Ostrawa, Czech Republic: 2001. [Google Scholar]

15. Chen J., Zhao J. Light Metals. TMS; Warrendale, PA, USA: 1995. Bubble distribution in a melt treatment water model; pp. 1227–1231. [Google Scholar]

16. Saternus M. Model Matematyczny do Sterowania Procesem Rafinacji Ciekłych Stopów Aluminium Przy Zastosowaniu URO-200. Katowice, Poland: 2004. Research Project Nr 7 T08B 019 21. [Google Scholar]

17. Pietrewicz L., Wężyk W. Urządzenia do rafinacji gazowej typu URO-200 sześć lat produkcji i doświadczeń; Proceedings of the Aluminum Conference; Zakopane, Poland. 12–16 October 1998. [Google Scholar]

18. Flow3d User’s Guide. Flow Science, Inc.; Santa Fe, NM, USA: 2020. [Google Scholar]

19. Sinelnikov V., Szucki M., Merder T., Pieprzyca J., Kalisz D. Physical and numerical modeling of the slag splashing process. Materials. 2021;14:2289. doi: 10.3390/ma14092289. [PMC free article] [PubMed] [CrossRef] [Google Scholar]

20. White F. Fluid Mechanics. McGraw-Hill; New York, NY, USA: 2010. (McGraw-Hill Series in Mechanical Engineering). [Google Scholar]

21. Yang Z., Yang L., Cheng T., Chen F., Zheng F., Wang S., Guo Y. Fluid Flow Characteristic of EAF Molten Steel with Different Bottom-Blowing Gas Flow Rate Distributions. ISIJ. 2020;60:1957–1967. doi: 10.2355/isijinternational.ISIJINT-2019-794. [CrossRef] [Google Scholar]

22. Nichols B.D., Hirt C.W. Methods for calculating multi-dimensional, transient free surface flows past bodies; Proceedings of the First International Conference on Numerical Ship Hydrodynamics; Gaithersburg, MD, USA. 20–22 October 1975. [Google Scholar]

23. Hirt C.W., Nichols B.D. Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries. J. Comput. Phys. 1981;39:201–255. doi: 10.1016/0021-9991(81)90145-5. [CrossRef] [Google Scholar]

24. Szucki M., Suchy J.S., Lelito J., Malinowski P., Sobczyk J. Application of the lattice Boltzmann method for simulation of the mold filling process in the casting industry. Heat Mass Transf. 2017;53:3421–3431. doi: 10.1007/s00231-017-2069-5. [CrossRef] [Google Scholar]

25. Themelis N.J., Goyal P. Gas injection in steelmaking. Candian Metall. Trans. 1983;22:313–320. [Google Scholar]

26. Zhang L., Jing X., Li Y., Xu Z., Cai K. Mathematical model of decarburization of ultralow carbon steel during RH treatment. J. Univ. Sci. Technol. Beijing. 1997;4:19–23. [Google Scholar]

27. Chiti F., Paglianti A., Bujalshi W. A mechanistic model to estimate powder consumption and mixing time in aluminium industries. Chem. Eng. Res. Des. 2004;82:1105–1111. doi: 10.1205/cerd.82.9.1105.44156. [CrossRef] [Google Scholar]

28. Bouaifi M., Roustan M. Power consumption, mixing time and homogenization energy in dual-impeller agitated gas-liquid reactors. Chem. Eng. Process. 2011;40:87–95. doi: 10.1016/S0255-2701(00)00128-8. [CrossRef] [Google Scholar]

29. Kang J., Lee C.H., Haam S., Koo K.K., Kim W.S. Studies on the overall oxygen transfer rate and mixing time in pilot-scale surface aeration vessel. Environ. Technol. 2001;22:1055–1068. doi: 10.1080/09593332208618215. [PubMed] [CrossRef] [Google Scholar]

30. Moucha T., Linek V., Prokopov E. Gas hold-up, mixing time and gas-liquid volumetric mass transfer coefficient of various multiple-impeller configurations: Rushton turbine, pitched blade and techmix impeller and their combinations. Chem. Eng. Sci. 2003;58:1839–1846. doi: 10.1016/S0009-2509(02)00682-6. [CrossRef] [Google Scholar]

31. Szekely J. Flow phenomena, mixing and mass transfer in argon-stirred ladles. Ironmak. Steelmak. 1979;6:285–293. [Google Scholar]

32. Iguchi M., Nakamura K., Tsujino R. Mixing time and fluid flow phenomena in liquids of varying kinematic viscosities agitated by bottom gas injection. Metall. Mat. Trans. 1998;29:569–575. doi: 10.1007/s11663-998-0091-1. [CrossRef] [Google Scholar]

33. Hjelle O., Engh T.A., Rasch B. Removal of Sodium from Aluminiummagnesium Alloys by Purging with Cl₂. Aluminium-Verlag GmbH; Dusseldorf, Germany: 1985. pp. 343–360. [Google Scholar]

34. Zhang L., Taniguchi S. Fundamentals of inclusion removal from liquid steel by bubble flotation. Int. Mat. Rev. 2000;45:59–82. doi: 10.1179/095066000101528313. [CrossRef] [Google Scholar]

원문 보기

Figura 7. Influencia del modelo de turbulencia. Qmodelo=27.95l/s.

Flow-3D를 사용하여 전산유체역학(CFD)을 적용한 빠른 단계의 플러시 유동 수치 모델링

Numerical Modeling of Flush Flow in a Rapid Step Applying Computational Fluid Dynamics (CFD) Using Flow-3D.

레브 폴리텍. (Quito) [온라인]. 2018, vol.41, n.2, pp.53-64. ISSN 2477-8990.

이 프로젝트의 주요 목표는 FLOW-3D를 사용하여 계단식 여수로에서 스키밍 흐름의 수치 모델링을 개발하는 것입니다. 이러한 구조의 설계는 물리적 모델링에서 얻은 경험적 표현과 CFD 코드를 지원하는 계단식 여수로를 통한 흐름의 수치 모델링에서 보완 연구를 기반으로 합니다. 수치 모델은 균일한 영역의 유속과 계단 여수로의 마찰 계수를 추정하는 데 사용됩니다(ϴ = 45º, Hd=4.61m). 흐름에 대한 자동 통기의 표현은 복잡하므로 프로그램은 공기 연행 모델을 사용하여 특정 제한이 있는 솔루션에 근접합니다.

The main objective of this project is to develop the numerical modeling of the skimming flow in a stepped spillway using FLOW-3D. The design of these structures is based on the use of empirical expressions obtained from physical modeling and complementary studies in the numerical modeling of flow over the stepped spillway with support of CFD code. The numerical model is used to estimate the flow velocity in the uniform region and the friction coefficient of the stepped spillway (ϴ = 45º, Hd=4.61m). The representation of auto aeration a flow is complex, so the program approximates the solution with certain limitations, using an air entrainment model; drift flux model and turbulence model k-ԑ RNG. The results obtained with numerical modeling and physical modeling at the beginning of natural auto aeration of flow and depth of the biphasic flow in the uniform region presents deviations above to 10% perhaps the flow is highly turbulent.

Keywords : Stepped spillway; skimming flow; air entrainment; drift flux; numerical modeling; FLOW-3D.

Keywords : 계단식 여수로; 스키밍 흐름; 공기 연행; 드리프트 플럭스; 수치 모델링; 흐름-3D.·

스페인어로 된 초록 · 스페인어 로 된 텍스트 · 스페인어로 된 텍스트( pdf )

Figure 1. Grazing flow over a rapid step.

Figura 2. Principales regiones existentes en un flujo rasante.

Figure 3. Dimensions of the El Batán stepped rapid.

Figure 4. 3D physical model of the El Batán stepped rapid

REFERENCIAS

ARAGUA. (2013). “Modelación numérica y experimental de flujos aire-agua
en caídas en colectores.”, Laboratório Nacional de Engenharia Civil, I.
P. Av do Brasil 101 • 1700-066 Lisboa.
Bombardelli, F.A., Meireles, I. and Matos, J., (2010), “Laboratory
measurement and multi-block numerical simulations of the mean flow
and turbulence in the non-aerated skimming flow region of steep stepped
spillways”, Environ Fluid Mechanics.
Castro M. (2015) “Análisis Dimensional y Modelación física en Hidráulica”.
Escuela Politécnica Nacional. Quito Ecuador. 50 p.
Chanson H., D. B. Bung., J. Matos (2015). “Stepped spillways and cascades”.
IAHR Monograph. School of Civil Engineering, University of
Queensland, Brisbane, Australia.
Chanson H. (1993). “Stepped Spillway Flows and Air Entrainment.” Can. Jl
of Civil Eng., Vol. 20, No. 3, June, pp. 422-435 (ISSN 0315-1468).
CIERHI, EPN TECH, (2016). “Estudio experimental en modelo físico de las
rápidas con perfil escalonado y liso de la quebrada el Batán Fase I y Fase
II”, Escuela Politécnica Nacional, Quito Ecuador.
Fernández Oro J. M. (2012)., “Técnicas Numéricas en Ingeniería de Fluidos:
Introducción a la Dinámica de Fluidos Computacional (CFD) por el
Método de Volúmenes Finitos”. Barcelona: Reverté.
Flow Science, Inc. (2012). “FLOW 3D 10.1.0 Documentation Release.
Manual de Usuario”, Los Alamos National Laboratory. Santa Fe, New
México
Khatsuria, R.M., (2005)., “Hydraulics of Spillways and Energy Dissipators”.
Department of Civil and Environmental Engineering Georgia Institute
of Technology Atlanta.
Lucio I., Matos J., Meireles I. (2015). “Stepped spillway flow over small
embankment dams: some computational experiments”. 15th FLOW-3D
European users conference.
Mohammad S., Jalal A. and Michael P., (2012). “Numerical Computation of
Inception Point Location for Steeply Sloping Stepped Spillways” 9th
International Congress on Civil Engineering. Isfahan University of
Technology (IUT), Isfahan, Iran
Pfister M., Chanson H., (2013), “Scale Effects in Modelling Two-phase Airwater Flows”, Proceedings of 2013 IAHR World Congress.
Sarfaraz, M. and Attari, J. (2011), “Numerical Simulation of Uniform Flow
Region over a Steeply Sloping Stepped Spillway”, 6th National
Congress on Civil Engineering, Semnan University, Semnan, Iran.
Valero, D., Bung, D., (2015), “Hybrid investigation of air transport processes
in moderately sloped stepped spillway flows”, E-proceedings of the 36th
IAHR World Congress 28 June – 3 July, 2015, The Hague, the Netherlands.

원문보기

Serife Yurdagul Kumcu−2−KSCE Journal of Civil Engineeringthe use of CFD for the assessment of a design, as well as screeningand optimizing of hydraulic structures and cofferdam layouts. Theyconclude that CFD has been successful in optimizing the finalconceptual configuration for the hydraulics design of the project,but recommend that physical modeling still be used as a finalconfirmation.This paper provides experimental studies performed on Kav akDam and analyses the stability of spillway design by usingFLOW-3D model. It compares the hydraulic model tests withFLOW-3D simulation results and gives information on howaccurately a commercially available Computational Fluid Dynamic(CFD) model can predict the spillway discharge capacity andpressure distribution along the spillway bottom surface. 2. Physical ModelA 1/50-scaled undistorted physical model of the Kavsak Damspillway and stilling basin was built and tested at the HydraulicModel Laboratory of State Hydraulic Works of Turkey (DSI).The model was constructed of plexiglas and was fabricated toconform to the distinctive shape of an ogee crest. The spillwayhas 45.8 m in width and 57 m long with a bottom slope of 125%.The length of the stilling basin is about 90 m. During model tests,flow velocities were measured with an ultrasonic flow meter.Pressures on the spillway were measured using a piezometerssçTable 1. Upstream and Downstream Operating Conditions of theKavsak DamRun Upstream reservoir elevation (m)Downstream tailwater elevation (m)1 306.55 168.002 311.35 174.503 314.00 178.904 316.50 182.55Fig. 1. (a) Original Project Design and Final Project Design after Experimental Investigations and Flow Measurement Sections at theApproach, (b) Top View Experimentally Modified Approach in the Laboratory, (c) Side View of the Experimentally Modified Approachin the Laboratory

Investigation of flow over spillway modeling and comparison between experimental data and CFD analysis

여수로 모델링 및 실험 데이터와 CFD 해석의 비교에 대한 조사

June 2016
KSCE Journal of Civil Engineering 21(3)

DOI:10.1007/s12205-016-1257-z

Authors:

Serife Yurdagul Kumcu

Necmettin Erbakan Üniversitesi

Abstract and Figures

As a part of design process for hydro-electric generating stations, hydraulic engineers typically conduct some form of model testing. The desired outcome from the testing can vary considerably depending on the specific situation, but often characteristics such as velocity patterns, discharge rating curves, water surface profiles, and pressures at various locations are measured. Due to recent advances in computational power and numerical techniques, it is now also possible to obtain much of this information through numerical modeling. In this paper, hydraulic characteristics of Kavsak Dam and Hydroelectric Power Plant (HEPP), which are under construction and built for producing energy in Turkey, were investigated experimentally by physical model studies. The 1/50-scaled physical model was used in conducting experiments. Flow depth, discharge and pressure data were recorded for different flow conditions. Serious modification was made on the original project with the experimental study. In order to evaluate the capability of the computational fluid dynamics on modeling spillway flow a comparative study was made by using results obtained from physical modeling and Computational Fluid Dynamics (CFD) simulation. A commercially available CFD program, which solves the Reynolds-averaged Navier-Stokes (RANS) equations, was used to model the numerical model setup by defining cells where the flow is partially or completely restricted in the computational space. Discharge rating curves, velocity patterns and pressures were used to compare the results of the physical model and the numerical model. It was shown that there is reasonably good agreement between the physical and numerical models in flow characteristics.

수력 발전소 설계 프로세스의 일부로 수력 엔지니어는 일반적으로 어떤 형태의 모델 테스트를 수행합니다. 테스트에서 원하는 결과는 특정 상황에 따라 상당히 다를 수 있지만 속도 패턴, 방전 등급 곡선, 수면 프로파일 및 다양한 위치에서의 압력과 같은 특성이 측정되는 경우가 많습니다. 최근 계산 능력과 수치 기법의 발전으로 인해 이제는 수치 모델링을 통해 이러한 정보의 대부분을 얻을 수도 있습니다.

본 논문에서는 터키에서 에너지 생산을 위해 건설 중인 Kavsak 댐과 수력발전소(HEPP)의 수력학적 특성을 물리적 모델 연구를 통해 실험적으로 조사하였다. 1/50 스케일의 물리적 모델이 실험 수행에 사용되었습니다. 다양한 흐름 조건에 대해 흐름 깊이, 배출 및 압력 데이터가 기록되었습니다. 실험 연구를 통해 원래 프로젝트에 대대적인 수정이 이루어졌습니다.

배수로 흐름 모델링에 대한 전산유체역학의 능력을 평가하기 위해 물리적 모델링과 전산유체역학(CFD) 시뮬레이션 결과를 이용하여 비교 연구를 수행하였습니다. RANS(Reynolds-averaged Navier-Stokes) 방정식을 푸는 상업적으로 이용 가능한 CFD 프로그램은 흐름이 계산 공간에서 부분적으로 또는 완전히 제한되는 셀을 정의하여 수치 모델 설정을 모델링하는 데 사용되었습니다.

물리적 모델과 수치 모델의 결과를 비교하기 위해 배출 등급 곡선, 속도 패턴 및 압력을 사용했습니다. 유동 특성에서 물리적 모델과 수치 모델 간에 상당히 좋은 일치가 있는 것으로 나타났습니다.

References

Bureau of Reclamation (1977). Design of small dams, U.S. Government Printing Office, Washington, D.C., U.S.

Bureau of Reclamation (1990). Cavitation in chute and spillways, Engineering Monograph, No.42, U.S. Chanel, P. G. (2008). An evaluation of computational fluid dynamics for

spillway modeling, MSc Thesis, University of Manitoba Winnipeg, Manitoba, Canada.

Chanson, H. (2002). The hydraulics of stepped chutes and spillways,Balkema, Lisse, The Netherlands.

Chanson, H. and Gonzalez, C. A. (2005). “Physical modeling and scale effects of air-water flows on stepped spillways.” Journal of Zhejiang University Science, Vol. 6A, No. 3, pp. 243-250.

Demiroz, E. (1986). “Specifications of aeration structures which are added to the spillways.” DSI Report, HI-754, DSI-TAKK Publications, Ankara, Turkey.

Erfanain-Azmoudeh, M. H. and Kamanbedast, A. A. (2013). “Determine the appropriate location of aerator system on gotvandoliadam’s spillway using Flow 3D.” American-Eurasian J. Agric. & Environ. Sci., Vol. 13, No. 3, pp. 378-383, DOI: 10.5829/idosi.aejaes.2013. 13.03. 458.

Falvey, H. T. (1990). Cavitation in chutes and spillways, Engineering Monograph 42 Water Resources Technical Publication US Printing Office, Bureau of Reclamation, Denver.

Flow-3D User ’s Manual (2012). Flow science, Inc., Santa Fe, N.M.

Hirt, C. W. (1992). “Volume-fraction techniques: Powerful tools for flow

modeling.” Flow Science Report, No. FSI-92-00-02, Flow Science, Inc., Santa Fe, N.M.

Hirt C. W. and Nichols B. D. (1981). “Volume of Fluid (VOF) method for the dynamics of free boundaries.”Jornal of Computational Physics, Vol. 39, pp. 201-225, DOI: 10.1016/0021-9991(81)90145-5.

Hirt, C. W. and Sicilian, J. M. (1985). “A Porosity technique for the definition of obstacles in rectangular cell meshes.” Proceedings of the 4th International Conference on Ship Hydro-dynamics, 24-27 September 1985, National Academic of Sciences, Washington DC.

Ho, D., Boyes, K., Donohoo, S., and Cooper, B. (2003). “Numerical flow analysis for spillways.” 43rd ANCOLD Conference, Hobart, Tas m a nia .

Johnson, M. C. and Savage, B. M. (2006). “Physical and numerical comparison of flow over ogee spillway in the presence of tailwater.”

Journal of Hydraulic Engineering, Vol. 132, No. 12, pp. 1353-135, DOI: 10.1061/(ASCE)0733-9429.

Kim, S. D., Lee, H. J., and An, S. D. (2010). “Improvement of hydraulic stability for spillway using CFD model.” Int. Journal of the Physical Sciences, Vol. 5, No. 6, pp. 774-780.

Kokpinar, M. A. and Gogus, M. (2002). “High speed jet flows over spillway aerators.” Canadian Journal of Civil Engineering, Vol. 29, No. 6, pp. 885-898, DOI: 10.1139/l02-088.

Kumcu, S. Y. (2010). Hydraulic model studies of Kavsak Dam and HEPP, DSI Report, HI-1005, DSI-TAKK Publications, Ankara, Turkey.

Margeirsson, B. (2007). Computational modeling of flow over a spillway, MSc Thesis, Chalmers University of Technology, Gothenburg, Sweden.

Nichols, B. D. and Hirt, C. W. (1975). “Methods for calculating multi-dimensional, transient free surface flows past bodies.” Proc. First Intern. Conf. Num., Ship Hydrodynamics, Gaithersburg, ML.

Savage, B. M. and Johnson, M. C. (2001). “Flow over ogee spillway: Physical and numerical model case study.” Journal of Hydraulic Engineering, ASCE, Vol. 127, No. 8, pp. 640-649, DOI: 10.1061/(ASCE)0733-9429.

Souders, D. T. and Hirt, C. W. (2004). “Modeling entrainment of air at turbulent free surfaces.” Critical Transitions in Water and Environmental resources Management, pp. 1-10.

entürk, F. (1994). Hydraulics of dams and reservoirs, Water Resources Publication Colorado, USA.

Teklemariam, E., Korbaylo, B, Groeneveld, J., Sydor, K., and Fuchs, D. (2001). Optimization of hydraulic design using computational fluid dynamics, Waterpower XII, Salt Lake City, Utah.

Teklemariam, E., Shumilak, B., Sydor, K., Murray, D., Fuchs, D., and Holder, G. (2008). “An integral approach using both physical and computational modeling can be beneficial in addressing the full range of hydraulic design issues.” CDA Annual Conference, Winnipeg, Canada.

Usta, E. (2014). Numerical investigation of hydraulic characteristics of Laleli Dam spillway and comparison with physical model study, Master Thesis, Middle East Technical University, Ankara, Turkey.

Versteeg, H. K. and Malalasekera, W. (1996). An introduction to computational fluid dynamics, Longman Scientific and Technical, Longman Group Limited, Harlow, England.

Vischer, D. L. and Hager, W. H. (1997). Dam hydraulics, J. Wiley & Sons Ltd., England.

Wagner, W. E. (1967). “Glen Canyon diversion tunnel outlets.” J. Hydraulic Division, ASCE, Vol. 93, No. HY6, pp. 113-134.

Willey, J., Ewing, T., Wark, B., and Lesleighter, E. (2012). Comple-mentary use of physical and numerical modeling techniques in spillway design refinement, Commission Internationale Des Grands Barrages, Kyoto, June 2012.

원문 보기

Assessment of spillway modeling using computational fluid dynamics

전산유체역학을 이용한 여수로 모델링 평가

Authors: Paul G. Chanel and John C. Doering AUTHORS INFO & AFFILIATIONS

Publication: Canadian Journal of Civil Engineering

3 December 2008

Abstract

Throughout the design and planning period for future hydroelectric generating stations, hydraulic engineers are increasingly integrating computational fluid dynamics (CFD) into the process. As a result, hydraulic engineers are interested in the reliability of CFD software to provide accurate flow data for a wide range of structures, including a variety of different spillways. In the literature, CFD results have generally been in agreement with physical model experimental data. Despite past success, there has not been a comprehensive assessment that looks at the ability of CFD to model a range of different spillway configurations, including flows with various gate openings. In this article, Flow-3D is used to model the discharge over ogee-crested spillways. The numerical model results are compared with physical model studies for three case study evaluations. The comparison indicates that the accuracy of Flow-3D is related to the parameter P/H_d.

미래의 수력 발전소를 위한 설계 및 계획 기간 동안 유압 엔지니어는 전산유체역학(CFD)을 프로세스에 점점 더 많이 통합하고 있습니다. 결과적으로 유압 엔지니어는 다양한 여수로를 포함하여 광범위한 구조에 대한 정확한 흐름 데이터를 제공하는 CFD 소프트웨어의 신뢰성에 관심을 갖고 있습니다. 문헌에서 CFD 결과는 일반적으로 물리적 모델 실험 데이터와 일치했습니다. 과거의 성공에도 불구하고 다양한 게이트 개구부가 있는 흐름을 포함하여 다양한 여수로 구성을 모델링하는 CFD의 기능을 살펴보는 포괄적인 평가는 없었습니다. 이 기사에서는 Flow-3D를 사용하여 ogee-crested 방수로의 배출을 모델링합니다. 세 가지 사례 연구 평가를 위해 수치 모델 결과를 물리적 모델 연구와 비교합니다. 비교는 Flow-3D의 정확도가 매개변수 P/Hd와 관련되어 있음을 나타냅니다.

Résumé

Les ingénieurs en hydraulique intègrent de plus en plus la dynamique des fluides numérique (« CFD ») dans le processus de conception et de planification des futures centrales. Ainsi, les ingénieurs en hydraulique s’intéressent à la fiabilité du logiciel de « CFD » afin de fournir des données précises sur le débit pour une large gamme de structures, incluant différents types d’évacuateurs. Les résultats de « CFD » dans la littérature ont été globalement sont généralement en accord avec les données expérimentales des essais physiques. Malgré les succès antérieurs, il n’y avait aucune évaluation complète de la capacité des « CFD » à modéliser une plage de configuration des évacuateurs, incluant les débits à diverses ouvertures de vannes. Dans le présent article, le logiciel Flow-3D est utilisé pour modéliser le débit par des évacuateurs en doucine. Les résultats du modèle de calcul sont comparés à ceux des essais physiques pour trois études de cas. La comparaison montre que la précision du logiciel Flow-3D est associée au paramètre P/H_d.

Get full access to this article

View all available purchase options and get full access to this article.

GET ACCESSALREADY A SUBSCRIBER? SIGN IN AS AN INDIVIDUAL OR VIA YOUR INSTITUTION

References

Chanel, P.G., and Doering, J.C. 2007. An evaluation of computational fluid dynamics for spillway modelling. In Proceedings of the 16th Australasian Fluid Mechanics Conference (AFMC), Gold Coast, Queensland, Australia, 3–7 December 2007. pp. 1201–1206.

Flow Science, Inc. 2007. Flow-3D user’s manuals. Version 9.2. Flow Science, Inc., Santa Fe, N.M.

Gessler, D. 2005. CFD modeling of spillway performance, EWRI 2005: Impacts of global climate change. In Proceedings of the World Water and Environmental Resources Congress, Anchorage, Alaska, 15–19 May 2005. Edited by R. Walton. American Society of Civil Engineers, Reston, Va.

Hirt, C.W., and Nichols, B.D. 1981. Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of Computational Physics, 39(1): 201–225.

Crossref

ISI

Hirt, C.W., and Sicilian, J.M. 1985. A porosity technique for the definition of obstacles in rectangular cell meshes. In Proceedings of the 4th International Conference on Ship Hydro-dynamics, Washington, D.C., 24–27 September 1985. National Academy of Sciences, Washington, D.C.

Ho, D., Cooper, B., Riddette, K., and Donohoo, S. 2006. Application of numerical modelling to spillways in Australia. In Dams and Reservoirs, Societies and Environment in the 21st Century. Edited by Berga et al. Taylor and Francis Group, London.

LaSalle Consulting Group Inc. 1992. Conawapa generating station. Sectional model study of the spillway. LaSalle Consulting Group Inc., Montréal, Que.

Lemke, D.E. 1989. A comparison of the hydraulic performance of an orifice and an overflow spillway in a northern application using physical modeling. M.Sc. thesis, University of Manitoba, Winnipeg, Man.

Savage, B.M., and Johnson, M.C. 2001. Flow over ogee spillway: Physical and numerical model case study. Journal of Hydraulic Engineering, 127(8): 640–649.

Crossref

ISI

Teklemariam, E., Korbaylo, B., Groeneveld, J., Sydor, K., and Fuchs, D. 2001. Optimization of hydraulic design using computational fluid dynamics. In Proceedings of Waterpower XII, Salt Lake City, Utah, 9–11 July 2001.

Teklemariam, E., Korbaylo, B., Groeneveld, J., and Fuchs, D. 2002. Computational fluid dynamics: Diverse applications in hydropower project’s design and analysis. In Proceedings of the CWRA 55th Annual Conference, Winnipeg, Man., 11–14 June 2002. Canadian Water Resources Association, Cambridge, Ontario.

Western Canadian Hydraulic Laboratories Inc. 1980. Hydraulics model studies limestone generating station spillway/diversion structure flume study. Final report. Western Canadian Hydraulic Laboratories Inc., Port Coquitlam, B.C.