Monitoring the convergence status of a model is important but often difficult due the complexity of the sediment transport models. The sediment log output may be grouped similar to the HDF5 spatial output as: (1) Sediment Transport (sediment concentrations), (2) Sediment Bed (grain class fractions and dry densities in the active layer). In addition, the logging output is written out to two different places: (1) the Computation Log File (.bco), and (2) the HDF5 Output File (.p.hdf).
The amount of information which is written to the Computation Log File depends on the sediment Output Level (see Section 2.4 Sediment Output Options) and the convergence status. The log output in the Computation Log File is output only when a convergence tolerance is exceeded for Output Levels less than 5. If the output level is 5 or 6, a log record is written at a constant computation time step interval. For an Output Level of 5, the log record only contains convergence information for the last outer-loop iteration. For an Output Level of 6 the log record contains convergence information for all iterations.
The log output in the HDF5 Output File is written to the same location as the Sediment Transport and the Sediment Bed output blocks and is therefore written at the corresponding output intervals. The output interval for the Sediment Transport and the Sediment Bed output blocks can be controlled by the user by specifying a multiple of the Mapping Output Interval (see Section 2.4 Sediment Output Options for details). The log output is written as floating point datasets instead of compound datasets or other forms for computational speeds.
The log output for Sediment Transport is written in the HDF5 Output File under
Dataset Name: Log Sediment Transport
Dataset Path: /Results/Unsteady/Output/Output Blocks/Sediment Transport/Unsteady Time Series/2D Flows Areas/<2D Flow Area Name>/
where <2D Flow Area Name> is the name of each 2D flow area.
Table 2 16. Sediment Transport Log Output. Note: HDF5 columns begin at 0.

Column

Log Output Variable

Units

0

Overall convergence status flag

-

1

Transport convergence status flag

-

2

Number of outer-loop iterations

-

3

Maximum total total-load concentration

mg/L

4

Cell location of maximum total total-load concentration

-

5

Maximum fractional total-load concentration

mg/L

6

Cell location of maximum fractional total-load concentration

-

7

Grain class corresponding to maximum fractional total-load concentration

-

8

Concentration maximum absolute correction for all grain classes

mg/L

9

Cell location of maximum absolute concentration correction for all grain classes

-

10

Grain class corresponding to maximum absolute concentration correction for all grain classes

-

11

Percent of active grid

%

12

Minimum total inner loop iterations (only for iterative matrix solvers)

-

13

Grain class corresponding to minimum inner loop iterations (only for iterative matrix solvers)

-

14

Maximum total inner loop iterations (only for iterative matrix solvers)

-

15

Grain class corresponding to maximum inner loop iterations (only for iterative matrix solvers)

-

16

Maximum root-mean-squared residual of inner-loop iterations (only for iterative matrix solvers)

-

17

Grain class corresponding to maximum root-mean-squared residual of inner-loop iterations (only for iterative matrix solvers)

-


Currently when the model diverges, it stops the simulation. In future versions, diverged time steps will be repeated with subcycles in order to try to obtain a converged solution. Usually divergence indicates a NaN (not a number for bed change, bed fractions, or sediment concentrations). However, it could also mean an internal divergence tolerance was exceeded. Future versions of the model will include more information on how the model diverged.
Similarly, the logging output for the Sediment Bed is written in the HDF5 Output File under
Dataset Name: Log Sediment Bed
Dataset Path: /Results/Unsteady/Output/Output Blocks/Sediment Bed/Unsteady Time Series/2D Flows Areas/<2D Flow Area Name>/
where again <2D Flow Area Name> is the name of each 2D flow area.
The convergence status flags indicate different model states (see Table 2-17). A status is computed each outer iteration for the sediment transport and bed calculations. The overall convergence status takes into account transport and bed sorting.
Table 2 17. Convergence Status Values

Status Value

Status Name

Description

-1

Diverged

The solution has diverged. Divergence is determined with a variety of tolerances and checks for many variables.

1

Iterating

Either the convergence tolerances or minimum number of iterations have not been satisfied.

0

Converged

Both the convergence tolerances and minimum number of iterations have been satisfied.

2

Continue

Either the convergence has stalled, or the maximum number of iterations has been reached. The solution will continue to the next time step.



Table 2 18. Sediment Bed Log Output. Note: HDF5 begin start at 0.

Column

Log Output Variable

Units

0

Overall status flag

-

1

Sediment bed status flag

-

2

Number of outer-loop iterations

-

3

Maximum absolute value of active layer fraction correction (difference between outer-loop iteration values)

%

4

Cell location of maximum absolute value of active layer fraction correction

-

5

Subarea corresponding to maximum absolute value of active layer fraction correction

-

6

Grain class corresponding to maximum absolute value of active layer fraction correction

-

7

Number of avalanching iterations

-

8

Number of avalanches at last iteration

-

9

Maximum bed change due to avalanching

ft or m

10

Cell location of maximum bed change to due avalanching

-











References
Ackers, P., and White, W.R. 1973. Sediment transport: New approach and analysis, Proceedings. ASCE, Journal of the Hydraulic Division, 99, HY11, pp. 2041-2060.
Ashida, K. and Michiue, M. 1971. An investigation of river bed degradation downstream of a dam, Proc. 14th Congress of IAHR, Paris, France, 3, 247–256.
Ariathurai, C.R. 1974. A finite element model for sediment transport in estuaries, PhD Thesis, University of California, Berkley, USA.
Armanini, A., and G. di Silvio. 1986. Discussion on the paper 'A depth-averaged model for suspended sediment transport' by Galappattti, G., and Vreugdenhil, C.B. Journal of Hydraulic Research. 24(5), 437–441.
Chazan, D., and Miranker, W. 1969. Chaotic relaxation. Linear Algebra and Its Applications, 2(2), 199–222.
Colby, B.R. 1963. Discussion of "Sediment transportation mechanics: Introduction and properties of sediment", Progress report by the Task Committee on Preparation of Sediment Manual of the Committee on Sedimentation of the Hydraulics Division, V. A. Vanoni, Chmn., Journal of Hydraulics Division., ASCE, 89(1): 266–268.
Costa, R.C.F.G. 1989. Flow-fine sediment hysteresis in sediment-stratified coastal waters, MS Thesis, University of Florida, Gainesville, FL.
Dahl, T.A, Heath, R.E., Gibson, S.A., and Nygaard, C.J. 2018. HEC-RAS unsteady flow and sediment model of the Mississippi River: Tarbet Landing to the Gulf, Mississippi River Geomorphology & Potomology Program, Report No. 25, Engineer Research and Development Center, Vicksburg, MS, USA.
Day, T.J. 1980. A study of the transport of graded sediments, Report N. IT 190, HR Wallingford, U.K.
DeLong, M.A. 1997. SOR as a Preconditioner. PhD Dissertation. School of Engineering and Applied Science, University of Virginia, 102 p.
Dietrich, W.E. 1982. Settling velocity of natural particles, Water Resources Research, 18(6), 1615–1626.
DHL 1980. Computation of Siltation in Dredged Trenches, Delft Hydraulics Laboratory, Storm Surge Barrier Oosterschelde, Report 1267-V/M 1570.
Dixit, J.G. 1982. Resuspension potential of deposited kaolinite beds. M.S. thesis, University of Florida.
Egiazaroff, I.V. 1965. Calculation of nonuniform sediment concentration, J. Hydr. Div., ASCE, 91(HY4), 225–247.
Engelund, F., and Hansen, E. 1967. A monograph on sediment transport in alluvial streams, Teknisk Forlag, Copenhagen, 65.
Gallappatti, G., and Vreugdenhil, C.B. 1985. A depth-integrated model for suspended sediment transport. J. Hydr. Res., 23(4), 359–377.
Ganju, N.K. 2001. Trapping organic-rich fine sediment in an estuary. MS Thesis, University of Florida, Gaiesville, FL.
Harten, A. 1983. High resolution schemes for hyperbolic conservation laws, Journal of Computational Physics. 49, 357–393.
Hayashi, T.S., Ozaki and Ichibashi, T. 1980. Study on bed load transport of sediment mixture, Proc. 24th Japanese Conf. on Hydraulics, Japan.
Hwang, K.-N. 1989. Erodibility of fine sediment in wave dominated environments. MS Thesis, University of Florida, Gainesville, FL.
Krone, R.B, 1962. Flume studies of transport of sediment in estuarial shoaling processes. Final Report, Hydraulic Engineering Laboratory and Sanitary Engineering Research Laboratory, University of California, Berkeley, CA.
Kahan, W. 1958. Gauss-Seidel methods of solving large systems of linear equations, PhD thesis, University of Toronto.
Kuijper, C., Cornelisse, J.M., and Winterwerp, J.C. 1990. Erosion and deposition of natural muds–sediments from the Western Scheldt (Breskens), Rep. 29, Rijkswaterstaat Delft Hydraul. Cohesive Sed., Delft, Netherlands.
Leone, R.De., and Mangasarian, O.L. 1988. Asynchronous parallel successive overrelaxation for the symmetric linear complementarity problem. Mathematical Programming, 42(1–3), 347–361.
Liu, W.-C., Hsu, M.-H., and Kuo, A.Y. 2002. Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system, Taiwan, Marine Pollution Bulletin, 44. 1076–1088.
Marván, F.G. 2001. A two-dimensional numerical transport model for organic-rich cohesive sediments in estuarine waters. PhD Thesis, Heriot-Watt University, Edinburgh, UK.
Nichols, M. 1984. Fluid mud accumulation processes in an estuary. Geo-Marine Letters, 4–5, 192–198.
Owen, M.W. 1970. A detailed study of the settling of an estuary mud. Report No IT 78, HR Wallingford, Wallingford, UK.
Patankar, S.V. 1980. Numerical Heat Transfer and Fluid Flow. Hemisphere: Washington, DC, 197 p.
Partheniades, E. 1962. A study of erosion and deposition of cohesive soils in salt water. Ph.D. thesis, University of California, Berkeley, 182 p
Richardson, J.F. and Zaki, W.N. 1954. Sedimentation and fluidization, Part I. Transaction of the Institute of Chemical Engineers., 32(1), 35–53.
Rubey, W.W. 1933. Settling velocities of gravel, sand, and silt particles, American Journal of Science, 5th Series, 25(148), 325–338.
Saad, Y. 1993. A flexible inner-outer preconditioned GMRES algorithm. SIAM Journal Scientific Computing (14): 461–469.
Schenk, O. and Gartner, K. 2011. Parallel Sparse Direct and Multi-Recursive Iterative Linear Solvers: Pardiso User Guide 4.1.2.
Schenk, O. and Gartner, K. 2004. Solving unsymmetric sparse systems of linear equations with PARDISO. Journal of Future Generation Computer Systems, 20(3):475–487.
Soulsby, R.L. 1997. Dynamics of marine sands. A manual for practical applications. Thomas Telford Publications, London, England, 249 p.
Toffaleti, F.B. 1968. A procedure for computation of total river sand discharge and detailed distribution, bed to surface, Technical Report No. 5, Committee on Channel Stabilitzation, U.S. Army Corps of Engineers, November.
van Rijn, L.C. 1993. Principles of sediment transport in river, estuaries, coastal seas, and oceans. International Institute for Infrastructural, Hydraulic, and Environmental Engineering, Delft, The Netherlands.
Wooster, J.K., Dusterhoff, S.R., Cui, Y., Sklar, L.S., Dietrich, W.E., and Malko, M. 2008. Sediment supply and relative size distribution effects on fine sediment infiltration into immobile gravels, Water Resources Research, 44, W03424.
Wei, H., Nearing, M.A., Stone, J.J., Guertin, D.P., Spaeth, K.E., Pierson, F.B., Nichols, M.H., and Moffett, C.A. 2009. A new splash and sheet erosion equation for rangelands. Soil Science Society of America Journal 73(4): 1386–1392.
Wu. W., Altinakar, M. and Wang, S.S.Y. 2006. Depth-average analysis of hysteresis between flow and sediment transport under unsteady conditions, Int. J. Sediment Res., 21(2), 101–112.
Wu, W., Wang, S.S.Y. and Jia, Y. 2000. Nonuniform sediment transport in alluvial rivers, J. Hydr. Res., IAHR, 38(6), 427–434.
Wu, W. and Wang, S.S.Y. 2006. Formulas for sediment porosity and settling velocity, Journal of Hydraulic Engineering, ASCE, 132(8), 858–862.
Zhu, J. 1992. On the higher-order bounded discretization schemes for finite volume computations of incompressible flows, Computational Methods in Applied Mechanical Engineering, 98, 345.