flow
¶
Note
Either boun
or flow
are required for a flow problem.
Flow data. Source and sink parameters are input and may be used to apply boundary conditions. Note that the alternative definitions for isothermal models apply when flow is used in conjunction with control statement **airwater **(See Control statement airwater or air (optional)).
Input Variable | Format | Default | Description |
---|---|---|---|
Non-Isothermal model | |||
SKD | real | Heat and mass source strength (kg/s), heat only (MJ/s). Negative value indicates injection into the rock mass. | |
EFLOW | real | Enthalpy of fluid injected (MJ/kg). If the fluid is flowing from the rock mass, then the in-place enthalpy is used. | |
If EFLOW < 0, then ABS(EFLOW) is interpreted as a temperature (oC) and the enthalpy (assuming water only) calculated accordingly. In heat only problems with EFLOW < 0, the node is in contact with a large heat pipe that supplies heat to the node through an impedance AIPED so as to maintain its temperature near ABS (EFLOW). Large values (approximately 1000) of AIPED are recommended. | |||
AIPED | real | Impedance parameter. If AIPED is nonzero, the code interprets SKD as a flowing wellbore pressure (MPa) with an impedance ABS(AIPED). If AIPED < 0, flow is only allowed out of the well. For heat only, AIPED is the thermal resistance. If AIPED = 0, SKD is flow rate. If AIPED ≠ 0 and SKD = 0 the initial value of pressure will be used for the flowing pressure. |
Note
If the porosity of the node is zero, then there is only a temperature solution, and the code forms a source proportional to the enthalpy difference. The source term is given by \(Q = AIPED \cdot (E - EFLOW)\), where E is the in-place enthalpy and EFLOW is a specified enthalpy.
Isothermal model | |||
Case 1: AIPED = 0 (Constant Mass Rate, 1- or 2-Phase Source or Sink) | |||
SKD | real | Mass source strength (kg/s). Negative value indicates injection into the rock mass. | |
EFLOW | real | a) EFLOW ≥ 0, EFLOW is the source liquid saturation,
* \(Q_w = SKD \cdot EFLOW\) (kg/s)
* \(Q_a = SKD \cdot (1 - EFLOW)\) (kg/s)
|
|
b) EFLOW < 0, ABS(EFLOW) is the source air pressure (MPa)
* \(Q_w = SKD\) (kg/s)
* \(Q_a = 1.0 \cdot (P_a - ABS(EEFLOW))\) (kg/s)
|
|||
AIPED | Used only as flag to indicatee Constant Mass Rate |
In the above and following relations, \(Q_w\) is the source term for water, \(Q_a\) is the source term for air, and \(P_a\) is the in-place air pressure. The second case works well in situations where inflow is specified and it is desired to hold the air pressure at a constant value.
Case 2: AIPED > 0 (Constant Pressure, Constant Liquid Saturation Source or Sink) | |||
SKD | real | Specified source pressure (MPa). | |
EFLOW | real | a) EFLOW < 0, air only source.
* \(Q_a = AIPED \cdot (P_a - SKD)\) (kg/s)
b) 0 < EFLOW < 1, EFLOW is specified source liquid saturation,
for SKD ≥ 0, 2-phase source,
* \(Q_a = AIPED \cdot (P_a - SKD)\) (kg/s)
* \(Q_w = AIPED \cdot (S_l - EEFLOW) \cdot P_0\) (kg/s)
when SKD < 0, water only source \(Q_a = 0\).
c) EFLOW = 1, water only source.
* \(Q_w = AIPED \cdot (P_l - SKD)\) kg/s
|
|
AIPED | real | Impedance parameter. A large value is recommended (102 - 106) in order to create a flow
term large enough to maintain constant pressure.
|
Case 3: AIPED < 0 (Outflow only, if \(P_l > SKD\)) | |||
SKD | real | Pressure above which outflow occurs (MPa) | |
EFLOW | real | Not used. | |
AIPED | real | Impedance parameter. | |
\(Q_w = ABS(AIPED) \cdot R_l / \mu_l(P_l - SKD)\) (kg/s) |
where \(R_l\) is the water relative peremeability and \(\mu_l\) is the water viscosity.
PARAMETER | SKD | EFLOW | |||
Physics Model | AIPED* | \(\ge 0\) | \(< 0\) | \(\ge 0\) | \(< 0\) |
Non-isothermal | \(> 0\) | Flowing wellbore pressure (MPa) | Flowing wellbore pressure (MPa) - injection into rock mass | Enthalpy (MJ/kg) | Temperature (C) |
\(= 0\) | Flow rate (kg/s)
Heat only (MJ/s)
|
Flow rate (kg/s)
Heat only (MJ/s) - injection into rock mass
|
|||
\(< 0\) | Flowing wellbore pressure (MPa) - outflow only | N/A | |||
Isothermal | \(> 0\) | Specified source pressure (MPa) | Specified source pressure (MPa) - Water only source | Source liquid saturation | Air only source (used as flag) |
\(= 0\) | Mass source strength (kg/s) | Mass source strength (kg/s) - injection into rock mass | Source liquid saturation | Source air pressure (MPa) | |
\(< 0\) | Pressure above which outflow occurs (MPa) | N/A | N/A | N/A | |
* Impedance parametere |
The following are examples of flow. In the first example, at node 88, a mass flow of 0.05 kg/s at 25 C is being withdrawn from the model. Because fluid is being withdrawn the in-place temperature will actually be used. For every 14th node from node 14 to 140, the pressure and temperature are being held constant at 3.6 MPa and 160 C, respectively. This represents a constant temperature recharge along one of the problem boundaries.
flow | ||||
88 | 88 | 1 | 0.050 | -25.0 |
14 | 140 | 14 | 3.600 | -160.0 |
In the second example, the corresponding input for airwater, is included, indicating an isothermal air-water two-phase simulation is being run with a reference temperature for property evaluation of 20 C and a reference pressure of 0.1 MPa. At nodes 26 and 52, water saturation is 100% and water is being injected at 2.e-3 kg/s. At nodes 1 and 27, there is an air only source, with a specified pressure of 0.1 MPa, and the air is being injected at the rate of 100*(Pa - 0.1) kg/s.
airwater | ||||
2 | ||||
20.0 | 0.1 | |||
flow | ||||
26 | 52 | 26 | -2.e-3 | 1.0 |
1 | 27 | 26 | 0.1 | -0.2 |
© Copyright 2018, Los Alamos National Laboratory