Refine
Year of publication
- 2016 (14) (remove)
Document Type
- Article (6)
- ZIB-Report (4)
- In Proceedings (3)
- Poster (1)
Language
- English (14)
Is part of the Bibliography
- no (14) (remove)
Keywords
Institute
ORBKIT is a toolbox for postprocessing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wave function, as well as several grid-independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemistry programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user-written modules to determine any other derived quantity. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT offers routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open-source under GNU-LGPLv3 license and freely available at https://github.com/orbkit/orbkit/.
This article provides an overview of ORBKIT with particular focus
on its capabilities and applicability, and includes several example
calculations.
Purpose: To account for the impact of turbulence in blood damage modeling, a novel approach based on the generation of instantaneous flow fields from RANS simulations is proposed.
Methods: Turbulent flow in a bileaflet mechanical heart valve was simulated using RANS-based (SST k-ω) flow solver using FLUENT 14.5. The calculated Reynolds shear stress (RSS) field is transformed into a set of divergence-free random vector fields representing turbulent velocity fluctuations using procedural noise functions. To consider the random path of the blood cells, instantaneous flow fields were computed for each time step by summation of RSS-based divergence-free random and mean velocity fields. Using those instantaneous flow fields, instantaneous pathlines and corresponding point-wise instantaneous shear stresses were calculated. For a comparison, averaged pathlines based on mean velocity field and respective viscous shear stresses together with
RSS values were calculated. Finally, the blood damage index (hemolysis) was integrated along the averaged and instantaneous pathlines using a power law approach and then compared.
Results: Using RSS in blood damage modeling without a correction factor overestimates damaging stress and thus the blood damage (hemolysis). Blood damage histograms based on both presented approaches differ.
Conclusions: A novel approach to calculate blood damage without using RSS as a damaging parameter is established. The results of our numerical experiment support the hypothesis that the use of RSS as a damaging parameter should be avoided.
Traditionally, Lagrangian fields such as finite-time Lyapunov exponents (FTLE)
are precomputed on a discrete grid and are ray casted afterwards. This, however,
introduces both grid discretization errors and sampling errors during ray marching.
In this work, we apply a progressive, view-dependent Monte Carlo-based approach
for the visualization of such Lagrangian fields in time-dependent flows. Our ap-
proach avoids grid discretization and ray marching errors completely, is consistent,
and has a low memory consumption. The system provides noisy previews that con-
verge over time to an accurate high-quality visualization. Compared to traditional
approaches, the proposed system avoids explicitly predefined fieldline seeding
structures, and uses a Monte Carlo sampling strategy named Woodcock tracking
to distribute samples along the view ray. An acceleration of this sampling strategy
requires local upper bounds for the FTLE values, which we progressively acquire
during the rendering. Our approach is tailored for high-quality visualizations of
complex FTLE fields and is guaranteed to faithfully represent detailed ridge surface
structures as indicators for Lagrangian coherent structures (LCS). We demonstrate
the effectiveness of our approach by using a set of analytic test cases and real-world numerical simulations.
Many scientific applications deal with data from a multitude of different sources, e.g., measurements, imaging and simulations. Each source provides an additional perspective on the phenomenon of interest, but also comes with specific limitations, e.g. regarding accuracy, spatial and temporal availability. Effectively combining and analyzing such multimodal and partially incomplete data of limited accuracy in an integrated way is challenging. In this work, we outline an approach for an integrated analysis and visualization of the atmospheric impact of volcano eruptions. The data sets comprise observation and imaging data from satellites as well as results from numerical particle simulations. To analyze the clouds from the volcano eruption in the spatiotemporal domain we apply topological methods. Extremal structures reveal structures in the data that support clustering and comparison. We further discuss the robustness of those methods with respect to different properties of the data and different parameter setups. Finally we outline open challenges for the effective integrated visualization using topological methods.