Large-Eddy Simulation von bewölkten Grenzschichten zur Untersuchung von Bilanzen der statistischen Momente zweiter Ordnung und zur Überprüfung von Turbulenzmodellen

[no abstract

The Parallelized Large-Eddy Simulation Model (PALM) version 4.0 for atmospheric and oceanic flows: model formulation, recent developments, and future perspectives

In this paper we present the current version of the Parallelized Large-Eddy Simulation Model (PALM) whose core has been developed at the Institute of Meteorology and Climatology at Leibniz Universität Hannover (Germany). PALM is a Fortran 95-based 5 code with some Fortran 2003 extensions and has been applied for the simulation of a variety of atmospheric and oceanic boundary layers for more than 15 years. PALM is optimized for use on massively parallel computer architectures and was recently ported to general-purpose graphics processing units. In the present paper we give a detailed description of the current version of the model and its features, such as an embedded 10 Lagrangian cloud model and the possibility to use Cartesian topography. Moreover, we discuss recent model developments and future perspectives for LES applications.DFG/RA/617/3DFG/RA/617/6DFG/RA/617/16DFG/RA/617/27-

Second-moment budgets in cloud topped boundary layers: A large-eddy simulation study

A detailed analysis of second-order moment budgets for cloud topped boundary layers (CTBLs) is performed using high-resolution large-eddy simulation (LES). Two CTBLs are simulated—one with trade wind shallow cumuli, and the other with nocturnal marine stratocumuli. Approximations to the ensemble-mean budgets of the Reynolds-stress components, of the fluxes of two quasi-conservative scalars, and of the scalar variances and covariance are computed by averaging the LES data over horizontal planes and over several hundred time steps. Importantly, the subgrid scale contributions to the budget terms are accounted for. Analysis of the LES-based second-moment budgets reveals, among other things, a paramount importance of the pressure scrambling terms in the Reynolds-stress and scalar-flux budgets. The pressure-strain correlation tends to evenly redistribute kinetic energy between the components, leading to the growth of horizontal-velocity variances at the expense of the vertical-velocity variance which is produced by buoyancy over most of both CTBLs. The pressure gradient-scalar covariances are the major sink terms in the budgets of scalar fluxes. The third-order transport proves to be of secondary importance in the scalar-flux budgets. However, it plays a key role in maintaining budgets of TKE and of the scalar variances and covariance. Results from the second-moment budget analysis suggest that the accuracy of description of the CTBL structure within the second-order closure framework strongly depends on the fidelity of parameterizations of the pressure scrambling terms in the flux budgets and of the third-order transport terms in the variance budgets.COST/ES/0905DF

Detection and attribution of aerosol-cloud interactions in large-domain large-eddy simulations with the ICOsahedral Non-hydrostatic model

Clouds and aerosols contribute the largest uncertainty to current estimates and interpretations of the Earth’s changing energy budget. Here we use a new-generation large-domain large-eddy model, ICON-LEM (ICOsahedral Non-hydrostatic Large Eddy Model), to simulate the response of clouds to realistic anthropogenic perturbations in aerosols serving as cloud condensation nuclei (CCN). The novelty compared to previous studies is that (i) the LEM is run in weather prediction mode and with fully interactive land surface over a large domain and (ii) a large range of data from various sources are used for the detection and attribution. The aerosol perturbation was chosen as peak-aerosol conditions over Europe in 1985, with more than fivefold more sulfate than in 2013. Observational data from various satellite and ground-based remote sensing instruments are used, aiming at the detection and attribution of this response. The simulation was run for a selected day (2 May 2013) in which a large variety of cloud regimes was present over the selected domain of central Europe. It is first demonstrated that the aerosol fields used in the model are consistent with corresponding satellite aerosol optical depth retrievals for both 1985 (perturbed) and 2013 (reference) conditions. In comparison to retrievals from ground-based lidar for 2013, CCN profiles for the reference conditions were consistent with the observations, while the ones for the 1985 conditions were not. Similarly, the detection and attribution process was successful for droplet number concentrations: the ones simulated for the 2013 conditions were consistent with satellite as well as new ground-based lidar retrievals, while the ones for the 1985 conditions were outside the observational range. For other cloud quantities, including cloud fraction, liquid water path, cloud base altitude and cloud lifetime, the aerosol response was small compared to their natural variability. Also, large uncertainties in satellite and ground-based observations make the detection and attribution difficult for these quantities. An exception to this is the fact that at a large liquid water path value (LWP > 200 g m−2), the control simulation matches the observations, while the perturbed one shows an LWP which is too large. The model simulations allowed for quantifying the radiative forcing due to aerosol–cloud interactions, as well as the adjustments to this forcing. The latter were small compared to the variability and showed overall a small positive radiative effect. The overall effective radiative forcing (ERF) due to aerosol–cloud interactions (ERFaci) in the simulation was dominated thus by the Twomey effect and yielded for this day, region and aerosol perturbation −2.6 W m$^{-2}$. Using general circulation models to scale this to a global-mean present-day vs. pre-industrial ERFaci yields a global ERFaci of −0.8 W m$^{-2}$

Large eddy simulation using the general circulation model ICON

ICON (ICOsahedral Nonhydrostatic) is a unified modeling system for global numerical weather prediction (NWP) and climate studies. Validation of its dynamical core against a test suite for numerical weather forecasting has been recently published by Zängl et al. (2014). In the present work, an extension of ICON is presented that enables it to perform as a large eddy simulation (LES) model. The details of the implementation of the LES turbulence scheme in ICON are explained and test cases are performed to validate it against two standard LES models. Despite the limitations that ICON inherits from being a unified modeling system, it performs well in capturing the mean flow characteristics and the turbulent statistics of two simulated flow configurations - one being a dry convective boundary layer and the other a cumulus-topped planetary boundary layer.BMBF/01LK1202

The HD(CP)² Observational Prototype Experiment (HOPE) – an overview

The HD(CP)2 Observational Prototype Experiment (HOPE) was performed as a major 2-month field experiment in Jülich, Germany, in April and May 2013, followed by a smaller campaign in Melpitz, Germany, in September 2013. HOPE has been designed to provide an observational dataset for a critical evaluation of the new German community atmospheric icosahedral non-hydrostatic (ICON) model at the scale of the model simulations and further to provide information on land-surface–atmospheric boundary layer exchange, cloud and precipitation processes, as well as sub-grid variability and microphysical properties that are subject to parameterizations. HOPE focuses on the onset of clouds and precipitation in the convective atmospheric boundary layer. This paper summarizes the instrument set-ups, the intensive observation periods, and example results from both campaigns. HOPE-Jülich instrumentation included a radio sounding station, 4 Doppler lidars, 4 Raman lidars (3 of them provide temperature, 3 of them water vapour, and all of them particle backscatter data), 1 water vapour differential absorption lidar, 3 cloud radars, 5 microwave radiometers, 3 rain radars, 6 sky imagers, 99 pyranometers, and 5 sun photometers operated at different sites, some of them in synergy. The HOPE-Melpitz campaign combined ground-based remote sensing of aerosols and clouds with helicopter- and balloon-based in situ observations in the atmospheric column and at the surface. HOPE provided an unprecedented collection of atmospheric dynamical, thermodynamical, and micro- and macrophysical properties of aerosols, clouds, and precipitation with high spatial and temporal resolution within a cube of approximately 10  ×  10  ×  10 km3. HOPE data will significantly contribute to our understanding of boundary layer dynamics and the formation of clouds and precipitation. The datasets have been made available through a dedicated data portal. First applications of HOPE data for model evaluation have shown a general agreement between observed and modelled boundary layer height, turbulence characteristics, and cloud coverage, but they also point to significant differences that deserve further investigations from both the observational and the modelling perspective

Discovery and Characterization of a Faint Stellar Companion to the A3V Star Zeta Virginis

Through the combination of high-order Adaptive Optics and coronagraphy, we report the discovery of a faint stellar companion to the A3V star zeta Virginis. This companion is ~7 magnitudes fainter than its host star in the H-band, and infrared imaging spanning 4.75 years over five epochs indicates this companion has common proper motion with its host star. Using evolutionary models, we estimate its mass to be 0.168+/-.016 solar masses, giving a mass ratio for this system q = 0.082. Assuming the two objects are coeval, this mass suggests a M4V-M7V spectral type for the companion, which is confirmed through integral field spectroscopic measurements. We see clear evidence for orbital motion from this companion and are able to constrain the semi-major axis to be greater than 24.9 AU, the period > 124$ yrs, and eccentricity > 0.16. Multiplicity studies of higher mass stars are relatively rare, and binary companions such as this one at the extreme low end of the mass ratio distribution are useful additions to surveys incomplete at such a low mass ratio. Moreover, the frequency of binary companions can help to discriminate between binary formation scenarios that predict an abundance of low-mass companions forming from the early fragmentation of a massive circumstellar disk. A system such as this may provide insight into the anomalous X-ray emission from A stars, hypothesized to be from unseen late-type stellar companions. Indeed, we calculate that the presence of this M-dwarf companion easily accounts for the X-ray emission from this star detected by ROSAT.Comment: 9 pages, 6 figures, Accepted to Ap

The Added Value of Large-Eddy and Storm-Resolving Models for Simulating Clouds and Precipitation

More than one hundred days were simulated over very large domains with fine (0.156 km to 2.5 km) grid spacing for realistic conditions to test the hypothesis that storm (kilometer) and large-eddy (hectometer) resolving simulations would provide an improved representation of clouds and precipitation in atmospheric simulations. At scales that resolve convective storms (storm-resolving for short), the vertical velocity variance becomes resolved and a better physical basis is achieved for representing clouds and precipitation. Similarly to past studies we found an improved representation of precipitation at kilometer scales, as compared to models with parameterized convection. The main precipitation features (location, diurnal cycle and spatial propagation) are well captured already at kilometer scales, and refining resolution to hectometer scales does not substantially change the simulations in these respects. It does, however, lead to a reduction in the precipitation on the time-scales considered – most notably over the ocean in the tropics. Changes in the distribution of precipitation, with less frequent extremes are also found in simulations incorporating hectometer scales. Hectometer scales appear to be more important for the representation of clouds, and make it possible to capture many important aspects of the cloud field, from the vertical distribution of cloud cover, to the distribution of cloud sizes, and to the diel (daily) cycle. Qualitative improvements, particularly in the ability to differentiate cumulus from stratiform clouds, are seen when one reduces the grid spacing from kilometer to hectometer scales. At the hectometer scale new challenges arise, but the similarity of observed and simulated scales, and the more direct connection between the circulation and the unconstrained degrees of freedom make these challenges less daunting. This quality, combined with already improved simulation as compared to more parameterized models, underpins our conviction that the use and further development of storm-resolving models offers exciting opportunities for advancing understanding of climate and climate change

Evaluation of Flutter Analyses in the Context of New Concepts of Preliminary Aircraft Design

The process chain of the aeroelasic flutter analysis for new aircraft configurations is presented in the context of the preliminary aircraft design and optimisation program PrADO of the IFL. Therefore a recalculated Fairchild Dornier 728JET is used as an example. The different software modules used in this process chain are described together with the models, which are needed for the computation: PrADO is used to calculate a fully converged aircraft design by an iterative process. As eigen frequencies and eigen modes of the aircraft structure are used by flutter analysis, a structural model for modal analyses with MSC.NASTRAN is generated by a PrADO sub module and transformed into the NASTRAN bulkdata format. With geometry information of PrADO database a ZAERO model is created to perform automated evaluation of flutter proneness. The results of the flutter analysis of the FD 728JET, performed at Fairchild Dornier, are compared to the results of the PrADO design. The differences between both are detailed. Further a modified PrADO 728JET with winglets is compared to the baseline design

Analysis of pressure-strain and pressure gradient-scalar covariances in cloud-topped boundary layers: A large-eddy simulation study

A detailed analysis of the pressure-scrambling terms (i.e., the pressure-strain and pressure gradient-scalar covariances) in the Reynolds-stress and scalar-flux budgets for cloud-topped boundary layers (CTBLs) is performed using high-resolution large-eddy simulation (LES). Two CTBLs are simulated — one with trade wind shallow cumuli, and the other with nocturnal marine stratocumuli. The pressure-scrambling terms are decomposed into contributions due to turbulence-turbulence interactions, mean velocity shear, buoyancy, and Coriolis effects. Commonly used models of these contributions, including a simple linear model most often used in geophysical applications and a more sophisticated two-component-limit (TCL) nonlinear model, are tested against the LES data. The decomposition of the pressure-scrambling terms shows that the turbulence-turbulence and buoyancy contributions are most significant for cloud-topped boundary layers. The Coriolis contribution is negligible. The shear contribution is generally of minor importance inside the cloudy layers, but it is the leading-order contribution near the surface. A comparison of models of the pressure-scrambling terms with the LES data suggests that the more complex TCL model is superior to the simple linear model only for a few contributions. The linear model is able to reproduce the principal features of the pressure-scrambling terms reasonably well. It can be applied in the second-order turbulence modeling of cloud-topped boundary layer flows, provided some uncertainties are tolerated.Deutscher WetterdienstEuropean Commission/COST Action ES090