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)
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.
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.
Long-lived flow patterns in the atmosphere such as weather fronts, mid-latitude blockings or tropical cyclones often induce extreme weather conditions. As a consequence, their description, detection, and tracking has received increasing attention in recent years. Similar objectives also arise in diverse fields such as turbulence and combustion research, image analysis, and medical diagnostics under the headlines of "feature tracking", "coherent structure detection" or "image registration" - to name just a few. A host of different approaches to addressing the underlying, often very similar, tasks have been developed and successfully used. Here, several typical examples of such approaches are summarized, further developed and applied to meteorological data sets. Common abstract operational steps form the basis for a unifying framework for the specification of "persistent structures" involving the definition of the physical state of a system, the features of interest, and means of measuring their persistence.
Thermodynamic relations at the coupling boundary in adaptive resolution simulations for open systems
(2021)
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.
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.
Power-User und Supercomputer
(1999)
WaveTrain is an open-source software for numerical simulations of chain-like quantum systems with nearest-neighbor (NN) interactions only. The Python package is centered around tensor train (TT, or matrix product) format representations of Hamiltonian operators and (stationary or time-evolving) state vectors. It builds on the Python tensor train toolbox Scikit_tt, which provides efficient construction methods and storage schemes for the TT format. Its solvers for eigenvalue problems and linear differential equations are used in WaveTrain for the time-independent and time-dependent Schrödinger equations, respectively. Employing efficient decompositions to construct low-rank representations, the tensor-train ranks of state vectors are often found to depend only marginally on the chain length N. This results in the computational effort growing only slightly more than linearly with N, thus mitigating the curse of dimensionality. As a complement to the classes for full quantum mechanics, WaveTrain also contains classes for fully classical and mixed quantum–classical (Ehrenfest or mean field) dynamics of bipartite systems. The graphical capabilities allow visualization of quantum dynamics “on the fly,” with a choice of several different representations based on reduced density matrices. Even though developed for treating quasi-one-dimensional excitonic energy transport in molecular solids or conjugated organic polymers, including coupling to phonons, WaveTrain can be used for any kind of chain-like quantum systems, with or without periodic boundary conditions and with NN interactions only. The present work describes version 1.0 of our WaveTrain software, based on version 1.2 of scikit_tt, both of which are freely available from the GitHub platform where they will also be further developed. Moreover, WaveTrain is mirrored at SourceForge, within the framework of the WavePacket project for numerical quantum dynamics. Worked-out demonstration examples with complete input and output, including animated graphics, are available.