TY - GEN A1 - Ruprecht, Daniel A1 - Klein, Rupert A1 - Majda, Andrew J. T1 - Moisture-Gravity Wave Interactions in a Multiscale Environment N2 - Starting from the conservation laws for mass, momentum and energy together with a three species, bulk microphysic model, a model for the interaction of internal gravity waves and deep convective hot towers is derived by using multiscale asymptotic techniques. From the resulting leading order equations, a closed model is obtained by applying weighted averages to the smallscale hot towers without requiring further closure approximations. The resulting model is an extension of the linear, anelastic equations, into which moisture enters as the area fraction of saturated regions on the microscale with two way coupling between the large and small scale. Moisture reduces the effective stability in the model and defines a potential temperature sourceterm related to the net effect of latent heat release or consumption by microscale up- and downdrafts. The dispersion relation and group velocity of the system is analyzed and moisture is found to have several effects: It reduces energy transport by waves, increases the vertical wavenumber but decreases the slope at which wave packets travel and it introduces a lower horizontal cutoff wavenumber, below which modes turn into evanescent. Further, moisture can cause critical layers. Numerical examples for steadystate and timedependent mountain waves are shown and the effects of moisture on these waves are investigated. T3 - ZIB-Report - 09-21 Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11358 SN - 1438-0064 ER - TY - GEN A1 - Carqué, Gunter A1 - Schmidt, Heiko A1 - Stevens, Bjorn A1 - Klein, Rupert T1 - Plausibility Check of an Asymptotic Column Model for Deep Convective Clouds N2 - By use of asymptotic analysis Carqué et al. [ZIB-Report 08-03] derived an asymptotic column model for deep convective clouds based on the three dimensional compressible flow equations and a bulk microphysics parameterization. In the present study we check the plausibility of the reduced model equations by comparing implications of the model for the scaling of various terms in the governing equations with those extracted from large eddy simulation data based on the computational model UCLA-LES1.1. This code solves an anelastic system of equations with complete droplet based microphysics and LES closures. We observe that the simulation data corroborate the basic assumptions of the asymptotic analysis and the main conclusions implied by the asymptotically reduced model. The code output reflects the scales of space and time: The deep convective clouds show an anisotropic structure where the horizontal scale is considerably narrower than the vertical scale; with a period of about 20 min, from emergence to breakup, the life cycle of one particular deep convective cloud corresponds exactly to the reference time of the reduced model. The characteristic properties of dynamics as predicted by the reduced model are also reflected in the simulation data: The horizontal flow is controlled by the pressure field; the vertical velocity develops freely independent of pressure over the depth of the convective column; the vertical velocity is directly determined by the buoyancy induced by the potential temperature deviation relative to the background stratification. With respect to grid resolution we observe that refining the spatial step size of the equidistant computational grid from 125 m to 62.5 m does not influence the results: Even with the coarser grid the relevant physical phenomena are sufficiently resolved. Somewhat surprisingly, the Coriolis term involving vertical velocity and acting on the horizontal (east-west) velocity component appears at leading order in the asymptotics. Accordingly, we expected to find a nontrivial impact of this Coriolis effect on the horizontal flow velocity components within columns of updrafts. However, switching the term on and off in subsequent simulations did not sizeably affect the results. T3 - ZIB-Report - 08-44 KW - atmospheric moist convection KW - deep convection KW - precipitating clouds KW - asymptotic analysis KW - Large Eddy Simulation Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10967 SN - 1438-0064 ER - TY - GEN A1 - Oevermann, Michael A1 - Scharfenberg, Carsten A1 - Klein, Rupert T1 - A sharp interface finite volume method for elliptic equations on Cartesian grids N2 - We present a second order sharp interface finite volume method for the solution of the three-dimensional poisson equation with variable coefficients on Cartesian grids. In particular, we focus on interface problems with discontinuities in the coefficient, the source term, the solution, and the fluxes across the interface. The method uses standard piecewiese trilinear finite elements for normal cells and a double piecewise trilinear ansatz for the solution on cells intersected by the interface resulting always in a compact 27-point stencil. Singularities associated with vanishing partial volumes of intersected grid cells are removed by a two-term asymptotic approach. In contrast to the 2D method presented by two of the authors in [M.~Oevermann, R.~Klein: A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, J.~Comp.~Phys.~219 (2006)] we use a minimization technique to determine the unknown coefficients of the double trilinear ansatz. This simplifies the treatment of the different cut-cell types and avoids additional special operations for degenerated interface topologies. The resulting set of linear equations has been solved with a BiCGSTAB solver preconditioned with an algebraic multigrid. In various testcases -- including large coefficient ratios and non-smooth interfaces -- the method achieves second order of accuracy in the L_inf and L_2 norm. T3 - ZIB-Report - 08-38 KW - elliptische Gleichungen KW - Finite Volumen KW - eingebettete Ränder KW - variable und diskontinuierliche Koeffizienten KW - diskontinuierliche Lösung KW - Elliptic equations KW - finite volume methods KW - embedded interface KW - variable and discontinuous coefficients KW - discontinuous solution Y1 - 2008 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-10900 SN - 1438-0064 ER - TY - GEN A1 - Oevermann, Michael A1 - Klein, Rupert T1 - A cartesian grid finite volume method for the solution of the Poisson equation with variable coefficients and embedded interfaces N2 - We present a finite volume method for the solution of the two-dimensional Poisson equation $ \nabla\cdot( \beta( {\mbox{\boldmath $x$}}) \nabla u({\mbox{\boldmath $x$}})) = f(\mbox{\boldmath $x$}) $ with variable, discontinuous coefficients and solution discontinuities on irregular domains. The method uses bilinear ansatz functions on Cartesian grids for the solution $u({\mbox{\boldmath $x$})$ resulting in a compact nine-point stencil. The resulting linear problem has been solved with a standard multigrid solver. Singularities associated with vanishing partial volumes of intersected grid cells or the dual bilinear ansatz itself are removed by a two-step asymptotic approach. The method achieves second order of accuracy in the $L^\infty$ and $L^2$ norm. T3 - ZIB-Report - 06-05 KW - Poisson equation KW - finite volume methods KW - embedded interface KW - ariable and discontinuous coefficients KW - discontinuous solution Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-8996 ER - TY - JOUR A1 - Delle Site, Luigi A1 - Krekeler, Christian A1 - Whittaker, John A1 - Agarwal, Animesh A1 - Klein, Rupert A1 - Höfling, Felix T1 - Molecular Dynamics of Open Systems: Construction of a Mean‐Field Particle Reservoir JF - Advanced Theory and Simulations N2 - The simulation of open molecular systems requires explicit or implicit reservoirs of energy and particles. Whereas full atomistic resolution is desired in the region of interest, there is some freedom in the implementation of the reservoirs. Here, a combined, explicit reservoir is constructed by interfacing the atomistic region with regions of point-like, non-interacting particles (tracers) embedded in a thermodynamic mean field. The tracer molecules acquire atomistic resolution upon entering the atomistic region and equilibrate with this environment, while atomistic molecules become tracers governed by an effective mean-field potential after crossing the atomistic boundary. The approach is extensively tested on thermodynamic, structural, and dynamic properties of liquid water. Conceptual and numerical advantages of the procedure as well as new perspectives are highlighted and discussed. Y1 - 2019 U6 - https://doi.org/10.1002/adts.201900014 VL - 2 SP - 1900014 ER - TY - JOUR A1 - Horenko, Illia A1 - Klein, Rupert A1 - Dolaptchiev, S. A1 - Schütte, Christof T1 - Automated Generation of Reduced Stochastic Weather Models I JF - Mult. Mod. Sim. Y1 - 2008 UR - http://publications.imp.fu-berlin.de/15/ U6 - https://doi.org/10.1137/060670535 VL - 6 IS - 4 SP - 1125 EP - 1145 ER - TY - GEN A1 - Schneider, Thomas A1 - Botta, Nicola A1 - Geratz, Karl Josef A1 - Klein, Rupert T1 - Extension of finite volume compressible flow solvers to multi-dimensional, variable density zero Mach number flow N2 - When attempting to compute unsteady, variable density flows at very small or zero Mach number using a standard finite volume compressible flow solver one faces at least the following difficulties: (i) Spatial pressure variations vanish as the Mach number $M \rightarrow 0$, but they do affect the velocity field at leading order; (ii) the resulting spatial homogeneity of the leading order pressure implies an elliptic divergence constraint for the energy flux; (iii) violation of this constraint would crucially affect the transport of mass, thereby disabling a code to properly advect even a constant density distribution. A previous companion paper derived the above observations from a single time - multiple length scale asymptotic analysis for $M \ll 1$, applied to the conservation form of the governing equations and assuming an ideal gas with constant specific heats. The paper then restricted to weakly compressible one-dimensional flows and introduced a semi-implicit extension of a compressible flow solver, designed to handle the interaction of long wavelength acoustics with small scale, large amplitude density fluctuations. In the present paper we concentrate on the limit of zero Mach number for multi-dimensional, variable density flows. The construction of numerical fluxes for all conserved quantities involves: An explicit upwind step (1) yielding predictions for the nonlinear convective flux components. This procedure still neglects the influence of pressure gradients on the convective fluxes during the time step. Suitable corrections are applied in step (2), which guarantees compliance of the convective fluxes with the divergence constraint. This step requires the solution of a Poisson-type equation to obtain the relevant pressure gradients. Step (3), which requires the solution of a second Poisson-type equation, yields the yet unknown (non-convective) pressure contribution to the total flux of momentum. The final, cell centered velocity field exactly satisfies a discrete divergence constraint consistent with the asymptotic limit. Notice that step (1) can be done by any standard finite volume compressible flow solver and that the input to steps (2) and (3) involves solely the fluxes from step (1), but is independent on how these were obtained. Thus, we claim that our approach allows any such solver to be extended to simulate incompressible flows. Extensions to the weakly compressible regime $0 < M \ll 1$, reactive flows and more complex equations of state will be addressed in follow-up publications. T3 - ZIB-Report - SC-98-31 KW - Incompressible variable density flows KW - conservation laws Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3749 ER - TY - GEN A1 - Helling, Christiane A1 - Klein, Rupert A1 - Woitke, Peter A1 - Nowak, Ulrich A1 - Sedlmayr, Erwin T1 - Dust in Brown Dwarfs IV. Dust formation and driven turbulence on mesoscopic scales N2 - Dust formation {in brown dwarf atmospheres} is studied by utilizing a model for driven turbulence in the mesoscopic scale regime. We apply a pseudo-spectral method where waves are created and superimposed {within} a {limited} wavenumber interval. The turbulent kinetic energy distribution follows the Kolmogoroff spectrum which is assumed to be the most likely value. Such superimposed, stochastic waves may occur in a convectively active environment. They cause nucleation fronts and nucleation events and thereby initiate the dust formation process which { continues until} all condensible material is consumed. Small disturbances {are found to} have a large impact on the dust forming system. An initially dust-hostile region, which may originally be optically thin, becomes optically thick in a patchy way showing considerable variations in the dust properties during the formation process. The dust appears in lanes and curls as a result of the interaction with waves, i.e. turbulence, which form larger and larger structures with time. Aiming on a physical understanding of the variability of brown dwarfs, related to structure formation in substellar atmospheres, we work out first necessary criteria for small-scale closure models to be applied in macroscopic simulations of dust forming astrophysical systems. T3 - ZIB-Report - 03-33 KW - substellar atmospheres KW - dust formation KW - turbulence Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7555 ER - TY - GEN A1 - Helling, Christiane A1 - Klein, Rupert A1 - Sedlmayr, Erwin T1 - The multi-scale dust formation in substellar atmospheres N2 - Substellar atmospheres are observed to be irregularly variable for which the formation of dust clouds is the most promising candidate explanation. The atmospheric gas is convectively unstable and, last but not least, colliding convective cells are seen as cause for a turbulent fluid field. Since dust formation depends on the local properties of the fluid, turbulence influences the dust formation process and may even allow the dust formation in an initially dust-hostile gas. A regime-wise investigation of dust forming substellar atmospheric situations reveals that the largest scales are determined by the interplay between gravitational settling and convective replenishment which results in a dust-stratified atmosphere. The regime of small scales is determined by the interaction of turbulent fluctuations. Resulting lane-like and curled dust distributions combine to larger and larger structures. We compile necessary criteria for a subgrid model in the frame of large scale simulations as result of our study on small scale turbulence in dust forming gases. T3 - ZIB-Report - 03-34 KW - dust formation KW - turbuelnce KW - multi-scale problem KW - substellar atmosphere Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7567 ER - TY - GEN A1 - Helling, Christiane A1 - Klein, Rupert A1 - Woitke, Peter A1 - Sedlmayr, Erwin T1 - Dust formation in brown dwarf atmospheres under conditions of driven turbulence N2 - Based on the knowledge gained from direct numerical simulations which are only possible in the microscale regime, a concept of driven turbulence is presented which allows to enter the mesoscopic scale regime. Here, dust formation under stochastic hydro- and thermodynamic conditions is studied: constructively superimposed stochastic waves initiate dust formation by the creation of singular nucleation events. It, hence, results a varying mean grain size and dust density in space and time. The newly formed dust changes the thermodynamic behavior from almost isotherm to adiabatic and chemically depletes the gas phase. T3 - ZIB-Report - 03-07 KW - dust formation KW - turbulence modelling KW - stochastic boundaries Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7294 ER -