NPG Nonlinear Processes in Geophysics NPG Nonlin. Processes Geophys. 1607-7946 Copernicus Publications Göttingen, Germany 10.5194/npg-20-1095-2013 Large eddy simulation model for wind-driven sea circulation in coastal areas Petronio A. 1 Roman F. 1 Nasello C. 2 Armenio V. 1 Dipartimento di Ingegneria e Architettura, Università di Trieste, Piazzale Europa 1, 34127 Trieste, Italy Dipartimento di Ingegneria Civile, Ambientale e Aerospaziale, Università di Palermo, Viale delle Scienze, 90128 Palermo, Italy 10 12 2013 20 6 1095 1112 Copyright: © 2013 A. Petronio et al. 2013 This work is licensed under the Creative Commons Attribution 3.0 Unported License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/3.0/ This article is available from https://npg.copernicus.org/articles/20/1095/2013/npg-20-1095-2013.html The full text article is available as a PDF file from https://npg.copernicus.org/articles/20/1095/2013/npg-20-1095-2013.pdf

In the present paper a state-of-the-art large eddy simulation model (LES-COAST), suited for the analysis of water circulation and mixing in closed or semi-closed areas, is presented and applied to the study of the hydrodynamic characteristics of the Muggia bay, the industrial harbor of the city of Trieste, Italy. The model solves the non-hydrostatic, unsteady Navier–Stokes equations, under the Boussinesq approximation for temperature and salinity buoyancy effects, using a novel, two-eddy viscosity Smagorinsky model for the closure of the subgrid-scale momentum fluxes. The model employs: a simple and effective technique to take into account wind-stress inhomogeneity related to the blocking effect of emerged structures, which, in turn, can drive local-scale, short-term pollutant dispersion; a new nesting procedure to reconstruct instantaneous, turbulent velocity components, temperature and salinity at the open boundaries of the domain using data coming from large-scale circulation models (LCM). Validation tests have shown that the model reproduces field measurement satisfactorily. The analysis of water circulation and mixing in the Muggia bay has been carried out under three typical breeze conditions. Water circulation has been shown to behave as in typical semi-closed basins, with an upper layer moving along the wind direction (apart from the anti-cyclonic veering associated with the Coriolis force) and a bottom layer, thicker and slower than the upper one, moving along the opposite direction. The study has shown that water vertical mixing in the bay is inhibited by a large level of stable stratification, mainly associated with vertical variation in salinity and, to a minor extent, with temperature variation along the water column. More intense mixing, quantified by sub-critical values of the gradient Richardson number, is present in near-coastal regions where upwelling/downwelling phenomena occur. The analysis of instantaneous fields has detected the presence of large cross-sectional eddies spanning the whole water column and contributing to vertical mixing, associated with the presence of sub-surface horizontal turbulent structures. Analysis of water renewal within the bay shows that, under the typical breeze regimes considered in the study, the residence time of water in the bay is of the order of a few days. Finally, vertical eddy viscosity has been calculated and shown to vary by a couple of orders of magnitude along the water column, with larger values near the bottom surface where density stratification is smaller.

 References Armenio, V. and Roman, F.: Large Eddy Simulation of environmental shallow water coastal flow, Ercoftac Bulletin, 78, 46–53, 2009. Balaras, E.: Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations, Comput. Fluids, 33, 375–404, 2004. Burchard, H., Craig, P., Gemmrich, J., Haren, H. V., Mathieu, P., Meier, H. M., Smith, W. N., Prandke, H., Rippeth, T., Skyllingstad, E., Smyth, W., Welsh, D., and Wijesekera, H.: Observational and numerical modeling methods for quantifying coastal ocean turbulence and mixing, Prog. Oceanogr., 76, 399–442, 2008. Cabot, W. and Moin, P.: Approximate wall boundary conditions in the Large-Eddy simulation of high Reynolds number flow, Flow, Turbulence Combustion, 63, 269–291, 1999. Coleman, G., Ferziger, J., and Spalart, P.: A numerical study of the turbulent Ekman layer, J. Fluid Mech., 213, 313–348, 1990. Fadlun, E. A., Verzicco, R., Orlandi, P., and Mohd-Yusof, J.: Combined Immersed-Boundary Finite-Difference methods for three-dimensional complex flow simulations, J. Comput. Phys., 61, 35–60, 2000. Fischer, H., List, E., Koh, R., Imberger, J., and Brooks, N.: Mixing in inland and coastal waters, Academic Press USA, 1979. Galea, A., Grifoll, M., Roman, F., Mestres, M., Armenio, V., and Sanchez-Arcilla, A.: Numerical Simulation of water mixing and renewal in the Barcelona harbour area, in: Book of Abstracts of the Fourth International Conference on the Application of Physical Modelling to Port and Coastal Protection, CoastLab12, 2012. Geyer, W., Morris, J., Pahl, F., and Jay, D.: Interaction between physical processes and ecosystem structure: A comparative approach, edited by: Hobbie, J. E., 2000. Hanna, S. R., Briggs, G. A., and Hosker, R.: Handbook on Atmospheric Diffusion, 21–22, United States Department of Energy, 1982. Hosker, R.: Atmospheric Science and Power Production, chap. Flow and diffusion near obstacles, 241–326, United States Department of Energy, 1984. Keating, A., Piomelli, U., Balaras, E., and Kaltenbach, H.: A priori and a posteriori tests of inflow conditions for large-eddy simulation, Phys. Fluids, 16, 4696–4712, 2004. Kuehn, T. and Coleman, J.: ASHRAE Handbook: Fundamentals, ASHRAE, 2005. Piomelli, U.: Wall-layer models for large-eddy simulations, Progr. Aerospace Sci., 44, 437–446, 2008. Pope, S.: Turbulent flows, Cambridge University Press, New York, USA, 771 pp., 2000. Ramachandran, S., Tandon, A., and Mahadevan, A.: Effect of subgrid-scale mixing on the evolution of forced submesoscale instabilities, Ocean Modell., 66, 45–65, 2013. Roman, F., Armenio, V., and Frölich, J.: A simple wall layer model for LES with IBM, Phys. Fluids, 21, 2009a. Roman, F., Napoli, E., Milici, B., and Armenio, V.: An improved Immersed Boundary Method for curvilinear grids, Compute. Fluids, 38, 1510–1527, 2009b. Roman, F., Stipcich, G., Armenio, V., Inghilesi, R., and Corsini, S.: Large Eddy Simulation of mixing in coastal areas, in: Proc. Turbulent Shear Flow Phenomena 6, Seoul, 2009c. Roman, F., Stipcich, G., Armenio, V., Inghilesi, R., and Corsini, S.: Large eddy simulation of mixing in coastal areas, Int. J. Heat Fluid Flow, 31, 327–341, 2010. Scire, J. S., Strimaitis, D. G., and Yamartino, R. J.: A User\rqs Guide for the CALPUFF Dispersion Model, 196 Baker Avenue Concord, MA 01742, 2000. Scotti, A.: Large eddy simulation in the ocean, Int. J. Comput. Fluid Dynam., 24, 393–406, 2010. Scotti, A., Meneveau, C., and Lilly, D.: Generalized Smagorinsky model for anisotropic grids, Phys. Fluids A: Fluid Dynamics, 5, 2306–2308, 1993. Skyllingstad, E. and Denbo, D.: An ocean large-eddy simulation of Langmuir circulation in the surface mixed layer, J. Geophys. Res., 100, 8501–8522, 1995. Stravisi, F.: Thalassia II, seawater temperature, salinity, current velocity, secchi disc and underwater irradiance data in the gulf of Trieste (1979–1984), Tech. rep., University of Trieste, Physics department, 1990. Stravisi, F., Specchi, M., and Princi, M.: Some morphological, physical and chemical observations on Muggia bay (gulf of Trieste – northern Adriatic), Mem. Biol. Marina Oceanogr., Suppl. X, 275–284, 1980. Wang, D., Large, W. G., and McWilliams, J. C.: Large-eddy simulation of the equatorial ocean boundary layer: Diurnal cycling, eddy viscosity, and horizontal rotation, J. Geophys. Res. Oceans, 101, 3649–3662, 1996. Wu, J.: Wind-stress coefficient over sea surface from breeze to hurricane, J. Geophys. Res., 87, 9704–9706, 1982. Zang, J., Street, R., and Koseff, J.: A non-staggered grid, fractional step method for time-dependent incompressible Navier-Stokes equations in curvilinear coordinates, J. Comput. Phys., 114, 18–33, 1994. Zikanov, O., Slinn, D. N., and Dhanak, M.: Large-eddy simulations of the wind-induced turbulent Ekman layer, J. Fluid Mech., 495, 343–368, 2003.