Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamical System Scaling”

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.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-scale dynamics of atmospheric-pressure discharges ignited over liquid electrodes

A multiphase computational model of atmospheric-pressure gas discharges ignited over liquid electrodes is developed. The model takes into account both the liquid electrode motion under applied electric fields and plasma generation in the gas phase. The influence of the applied voltage and liquid properties (density, viscosity, and surface tension) on the liquid surface deformation is analyzed. It is shown that the shape of dynamic cones formed on the liquid surface differs from the static Taylor's cones. The influence of the liquid surface protrusions on gas breakdown dynamics is demonstrated. It is shown that the breakdown develops in two stages: first, a fast ionization wave propagates from the cathode to the anode, and then the cathode sheath collapses due to secondary electron emission from the liquid surface. It is shown that dynamics of the multiphase system containing gas, plasma, and liquid states is characterized by three disparate time scales: the fast electron time scale, the intermediate ion time scale, and the slow liquid dynamics time scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗
DOE OSTI2022

Long-term 3D MHD simulations of black hole accretion discs formed in neutron star mergers

ABSTRACT We examine the long-term evolution of accretion tori around black hole (BH) remnants of compact object mergers involving at least one neutron star, to better understand their contribution to kilonovae and the synthesis of r-process elements. To this end, we modify the unsplit magnetohydrodynamic (MHD) solver in flash 4.5 to work in non-uniform three-dimensional spherical coordinates, enabling more efficient coverage of a large dynamic range in length scales while exploiting symmetries in the system. This modified code is used to perform BH accretion disc simulations that vary the initial magnetic field geometry and disc compactness, utilizing a physical equation of state, a neutrino leakage scheme for emission and absorption, and modelling the BH’s gravity with a pseudo-Newtonian potential. Simulations run for long enough to achieve a radiatively inefficient state in the disc. We find robust mass ejection with both poloidal and toroidal initial field geometries, and suppressed outflow at high disc compactness. With the included physics, we obtain bimodal velocity distributions that trace back to mass ejection by magnetic stresses at early times, and to thermal processes in the radiatively inefficient state at late times. The electron fraction distribution of the disc outflow is broad in all models, and the ejecta geometry follows a characteristic hourglass shape. We test the effect of removing neutrino absorption or nuclear recombination with axisymmetric models, finding ${\sim} 50{{\ \rm per\ cent}}$ less mass ejection and more neutron-rich composition without neutrino absorption, and a subdominant contribution from nuclear recombination. Tests of the MHD and neutrino leakage implementations are included.

Fahlman, Steven (ORCID:0000000348759940)↗
DOE OSTI2023

Determining the biogeochemical transformations of organic matter composition in rivers using molecular signatures

Inland waters are hotspots for biogeochemical activity, but the environmental and biological factors that govern the transformation of organic matter (OM) flowing through them are still poorly constrained. Here we evaluate data from a crowdsourced sampling campaign led by the Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems (WHONDRS) consortium to investigate broad continental-scale trends in OM composition compared to localized events that influence biogeochemical transformations. Samples from two different OM compartments, sediments and surface water, were collected from 97 streams throughout the Northern Hemisphere and analyzed to identify differences in biogeochemical processes involved in OM transformations. By using dimensional reduction techniques, we identified that putative biogeochemical transformations and microbial respiration rates vary across sediment and surface water along river continua independent of latitude (18°N–68°N). In contrast, we reveal small- and large-scale patterns in OM composition related to local (sediment vs. water column) and reach (stream order, latitude) characteristics. These patterns lay the foundation to modeling the linkage between ecological processes and biogeochemical signals. We further showed how spatial, physical, and biogeochemical factors influence the reactivity of the two OM pools in local reaches yet find emergent broad-scale patterns between OM concentrations and stream order. OM processing will likely change as hydrologic flow regimes shift and vertical mixing occurs on different spatial and temporal scales. As our planet continues to warm and the timing and magnitude of surface and subsurface flows shift, understanding changes in OM cycling across hydrologic systems is critical, given the unknown broad-scale responses and consequences for riverine OM.

rivers↗
NASA NTRS1992

Molecular structure of the coalescence of liquid interfaces

When two bodies of liquid merge, their interfaces must also rupture and rearrange into one. Virtually no information is available concerning the small-scale dynamics of this process. Molecular dynamics simulations of coalescence in systems of about 10,000 Lennard-Jones particles have been performed, arranged so as to mimic laboratory experiments on dense liquids. The coalescence event begins when molecules near the boundary of one liquid body thermally fluctuate into the range of attraction of the other, forming a string of mutually attracting molecules. These molecules gradually thicken into a tendril, which continues to thicken as the bodies smoothly combine in a zipper-like merger.

Koplik, Joel↗
NASA NTRS1995

Collisional Time Scales in the Kuiper Disk and Their Implications

We explore the rate of collisions among bodies in the present-day Kuiper Disk as a function of the total mass and population size structure of the disk. We find that collisional evolution is an important evolutionary process in the disk as a whole, and indeed, that it is likely the dominant evolutionary process beyond approx. 42 AU, where dynamical instability time scales exceed the age of the solar system. Two key findings we report from this modeling work are: that unless the disk's population structure is sharply truncated for radii smaller than approx. 1-2 km, collisions between comets and smaller debris are occurring so frequently in the disk, and with high enough velocities, that the small body (i.e., KM-class object) population in the disk has probably developed into a collisional cascade, thereby implying that the Kuiper Disk comets may not all be primordial, and that the rate of collisions of smaller bodies with larger 100 less R less 400 km objects (like 1992QB(sub 1) and its cohorts) is so low that there appears to be a dilemma in explaining how QB(sub 1)s could have grown by binary accretion in the disk as we know it. Given these findings, it appears that either the present-day paradigm for the formation of Kuiper Disk is failed in some fundamental respect, or that the present-day disk is no longer representative of the ancient structure from which it evolved. This in turn suggests the intriguing possibility that the present-day Kuiper Disk evolved through a more erosional stage reminiscent of the disks around the stars Beta Pictorus, alpha PsA, and alpha Lyr.

Stern, S. Alan↗
NASA NTRS2000

High-Level Clouds and Relation to Sea Surface Temperature as Inferred from Japan's GMS Measurements

High-level clouds have a significant impact on the radiation energy budgets and, hence, the climate of the Earth. Convective cloud systems, which are controlled by large-scale thermal and dynamical conditions, propagate rapidly within days. At this time scale, changes of sea surface temperature (SST) are small. Radiances measured by Japan's Geostationary Meteorological Satellite (GMS) are used to study the relation between high-level clouds and SST in the tropical western and central Pacific (30 S-30 N; 130 E-170 W), where the ocean is warm and deep convection is intensive. Twenty months (January 1998 - August, 1999) of GMS data are used, which cover the second half of the strong 1997-1998 El Nino. Brightness temperature at the 11-micron channel is used to identify high-level clouds. The core of convection is identified based on the difference in the brightness temperatures of the 11- and 12-micron channels. Because of the rapid movement of clouds, there is little correlation between clouds six hours apart. When most of deep convection moves to regions of high SST, the domain averaged high-level cloud amount decreases. A +2C change of SST in cloudy regions results in a relative change of -30% in high-level cloud amount. This large change in cloud amount is due to clouds moving from cool regions to warm regions but not the change in SST itself. A reduction in high-level cloud amount in the equatorial region implies an expanded dry upper troposphere in the off-equatorial region, and the greenhouse warming of high clouds and water vapor is reduced through enhanced longwave cooling to space. The results are important for understanding the physical processes relating SST, convection, and water vapor in the tropics. They are also important for validating climate simulations using global general circulation models.

Chou, Ming-Dah↗
NASA NTRS1999

Directional Solidification of a Binary Alloy into a Cellular Convective Flow: Localized Morphologies

A steady, two dimensional cellular convection modifies the morphological instability of a binary alloy that undergoes directional solidification. When the convection wavelength is far longer than that of the morphological cells, the behavior of the moving front is described by a slow, spatial-temporal dynamics obtained through a multiple-scale analysis. The resulting system has a "parametric-excitation" structure in space, with complex parameters characterizing the interactions between flow, solute diffusion, and rejection. The convection stabilizes two dimensional disturbances oriented with the flow, but destabilizes three dimensional disturbances in general. When the flow is weak, the morphological instability behaves incommensurably to the flow wavelength, but becomes quantized and forced to fit into the flow-box as the flow gets stronger. At large flow magnitudes the instability is localized, confined in narrow envelopes with cells traveling with the flow. In this case the solutions are discrete eigenstates in an unbounded space. Their stability boundary and asymptotics are obtained by the WKB analysis.

Chen, Y.- J.↗
NASA NTRS2009

Use of Patterned CNT Arrays for Display Purposes

Method and system for providing a dynamically reconfigurable display having nanometer-scale resolution, using a patterned array of multi-wall carbon nanotube (MWCNT) clusters. A diode, phosphor or other light source on each MWCNT cluster is independently activated, and different color light sources (e.g., red, green, blue, grey scale, infrared) can be mixed if desired. Resolution is estimated to be 40-100 nm, and reconfiguration time for each MWCNT cluster is no greater than one microsecond.

Delzeit, Lance D.↗
NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗
DOE OSTI2023

Understanding the role of morphology and environment in the dynamical evolution of isolated galaxy triplets

The environments where galaxies reside affects their evolutionary histories. Galaxy triplets (systems composed of three physically bound galaxies) are one of simplest groups of galaxies, and are therefore excellent laboratories to study evolutionary mechanisms where effects of the environment are minimal. We present a statistical study of the dynamical properties of isolated galaxy triplets as a function of their local and large-scale environments. To explore the connection of the dynamical evolution on the systems with the evolution of the galaxies composing the triplets, we consider observational properties such as morphology and star formation rate (SFR). We used the SDSS-based catalogue of Isolated Triplets (SIT), which contains 315 triplets. We classified each triplet according to galaxy morphologies and defined a parameter Q trip to quantify the total local tidal strengths in the systems. To quantify the dynamical stage of the system we used the parameters of harmonic radius R H , velocity dispersion σvr, crossing time H 0 t c , and virial mass M vir . Triplets composed of three early-type galaxies present the smallest R H , indicating that they are in general more compact than triplets with one or more late-type galaxies. Among triplets with low values of R H and H 0 t c , SIT triplets with Q trip < –2 are relaxed systems that are more dynamically evolved, while triplets with Q trip > –2 show compact configurations due to interactions within the system, such as ongoing mergers. We found that there is no dominant galaxy in triplets in terms of properties of stellar populations such as global colour and SFR. Moreover, the global SFR in isolated triplets composed of two or more early-type galaxies increases with the stellar mass ratio of the galaxies with respect to the central galaxy, therefore the system is globally ‘rejuvenated’.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Improving collaborations between empiricists and modelers to advance grassland community dynamics in ecosystem models

Climate change, increasing atmospheric CO 2 , and land use change have altered biogeochemical and hydrologic cycles world-wide, with grassland systems being particularly vulnerable to resulting vegetation shifts (Komatsu et al., 2019). Therefore, incorporating plant community dynamics into ecosystem models is critical for accurate forecasting of ecosystem responses to global change (Levine, 2016). Process-based ecosystem models, which simulate the biogeochemical transfers of mass and energy among biota, the subsurface, and atmosphere, require representation of dynamic composition of organisms within ecosystems. For example, these models simulate leaf and plant-level characteristics, such as electron transport rate and allometry of carbon (C) allocation, to predict how net primary productivity and other ecosystem processes respond to abiotic drivers. These models are particularly useful in scaling from organismal to ecosystem levels but are still underdeveloped in their ability to capture community change, especially in grassland ecosystems. To represent compositional changes, these models must simulate competition, mortality, establishment, and reproduction of plant populations within communities. Yet, current ecosystem modeling approaches to forecast plant community change have derived from studies of forested systems and are either too coarse to capture fine-scale community dynamics (e.g. dynamic global vegetation models (DGVMs)) or too complex to be used at large spatial scales (e.g. forest gap models).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Measurement-informed Dynamic Aggregation of Distribution Systems

This paper proposes a measurement-informed dynamic aggregation methodology in order to create equivalent representations of distribution systems that are compatible with large-scale transmission analysis. By optimizing an equivalent feeder parameters using time-series measurements of active power, reactive power, and voltage at the Point of Interconnection (POI), the approach yields simplified yet dynamically accurate equivalents. Implemented in PSCAD with models of photovoltaic–battery systems, three-phase motors, and static loads, the method employs hybrid differential evolution and bounded least-squares optimization laying the foundation for for real-time state estimation and optimized sensor placement in distribution networks.

Ahmed, Kazi Ishrak [University of Tennessee, Knoxv↗
DOE OSTI2020

CP2K: An Electronic Structure and Molecular Dynamics Software Package - Quickstep: Efficient and Accurate Electronic Structure Calculations

CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2K to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension. TDK has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No. 716142). VRR has been supported by the Swiss National Science Foundation in the form of Ambizione grant No. PZ00P2 174227 and RZK by the Natural Sciences and Engineering Research Council of Canada (NSERC) through Discovery Grants (RGPIN-2016-0505). GKS and CJM are supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. UK based work was funded under the embedded CSE programme of the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), grants eCSE03-011, eCSE06-6, eCSE08-9, eCSE13-17 and the EPSRC (EP/P022235/1) grant “Surface and Interface Toolkit for the Materials Chemistry Community". Computational resources were provided by the Swiss National Supercomputing Centre (CSCS) and Compute Canada. The generous allocation of computing time on the FPGA-based supercomputer “Noctua" at PC2 is kindly acknowledged.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2014

Complexity and Pilot Workload Metrics for the Evaluation of Adaptive Flight Controls on a Full Scale Piloted Aircraft

Flight research has shown the effectiveness of adaptive flight controls for improving aircraft safety and performance in the presence of uncertainties. The National Aeronautics and Space Administration's (NASA)'s Integrated Resilient Aircraft Control (IRAC) project designed and conducted a series of flight experiments to study the impact of variations in adaptive controller design complexity on performance and handling qualities. A novel complexity metric was devised to compare the degrees of simplicity achieved in three variations of a model reference adaptive controller (MRAC) for NASA's F-18 (McDonnell Douglas, now The Boeing Company, Chicago, Illinois) Full-Scale Advanced Systems Testbed (Gen-2A) aircraft. The complexity measures of these controllers are also compared to that of an earlier MRAC design for NASA's Intelligent Flight Control System (IFCS) project and flown on a highly modified F-15 aircraft (McDonnell Douglas, now The Boeing Company, Chicago, Illinois). Pilot comments during the IRAC research flights pointed to the importance of workload on handling qualities ratings for failure and damage scenarios. Modifications to existing pilot aggressiveness and duty cycle metrics are presented and applied to the IRAC controllers. Finally, while adaptive controllers may alleviate the effects of failures or damage on an aircraft's handling qualities, they also have the potential to introduce annoying changes to the flight dynamics or to the operation of aircraft systems. A nuisance rating scale is presented for the categorization of nuisance side-effects of adaptive controllers.

pilot workload metrics↗
NASA NTRS1975

Comparison of human driver dynamics in simulators with complex and simple visual displays and in an automobile on the road

As part of a comprehensive program exploring driver/vehicle system response in lateral steering tasks, driver/vehicle system describing functions and other dynamic data have been gathered in several milieu. These include a simple fixed base simulator with an elementary roadway delineation only display; a fixed base statically operating automobile with a terrain model based, wide angle projection system display; and a full scale moving base automobile operating on the road. Dynamic data with the two fixed base simulators compared favorably, implying that the impoverished visual scene, lack of engine noise, and simplified steering wheel feel characteristics in the simple simulator did not induce significant driver dynamic behavior variations. The fixed base vs. moving base comparisons showed substantially greater crossover frequencies and phase margins on the road course.

Mcruer, D. T.↗