Refine
Year of publication
Document Type
- ZIB-Report (12)
- Article (10)
Is part of the Bibliography
- no (22)
Keywords
- dust formation (3)
- asymptotic analysis (2)
- atmospheric moist convection (2)
- deep convection (2)
- discontinuous solution (2)
- embedded interface (2)
- finite volume methods (2)
- precipitating clouds (2)
- projection method (2)
- Benutzertreffen (1)
Institute
- ZIB Allgemein (9)
- Distributed Algorithms and Supercomputing (4)
- Numerical Mathematics (4)
- AI in Society, Science, and Technology (2)
- Visual and Data-centric Computing (2)
- Computational Systems Biology (1)
- Modeling and Simulation of Complex Processes (1)
- Parallel and Distributed Computing (1)
- Visual Data Analysis (1)
Asymptotic analyses of the three dimensional compressible flow equations coupled with transport equations for the mixing ratios of water vapour, cloud water and rain water are described. We obtain reduced systems of equations for two particular regimes of length and time scales: Models for the long time evolution of deep convective columns and for the short time evolution of shallow convective layers. The asymptotic deep convective column model is anelastic, yet the vertical motion is pressure free, i.e., it evolves freely in interaction with buoyancy while the horizontal divergence adjusts to fullfil the anelastic constraint. The perturbation pressure guaranteeing compliance with the horizontal divergence constraint obeys a Poisson-type equation. Surprisingly, the vertical velocity plays an important role in the horizontal dynamics through the Coriolis term. The vertical acceleration in a saturated column is directly determined by the buoyancy induced by potential temperature differences relative to the background stratification. This potential temperature deviation is a conserved quantity. Evaporation is the only important microphysical process in the undersaturated regime. The evaporation rate depends on the saturation deficit and the amount of rain water present and determines the (downward) vertical velocity and the distribution of water vapour. To connect the deep convective column solutions to top and bottom boundary conditions, a different flow regime needs to be accounted for. Within shallow layers whose depth is comparable to the column diameters, adjustment to physical boundary conditions can take place. This is the second regime considered in this report. The shallow convective layer regime is shown to be asymptotically described by Boussinesq-type equations. These equations are closed by evolution equations which show that, in the saturated regime, the distributions of potential temperature and cloud water are determined by a condensation rate that is directly proportional to the vertical velocity. In the undersaturated regime, the potential temperature distribution is determined by the amount of rain present, since the water vapour in this case is shown to be a conserved quantity. In both regimes the distribution of rain water depends on the rain water flux.
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.
The chemical diffusion master equation (CDME) describes the probabilistic dynamics of reaction--diffusion systems at the molecular level [del Razo et al., Lett. Math. Phys. 112:49, 2022]; it can be considered the master equation for reaction--diffusion processes. The CDME consists of an infinite ordered family of Fokker--Planck equations, where each level of the ordered family corresponds to a certain number of particles and each particle represents a molecule. The equations at each level describe the spatial diffusion of the corresponding set of particles, and they are coupled to each other via reaction operators --linear operators representing chemical reactions. These operators change the number of particles in the system, and thus transport probability between different levels in the family. In this work, we present three approaches to formulate the CDME and show the relations between them. We further deduce the non-trivial combinatorial factors contained in the reaction operators, and we elucidate the relation to the original formulation of the CDME, which is based on creation and annihilation operators acting on many-particle probability density functions. Finally we discuss applications to multiscale simulations of biochemical systems among other future prospects.
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.
Power-User und Supercomputer
(1999)
Päschke et al. (J Fluid Mech, 2012) studied the nonlinear dynamics of strongly tilted vortices subject to asymmetric diabatic heating by asymptotic methods. They found, inter alia, that an azimuthal Fourier mode 1 heating pattern can intensify or attenuate such a vortex depending on the relative orientation of the tilt and the heating asymmetries. The theory originally addressed the gradient wind regime which, asymptotically speaking, corresponds to vortex Rossby numbers of order unity in the limit. Formally, this restricts the applicability of the theory to rather weak vortices. It is shown below that said theory is, in contrast, uniformly valid for vanishing Coriolis parameter and thus applicable to vortices up to low hurricane strengths. An extended discussion of the asymptotics as regards their physical interpretation and their implications for the overall vortex dynamics is also provided in this context. The paper’s second contribution is a series of three-dimensional numerical simulations examining the effect of different orientations of dipolar diabatic heating on idealized tropical cyclones. Comparisons with numerical solutions of the asymptotic equations yield evidence that supports the original theoretical predictions of Päschke et al. In addition, the influence of asymmetric diabatic heating on the time evolution of the vortex centerline is further analyzed, and a steering mechanism that depends on the orientation of the heating dipole is revealed. Finally, the steering mechanism is traced back to the correlation of dipolar perturbations of potential temperature, induced by the vortex tilt, and vertical velocity, for which diabatic heating not necessarily needs to be responsible, but which may have other origins.
We consider the theoretical model of Bergmann and Lebowitz for open systems out of equilibrium and translate its principles in the adaptive resolution simulation molecular dynamics technique. We simulate Lennard-Jones fluids with open boundaries in a thermal gradient and find excellent agreement of the stationary responses with the results obtained from the simulation of a larger locally forced closed system. The encouraging results pave the way for a computational treatment of open systems far from equilibrium framed in a well-established theoretical model that avoids possible numerical artifacts and physical misinterpretations.
We demonstrate how to apply the tensor-train format to solve the time-independent Schrödinger equation for quasi-one-dimensional excitonic chain systems with and without periodic boundary conditions. The coupled excitons and phonons are modeled by Fröhlich–Holstein type Hamiltonians with on-site and nearest-neighbor interactions only. We reduce the memory consumption as well as the computational costs significantly by employing efficient decompositions to construct low-rank tensor-train representations, thus mitigating the curse of dimensionality. In order to compute also higher quantum states, we introduce an approach that directly incorporates the Wielandt deflation technique into the alternating linear scheme for the solution of eigenproblems. Besides systems with coupled excitons and phonons, we also investigate uncoupled problems for which (semi-)analytical results exist. There, we find that in the case of homogeneous systems, the tensor-train ranks of state vectors only marginally depend on the chain length, which results in a linear growth of the storage consumption. However, the central processing unit time increases slightly faster with the chain length than the storage consumption because the alternating linear scheme adopted in our work requires more iterations to achieve convergence for longer chains and a given rank. Finally, we demonstrate that the tensor-train approach to the quantum treatment of coupled excitons and phonons makes it possible to directly tackle the phenomenon of mutual self-trapping. We are able to confirm the main results of the Davydov theory, i.e., the dependence of the wave packet width and the corresponding stabilization energy on the exciton–phonon coupling strength, although only for a certain range of that parameter. In future work, our approach will allow calculations also beyond the validity regime of that theory and/or beyond the restrictions of the Fröhlich–Holstein type Hamiltonians.
Thermodynamic relations at the coupling boundary in adaptive resolution simulations for open systems
(2021)
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.