Reducing Numerical Diffusion Effects with Pycnocline Filter
Publication: Journal of Hydraulic Engineering
Volume 129, Issue 3
Abstract
Numerical or artificial diffusion is the unintentional smoothing of gradients associated with the discretization of the transport equations. In lakes and reservoirs where through-flow is small, the effects of numerical diffusion of mass are cumulative, leading to a progressive weakening of vertical density stratification. This density field misrepresentation precludes accurate, long-term, three-dimensional (3D), hydrodynamic simulations on fixed grids in closed basins with an active thermocline. An ad hoc technique to limit the destratifying effects of numerical diffusion of mass is presented and tested for a 3D, hydrostatic, Z-coordinate numerical model. The technique quantifies the domain-integrated numerical diffusion by assessing the change in the background potential energy At each time step, the change in associated with numerical diffusion is calculated, then removed using a sharpening filter applied to each water column. In idealized test cases, the filtering technique is effective in maintaining density stratification over one year while undergoing periodic, large-amplitude forcing by internal waves. Forty-day simulations of Lake Kinneret compared to field measurements demonstrate improved representation of density stratification using the filtering technique.
Get full access to this article
View all available purchase options and get full access to this article.
References
Ahsan, A. K. M. Q., and Blumberg, A. F.(1999). “Three-dimensional hydrothermal model of Onondaga Lake, New York.” J. Hydraul. Eng., 125(9), 912–923.
Antenucci, J. P., Imberger, J., and Saggio, A.(2000). “Seasonal evolution of the basin-scale internal wave field in a large stratified lake.” Limnol. Oceanogr., 45(7), 1621–1638.
Beletsky, D., and Schwab, D. J.(2001). “Modeling circulation and thermal structure in Lake Michigan: Annual cycle and interannual variability.” J. Geophys. Res., [Oceans], 106(C9), 19745–19771.
Casulli, V., and Cheng, R. T.(1992). “Semi-implicit finite difference methods for three-dimensional shallow water flow.” Int. J. Numer. Methods Fluids, 15, 629–648.
Farmer, D. M.(1978). “Observations of long nonlinear internal waves in a lake.” J. Phys. Oceanogr., 8(1), 63–73.
Farrow, D., and Stevens, D. P.(1995). “A new tracer advection scheme for Bryan and Cox type ocean general circulation models.” J. Phys. Oceanogr., 25(7), 1731–1741.
Fozdar, F. M., Parker, G. J., and Imberger, J.(1985). “Matching temperature and conductivity sensor response characteristics.” J. Phys. Oceanogr., 15(11), 1557–1569.
Griffies, S. M., Pacanowski, R. C., and Hallberg, R. W.(2000). “Spurious diapycnal mixing associated with advection in a z-coordinate ocean model.” Mon. Weather Rev., 128(3), 538–564.
Gross, E. S., Koseff, J. R., and Monismith, S. G.(1999a). “Evaluation of advective schemes for estuarine salinity simulations.” J. Hydraul. Eng., 125(1), 32–46.
Gross, E. S., Koseff, J. R., and Monismith, S. G.(1999b). “Three-dimensional salinity simulations of South San Francisco Bay.” J. Hydraul. Eng., 125(11), 1199–1209.
Hecht, M. W., Bryan, F. O., and Holland, W. R.(1998). “A consideration of tracer advection schemes in a primitive equation ocean model.” J. Geophys. Res., [Oceans], 103(C2), 3301–3321.
Hecht, M. W., Holland, W. R., and Rasch, P. J.(1995). “Upwind-weighted advection schemes for ocean tracer transport: An evaluation in a passive tracer context.” J. Geophys. Res., [Oceans], 100(C10), 20763–20778.
Hecht, M. W., Wingate, B. A., and Kassis, P.(2000). “A better, more discriminating test problem for ocean tracer transport.” Ocean Modelling, 2(January), 1–15.
Hodges, B. R. (2000). “Numerical techniques in CWR-ELCOM.” Technical Rep., Centre for Water Research, University of Western Australia, Nedlands, Western Australia, 6907. Reference No. WP 1422-BH.
Hodges, B. R., Imberger, J., Saggio, A., and Winters, K.(2000). “Modeling basin-scale internal waves in a stratified lake.” Limnol. Oceanogr., 45(7), 1603–1620.
Horn, D. A., Imberger, J., and Ivey, G. N.(2001). “The degeneration of large-scale interfacial waves in lakes.” J. Fluid Mech., 434, 181–207.
Imberger, J. (1994). “Transport processes in lakes: A review.” Limnology now: A paradigm of planetary problems, Margalef, R., ed., Elsevier Science, Amsterdam, 99–193.
Imberger, J., and Patterson, John C.(1990). “Physical limnology.” Adv. Fluid Mech., 27, 303–475.
Kundu, P. K. (1990). Fluid mechanics, Academic, New York, 638.
Laval, B., Imberger, J., and Hodges, B. R. (2002). “Modelling circulation in lakes: Spatial and temporal variations.” Limnol. Oceanogr., in press.
Leonard, B. P.(1991). “The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection.” Comput. Methods Appl. Mech. Eng., 88, 17–74.
Lin, B., and Falconer, R. A.(1997). “Tidal flow and transport modeling using ULTIMATE QUICKEST scheme.” J. Hydraul. Eng., 123(4), 303–314.
Lorenz, E. N.(1955). “Available potential energy and the maintenance of the general circulation.” Tellus, 7(2), 157–167.
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P. (1992). Numerical recipes in Fortran 77: The art of scientific computing, 2nd Ed., Vol. 1. Cambridge University, Cambridge, England.
Saggio, A., and Imberger, J.(1998). “Internal wave weather in a stratified lake.” Limnol. Oceanogr., 43(8), 1780–1795.
Saggio, A., and Imberger, J.(2001). “Mixing and turbulent fluxes in the metalimnion of a stratified lake.” Limnol. Oceanogr., 46(2), 392–409.
Sanderson, B. G.(1998). “Order and resolution for computational ocean dynamics.” J. Phys. Oceanogr., 28(6), 1271–1286.
Staniforth, A., and Co⁁té, J.(1991). “Semi-Lagrangian integration schemes for atmospheric models—A review.” Mon. Weather Rev., 119(9), 2206–2223.
Thorpe, S. A., Keen, J. M., Jiang, R., and Lemmin, U.(1996). “High-frequency internal waves in Lake Geneva.” Philos. Trans. R. Soc. London, Ser. A, 354, 237–257.
Winters, K. B., Lombard, P. N., Riley, J. J., and D’Asaro, E. A.(1995). “Available potential energy and mixing in density-stratified fluids.” J. Fluid Mech., 289, 115–128.
Wüest, A., Piepke, G., and Van Senden, D. C.(2000). “Turbulent kinetic energy balance as a tool for estimating vertical diffusivity in wind-forced stratified waters.” Limnol. Oceanogr., 45(6), 1388–1400.
Wüest, A., Van Senden, D. C., Imberger, J., Piepke, G., and Gloor, M.(1996). “Comparison of diapycnal diffusivity measured by tracer and microstructure techniques.” Dyn. Atmos. Oceans, 24(1), 27–39.
Information & Authors
Information
Published In
Copyright
Copyright © 2003 American Society of Civil Engineers.
History
Received: Mar 20, 2001
Accepted: Sep 3, 2002
Published online: Feb 14, 2003
Published in print: Mar 2003
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.