Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

The potential impact of scatterometry on oceanography - A wave forecasting case

A series of observing system simulation experiments have been performed in order to assess the potential impact of marine surface wind data on numerical weather prediction. In addition to conventional data, the experiments simulated the time-continuous assimilation of remotely sensed marine surface wind or temperature sounding data. The wind data were fabricated directly for model grid points intercepted by a Seasat-1 scatterometer swath and were assimilated into the lowest active level (945 mb) of the model using a localized successive correction method. It is shown that Seasat wind data can greatly improve numerical weather forecasts due to better definition of specific features. The case of the QE II storm is examined.

Cane, M. A.↗
DOE OSTI2020

An intra-node HPC network architecture with nanosecond-scale photonic switches

We propose a single-stage network architecture for intra-node connectivity that makes use of nanosecond-scale photonic switches. Although buffering at the switch points is of vital importance for complex multi-stage networks, this is not the case for smaller-scale single-stage networks where the end nodes are located only one hop apart. By limiting the buffering to the end points, the proposed architecture manages to minimize the required electro-optic and opto-electronic conversions, leading in this way to both low end-to-end latency and better energy efficiency. Combining these advantages with nanosecond-scale switching times can allow for high-throughput operation even for frequent switch reconfigurations. The performance of the proposed architecture is evaluated via discrete-event simulations for a wide range of synthetic-traffic cases. The simulation results show that high-throughput operation of ≥90% can be achieved even for small message sizes, i.e., 32 KB for all-to-all communication and 2 KB for uniform random traffic, at a data rate of 400 Gb/s and a switch reconfiguration time of ≤72 ns. Moreover, if 100-ns reconfiguration times are achievable as opposed to 150-ns, then for the all-to-all traffic case a 16% and 28% reduction in completion time can be achieved for message sizes of 8 KB and 1 KB, respectively. In a forthcoming era of optically interfaced processors and accelerators, nanosecond-scale photonic switches appear as a highly promising solution for keeping up with the intra-node bandwidth scaling due to their high-bandwidth, low-latency and fast-switching capabilities.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗
DOE OSTI2020

Thermodynamic analysis and optimization of supercritical carbon dioxide Brayton cycles for use with low-grade geothermal heat sources

The transcritical carbon dioxide Rankine cycle has emerged as an alternative for power generation in low-grade heat applications. However, the low heat rejection temperatures required for condensing the carbon dioxide is almost prohibitive for many locations. A simple solution is the adoption of a supercritical carbon dioxide Brayton cycle which does not need such low temperature heat sinks. This power cycle has been poorly studied previously in low-grade heat applications. Therefore, with the aim of determining the performance and feasibility of the supercritical carbon dioxide Brayton cycle from thermodynamic viewpoint, using a low-grade geothermal heat source, this paper presents a comparative thermodynamic analysis between four different supercritical carbon dioxide Brayton cycles. For this purpose, detailed models were developed and solved by coupling the thermodynamic model with the thermal-hydraulic model of a Printed Circuit Heat Exchanger. Subsequently, simulations of a base case, a parametric analysis, and an optimization using Genetic Algorithms, were carried out. Results show that there is a combination of operating conditions that maximizes the electrical output of the system. Besides, some values of the minimum pressure of the cycles can increase substantially the precooler cooling water pumping power. At the optimum design point, using a 20 kg/s geothermal brine stream at 150 °C as heat source with a minimum allowable reinjection temperature of 70 °C, the Intercooled Recuperated Brayton Cycle achieved the highest electric power output, energy and exergy efficiencies, obtaining values of 779.99 kW, 11.51%, and 52.49%, respectively. The other feasible alternatives were the Recuperated Brayton Cycle, the Simple Brayton Cycle and the Intercooled Brayton Cycle, ranked, in that order.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Stochastic AC optimal power flow: A data-driven approach

There is an emerging need for efficient solutions to stochastic AC Optimal Power Flow (AC-OPF) to ensure optimal and reliable grid operations in the presence of increasing demand and generation uncertainty. Herein this paper presents a highly scalable data-driven algorithm for stochastic AC-OPF that has extremely low sample requirement. The novelty behind the algorithm’s performance involves an iterative scenario design approach that merges information regarding constraint violations in the system with data-driven sparse regression. Compared to conventional methods with random scenario sampling, our approach is able to provide feasible operating points for realistic systems with much lower sample requirements. Furthermore, multiple sub-tasks in our approach can be easily paralleled and based on historical data to enhance its performance and application. We demonstrate the computational improvements of our approach through simulations on different test cases in the IEEE PES PGLib-OPF benchmark library.

42 ENGINEERING↗
DOE OSTI2024

A machine-learning-aided data recovery approach for predicting multi-material thermal behaviors in advanced test reactor capsules

Instrumented experiments conducted at test reactors are essential to the deployment of new advanced reactor systems. Designing new experiments and generating data on specific reactor conditions require significant investments in terms of both time and cost. Finite element analysis software can be used to create high-fidelity models of experiment environments in order to support the actual experiments, but computation time remains a concern in terms of applying outcomes to real-time usage of data (e.g., a digital twin [DT]). Here, the present research proposes a machine-learning (ML) aided approach to making temperature and displacement predictions based on the thickness of the outer gas gap on the experimental capsule used for in-pile demonstration of a novel new thermal conductivity probe in the Advanced Test Reactor (ATR). This capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. Gas gaps existed between the fuel and the rodlet, and between the inner and the outer capsule. The learning data pertained to an experimental capsule's radial distributions of temperature and displacement, as obtained based on Abaqus and the physical features. For the first step of ML sequence, the temperature was predicted using three positional parameters. Next, the displacement was predicted using seven additional parameters. Each physical feature was normalized in order to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement with the simulation results in all cases involving interpolation and extrapolation. Furthermore, data similarity enhancement increased the similarity between the training and the target data, thereby increasing the predictive accuracy of the ML models. In certain extrapolation cases involving limited original ML model accuracy, data similarity enhancement and data recovery was able to somewhat improve this accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Testing and signal identification for two-sample high-dimensional covariances via multi-level thresholding

The paper considers testing and signal identification for covariance matrices from two populations of marginally sub-Gaussian distributed. A multi-level thresholding procedure is proposed for testing the equality of two high-dimensional covariance matrices, which is designed to detect sparse and faint differences between the covariances. A novel -statistic composition is developed to establish the asymptotic distribution of the thresholding statistics in conjunction with the matrix blocking and the coupling techniques. It is shown that the proposed test is more powerful than the existing tests in detecting sparse and weak signals in covariances. Multiple testing procedures are constructed to discover different covariances and the sub-groups of variables with different covariance structures between the two populations. The proposed procedures are based on the multi-level thresholding test, which are able to control the false discovery proportion () with high power. Simulation experiments and a case study on the returns of the S&P 500 stocks before and after the COVID-19 pandemic are conducted to demonstrate and compare the utilities of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An Experimental and Numerical Investigation to Characterize the Low-Temperature Heat Release in Stoichiometric and Lean Combustion

This work reports on an experimental and modeling study on the low-temperature heat release (LTHR) characteristics for three RON 90 binary blends (n-heptane blended with isooctane, toluene and ethanol) in a Cooperative Fuel Research (CFR) engine at lean and stoichiometric conditions that are representative of homogeneous charge compression ignition (HCCI) and spark-ignition (SI) end-gas combustion conditions, respectively. An analysis of the end-gas temperature-pressure (T-P) trajectories was performed to identify the intake conditions leading to similar T-P trajectories between the two lambdas for each fuel blend. A heat release analysis was then conducted for the identified cases, where fuel-to-fuel differences in LTHR were identified and found to be sensitive to the operating condition. Simulations were conducted for these cases using a recently updated chemical kinetic model and a 0-D engine model, where good qualitative and reasonable quantitative agreements in LTHR were obtained. Sensitivity analysis was also performed directly on the rates of LTHR, to understand the controlling chemical reactions of LTHR, providing further insights into the fuel-to-fuel differences. The results demonstrate the significant promoting effect of n-heptane on LTHR rates, while inhibiting effects were seen for ethanol and toluene. Also highlighted was the importance of H atom abstraction reactions from the chemistry of each fuel component, which could lead to contradictory fuel behavior depending on the locations of the H site of the abstraction reaction due to the different ensuing pathways for the primary fuel radicals.

LTHR, RON 90 binary blends, sensitivity analysis↗
DOE OSTI2020

Xcompact3D: An open-source framework for solving turbulence problems on a Cartesian mesh

Xcompact3D is a Fortran 90–95 open-source framework designed for fast and accurate simulations of turbulent flows, targeting CPU-based supercomputers. It is an evolution of the flow solver Incompact3D which was initially designed in France in the mid-90’s for serial processors to solve the incompressible Navier–Stokes equations. Incompact3D was then ported to parallel High Performance Computing (HPC) systems in the early 2010’s. Very recently the capabilities of Incompact3D have been extended so that it can now tackle more flow regimes (from incompressible flows to compressible flows at low Mach numbers), resulting in the design of a new user-friendly framework called Xcompact3D. The present manuscript presents an overview of Xcompact3D with a particular focus on its functionalities, its ready-to-run simulations and a few case studies to demonstrate its impact.

17 WIND ENERGY↗
DOE OSTI2022

Simulating a solar parabolic trough collector plant used for industrial process heat using an optimized operating scheme that utilizes flexible heat integration

This study focuses on the modeling and simulation of a novel design and operation of a solar industrial process heat (SIPH) plant that uses a parabolic trough collector system for generating process heat. The SIPH plant incorporates flexible heat integration (FHI) by having two options for heat sinks as well as a flexible collection temperature for the heat transfer fluid. Leveraging the degrees of freedom created by FHI allows an optimized operating scheme to be created which maximizes the solar heat delivered to the industrial processes. Yearly results for a 27 MW-t plant located in Salt Lake City, UT are simulated for a base case and a case utilizing optimized FHI. For the main process temperature case presented, FHI increases the total solar share by 13.0% and reduces the levelized cost of heat of the SIPH plant by 10.4% relative to the base case. Furthermore, three case study days are presented which highlight how FHI maximizes solar share, specifically during mornings, afternoons, and days with intermittent solar conditions. An environmental analysis is also performed showing that the SIPH plant can reduce harmful emissions by 15.4% for the base case and 17.4% for the FHI case compared to a plant that only utilizes a natural gas furnace for process heat. Overall, the work thoroughly explores an effective operating scheme for a SIPH plant, helping to make it more economical and environmentally beneficial.

14 SOLAR ENERGY↗
DOE OSTI2026

Fe- and Ru-H-Mordenites for Polyethylene Upcycling: Insights from Thermo-Catalytic Pyrolysis and DFT Studies

Catalytic pyrolysis of polyethylene (PE) was performed using Fe- and Ru-impregnated H-mordenite (HM) catalysts to produce lighter hydrocarbons. Catalyst surface areas were analyzed by N 2 adsorption−desorption and acidity by NH 3 - TPD studies. Reducibility of metal from H 2 -TPR, experiments indicate multiple oxidation states of Fe due to higher H 2 consumption. The catalyst activity studies monitored by GC− MS analysis highlight the importance of temperature and PE: catalyst ratio in optimizing hydrocarbon selectivity to C 1 −C 4 alkanes and C 2 −C 5 olefins. Studies with polymer-to-catalyst ratios of 1:1 and 1:2 at higher temperature increased PE conversion with stronger acid sites favoring conversion and moderate acid sites improving selectivity toward lighter olefins. While propane dominated at lower temperatures, propene was the main product at 500 °C. Both catalysts exhibited overall similar conversion, while Ru-HM yielded >40% propene selectivity compared to <40% for Fe-HM. Ethylene selectivity exceeded 10% in both cases. Density functional theory simulations using a C 4 surrogate on Ru-HM confirmed that late-stage cracking dominates dehydrogenation, validating the proposed reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Stratospheric Aerosol Injection Can Reduce Risks to Antarctic Ice Loss Depending on Injection Location and Amount

Owing to increasing greenhouse gas emissions, the Antarctic Ice Sheet is vulnerable to rapid ice loss in the upcoming decades and centuries. Here this study examines the effectiveness of using stratospheric aerosol injection (SAI) that minimizes global mean temperature (GMT) change to slow projected 21st century Antarctic ice loss. We simulate 11 different SAI cases which vary by the latitudinal location(s) and the amount(s) of the injection(s) to examine the climatic response near Antarctica in each case as compared to the reference climate at the turn of the last century. We demonstrate that injecting at a single latitude in the northern hemisphere or at the Equator increases Antarctic shelf ocean temperatures pertinent to ice shelf basal melt, while injecting only in the southern hemisphere minimizes this temperature change. We use these results to analyze the results of more complex multi-latitude injection strategies that maintain GMT at or below 1.5°C above the pre-industrial. All these multi-latitude cases will slow Antarctic ice loss relative to the mid-to-late 21st century SSP2-4.5 emissions pathway. Yet, to avoid a GMT threshold estimated by previous studies pertaining to rapid West Antarctic ice loss (1.5°C above the pre-industrial GMT, though large uncertainty), our study suggests SAI would need to cool about 1.0°C below this threshold and predominately inject at low southern hemisphere latitudes (~15°S - 30°S). These results highlight the complexity of factors impacting the Antarctic response to SAI and the critical role of the injection strategy in preventing future ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Statistical Modeling and Analysis of k -Layer Coverage of Two-Dimensional Materials in Inkjet Printing Processes

Two-dimensional layered materials/flakes, also known as crystalline atom-thick layer nanosheets, have recently been receiving great attention in electronics fabrication due to their unique and intriguing properties. The k-layer coverage area (i.e., the area covered by k number of overlapping layers) of the printed flake pattern significantly impacts on the properties of the printed electronics. In this work, we constructed a statistical model to describe the k-layer coverage of randomly distributed two-dimensional materials. A series of results are obtained to provide not only the expectation but also the variance of the coverage area. The boundary effects on the random flakes coverage are also studied. In addition, an approximated statistical testing approach is also developed in this work to detect abnormal coverage patterns. Furthermore, the case studies based on simulated data and real flakes images obtained from the inkjet printing process demonstrate the accuracy and effectiveness of the proposed model and analysis methods.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constructing a Simulation Surrogate with Partially Observed Output

Gaussian process surrogates are a popular alternative to directly using computationally expensive simulation models. When the simulation output consists of many responses, dimension-reduction techniques are often employed to construct these surrogates. However, surrogate methods with dimension reduction generally rely on complete output training data. This article proposes a new Gaussian process surrogate method that permits the use of partially observed output while remaining computationally efficient. The new method involves the imputation of missing values and the adjustment of the covariance matrix used for Gaussian process inference. The resulting surrogate represents the available responses, disregards the missing responses, and provides meaningful uncertainty quantification. In conclusion, the proposed approach is shown to offer sharper inference than alternatives in a simulation study and a case study where an energy density functional model that frequently returns incomplete output is calibrated.

42 ENGINEERING↗
DOE OSTI2025

Personalized Tucker Decomposition: Modeling Commonality and Peculiarity on Tensor Data

In this paper, we propose a personalized Tucker decomposition (perTucker) to address the limitations of traditional tensor decomposition methods in capturing heterogeneity across different datasets. perTucker decomposes tensor data into shared global components and personalized local components. We introduce an order orthogonality assumption and develop a proximal gradient regularized block coordinate descent algorithm guaranteed to converge to a stationary point. The unique and common representations learned by perTucker reveal intrinsic statistical patterns in data and provide valuable information for a wide range of downstream analytics, including anomaly detection, source classification, and clustering. We demonstrate perTucker’s effectiveness through a simulation study and two case studies on solar flare detection and tonnage signal classification.

14 SOLAR ENERGY↗
DOE OSTI2023

Quantum algorithm for the linear Vlasov equation with collisions

The Vlasov equation is a nonlinear partial differential equation that provides a first-principles description of the dynamics of plasmas. Its linear limit is routinely used in plasma physics to investigate plasma oscillations and stability. In this paper, we present a quantum algorithm that simulates the linearized Vlasov equation with and without collisions, in the one-dimensional electrostatic limit. Rather than solving this equation in its native spatial and velocity phase space, we adopt an efficient representation in the dual space yielded by a Fourier-Hermite expansion. For a given simulation time, the Fourier-Hermite representation is exponentially more compact, thus yielding a classical algorithm that can match the performance of a previously proposed quantum algorithm for this problem. Further, this representation results in a system of linear ordinary differential equations (ODEs) which can be solved with well-developed quantum algorithms: a Hamiltonian simulation in the collisionless case, and quantum ODE solvers in the collisional case. In particular, we demonstrate that a quadratic speedup in system size is attainable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Towards a classification of bifurcations in Vlasov equations

Here, we propose a classification of bifurcations of Vlasov equations, based on the strength of the resonance between the unstable mode and the continuous spectrum on the imaginary axis. We then identify and characterize a new type of generic bifurcation where this resonance is weak, but the unstable mode couples with a stable mode and a Casimir invariant of the system to form a size-3 Jordan block. We derive a three-dimensional reduced noncanonical Hamiltonian system describing this bifurcation. Comparison of the reduced dynamics with direct numerical simulations on a test case gives excellent agreement. We finally discuss the relevance of this bifurcation to specific physical situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗