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 1,045 records · Page 58

SOC Microstructural Property Estimator

This pre-trained ML model is a tool that uses basic compositional parameters for porous solid oxide cell (SOC) electrodes - the phase fractions and mean particle/pore diameters – as inputs and uses them to estimate additional electrochemical performance parameters: active (i.e., connected) TPB density, all tortuosity factors, and phase pair specific interfacial areas. The electrode is assumed to be composed of two solid phases and a pore phase. The property calculations are performed using neural network regression models trained on a large bank of synthetic electrode microstructural data that NETL has generated using the program DREAM3D (that bank is also hosted on EDX: https://edx.netl.doe.gov/dataset/soc-synthetic-microstructure-bank). This means the generated parameters are based on training from actual measured properties from 3D microstructures, not estimated from geometric simplifications. This tool was developed and is intended to replace percolation theory calculations in models that use hypothetical electrode properties. An example use case would be running SOC performance simulations across a parametric sweep of electrode designs (e.g., varying phase fractions and particle sizes) and assessing how it impacts the electrochemical performance of the SOC. Within the parameter space of the training data (statistics of that parameter space is provided in the readme file), this model achieves sub-5% mean absolute percent errors, an order of magnitude less error than percolation theory across the same parameter space. However, be aware that this tool was developed with parametric simulations in mind, and users are encouraged to assess accuracy for their own specific use case rather than taking accuracy metrics at face value. More info, including a usage guide, is in the included readme file. This tool should be cited with the DOI number provided.

Electrode Microstructure↗

Ab initio Molecular Dynamics Beyond Density Functional Theory

The funding from this award supported the development of new theoretical methods that both increase the accuracy of quantum simulations beyond that of density functional theory, as well as reduce the cost of such simulations, for example, through machine learning and algorithmic improvements. The improvements lead to, in some cases, orders of magnitude increases in simulation speed as well as orders of magnitude increase in accuracy, compared to before this work was started. This translates to both greater certainty in making predictions about chemical phenomena that are currently studied computationally, as well as laying the ground-work for the study of new areas of chemistry that have not traditionally been simulated, because of the new time-scales and accuracies that can be reached.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Application of Modeling and Simulation in Capacity Management within the ITIL Framework

Tightly integrating modeling and simulation techniques into Information Technology Infrastructure Library (ITIL) practices can be one of the driving factors behind a successful and cost-effective capacity management effort for any Information Technology (IT) system. ITIL is a best practices framework for managing IT infrastructure, development and operations. Translating ITIL theory into operational reality can be a challenge. This paper aims to highlight how to best integrate modeling and simulation into an ITIL implementation. For cases where the project team initially has difficulty gaining consensus on investing in modeling and simulation resources, a clear definition for M&S implementation into the ITIL framework, specifically its role in supporting Capacity Management, is critical to gaining the support required to garner these resources. This implementation should also help to clearly define M&S support to the overall system mission. This paper will describe the development of an integrated modeling approach and how best to tie M&S to definitive goals for evaluating system capacity and performance requirements. Specifically the paper will discuss best practices for implementing modeling and simulation into ITIL. These practices hinge on implementing integrated M&S methods that 1) encompass at least two or more predictive modeling techniques, 2) complement each one's respective strengths and weaknesses to support the validation of predicted results, and 3) are tied to the system's performance and workload monitoring efforts. How to structure two forms of modeling: statistical and simUlation in the development of "As Is" and "To Be" efforts will be used to exemplify the integrated M&S methods. The paper will show how these methods can better support the project's overall capacity management efforts.

Rahmani, Sonya↗

Application of a Comprehensive Lagrangian–Eulerian Spark-Ignition Model to Different Operating Conditions

Increasing engine efficiency is essential to reducing emissions, which is a priority for automakers. Unconventional modes such as boosted and highly dilute operation have the potential to increase engine efficiency but suffer from stability concerns and cyclic variability. To aid engineers in designing ignition systems that reduce cyclic variability in such engine operation modes, reliable and accurate spark-ignition models are necessary. Here, in this article, a Lagrangian–Eulerian spark-ignition (LESI) model is used to simulate electrical discharge, spark channel elongation, and ignition in inert or reacting crossflow within a combustion vessel, at different pressures, flow speeds, and dilution rates. First the model formulation is briefly revisited. Then, the experimental and simulations setups are presented. The results showcase the model’s ability to predict the secondary circuit voltage, current, and power signals, in addition to the spark channel elongation, for the inert cases, or flame front growth, for the reacting cases. The results also compare simulation spark channel and flame growth plots to experimental Schlieren images at different instants in time. This work serves to highlight LESI’s ability to predict the characteristics of discharge and ignition across a variety of operating conditions.

42 ENGINEERING↗

Equations of motion of slung load systems with results for dual lift

General simulation equations are derived for the rigid body motion of slung load systems. These systems are viewed as consisting of several rigid bodies connected by straight-line cables or links. The suspension can be assumed to be elastic or inelastic, both cases being of interest in simulation and control studies. Equations for the general system are obtained via D'Alembert's principle and the introduction of generalized velocity coordinates. Three forms are obtained. Two of these generalize previous case-specific results for single helicopter systems with elastic or inelastic suspensions. The third is a new formulation for inelastic suspensions. It is derived from the elastic suspension equations by choosing the generalized coordinates so as to separate motion due to cable stretching from motion with invariant cable lengths. The result is computationally more efficient than the conventional formulation, and is readily integrated with the elastic suspension formulation and readily applied to the complex dual lift and multilift systems. Equations are derived for dual lift systems. Three proposed suspension arrangements can be integrated in a single equation set. The equations are given in terms of the natural vectors and matrices of three-dimensional rigid body mechanics and are tractable for both analysis and programming.

Cicolani, Luigi S.↗

Simulation of Range Safety for the NASA Space Shuttle

This paper describes a simulation environment that seamlessly combines a number of safety and environmental models for the launch phase of a NASA Space Shuttle mission. The components of this simulation environment represent the different systems that must interact in order to determine the Expectation of casualties (E(sub c)) resulting from the toxic effects of the gas dispersion that occurs after a disaster affecting a Space Shuttle within 120 seconds of lift-off. The utilization of the Space Shuttle reliability models, trajectory models, weather dissemination systems, population models, amount and type of toxicants, gas dispersion models, human response functions to toxicants, and a geographical information system are all integrated to create this environment. This simulation environment can help safety managers estimate the population at risk in order to plan evacuation, make sheltering decisions, determine the resources required to provide aid and comfort, and mitigate damages in case of a disaster. This simulation environment may also be modified and used for the landing phase of a space vehicle but will not be discussed in this paper.

Rabelo, Luis↗

Numerical Simulation of Pressure-induced Separation of Turbulent Flat-plate Boundary Layers: Definition and Overview of New Cases with Suction-only Transpiration and a Step in Reynolds Number

Results from two new direct numerical simulations (DNSs) are added to the database described in Coleman, Rumsey & Spalart (2018) (henceforth CRS18), and available at the NASA TMR website (https://turbmodels.larc.nasa.gov). The new flows, Cases D and E, are similar to the earlier ones (Cases A–C) in that a transpiration profile above a flat plate creates a prolonged adverse pressure gradient (APG) that drives a canonical turbulent zero-pressure-gradient (ZPG), flat-plate boundary layer towards separation; they differ in that for the new cases (1) the APG is not followed by a favorable gradient (FPG), since only suction transpiration is employed, and (2) the mean shear stress does not cross zero, so that the width of the mean bubble, strictly defined, is essentially zero. However, the probability of reversed skin friction exceeds 49%. The motivation was to produce a flow that is more akin to technological flows, including shock-boundary-layer interaction, with a gradual turbulence-controlled recovery. One of the new simulations is also at a significantly higher Reynolds number than the earlier flows.

Turbulence↗

HERO WEC V1.0 - WEC-Sim Detailed Simulation Runs and Summary Data

Data sets from simulation runs include 144 wave cases that were run based on the WEC-Sim model (https://mhkdr.openei.org/submissions/483). The 144 wave cases represent waves with the following wave height and wave period ranges: - Significant Wave Height: 0.25 - 4.0m in 0.25m increments - Wave Period: 5 - 13 sec in 1 sec increments Each run was simulated using a Pierson-Moskowitz irregular wave spectrum with a 200 second ramp time and a total simulation time of 2,200 seconds. Summary data set includes a spreadsheet and image files with matrices that are associated with data from simulation runs. All matrices cover the same significant wave height and wave periods from the sim runs, in the same increments. The following matrices are included: - Mechanical Power Input: Absorbed power calculated using linear input velocity at PTO (power take-off) multiplied by winch tension. - Mechanical Power Output: The hydraulic power calculated at the pump output prior to the pressure relief valve dump. - Feed Pressure: Average pressure at RO (reverse osmosis) input prior to pre-filters - RO Power Input: The hydraulic power calculated at the RO system inlet, accounts for flow loss due to pressure relief valve - Water Production: Average clean water flow at RO output

16 TIDAL AND WAVE POWER↗

A nonparametric approach to modeling dependence among production parameters for energy cost estimation

Dependence existing among production parameters can significantly affect energy costs. For material removal processes, processing times and material amounts to be cut have positive dependence, and the dependence can affect energy costs by varying a peak electricity load. Existing methods to model the dependence of production parameters are based on parametric approaches that cannot accurately represent the dependence due to the rigid nature of parametric models. To address this issue, this study proposes a method to quantify dependence among manufacturing parameters through the application of the empirical copula to the dependence between the milling processing time and the amount of volume to be cut. The proposed method is illustrated by a case study of a manufacturing facility consisting of milling machines; a total of 27 scenarios are simulated for energy cost estimations. The case study clearly shows that the proposed method can capture the dependence between the milling parameters more accurately than a conventional dependence measure from parametric correlation models. Furthermore, the findings from this study would be useful to estimate energy costs in that the proposed method provides a better fit to real production data than conventional parametric approaches can.

42 ENGINEERING↗

Effective mie-scattering and CO2 absorption in the dust-laden Martian atmosphere and its impact on radiative-convective temperature changes in the lower scale heights

A time dependent computer model of radiative-convective-conductive heat transfer in the Martian ground-atmosphere system was refined by incorporating an intermediate line strength CO2 band absorption which together with the strong-and weak-line approximation closely simulated the radiative transmission through a vertically inhomogeneous stratification. About 33,000 CO2 lines were processed to cover the spectral range of solar and planetary radiation. Absorption by silicate dust particulates, was taken into consideration to study its impact on the ground-atmosphere temperature field as a function of time. This model was subsequently attuned to IRIS, IR-radiometric and S-band occultation data. Satisfactory simulations of the measured IRIS spectra were accomplished for the dust-free condition. In the case of variable dust loads, the simulations were sufficiently fair so that some inferences into the effect of dust on temperature were justified.

Pallmann, A. J.↗

How equation of state selection impacts accuracy near the critical point: Forced convection supercritical CO 2 flow over a cylinder

Increasing application of supercritical fluids leads to increased use of numerical simulations to guide designs. However, the impact of Equations of State (EoS) selection on computational cost and accuracy is not well understood. Laminar flow around a heated cylinder is used as a test case. A robust matrix of 99-simulations is used to analyze the flow of carbon dioxide just above the critical point at the isobaric condition of 8 MPa and temperatures between 305 and 390 K. Four EoS were studied: Ideal Gas Law, cubic EoS’s Peng-Robinson and Soave-Redlich-Kwong, and the Span-Wagner EoS. Results were identical for isothermal cases, but for heated or cooled cases, the high-fidelity Span-Wagner EoS simulations presented drastically different wake structures. Compared to the high-fidelity SW EoS, cEoS’s underestimated Nu values by up to 38% and cEoS’s underestimated CD values by 18–42%. In conclusion, the Span-Wagner EoS required only 17% more computation time then the cubic EoS simulations, making the Span-Wagner EoS recommended for engineering design.

42 ENGINEERING↗

Numerical Investigation of Near-Field Plasma Flows in Magnetic Nozzles

The development and application of a multidimensional numerical simulation code for investigating near-field plasma processes in magnetic nozzles are presented. The code calculates the time-dependent evolution of all three spatial components of both the magnetic field and velocity in a plasma flow, and includes physical models of relevant transport phenomena. It has been applied to an investigation of the behavior of plasma flows found in high-power thrusters, employing a realistic magnetic nozzle configuration. Simulation of a channel-flow case where the flow was super-Alfvenic has demonstrated that such a flow produces adequate back-emf to significantly alter the shape of the total magnetic field, preventing the flow from curving back to the magnetic field coil in the near-field region. Results from this simulation can be insightful in predicting far-field behavior and can be used as a set of self-consistent boundary conditions for far-field simulations. Future investigations will focus on cases where the inlet flow is sub-Alfvenic and where the flow is allowed to freely expand in the radial direction once it is downstream of the coil.

Sankaran, Kamesh↗

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗

Factors Controlling Precipitation Onset and Maintenance: Inferences From Large Eddy Simulations of Two Afternoon Deep Convective Regimes Over Amazon

To investigate the environmental factors controlling the onset and maintenance of afternoon precipitation over tropical rainforests, two contrasting cases are created for large eddy simulations (LES) using GoAmazon observations: one features a single-pulse rain dissipating quickly in earlier afternoon, while the other shows double pulses lasting until evening. For these two specific cases under consideration, LES confirms that early-morning relative humidity dominates afternoon rain patterns. The attributed impacts are distinct: moisture determines the rain onset timing while temperature affects the peak rain intensity. Single-pulse day observes one round of strong precipitation and cold pools, which further suppresses convection. On double-pulse day, the first precipitation peak results from an intermediate development of congestus clouds, whose detrainment leads to a gradual moistening of the lower-to-middle troposphere. This behavior of convection favors a second pulse of stronger precipitation with more convective organization, whose development decouples from surface fluxes and sustains until evening.

Amazon↗

Check-Cases for Verification of 6-Degree-of-Freedom Flight Vehicle Simulations: Appendices - Volume 2

This NASA Engineering and Safety Center (NESC) assessment was established to develop a set of time histories for the flight behavior of increasingly complex example aerospacecraft that could be used to partially validate various simulation frameworks. The assessment was conducted by representatives from several NASA Centers and an open-source simulation project. This document contains details on models, implementation, and results.

Murri, Daniel G.↗

DSO+T: Integrated System Simulation (DSO+T Study: Volume 2)

This report summarizes an integrated co-simulation model used by the Distribution System Operator with Transactive (DSO+T) study to represent an electrical generation, delivery, and end-load systems for the purposes of assessing the viability and value proposition of transactive energy coordination of flexible assets versus a business-as-usual case. The integrated co-simulation model includes the bulk generation and transmission system, including the day-ahead and real-time scheduling and dispatch of thermal generators. Forty distribution system operators were modelled in detail, including tens of thousands of residential and commercial buildings and their flexible end-loads. These included HVAC systems, residential water heaters, electric vehicles, and stationary, behind-the-meter, batteries. Both wholesale market and end-load results for the business-as-usual case are presented and compared to actual ERCOT system data to assess the accuracy and representativeness of the resulting model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Detached-Eddy Simulation Based on the V2-F Model

Detached-eddy simulation (DES) based on the v(sup 2)-f Reynolds-averaged Navier-Stokes (RANS) model is developed and tested. The v(sup 2)-f model incorporates the anisotropy of near-wall turbulence which is absent in other RANS models commonly used in the DES community. The v(sup 2)-f RANS model is modified in order the proposed v(sup 2)-f-based DES formulation reduces to a transport equation for the subgrid-scale kinetic energy isotropic turbulence. First, three coefficients in the elliptic relaxation equation are modified, which is tested in channel flows with friction Reynolds number up to 2000. Then, the proposed v(sup 2)-f DES model formulation is derived. The constant, C(sub DES), required in the DES formulation was calibrated by simulating both decaying and statistically-steady isotropic turbulence. After C(sub DES) was calibrated, the v(sub 2)-f DES formulation is tested for flow around a circular cylinder at a Reynolds number of 3900, in which case turbulence develops after separation. Simulations indicate that this model represents the turbulent wake nearly as accurately as the dynamic Smagorinsky model. Spalart-Allmaras-based DES is also included in the cylinder flow simulation for comparison.

Jee, Sol Keun↗

Composite Analysis Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Sources.

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP1) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗