Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Sensitivity to Snapshot Frequency in the POD-Based Reduced-Order Modeling of Flow over a Gaussian Bump

The flow over a Gaussian bump at a Reynolds number of 10^6, based on the bump length and freestream conditions, and a Mach number of 0.2, was simulated using wall-modeled large eddy simulation (WMLES). The incoming turbulent flow exhibits relaminarization upstream of the bump apex, due to the favourable pressure gradient acceleration induced by the bump curvature, as well as separation, reattachment and a detached shear layer downstream of the apex. A Galerkin reduced-order model (ROM) was then constructed using a proper orthogonal decomposition (POD) modal basis. The complex, three-dimensional flow features are captured by spatial POD modes differently based on the frequency of the high-fidelity snapshots originally provided. POD modes constructed from differently sampled snapshot matrices were then truncated to retain 90% of the turbulence kinetic energy in the domain of interest. A parametric study is conducted to evaluate the effect of the frequency of snapshots used to construct the POD-ROM.

turbulence↗
NASA NTRS2002

Development of an Efficient Binaural Simulation for the Analysis of Structural Acoustic Data

Binaural or "virtual acoustic" representation has been proposed as a method of analyzing acoustic and vibroacoustic data. Unfortunately, this binaural representation can require extensive computer power to apply the Head Related Transfer Functions (HRTFs) to a large number of sources, as with a vibrating structure. This work focuses on reducing the number of real-time computations required in this binaural analysis through the use of Singular Value Decomposition (SVD) and Equivalent Source Reduction (ESR). The SVD method reduces the complexity of the HRTF computations by breaking the HRTFs into dominant singular values (and vectors). The ESR method reduces the number of sources to be analyzed in real-time computation by replacing sources on the scale of a structural wavelength with sources on the scale of an acoustic wavelength. It is shown that the effectiveness of the SVD and ESR methods improves as the complexity of the source increases. In addition, preliminary auralization tests have shown that the results from both the SVD and ESR methods are indistinguishable from the results found with the exhaustive method.

Lalime, Aimee L.↗
DOE OSTI2023

Model-parallel Fourier neural operators as learned surrogates for large-scale parametric PDEs

Fourier neural operators (FNOs) are a recently introduced neural network architecture for learning solution operators of partial differential equations (PDEs), which have been shown to perform significantly better than comparable deep learning approaches. Once trained, FNOs can achieve speed-ups of multiple orders of magnitude over conventional numerical PDE solvers. However, due to the high dimensionality of their input data and network weights, FNOs have so far only been applied to two-dimensional or small three-dimensional problems. To remove this limited problem-size barrier, we propose a model-parallel version of FNOs based on domain-decomposition of both the input data and network weights. Here, we demonstrate that our model-parallel FNO is able to predict time-varying PDE solutions of over 2.6 billion variables on Perlmutter using up to 512 A100 GPUs and show an example of training a distributed FNO on the Azure cloud for simulating multiphase CO 2 dynamics in the Earth’s subsurface.

58 GEOSCIENCES↗
DOE OSTI2021

Electric vehicle fast charging infrastructure planning in urban networks considering daily travel and charging behavior

Electric vehicles are a sustainable substitution to conventional vehicles. This work introduces an integrated framework for urban fast charging infrastructure to address the range anxiety issue. A mesoscopic simulation tool is developed to generate trip trajectories, and simulate charging behavior based on various trip attributes. The resulting charging demand is the key input to a mixed-integer nonlinear program that seeks charging station configuration. The model minimizes the total system cost including charging station and charger installation costs, and charging, queuing, and detouring delays. The problem is solved using a decomposition technique incorporating a commercial solver for small networks, and a heuristic algorithm for large-scale networks, in addition to the Golden Section method. The solution quality and significant superiority in the computational efficiency of the decomposition approach are confirmed in comparison with the implicit enumeration approach. Furthermore, the required infrastructure to support urban trips is explored for future market shares and technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Scalable freeform optimization of wide-aperture 3D metalenses by zoned discrete axisymmetry

We introduce a novel framework for design and optimization of 3D freeform metalenses that attains nearly linear scaling of computational cost with diameter, by breaking the lens into a sequence of radial “zones” with 𝑛-fold discrete axisymmetry, where 𝑛 increases with radius. This allows vastly more design freedom than imposing continuous axisymmetry, while avoiding the compromises of the locally periodic approximation (LPA) or scalar diffraction theory. Using a GPU-accelerated finite-difference time-domain (FDTD) solver in cylindrical coordinates, we perform full-wave simulation and topology optimization within each supra-wavelength zone. We validate our approach by designing millimeter and centimeter-scale, poly-achromatic, 3D freeform metalenses which outperform the state of the art. By demonstrating the scalability and resulting optical performance enabled by our “zoned discrete axisymmetry” (ZDA) and supra-wavelength domain decomposition, we highlight the potential of our framework to advance large-scale meta-optics and next-generation photonic technologies.

Sun, Mengdi [Wesleyan University]↗
NASA NTRS1987

Evaluation of constraint stabilization procedures for multibody dynamical systems

Comparative numerical studies of four constraint treatment techniques for the simulation of general multibody dynamic systems are presented, and results are presented for the example of a classical crank mechanism and for a simplified version of the seven-link manipulator deployment problem. The staggered stabilization technique (Park, 1986) is found to yield improved accuracy and robustness over Baumgarte's (1972) technique, the singular decomposition technique (Walton and Steeves, 1969), and the penalty technique (Lotstedt, 1979). Furthermore, the staggered stabilization technique offers software modularity, and the only data each solution module needs to exchange with the other is a set of vectors plus a common module to generate the gradient matrix of the constraints, B.

Park, K. C.↗
DOE OSTI2023

Photochemical Decomposition of Y‐Series Non‐Fullerene Acceptors Is Responsible for Degradation of High‐Efficiency Organic Solar Cells

Abstract Organic photovoltaic cells that employ Y‐series non‐fullerene acceptors (NFAs) have recently achieved impressive power‐conversion efficiencies (>18%). To fulfill their commercial promise, it is important to quantify their operational lifetimes and understand their degradation mechanisms. In this work, the spectral‐dependent photostability of films and solar cells comprising several Y‐series acceptors and the donor polymer PM6 is investigated systematically. By applying longpass filters during aging, it is shown that UV/near‐UV photons are responsible for the photochemical decomposition of Y‐series acceptors; this degradation is the primary driver of early solar cell performance losses. Using mass spectrometry, the vinylene linkage between the core and electron‐accepting moieties of Y‐series acceptors is identified as the weak point susceptible to cleavage under UV‐illumination. Employing a series of device characterization, along with numerical simulations, the efficiency losses in organic photovoltaic cells are attributed to the formation of traps, which reduces charge extraction efficiency and facilitates non‐radiative recombination as the Y‐series acceptors degrade. This study provides new insights for molecular degradation of organic photovoltaic absorber materials and highlights the importance of future molecular design and strategies for improved solar cell stability.

14 SOLAR ENERGY↗
DOE OSTI2021

Stochastic pre-event preparation for enhancing resilience of distribution systems

Extreme weather events are the common causes for power supply interruptions and power outages in electrical distribution systems. Improving the distribution system and enhancing its resilience is becoming crucial due to the increased frequency of extreme weather events. Preparation and allocation of multiple flexible resources, such as mobile resources, fuel resources, and labor resources before extreme weather events can mitigate the effects of extreme weather events and enhance the resilience of power distribution systems. Here, in this paper, a two-stage stochastic mixed-integer linear programming (SMILP) is proposed to optimize the preparation and resource allocation process for upcoming extreme weather events, which leads to faster and more efficient post-event restoration. The objective of the proposed two-stage SMILP is to maximize the served load and minimize the operating cost of flexible resources. The first stage in the optimization problem selects the amounts and locations of different resources. The second stage considers the operational constraints of the distribution system and repair crew scheduling constraints. The proposed stochastic pre-event preparation model is solved by a scenario decomposition method, Progressive Hedging (PH), to ease the computational complexity introduced by a large number of scenarios. Furthermore, to show the impact of solar photovoltaic (PV) generation on system resilience, three types of PV systems are considered during a power outage and the resilience improvements with different PV penetration levels are compared. Numerical results from simulations on a large-scale (more than 10,000 nodes) distribution feeder have been used to validate the effectiveness and scalability of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH 4 ) 2

Complex borohydrides such as Mg(BH 4 ) 2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH 4 ) 2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH 4 – –B 3 H 8 – conversion is the formation of the B 2 H 7 – intermediate. We further emphasize and clarify that the B 3 H 8 – and H – intermediates, formed during initial Mg(BH 4 ) 2 decomposition, appear as molecular species that are embedded in the Mg–BH 4 –Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH 2 and Mg(B 3 H 8 ) 2 as previously assumed in theoretical predictions of the thermodynamics.

08 HYDROGEN↗
DOE OSTI2025

High resolution characterization of soil dissolved organic matter with FTICR-MS (Fourier-transform ion cyclotron resonance mass spectrometry) from soil samples in control and warming plots in Blodgett Forest, CA (2014 and 2018)

The pathways of carbon transport and loss through and from soils—soil organic matter (SOM) depolymerization to dissolved organic carbon and mineralization to carbon dioxide (CO2)—are fundamentally driven by microbial activity, which is strongly regulated by environmental conditions. As part of Lawrence Berkeley National Laboratory (LBNL) Terrestrial Ecosystem Science (TES) Belowground Biogeochemistry Science Focus Area (SFA), we have established a novel whole-soil long-term warming experiment at the University of California (UC) Blodgett Forest Research Station (Sierra Nevada) in 2014, where we study the role of biogeochemical, microbial and geochemical process interactions in SOM decomposition and stabilization.This package contains Fourier transform ion cyclotron resonance mass spectrometry (21 Tesla FTICR-MS) data measured in negative and positive ionization mode from water and methanol soil extracts. Soil samples were collected in 2014/06/03 and 2018/06/04 from 3 replicated paired plots that had been subjected to experimental warming since June 2014 to simulate a predicted climate change scenario for northern California. The following files are included: (1) fticr_neg_h2oMeoh_data_raw.csv: raw data from combined water (H2O) and methanol (MeOH) extracts in negative ion mode, (2) fticr_neg_h2oMeoh_data_processed.csv: processed data from combined water (H2O) and methanol (MeOH) extracts in negative ion mode, (3) fticr_neg_metadata.csv: metadata for samples/measurements in negative ion mode, (4) fticr_pos_h2oMeoh_data_raw.csv: raw data from combined water (H2O) and methanol (MeOH) extracts in positive ion mode, (5) fticr_pos_h2oMeoh_data_processed.csv: processed data from combined water (H2O) and methanol (MeOH) extracts in positive ion mode, (6) fticr_pos_metadata.csv: metadata for samples/measurements in positive ion mode.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Error-Bounded Learned Scientific Data Compression with Preservation of Derived Quantities

Scientific applications continue to grow and produce extremely large amounts of data, which require efficient compression algorithms for long-term storage. Compression errors in scientific applications can have a deleterious impact on downstream processing. Thus, it is crucial to preserve all the “known” Quantities of Interest (QoI) during compression. To address this issue, most existing approaches guarantee the reconstruction error of the original data or primary data (PD), but cannot directly control the problem of preserving the QoI. In this work, we propose a physics-informed compression technique that is composed of two parts: (i) reduction of the PD with bounded errors and (ii) preservation of the QoI. In the first step, we combine tensor decompositions, autoencoders, product quantizers, and error-bounded lossy compressors to bound the reconstruction error at high levels of compression. In the second step, we use constraint satisfaction post-processing followed by quantization to preserve the QoI. To illustrate the challenges of reducing the reconstruction errors of the PD and QoI, we focus on simulation data generated by a large-scale fusion code, XGC, which can produce tens of petabytes in a single day. The results show that our approach can achieve a high compression amount while accurately preserving the QoI within scientifically acceptable bounds.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2006

Role of Climate Change in Global Predictions of Future Tropospheric Ozone and Aerosols

A unified tropospheric chemistry-aerosol model within the Goddard Institute for Space Studies general circulation model II is applied to simulate an equilibrium CO2-forced climate in the year 2100 to examine the effects of climate change on global distributions of tropospheric ozone and sulfate, nitrate, ammonium, black carbon, primary organic carbon, secondary organic carbon, sea salt, and mineral dust aerosols. The year 2100 CO2 concentration as well as the anthropogenic emissions of ozone precursors and aerosols/aerosol precursors are based on the Intergovernmental Panel on Climate Change Special Report on Emissions Scenarios (SRES) A2. Year 2100 global O3 and aerosol burdens predicted with changes in both climate and emissions are generally 5-20% lower than those simulated with changes in emissions alone; as exceptions, the nitrate burden is 38% lower, and the secondary organic aerosol burden is 17% higher. Although the CO2-driven climate change alone is predicted to reduce the global O3 concentrations over or near populated and biomass burning areas because of slower transport, enhanced biogenic hydrocarbon emissions, decomposition of peroxyacetyl nitrate at higher temperatures, and the increase of O3 production by increased water vapor at high NOx levels. The warmer climate influences aerosol burdens by increasing aerosol wet deposition, altering climate-sensitive emissions, and shifting aerosol thermodynamic equilibrium. Climate change affects the estimates of the year 2100 direct radiative forcing as a result of the climate-induced changes in burdens and different climatological conditions; with full gas-aerosol coupling and accounting for ozone and direct radiative forcings by the O2, sulfate, nitrate, black carbon, and organic carbon are predicted to be +0.93, -0.72, -1.0, +1.26, and -0.56 W m(exp -2), respectively, using present-day climate and year 2100 emissions, while they are predicted to be +0.76, -0.72, 0.74, +0.97, and -0.58 W m(exp -2), respectively, with year 2100 climate and emissions.

Liao, Hong↗
NASA NTRS2006

Zero to Integration in Eight Months, the Dawn Ground Data System Engineering Challenge

The Dawn Project has presented the Ground Data System (GDS) with technical challenges driven by cost and schedule constraints commonly associated with National Aeronautics and Space Administration (NASA) Discovery Projects. The Dawn mission consists of a new and exciting Deep Space partnership among: the Jet Propulsion Laboratory (JPL), manages the project and is responsible for flight operation; Orbital Sciences Corporation (OSC), is the spacecraft builder and is responsible for flight system test and integration; and the University of California, at Los Angeles (UCLA), is responsible for science planning and operations. As a cost-capped mission, one of Dawn's implementation strategies is to leverage from both flight and ground heritage. OSC's ground data system is used for flight system test and integration as part of the flight heritage strategy. Mission operations, however, are to be conducted with JPL's ground system. The system engineering challenge of dealing with two heterogeneous ground systems emerged immediately. During the first technical interchange meeting between the JPL's GDS Team and OSC's Flight Software Team, August 2003, the need to integrate the ground system with the flight software was brought to the table. This need was driven by the project's commitment to enable instrument engineering model integration in a spacecraft simulator environment, for both demonstration and risk mitigation purposes, by April 2004. This paper will describe the system engineering approach that was undertaken by JPL's GDS Team in order to meet the technical challenge within a non-negotiable eight-month schedule. Key to the success was adherence to fundamental systems engineering practices: decomposition of the project request into manageable requirements; integration of multiple ground disciplines and experts into a focused team effort; definition of a structured yet flexible development process; definition of an in-process risk reduction plan; and aggregation of the intermediate products to an integrated final product. In addition, this paper will highlight the role of lessons learned from the integration experience. The lessons learned from an early GDS deployment have served as the foundation for the design and implementation of the Dawn Ground Data System.

Ground Data System (GDS)↗
DOE OSTI2025

Microbially mediated nitrification improves modeled temperate forest responses to declining nitrogen deposition

As nitrogen deposition declines across the US, uncertainty remains in whether temperate forests will continue to sequester carbon. This uncertainty is amplified by ecosystem models that inaccurately capture the microbial mechanisms that drive soil carbon sequestration and nitrogen loss. Further, even soil process models that represent decomposer microbes and better capture soil carbon retention under nitrogen deposition broadly underrepresent microbial nitrogen transformations. To address this limitation, we leveraged three decades of biogeochemical cycling data from a whole-watershed nitrogen fertilization experiment to incorporate microbially-driven nitrification in the FUN-CORPSE (Fixation and Uptake of Nitrogen-Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) model. Our objectives include: 1) Reproducing key ecosystem responses to fertilization, 2) Integrating microbially-explicit nitrification in FUN-CORPSE, and 3) Assessing modeled soil C and N under projected N deposition shifts. FUN-CORPSE accurately represented soil C and streamwater N losses under ambient N deposition and captured how fertilization induced a 25 % decline in plant C cost of N acquisition, reducing decomposition and increasing soil C. Furthermore, with microbially-explicit nitrification, FUN-CORPSE captured the 100 % increase in nitrification rates and the 50 % increase in streamwater nitrate loss under N fertilization. Specifically, incorporating microbial nitrification improved modeled streamwater N leaching from R 2 = 0.01 to R 2 = 0.57. Under declining N deposition, FUN-CORPSE simulates that N losses recover more quickly than soil C pools. However, the predicted return of soil C to pre-fertilized levels suggests that additional C sequestered due to N deposition may be vulnerable to loss over the next century.

Microbial nitrification model↗
DOE OSTI2024

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asynchronous and Load-Balanced Union-Find for Distributed and Parallel Scientific Data Visualization and Analysis

We present a novel distributed union-find algorithm that features asynchronous parallelism and k-d tree based load balancing for scalable visualization and analysis of scientific data. Applications of union-find include level set extraction and critical point tracking, but distributed union-find can suffer from high synchronization costs and imbalanced workloads across parallel processes. In this study, we prove that global synchronizations in existing distributed union-find can be eliminated without changing final results, allowing overlapped communications and computations for scalable processing. We also use a k-d tree decomposition to redistribute inputs, in order to improve workload balancing. We benchmark the scalability of our algorithm with up to 1,024 processes using both synthetic and application data. Here, we demonstrate the use of our algorithm in critical point tracking and super-level set extraction with high-speed imaging experiments and fusion plasma simulations, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling the effects of active wake mixing on wake behavior through large-scale coherent structures

The use of active wake mixing (AWM) to mitigate downstream turbine wakes has created new opportunities for reducing power losses in wind farms. However, many current analytical or semi-empirical wake models do not capture the flow instabilities that are excited through the blade pitch actuation. In this work, we develop a framework, which accounts for the impacts of the large-scale coherent structures and turbulence on the mean flow, for modeling AWM. The framework uses a triple-decomposition approach for the unsteady flow field and models the mean flow and fine-scale turbulence with a parabolized Reynolds-averaged Navier–Stokes (RANS) system. The wave components are modeled using a simplified spatial linear stability formulation that captures the growth and evolution of the coherent structures. Comparisons with high-fidelity large eddy simulations (LESs) of the turbine wakes showed that this framework was able to capture the additional wake mixing and faster wake recovery in the far-wake regions for both the pulse and helix AWM strategies with minimal computational expense. In the near-wake region, some differences are observed in both the RANS velocity profiles and initial growth of the large-scale structures, which may be due to some simplifying assumptions used in the model.

17 WIND ENERGY↗
NASA NTRS1995

Progress in incompressible Navier-Stokes computations for propulsion flows and its dual-use applications

Development of an incompressible Navier-Stokes solution procedure was performed for the analysis of a liquid rocket engine pump components and for the mechanical heart assist devices. The solution procedure for the propulsion systems is applicable to incompressible Navier-Stokes flows in a steadily rotating frame of reference for any general complex configurations. The computer codes were tested on different complex configurations such as liquid rocket engine inducer and impellers. As a spin-off technology from the turbopump component simulations, the flow analysis for an axial heart pump was conducted. The baseline Left Ventricular Assist Device (LVAD) design was improved by adding an inducer geometry by adapting from the liquid rocket engine pump. The time-accurate mode of the incompressible Navier-Stokes code was validated with flapping foil experiment by using different domain decomposition methods. In the flapping foil experiment, two upstream NACA 0025 foils perform high-frequency synchronized motion and generate unsteady flow conditions for a downstream larger stationary foil. Fairly good agreement was obtained between unsteady experimental data and numerical results from two different moving boundary procedures. Incompressible Navier-Stokes code (INS3D) has been extended for heat transfer applications. The temperature equation was written for both forced and natural convection phenomena. Flow in a square duct case was used for the validation of the code in both natural and forced convection.

Kiris, Cetin↗