Massively Parallel Constrained Force Algorithm for MD Simulation of Polymers and Dendrimers
Explore the source record for details and available documents.
Algorithm Parallel Constrained Force↗
SEARCH · Search NASA
Results for “Algorithm Parallel Constrained Force”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Radiative Characteristics of Clouds Embedded in Smoke Derived from Airborne Multiangular Measurements
Clouds in the presence of absorbing aerosols result in their apparent darkening, observed at the top of atmosphere (TOA), which is associated with the radiative effects of aerosol absorption. Owing to the large radiative effect and potential impacts on regional climate, above-cloud aerosols have recently been characterized in multiple satellite-based studies. While satellite data are particularly useful in showing the radiative impact of above-cloud aerosols at the TOA, recent literature indicates large uncertainties in satellite retrievals of above-cloud aerosol optical depth (AOD) and single scattering albedo (SSA), which are among the most important parameters in the assessment of associated radiative effects. In this study, we analyze radiative characteristics of clouds in the presence of wildfire smoke using airborne data primarily from NASA's Cloud Absorption Radiometer, collected during the Arctic Research of the Composition of the Troposphere from Aircraft and Satellites campaign in Canada during the 2008 summer season. We found a strong positive reflectance (R) gradient in the UV-visible (VIS)-near infrared (NIR) spectrum for clouds embedded in dense smoke, as opposed to an (expected) negative gradient for cloud-free smoke and a flat spectrum for smoke-free cloud cover. Several cases of clouds embedded in thick smoke were found, when the aircraft made circular/spiral measurements, which not only allowed the complete characterization of angular distribution of smoke scattering but also provided the vertical distribution of smoke and clouds (within 0.5-5 km). Specifically, the largest darkening by smoke was found in the UV/VIS, with R(sub 0.34 microns) reducing to 0.2 (or 20%), in contrast to 0.8 at NIR wavelengths (e.g., 1.27 microns). The observed darkening is associated with large AODs (0.5-3.0) and moderately low SSA (0.85-0.93 at 0.53 microns), resulting in a significantly large instantaneous aerosol forcing efficiency of 254 +/- 47 W/sq m/τ. Our observations of smoke-cloud radiative interactions were found to be physically consistent with theoretical plane-parallel 1-D and Monte Carlo 3-D radiative transfer calculations, capturing the observed gradient across UV-VIS-NIR. Results from this study offer insights into aerosol-cloud radiative interactions and may help in better constraining satellite retrieval algorithms.
absorbing aerosols↗
Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition
High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.
97 MATHEMATICS AND COMPUTING↗
kynema-fmb [SWR-23-07]
Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.
Sprague, MichaelA.↗