DATA SET 16: VELOCITY, MOISTURE CONTENT, AND LIQUID DENSITY

The first two data subsets (velocity and moisture content) are needed only if KVI < 0. If KVI > 0, the velocity and moisture data will be read via Logical Unit 11. (KVI is specified in Data Set 2).

The format for these two data subsets differs for KVI = -1 (steady-state velocities and moisture contents) and KVI = -2 (time series profile for velocity and moisture content).

I. KVI = -1: Steady-state hydrologic field

This option allows the input of velocity and moisture content values that are different for all nodes and elements, but always remain the same for the entire duration of the simulation.

  1. Velocity field: Normally, one line per node is needed. However, automatic generation can be used if there is a pattern. Format varies slightly among one-, two-, and three-dimensional problems in this group of data.

One-Dimensional Problems:

Each line is free-field format input containing five variables:

  1. NI = Node number of the first node in the sequence.
  2. NSEQ = NSEQ subsequent nodes will be automatically generated.
  3. NAD = Increment of node number in each of the NSEQ subsequent nodes.
  4. VXNI = x-velocity component of node NI (L/T).
  5. VXAD = Increment of x-velocity for each of the NSEQ subsequent nodes (L/T).

NOTE: A line with five zeros separated by spaces or commas signals the end of this data subset.

Two-Dimensional Problems:

Each line is free-field format input containing seven variables:

  1. NI = Node number of the first node in the sequence.
  2. NSEQ = NSEQ subsequent nodes will be automatically generated.
  3. NAD = Increment of node number in each of the NSEQ subsequent nodes.
  4. VXNI = x-velocity component of node NI (L/T).
  5. VZNI = z-velocity component of node NI (L/T).
  6. VXAD = Increment of x-velocity for each of the NSEQ subsequent nodes (L/T).
  7. VZAD = Increment of z-velocity for each of the NSEQ subsequent nodes (L/T).

NOTE: A line with seven zeros separated by spaces or commas signals the end of this data subset.

Three-Dimensional Problems:

Each line is free-field format input containing nine variables:

  1. NI = Node number of the first node in the sequence.
  2. NSEQ = NSEQ subsequent nodes will be automatically generated.
  3. NAD = Increment of node number in each of the NSEQ subsequent nodes.
  4. VXNI = x-velocity component of node NI (L/T).
  5. VYNI = y-velocity component of node NI (L/T).
  6. VZNI = z-velocity component of node NI (L/T).
  7. VXAD = Increment of x-velocity for each of the NSEQ subsequent nodes (L/T).
  8. VYAD = Increment of y-velocity for each of the NSEQ subsequent nodes (L/T).
  9. VZAD = Increment of z-velocity for each of the NSEQ subsequent nodes (L/T).

NOTE: A line with seven zeros separated by spaces or commas signals the end of this data subset.

  1. Moisture content field: Typically, one line per element is needed. However, automatic generation can be used if there is a pattern. Each line is free-field format input containing five variables:
  1. NI = Element number of the first element in the sequence.
  2. NSEQ = NSEQ subsequent elements will have the moisture content automatically generated.
  3. NAD = Increment of element number for each of the NSEQ subsequent elements.
  4. THNI = Moisture content of element NI (decimal point).
  5. THNIAD = Increment of moisture content for each of the NSEQ subsequent elements (decimal point).

NOTE: A line with five zeros separated by spaces or commas signals the end of this data subset.

  1. Liquid density field: Typically, one line per node is needed. However, automatic generation can be used. Each line is free-field format input containing six variables:
  1. NI = Node number of the first node in the sequence.
  2. NSEQ = NSEQ subsequent nodes will have liquid density automatically generated.
  3. NAD = Increment of the node number for each of the NSEQ subsequent nodes.
  4. RHOLNI = Liquid density at node NI (M/L3).
  5. RHOLAD = Incremental increase of the liquid density for each of the NSEQ subsequent nodes (M/L3).
  6. RHOLRD = Percent increase of the increment over its preceding increment (decimal fraction).

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

II. KVI = -2: Hydrologic time series profiles

This option allows the input of a time series profile for the velocity and the moisture content. The value of the velocity is the same for all nodes, the value for the moisture content is the same for all elements. These values may vary in time. This option is especially useful for one-dimensional simulations, e.g. column experiments with flow interrupts.

  1. Velocity and moisture content: The number of lines in this data subset is NPROF + 1, where NPROF is the number of data points in the hydrologic time series profile. The first line is:
1. NPROF = Number of data points in the hydrologic time series profile

Following this line, one set of data for each point in the profile is entered on one line. The number of variables per line depends on the dimensionality of the simulation.

One-Dimensional Problems:

Each of the following NPROF lines is free-field format input containing three variables:

For IPROF=1, NPROF

  1. TPROF(IPROF) = Time of IPROF-th data point in the hydrologic time series profile (T).
  2. VPROF(IPROF,1) = x-velocity component of all nodes for the IPROF-TH data point (L/T).
  3. THPROF(IPROF) = Moisture content of all elements for the IPROF-TH data point (decimal point).

Two-Dimensional Problems:

Each of the following NPROF lines is free-field format input containing four variables:

For IPROF=1, NPROF

  1. TPROF(IPROF) = Time of IPROF-th data point in the hydrologic time series profile (T).
  2. VPROF(IPROF,1) = x-velocity component of all nodes for the IPROF-TH data point (L/T).
  3. VPROF(IPROF,3) = z-velocity component of all nodes for the IPROF-TH data point (L/T).
  4. THPROF(IPROF) = Moisture content of all elements for the IPROF-TH data point (decimal point).

Three-Dimensional Problems:

Each of the following NPROF lines is free-field format input containing five variables:

For IPROF=1, NPROF

  1. TPROF(IPROF) = Time of IPROF-th data point in the hydrologic time series profile (T).
  2. VPROF(IPROF,1) = x-velocity component of all nodes for the IPROF-TH data point (L/T).
  3. VPROF(IPROF,2) = y-velocity component of all nodes for the IPROF-TH data point (L/T).
  4. VPROF(IPROF,3) = z-velocity component of all nodes for the IPROF-TH data point (L/T).
  5. THPROF(IPROF) = Moisture content of all elements for the IPROF-TH data point (decimal point).
  1. Density: The liquid density is constant throughout the simulation and is entered in the same way as in the case KVI = -1 (see above).

v1.1 - 6/8/1999

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