Test of the Volume-of-Fluid Method on Simulations of Marine Boundary Layer Clouds

C-Y. J. Kao Earth and Environment Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico

Search for other papers by C-Y. J. Kao in
Current site
Google Scholar
PubMed
Close
,
Y. H. Hang Earth and Environment Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico

Search for other papers by Y. H. Hang in
Current site
Google Scholar
PubMed
Close
,
J. M. Reisner Earth and Environment Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico

Search for other papers by J. M. Reisner in
Current site
Google Scholar
PubMed
Close
, and
W. S. Smith Earth and Environment Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico

Search for other papers by W. S. Smith in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

The impact of using grid-averaged thermodynamic properties (i.e., neglecting their subgrid variability due to partial cloudiness) to represent forcings for condensation or evaporation has long been recognized. In particular, numerical difficulties in terms of spurious oscillations and/or diffusion in vicinity of a cloud–environment interface have been encountered in most of the conventional finite-difference Eulerian advection schemes. This problem is equivalent to the inability of models to accurately track the cloud boundary within a grid cell, which eventually leads to spurious production or destruction of cloud water at leading or trailing edges of clouds. This paper employs a specialized technique called the “volume-of-fluid” (VOF) method to better parameterize the subgrid-scale advection process that accounts for the transport of material interfaces. VOF also determines the actual location of the partial cloudiness within a grid box. Consequently, relevant microphysical parameterizations in mixed cells can be consistently applied in “cloudy” and “clear” regions. The VOF technique is incorporated in a two-dimensional hydrodynamic model to simulate the diurnal cycle of the marine stratocumulus-capped boundary layer. The fidelity of VOF to advection–condensation processes under a diurnal radiative forcing is assessed by comparing the model simulation with data taken during the ISCCP FIRE observational period as well as with results from simulations without VOF. This study shows that the VOF method indeed suppresses the spurious cloud boundary instability and supports a multiday cloud evolution as observed. For the case without VOF, the spurious instability near the cloud top causes the dissipation of the entire cloud layer within a half of a diurnal cycle.

* On leave from Institute of Applied Physics and Computational Mathematics, Beijing, China.

Corresponding author address: Dr. C.-Y. J. Kao, Los Alamos National Laboratory, EES-8, MS C300, Los Alamos, NM 87545.

Abstract

The impact of using grid-averaged thermodynamic properties (i.e., neglecting their subgrid variability due to partial cloudiness) to represent forcings for condensation or evaporation has long been recognized. In particular, numerical difficulties in terms of spurious oscillations and/or diffusion in vicinity of a cloud–environment interface have been encountered in most of the conventional finite-difference Eulerian advection schemes. This problem is equivalent to the inability of models to accurately track the cloud boundary within a grid cell, which eventually leads to spurious production or destruction of cloud water at leading or trailing edges of clouds. This paper employs a specialized technique called the “volume-of-fluid” (VOF) method to better parameterize the subgrid-scale advection process that accounts for the transport of material interfaces. VOF also determines the actual location of the partial cloudiness within a grid box. Consequently, relevant microphysical parameterizations in mixed cells can be consistently applied in “cloudy” and “clear” regions. The VOF technique is incorporated in a two-dimensional hydrodynamic model to simulate the diurnal cycle of the marine stratocumulus-capped boundary layer. The fidelity of VOF to advection–condensation processes under a diurnal radiative forcing is assessed by comparing the model simulation with data taken during the ISCCP FIRE observational period as well as with results from simulations without VOF. This study shows that the VOF method indeed suppresses the spurious cloud boundary instability and supports a multiday cloud evolution as observed. For the case without VOF, the spurious instability near the cloud top causes the dissipation of the entire cloud layer within a half of a diurnal cycle.

* On leave from Institute of Applied Physics and Computational Mathematics, Beijing, China.

Corresponding author address: Dr. C.-Y. J. Kao, Los Alamos National Laboratory, EES-8, MS C300, Los Alamos, NM 87545.

Save
  • Albrecht, B. A., D. A. Randall, and S. Nicholls, 1988: Observations of marine stratocumulus during FIRE. Bull. Amer. Meteor. Soc.,69, 618–626.

  • Blaskovic, M., R. Davies, and J. B. Snider, 1991: Diurnal variation of marine stratocumulus over San Nicolas Island during July, 1987. Mon. Wea. Rev.,119, 1469–1478.

  • Duda, D. P., G. L. Stephens, and S. K. Cox, 1991: Microphysical and radiative properties of marine stratocumulus from tethered balloon measurements. J. Appl. Meteor.,30, 170–186.

  • Duynkerke, P. G., and P. Hignett, 1993: Simulations of diurnal variation in a stratocumulus-capped marine boundary layer during FIRE. Mon. Wea. Rev.,121, 3291–3330.

  • Grabowski, W. W., 1989: Numerical experiments on the dynamics of the cloud–environment interface: Small cumulus in a shear-free environment. J. Atmos. Sci.,46, 3513–3541.

  • ——, and P. K. Smolarkiewicz, 1990: Monotone finite-difference approximations to the advection–condensation problem. Mon. Wea. Rev.,118, 2082–2097.

  • ——, and ——, 1996: Two-time-level semi-Lagrangian modeling of precipitating clouds. Mon. Wea. Rev.,124, 487–497.

  • Hignett, P., 1991: Observations of diurnal variation in a cloud-capped marine boundary layer. J. Atmos. Sci.,48, 1472–1482.

  • Hirt, C. W., and B. D. Nichols, 1981: Volume of fluid (VOF) method for the dynamics of free boundaries. J. Comput. Phys.,39, 201–225.

  • Kloesel, K. A., B. A., Albrecht, and D. P. Wylie, 1988: FIRE marine stratocumulus observations—Summary of operations and synoptic conditions. FIRE Tech. Rep. 1, The Pennsylvania State University, 171 pp. [Available from Department of Meteorology, The Pennsylvania State University, University Park, PA 16802.].

  • Linden, P. F., J. M. Redondo, and D. L. Youngs, 1994: Molecular mixing in Rayleigh–Taylor instability. J. Fluid Mech.,265, 97–124.

  • Liou, K. N., and G. D. Wittman, 1979: Parameterization of radiative properties of clouds. J. Atmos. Sci.,36, 1261–1273.

  • ——, K. P. Freeman, and T. Sasamori, 1978: Cloud and aerosol effects on the solar heating rate of the atmosphere. Tellus,30, 62–70.

  • Manabe, S., and R. F. Strickler, 1964: Thermal equilibrium of the atmosphere with a convective adjustment. J. Atmos. Sci.,21, 361–385.

  • Margolin, L., J. M. Reisner, and P. K. Smolarkiewicz, 1997: Applications of the volume-of-fluid method to the advection–condensation problem. Mon. Wea. Rev.,125, 2265–2273.

  • ——, P. K. Smolarkiewicz, and Z. Sorbjan, 1999: Large-eddy simulations of convective boundary layers using nonoscillatory differencing. Physica D,133, 390–397.

  • Moeng, C.-H., and Coauthors, 1996: Simulation of a stratocumulus-topped planetary boundary layer: Intercomparison among different numerical codes. Bull. Amer. Meteor. Soc.,77, 261–278.

  • Oran, S. E., and J. P. Boris, 1993: Computing turbulent shear flows—A convenient conspiracy. Comput. Phys.,7, 523–533.

  • Parker, B. J., and D. L. Youngs, 1992: Two and three dimensional Eulerian simulation of fluid flow with material interfaces. Atomic Weapon Research Establishment (AWRE) Preprint 01/92, 29 pp.

  • Porter, D. H., A. Pouquet, and P. R. Woodward, 1994: Kolomogorov-like spectra in decaying three-dimensional supersonic flows. Phys. Fluids,6, 2133–2142.

  • Rawlins, F., and J. S. Foot, 1990: Remotely sensed measurements of stratocumulus properties during FIRE using the C130 aircraft multi-channel radiometer. J. Atmos. Sci.,47, 2488–2503.

  • Reisner, J., and P. K. Smolarkiewicz, 1994: Effects of surface thermal forcing on flow past isolated obstacles. J. Atmos. Sci.,51, 117–133.

  • Sasamori, T., 1968: The radiative cooling calculation for application to general circulation experiments. J. Appl. Meteor.,7, 721–729.

  • Slingo, A., and H. M. Schrecker, 1982: On the shortwave radiative properties of stratiform clouds. Quart. J. Roy. Meteor. Soc.,108, 407–426.

  • Smith, W. S., and C.-Y. J. Kao, 1996: Numerical simulations of stratiform clouds observed in 1987 FIRE IOP with the emphasis on diurnal variation. Mon. Wea. Rev.,124, 1803–1816.

  • Smolarkiewicz, P. K., and L. G. Margolin, 1993: On forward-in-time differencing for fluids: Extension to a curvilinear framework. Mon. Wea. Rev.,121, 1847–1859.

  • ——, and ——, 1994: Variational solver for elliptic problems in atmospheric flows. Appl. Math. Comput. Sci.,4, 527–551.

  • ——, and ——, 1997: On forward-in-time differencing for fluids: An Eulerian/semi-Lagrangian nonhydrostatic model for stratified flows. Atmos.–Ocean,35, 127–152.

  • ——, and ——, 1998: MPDATA: A finite-difference solver for geophysical flows. J. Comput. Phys.,140, 459–480.

  • Stevens, B., R. L. Walko, G. Feingold, and W. R. Cotton, 1996: The spurious production of cloud edge supersaturations by Eulerian models. Mon. Wea. Rev.,124, 1034–1041.

  • Syrett, W. J., 1988: Hourly wind, potential temperature, and Richardson number profiles at San Nicolas Island during project FIRE. FIRE Tech. Rep. 2, The Pennsylvania State University, 49 pp. [Available from Department of Meteorology, The Pennsylvania State University, University Park, PA 16802.].

  • Youngs, D. L., 1984: An interface tracking method for a 3D Eulerian hydrodynamics code. AWRE Preprint 44/92, 35 pp.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 162 67 10
PDF Downloads 77 25 3