DATA SET 13: ELEMENT (DISTRIBUTED) AND WELL (POINT) SOURCE/SINK

(1) Element Source/Sink:

This data set is needed if and only if NSEL > 0. When NSEL > 0, data subsets (b) and (c), together and in sequence as a supergroup, should be repeated (NOHA + NOKX + NOMB + 2) times for simulations without heat transfer or (NOHA + NOKX + NOMB + 3) times for simulations with heat transfer (i.e., ISHEAT > 0). (NSEL is specified in Data Set 12. NOHA, NOKX, and NOMB are specified in Data Set 10).

  1. Global element number of compressed source/sink element number: Free-field format input contains the following information:

LES(I), I = 1, NSEL, Global element number of the I-th distributed source/sink element.

  1. Element source/sink profile: Number of lines depends on NSPR, NSDP, and NON. Each line contains a number of data points = 2*NSDP numerical numbers. Each line is free-field format input. For a transport component, fluid density, flow rate, and temperature, the sequence of input goes by a (time, source strength) pair that is repeated until all the NSDP pairs are entered. This sequence of input is then repeated for the next component, flow rate, fluid density, and/or temperature.

For K = 1, NOHA + NOKX + NOMB + 3, I = 1, NSPR, J = 1, NSDP,

  1. TSOSF(J,I,K) = Time of J-th data point in I-th profile for K-th component (T).
  2. SOSF(J,I,K) = Source/sink value of J-th data point in I-th profile for K-th component:
  • concentration (moles/dm3 liquid) if K <= (NOHA + NOKX + NOMB),
  • flow rate of liquid (L**3/T) if K = (NOHA + NOKX + NOMB + 1),
  • density of incoming liquid (kg/dm3) if K = (NOHA + NOKX + NOMB + 2), and
  • temperature of incoming liquid (oK) if ISHEAT > 0 and K = (NOHA + NOKX + NOMB + 3).
  1. Source/sink profile type in each source/sink element: Typically, one line per element is needed. However, automatic generation can be used. Each line is free-field format input containing the following:
  1. MI = Compressed source/sink element number of first element in the sequence.
  2. NSEQ = NSEQ source/sink elements will have profile type MITYP.
  3. MAD = Increment of MI for each of the NSEQ elements.
  4. MITYP= Source/sink profile type for element MI.
  5. MTYPAD = Increment of source profile type in each of the NSEQ subsequent elements.

NOTE: A line containing five zeros separated by spaces or commas must be used to end this data set.

(2) Point source/sink:

This data set needed if and only if NWNP > 0. When NWNP > 0, data subsets (b) and (c), together and in sequence as a supergroup, should be repeated (NOHA + NOKX + NOMB + 2) times for simulations without heat transfer or (NOHA + NOKX + NOMB + 3) times for simulations with heat transfer (i.e., ISHEAT > 0). (NWNP is specified in Data Set 12. NOHA, NOKX, and NOMB are specified in Data Set 10).

  1. Global node number of compressed point source/sink number: Free-field format input contains the following information:

NPW(I), I = 1, NWNP, Global node number of the first source/sink node.

  1. Point source/sink profile: Number of lines depends on NWPR, NWDP, and NOHA + NOKX + NOMB. Each line contains a number of data points = 2*NWDP numerical numbers. Each line is free-field format input. For a transport component, fluid density, flow rate, and temperature, the sequence of input goes by a (time, source strength) pair that is repeated until all the NWDP pairs are entered. This sequence of input is then repeated for the next component, flow rate, fluid density, and/or temperature.

For K =1, NOHA + NOKX + NOMB + 3, I = 1, NWPR, J = 1, NWDP,

  1. TWSSF(J,I,K) = Time of J-th data point in I-th profile for K-th component (T).
  2. WSSF(J,I,K) = Source/sink value of J-th data point in I-th profile for K-th component:
  1. Type of point source/sink nodes: Typically, one line per well node is needed. However, automatic generation may be used. Each line is free-field format input containing the following:
  1. NI = Compressed well node number of the first node in a sequence.
  2. NSEQ = NSEQ subsequent well nodes will be generated automatically for the source profile type.
  3. NIAD = Increment of compressed well node number for each of the NSEQ subsequent nodes.
  4. NITYP = Type of well source/sink profile assigned to NI-th well node.
  5. NITYPA= Increment of NITYP for each of the NSEQ subsequent nodes.

NOTE: A line containing five zeros separated by spaces or commas must be used to signal the end of this data subset.

v1.0 - 1/13/1999

Jin-Ping Gwo, email: jgwo@umbc.edu