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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

HEC Applications on Columbia Project

NASA's Columbia system consists of a cluster of twenty 512 processor SGI Altix systems. Each of these systems is 3 TFLOP/s in peak performance - approximately the same as the entire compute capability at NAS just one year ago. Each 512p system is a single system image machine with one Linunx O5, one high performance file system, and one globally shared memory. The NAS Terascale Applications Group (TAG) is chartered to assist in scaling NASA's mission critical codes to at least 512p in order to significantly improve emergency response during flight operations, as well as provide significant improvements in the codes. and rate of scientific discovery across the scientifc disciplines within NASA's Missions. Recent accomplishments are 4x improvements to codes in the ocean modeling community, 10x performance improvements in a number of computational fluid dynamics codes used in aero-vehicle design, and 5x improvements in a number of space science codes dealing in extreme physics. The TAG group will continue its scaling work to 2048p and beyond (10240 cpus) as the Columbia system becomes fully operational and the upgrades to the SGI NUMAlink memory fabric are in place. The NUMlink uprades dramatically improve system scalability for a single application. These upgrades will allow a number of codes to execute faster at higher fidelity than ever before on any other system, thus increasing the rate of scientific discovery even further

Taft, Jim↗
DOE OSTI2025

Experimental analysis of gas dynamics in the reactor cavity section of a High-Temperature Gas-Cooled Reactor during an accident scenario

Here, this study investigates the behavior of a High-Temperature Gas-Cooled Reactor (HTGR) during a break of the helium pressure boundary event. Understanding the gas dynamics that occur during a break could help develop safety systems that reduce the probability of air entering the reactor core. To achieve this objective, experimental studies were conducted in a scaled-down HTGR experimental facility to provide information on the air-helium gas mixture within the confinement building and the impact of the ventilation system on the helium/air gas concentration following the depressurization event. The study evaluated different configurations, including active ventilation time scales, break sizes, and locations. The analysis of the active ventilation time revealed that a ventilation time of 22 seconds produced positive results, and smaller break sizes resulted in longer depressurization times, which improved the ventilation process. The orientation and elevation of the breaks had little effect on the oxygen concentration in the cavity but did impact the gas velocities.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration in Relativistic Turbulent Magnetic Reconnection

In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. Additionally, we find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. Additionally, the main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2025

A neural master equation framework for multiscale modeling of molecular processes: application to atomic-scale plasma processes

Plasma-surface interactions (PSI) play a crucial role in microelectronics fabrication; however, their multiscale nature and array of complex, often unknown interactions make computational modeling of PSIs extremely difficult. To this end, we propose a general neural master equation (NME) framework that uses master equations to describe the dynamics of a molecular process, wherein neural networks learned from atomistic simulations represent unknown transitions between different system states. By leveraging the physics-based structure of master equations and data-driven state transitions, the NME framework promotes generalizability and physics interpretability, and can bridge disparate length and time scales. The framework is demonstrated for multiscale modeling of Si atomic layer etching and reactive ion etching, where the learned NME-based surface kinetic models exhibit good predictive and extrapolative capabilities for predicting experimentally relevant observables as a function of process parameters. The NME-based surface kinetic models obey physical constraints, which are violated in models based on neural ordinary differential equations. The proposed NME framework for multiscale modeling of molecular processes can pave the way for the discovery of new chemistries and materials in atomic-scale plasma processes.

Chemical engineering↗
NASA NTRS2012

Dynamically Reconfigurable Systolic Array Accelerator

A polymorphic systolic array framework has been developed that works in conjunction with an embedded microprocessor on a field-programmable gate array (FPGA), which allows for dynamic and complimentary scaling of acceleration levels of two algorithms active concurrently on the FPGA. Use is made of systolic arrays and a hardware-software co-design to obtain an efficient multi-application acceleration system. The flexible and simple framework allows hosting of a broader range of algorithms, and is extendable to more complex applications in the area of aerospace embedded systems. FPGA chips can be responsive to realtime demands for changing applications needs, but only if the electronic fabric can respond fast enough. This systolic array framework allows for rapid partial and dynamic reconfiguration of the chip in response to the real-time needs of scalability, and adaptability of executables.

Dasu, Aravind↗
NASA NTRS2012

Coupling Mars' Dust and Water Cycles: Effects on Dust Lifting Vigor, Spatial Extent and Seasonality

The dust cycle is an important component of Mars' current climate system. Airborne dust affects the radiative balance of the atmosphere, thus greatly influencing the thermal and dynamical state of the atmosphere. Dust raising events on Mars occur at spatial scales ranging from meters to planet-wide. Although the occurrence and season of large regional and global dust storms are highly variable from one year to the next, there are many features of the dust cycle that occur year after year. Generally, a low-level dust haze is maintained during northern spring and summer, while elevated levels of atmospheric dust occur during northern autumn and winter. During years without global-scale dust storms, two peaks in total dust loading were observed by MGS/TES: one peak occurred before northern winter solstice at Ls 200-240, and one peak occurred after northern winter solstice at L(sub s) 305-340. These maxima in dust loading are thought to be associated with transient eddy activity in the northern hemisphere, which has been observed to maximize pre- and post-solstice. Interactive dust cycle studies with Mars General Circulation Models (MGCMs) have included the lifting, transport, and sedimentation of radiatively active dust. Although the predicted global dust loadings from these simulations capture some aspects of the observed dust cycle, there are marked differences between the simulated and observed dust cycles. Most notably, the maximum dust loading is robustly predicted by models to occur near northern winter solstice and is due to dust lifting associated with down slope flows on the flanks of the Hellas basin. Thus far, models have had difficulty simulating the observed pre- and post- solstice peaks in dust loading. Interactive dust cycle studies typically have not included the formation of water ice clouds or their radiative effects. Water ice clouds can influence the dust cycle by scavenging dust from atmosphere and by interacting with solar and infrared radiation, thereby modifying the thermal structure of the atmosphere and its circulation. Results presented in other papers at this workshop show that including the radiative effects of water ice clouds greatly influence the water cycle and the vigor of weather systems in both the northern and southern hemispheres. Our goal is to investigate the effects of fully coupling the dust and water cycles on the dust cycle. We show that including water ice clouds and their radiative effects greatly affect the magnitude, spatial extent and seasonality of dust lifting and the season of maximum atmospheric dust loading.

Kahre, M. A.↗
DOE OSTI2021

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Budget-Constrained Sizing of Renewable and Energy Storage Systems for Farm-Scale Ammonia Production Within the Food-Energy-Water Nexus

In the transition toward sustainable agriculture, farms have emerged as eco-friendly pioneers, harnessing clean hybrid wind and solar systems to improve farm performance. A concern in this paradigm is the effective sizing of renewable energy systems to ensure optimal energy use within budget considerations. This research focuses on optimizing renewable energy sizing in small-scale ammonia production to meet specific farm demands and enhance local resilience, emphasizing the interplay between environmental and economic factors. These findings promise increased energy efficiency and sustainability in this innovative agricultural sector. Additionally, our approach considers small-scale ammonia plant needs and the dynamic relationships between ammonia, water, and farm demands. Simulations demonstrate substantial cost savings in farm electricity consumption. Specifically, scenarios with renewable energy integration in the farm can reduce at least 13% electricity cost compared to a grid-dependent system in the 15-year simulation.

ammonia↗
NASA NTRS1982

Methods of sequential estimation for determining initial data in numerical weather prediction

Numerical weather prediction (NWP) is an initial-value problem for a system of nonlinear differential equations, in which initial values are known incompletely and inaccurately. Observational data available at the initial time must therefore be supplemented by data available prior to the initial time, a problem known as meteorological data assimilation. A further complication in NWP is that solutions of the governing equations evolve on two different time scales, a fast one and a slow one, whereas fast scale motions in the atmosphere are not reliably observed. This leads to the so called initialization problem: initial values must be constrained to result in a slowly evolving forecast. The theory of estimation of stochastic dynamic systems provides a natural approach to such problems. For linear stochastic dynamic models, the Kalman-Bucy (KB) sequential filter is the optimal data assimilation method, for linear models, the optimal combined data assimilation-initialization method is a modified version of the KB filter.

Cohn, S. E.↗
DOE OSTI2020

Fast Magnetic Reconnection with Turbulence in High Lundquist Number Limit

Here, we use extensive 3D resistive MHD simulations to study how large-scale current sheets will undergo fast reconnection in the high Lundquist number S limit (above ~104), when the system is subject to different externally driven turbulence levels and the self-generated turbulence produced by 3D reconnection dynamics. We find that the normalized global reconnection rate ~0.01–0.13 is weakly dependent on S. Global reconnection with the classic inflow/outflow configurations is observed, and 3D flux ropes are hierarchically formed and ejected from reconnection regions. A statistical separation of the reconnected magnetic field lines follows a superdiffusive behavior, from which the rate is measured to be very similar to that obtained from the mixing of tracer populations. We find that the reconnection rate scales roughly linearly with the turbulence level during the peak of reconnection. This scaling is consistent with the turbulence properties produced by both the externally driven and self-generation processes. These results imply that large-scale thin current sheets tend to undergo rigorous reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗
DOE OSTI2021

Hamiltonian reduction of Vlasov–Maxwell to a dark slow manifold

In this work, we show that non-relativistic scaling of the collisionless Vlasov–Maxwell system implies the existence of a formal invariant slow manifold in the infinite-dimensional Vlasov–Maxwell phase space. Vlasov–Maxwell dynamics restricted to the slow manifold recovers the Vlasov–Poisson and Vlasov–Darwin models as low-order approximations, and provides higher-order corrections to the Vlasov–Darwin model more generally. The slow manifold may be interpreted to all orders in perturbation theory as a collection of formal Vlasov–Maxwell solutions that do not excite light waves, and are therefore ‘dark’. We provide a heuristic lower bound for the time interval over which Vlasov–Maxwell solutions initialized optimally near the slow manifold remain dark. We also show how the dynamics on the slow manifold naturally inherits a Hamiltonian structure from the underlying system. After expressing this structure in a simple form, we use it to identify a manifestly Hamiltonian correction to the Vlasov–Darwin model. The derivation of higher-order terms is reduced to computing the corrections of the system Hamiltonian restricted to the slow manifold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Wind tunnel test of a fighter aircraft wing/store flutter suppression system: An international effort

A 30% scale, half span model of a lightweight fighter aircraft with an active wing/store flutter suppression system was tested in the NASA Langley Research Center sixteen foot transonic dynamics tunnel. The test featured a store configuration that was intentionally designed to exhibit a violent flutter condition. In addition to Northrop organized control laws, three European countries also contributed control laws to stabilize this condition. After the control laws were mechanized by Northrop, they were tested at the Langley facility. The model was tested up to 170% of the open loop flutter dynamic pressure in a number of cases, with the indication that a substantially greater improvement was achievable. Some special features of the test model are discussed and the design and implementation of the control laws as well as the test monitoring techniques and results are presented.

Hwang, C.↗
NASA NTRS1992

Linkages between terrestrial ecosystems and the atmosphere

The primary research issue in understanding the role of terrestrial ecosystems in global change is analyzing the coupling between processes with vastly differing rates of change, from photosynthesis to community change. Representing this coupling in models is the central challenge to modeling the terrestrial biosphere as part of the earth system. Terrestrial ecosystems participate in climate and in the biogeochemical cycles on several temporal scales. Some of the carbon fixed by photosynthesis is incorporated into plant tissue and is delayed from returning to the atmosphere until it is oxidized by decomposition or fire. This slower (i.e., days to months) carbon loop through the terrestrial component of the carbon cycle, which is matched by cycles of nutrients required by plants and decomposers, affects the increasing trend in atmospheric CO2 concentration and imposes a seasonal cycle on that trend. Moreover, this cycle includes key controls over biogenic trace gas production. The structure of terrestrial ecosystems, which responds on even longer time scales (annual to century), is the integrated response to the biogeochemical and environmental constraints that develop over the intermediate time scale. The loop is closed back to the climate system since it is the structure of ecosystems, including species composition, that sets the terrestrial boundary condition in the climate system through modification of surface roughness, albedo, and, to a great extent, latent heat exchange. These separate temporal scales contain explicit feedback loops which may modify ecosystem dynamics and linkages between ecosystems and the atmosphere. The long-term change in climate, resulting from increased atmospheric concentrations of greenhouse gases (e.g., CO2, CH4, and nitrous oxide (N2O)) will further modify the global environment and potentially induce further ecosystem change. Modeling these interactions requires coupling successional models to biogeochemical models to physiological models that describe the exchange of water, energy, and biogenic trace gases between the vegetation and the atmosphere at fine time scales. There does not appear to be any obvious way to allow direct reciprocal coupling of atmospheric general circulation models (GCM's), which inherently run with fine time steps, to ecosystem or successional models, which have coarse temporal resolution, without the interposition of physiological canopy models. This is equally true for biogeochemical models of the exchange of carbon dioxide and trace gases. This coupling across time scales is nontrivial and sets the focus for the modeling strategy.

Bretherton, Francis↗
NASA NTRSDate not supplied

Development of Buffet Forcing Functions for a Transonic Condition Exhibiting Bimodal Flow Behavior

A wind-tunnel test campaign was used to derive Buffet Forcing Functions (BFFs) of the Space Launch System Block 1B vehicle configuration by measuring and integrating unsteady surface pressures on a 3-percent scale rigid buffet model (RBM) tested at the NASA Langley Transonic Dynamics Tunnel. The model was tested at predetermined and repeatable pitch and side-slip angles and flow Mach numbers encompassing a full range of possible flight conditions. Although each data point was collected at a steady wind-tunnel condition, a transient, bimodal flow behavior was observed in some Mach 1.10 measurements. The pressure time series was alternating between two states, which differed either in mean or fluctuation amplitude values, or both. As a consequence of this behavior, segments of the resultant BFFs can differ based on the duration that each measurement spends in a particular state. In this paper, a methodology is proposed, which envelopes the range of BFFs magnitudes resulting from this random bimodal phenomenon in pressure time histories of certain sensors.

wind-tunnel testing↗
NASA NTRSDate not supplied

Development of Buffet Forcing Functions for a Transonic Condition Exhibiting Bimodal Flow Behavior

A wind-tunnel test campaign was used to derive Buffet Forcing Functions (BFFs) of the Space Launch System Block 1B vehicle configuration by measuring and integrating unsteady surface pressures on a 3-percent scale rigid buffet model (RBM) tested at the NASA Langley Transonic Dynamics Tunnel. The model was tested at predetermined and repeatable pitch and side-slip angles and flow Mach numbers encompassing a full range of possible flight conditions. Although each data point was collected at a steady wind-tunnel condition, a transient, bimodal flow behavior was observed in some Mach 1.10 measurements. The pressure time series was alternating between two states, which differed either in mean or fluctuation amplitude values, or both. As a consequence of this behavior, segments of the resultant BFFs can differ based on the duration that each measurement spends in a particular state. In this paper, a methodology is proposed, which envelopes the range of BFFs magnitudes resulting from this random bimodal phenomenon in pressure time histories of certain sensors.

wind-tunnel testing↗