Search NASA⌕ Search

SEARCH · Search NASA

Results for “computer simulation model development”

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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Towards a theory of economic value for modeling and simulation: incremental cost of parallel simulation (wip)

A useful perspective on the economic value of simulation in engineering comes from considering the role simulation plays in constructing a system. We propose a conceptual model that combines an embryonic concept of the value of an information bit with cost estimation tools used in software engineering. The model provides an economic justification for the common conception that models should be built to address specific questions rather than to comprehensively represent a referent system. Extending these ideas, we offer a model for the incremental value of parallel execution to explain why simulation applications that take advantage of sophisticated parallel discrete event simulation tools form a relatively narrow niche despite the ubiquity of parallel computing hardware. The proposed model enables us to bring value of information criteria into model development and to make predictions about the future of parallel simulation applications.

Nutaro, James J↗
NASA NTRS1981

The application of simulation modeling to the cost and performance ranking of solar thermal power plants

Small solar thermal power systems (up to 10 MWe in size) were tested. The solar thermal power plant ranking study was performed to aid in experiment activity and support decisions for the selection of the most appropriate technological approach. The cost and performance were determined for insolation conditions by utilizing the Solar Energy Simulation computer code (SESII). This model optimizes the size of the collector field and energy storage subsystem for given engine generator and energy transport characteristics. The development of the simulation tool, its operation, and the results achieved from the analysis are discussed.

Rosenberg, L. S.↗
NASA NTRS2019

Scenario Complexity for Unmanned Aircraft System Traffic

This work introduces an approach to estimate the complexity of a low-altitude air traffic scenario involving multiple UASs using mathematical programming. Given a set of multi-point UAS flight trajectories, vehicle dynamics, and a conflict resolution algorithm, an abstract model is developed such that it can be solved quickly using a mathematical programming optimization software without running high-fidelity simulations that can be computationally expensive and may not suit real-time applications. In the abstract model, each vehicle is represented by a time-varied vector associated with position, speed, and heading information. The total extra distance that aircraft need to divert from their original routes to avoid collisions is computed and used to setup a quadratic programming formula. The metrics including the number of conflicts and extra distances travelled by all vehicles are then utilized to estimate the complexity of a given UAS flight scenario. Results and verification against high-fidelity simulations will be provided in the final draft.

UTM↗
DOE OSTI2017

Limits on the Efficiency of Event-Based Algorithms for Monte Carlo Neutron Transport

The traditional form of parallelism in Monte Carlo particle transport simulations, wherein each individual particle history is considered a unit of work, does not lend itself well to data-level parallelism. Event-based algorithms, which were originally used for simulations on vector processors, may offer a path toward better utilizing data-level parallelism in modern computer architectures. In this study, a simple model is developed for estimating the efficiency of the event-based particle transport algorithm under two sets of assumptions. Data collected from simulations of four reactor problems using OpenMC was then used in conjunction with the models to calculate the speedup due to vectorization as a function of two parameters: the size of the particle bank and the vector width. When each event type is assumed to have constant execution time, the achievable speedup is directly related to the particle bank size. We observed that the bank size generally needs to be at least 20 times greater than vector size in order to achieve vector efficiency greater than 90%. When the execution times for events are allowed to vary, however, the vector speedup is also limited by differences in execution time for events being carried out in a single event-iteration. For some problems, this implies that vector effciencies over 50% may not be attainable. While there are many factors impacting performance of an event-based algorithm that are not captured by our model, it nevertheless provides insights into factors that may be limiting in a real implementation.

Romano, Paul K.↗
DOE OSTI2021

Quantifying Turbine-Level Risk to Golden Eagles Using a High-Fidelity Updraft Model and a Stochastic Behavioral Model

To minimize the effects of wind farms on Golden Eagle (Aquila chrysaetos) populations while enabling sustainable development of renewable energy resources, it is important to understand how eagles interact with atmospheric flows, terrain features, and anthropogenic structures. Models that predict migratory flight paths provide one tool that helps us grasp how the location of wind farms may influence interactions and impacts on migrating Golden Eagles. The current state-of-the-art in predicting migratory flight paths uses a deterministic fluid-flow analogy to predict eagle trajectory using only an orographic updraft potential computed from topographical features. This model does not take into account variables, such as thermal updrafts and time varying atmospheric conditions that are known to influence migratory behavior. In this work, we improve on the model with the objective of developing tools that advance our understanding of how atmospheric flows and terrain features affect migratory eagle behavior and their interactions with wind farms. Specifically, we 1) incorporate both orographic and thermal updraft information in simulating eagle flight paths; 2) incorporate stochasticity into eagle travel patterns to better capture the influence of exogenous factors on, and the inherent stochasticity of eagle behavior; 3) consider spatio-temporal atmospheric data at wind-farm-scale when computing updraft potential; and 4) account for how atmospheric conditions and the direction of migration change seasonally and how these changes affect eagle migratory flight behavior. We tested the model using a 50km by 50km region with 50 m resolution in the western United States. We simulated 900 independent, probabilistic eagle tracks during southerly and northerly migration, assuming eagles solely rely on orographic updrafts. The preliminary results indicate that the inclusion of finer resolution atmospheric data allows for the inclusion of realistic conditions that an eagle experiences. The stochasticity in eagle tracks provides a platform to include uncertainty in eagle decision making and help produce robust eagle presence maps. We will deploy updraft and downdraft velocities computed using a high-fidelity, wind farm scale, computational fluid dynamics solver under development at National Renewable Energy Laboratory. This work is a first step in the development of a predictive and generalizable eagle behavior model at the wind farm scale that does not rely on empirical data collection. Although the current model is intended for migratory eagles, we will extend and refine this model to inform the development of additional behavioral modes, including resident eagle behavior. This modeling approach improves our ability to understand eagle use of the landscape at a fine scale, and it is our hope that this work will ultimately help advance strategies that minimize the impact of wind development on Golden Eagle populations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2025

Insight into the deformation features and capacity loss mechanisms of lithium-ion pouch cells under spherical indentation conditions

Mechanical deformation under extreme conditions is one of the important reasons for the failure of lithium-ion batteries in automotive application. However, the deformation features and component failure of lithium-ion cells to external loading has never been a design consideration. Here, in this study, we conduct spherical indentation tests on a dozen of lithium-ion cells with different capacities under different control mode conditions to investigate their deformation features and capacity loss mechanisms. The experimental results show that, under mechanical deformation conditions, internal faults of cells occur in stages, and energy accumulation and sudden release are two key processes of cell's mechanical failure. The cells' state of charge is the main factor affecting their thermal runaway behaviors. In addition, a finite element model is developed to simulate the deformation features and the failure mechanism of key components of lithium-ion pouch cells; the 3D x-ray computed tomography is employed to demonstrate its internal configuration. With this model, the force-strain response, the deformation features as well as the size of the failure area of lithium-ion cells under spherical indentation conditions are accurately predicted. In 3D x-ray computed tomography images, unique mud cracks in cooper current collector are observed, and the influence mechanisms of the isolated fragments on the cell capacities are revealed. These results may provide useful information for the mechanical structure design of the components of lithium-ion pouch cells.

25 ENERGY STORAGE↗
DOE OSTI2019

Evaluation of a thermomechanical model for prediction of residual stress during laser powder bed fusion of Ti-6Al-4V

The build-up of residual stresses in a part during laser powder bed fusion provides a significant limitation to the adoption of this process. These residuals stresses may cause a part to fail during a build or fall outside the specified tolerances after fabrication. In the present work a thermomechanical model is used to simulate the build process and calculate the residual stress state for Ti–6Al–4V specimens built with continuous and island scan strategies. A layer agglomeration, or lumping, approach is used to speed up the computations. A material model is developed to naturally capture the strain-rate dependence and annealing behavior of Ti–6Al–4V at elevated temperatures. Conclusions from the thermomechanical simulations showed good agreement with synchrotron X-ray diffraction measurements used to determine the residual elastic strains in these parts. Yet, the experimental measurements showed higher residual strains for the specimen built with an island scan strategy; a trend not fully captured by the simulations. Parameter studies were performed to fully understand the advantages and limitations of the current simulation methodology. Reasons for both the computational and experimental findings are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An efficient hybrid downscaling framework to estimate high-resolution river hydrodynamics

Flow depth and velocity are the most important hydrodynamic variables that govern various river functions, including water resources, navigation, sediment transport, and biogeochemical cycling. Existing high-resolution flow depth simulations rely on either computationally expensive river hydrodynamic models (RHMs) or data-driven models with formidable training costs, whereas data-driven modeling of flow velocity has rarely been explored. Here, using the hybrid Low-fidelity, Spatial analysis, and Gaussian process learning (LSG) model, we developed a downscaling approach to construct high-resolution flow depth and velocity from a two-dimensional (2-D) RHM simulation at coarse resolution. The LSG models were trained and tested in an urban watershed in Houston using two different hurricane-driven flood events. The high-resolution (as fine as 30 m resolution) and low-resolution (mostly 1000 m resolution) meshes include 664 724 and 14 536 grid cells, respectively. The results showed that through downscaling, the simulation errors were reduced to less than one-fourth and one-third of the errors of the low-resolution 2-D RHM for flow depth and velocity, respectively. Our analysis further revealed that the dominant uncertainty sources of the downscaled hydrodynamics are different, with flow velocity dominated by the dimensionality reduction error, which we reduced by using a regionalized training procedure. The downscaling approach achieves an 84-fold acceleration in computational time compared to the high-resolution 2-D RHM, making high-fidelity ensemble flood modeling feasible. More importantly, the developed method provides an opportunity to couple large-scale hydrodynamical processes with local physical, chemical, and biological processes in river models.

Tan, Zeli [Pacific Northwest National Laboratory (↗
NASA NTRS1987

Integration of communications and tracking data processing simulation for space station

A simplified model of the communications network for the Communications and Tracking Data Processing System (CTDP) was developed. It was simulated by use of programs running on several on-site computers. These programs communicate with one another by means of both local area networks and direct serial connections. The domain of the model and its simulation is from Orbital Replaceable Unit (ORU) interface to Data Management Systems (DMS). The simulation was designed to allow status queries from remote entities across the DMS networks to be propagated through the model to several simulated ORU's. The ORU response is then propagated back to the remote entity which originated the request. Response times at the various levels were investigated in a multi-tasking, multi-user operating system environment. Results indicate that the effective bandwidth of the system may be too low to support expected data volume requirements under conventional operating systems. Instead, some form of embedded process control program may be required on the node computers.

Lacovara, Robert C.↗
NASA NTRS2021

Coupling of the FUN3D Unstructured Flow Solver and the LASTRAC Stability Code to Model Transition

We develop an iterative automated method to predict transition locations in boundary-layer flows by using the FUN3D solver to perform flow simulations and the LASTRAC code for linear stability computations. The coupling of FUN3D and LASTRAC allows for a robust physics-based approach to model boundary-layer transition by analyzing the growth of different instability waves and then using that information to iteratively update the resulting transition location. There is no user involvement during the iterative computations. We apply this automated method to subsonic flow over a flat plate with a sharp leading edge. The final solution has regions of laminar and turbulent flow with a transition onset location that agrees with experiments and stability-based correlations. This iterative automated method is also applied to an NLF(1)-0416 airfoil for conditions with and without a separation bubble. Along with predicting transition locations, we compare the streamwise distributions of surface-pressure and skin-friction coefficients to a transport-equation-based model. We consider a 6:1 prolate spheroid at three angles of attack, namely, five, ten, and fifteen degrees, where mixed-mode transition occurs due to both Tollmien-Schlichting and crossflow instabilities. The skin-friction contours and transition fronts at every angle of attack from our iterative automated method show good agreement with past experimental and computational results in the literature for the 6:1 prolate spheroid.

Transition↗
DOE OSTI2017

A numerical analysis of a magnetocaloric refrigerator with a 16-layer regenerator

A numerical analysis was conducted to study a room temperature magnetocaloric refrigerator with a 16-layer parallel plates active magnetic regenerator (AMR). Sixteen layers of LaFeMnSiH having different Curie temperatures were employed as magnetocaloric material (MCM) in the regenerator. Measured properties data was used. A transient one dimensional (1D) model was employed, in which a unique numerical method was developed to significantly accelerate the simulation speed of the multi-layer AMR system. As a result, the computation speed of a multi-layer AMR case was very close to the single-layer configuration. The performance of the 16-layer AMR system in different frequencies and utilizations has been investigated using this model. To optimize the layer length distribution of the 16-layer MCMs in the regenerator, a set of 137 simulations with different MCM distributions based on the Design of Experiments (DoE) method was conducted and the results were analyzed. The results show that the 16-layer AMR system can operate up to 84% of Carnot cycle COP at a temperature span of 41 K, which cannot be obtained using an AMR with fewer layers. Here, the DoE results indicate that for a 16-layer AMR system, the uniform distribution is very close to the optimized design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.

Computational Fluid Dynamics↗
DOE OSTI2019

A simplified parameterization of isoprene-epoxydiol-derived secondary organic aerosol (IEPOX-SOA) for global chemistry and climate models: a case study with GEOS-Chem v11-02-rc

Secondary organic aerosol derived from isoprene epoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction of total isoprene SOA, but the current volatility-based lumped SOA parameterizations are not appropriate to represent the reactive uptake of IEPOX onto acidified aerosols. A full explicit modeling of this chemistry is however computationally expensive owing to the many species and reactions tracked, which makes it difficult to include it in chemistry–climate models for long-term studies. Here we present three simplified parameterizations (version 1.0) for IEPOX-SOA simulation, based on an approximate analytical/fitting solution of the IEPOX-SOA yield and formation timescale. The yield and timescale can then be directly calculated using the global model fields of oxidants, NO, aerosol pH and other key properties, and dry deposition rates. The advantage of the proposed parameterizations is that they do not require the simulation of the intermediates while retaining the key physicochemical dependencies. We have implemented the new parameterizations into the GEOS-Chem v11-02-rc chemical transport model, which has two empirical treatments for isoprene SOA (the volatility-basis-set, VBS, approach and a fixed 3% yield parameterization), and compared all of them to the case with detailed fully explicit chemistry. The best parameterization (PAR3) captures the global tropospheric burden of IEPOX-SOA and its spatiotemporal distribution ( R 2 =0.94 ) vs. those simulated by the full chemistry, while being more computationally efficient ( ~5 times faster), and accurately captures the response to changes in NO x and SO 2 emissions. On the other hand, the constant 3 % yield that is now the default in GEOS-Chem deviates strongly ( R 2 =0.66 ), as does the VBS ( R 2 =0.47 , 49% underestimation), with neither parameterization capturing the response to emission changes. With the advent of new mass spectrometry instrumentation, many detailed SOA mechanisms are being developed, which will challenge global and especially climate models with their computational cost. The methods developed in this study can be applied to other SOA pathways, which can allow including accurate SOA simulations in climate and global modeling studies in the future.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

A CFD Framework for Film Boiling and Quenching in Cryogenic Chilldown

Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.

Computational Fluid Dynamics↗
DOE OSTI2021

Advancing Catalytic Fast Pyrolysis Through Integrated Experimentation and Multi-Scale Computational Modeling

This webinar will highlight recent results from a multi-disciplinary research effort in which integrated reaction testing was coupled with particle- and reactor-scale computational modeling to advance catalytic fast pyrolysis (CFP) for the production of renewable hydrocarbon fuels. Data will be presented from a series of ex situ CFP experiments in which a fixed bed of Pt/TiO2 was utilized with co-fed H2 to upgrade woody biomass pyrolysis vapors. Further discussion will include the application of these data towards the development of (1) a multiscale simulation framework to de-couple apparent kinetics from both intraparticle and reactor-scale transport phenomena and (2) a finite element computational model to understand and predict thermal excursions during catalyst regeneration. Throughout the presentation, the speakers will emphasize synergistic outcomes derived from the collaborative approach and highlight ongoing research efforts to accelerate technology maturation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

The Influence of Environment on Post-Detonation Chemistry and Debris Formation (Abbreviated Final Report: 20-SI-006)

Predicting, responding to, or interpreting the chemical record preserved in debris derived from nuclear events can be challenging due to chemical fractionation. Chemical fractionation is where different species of the evolving radionuclide inventory segregate and/or are lost from the system over the timescales of debris formation. Both historic data and recent research suggest that the interaction and character of the local environment may exert controls on chemical fractionation by influencing the cooling and evolution of the associated fireball as well as the composition of the vapor term and resultant speciation. Prior to this work, an integrated platform permitting dynamic and concurrent consideration of physical and chemical evolution of early time post-detonation event environments did not exist. Our work merged historic data and experimental approaches to support development of a computational framework able to simulate fundamental processes (e.g., entrainment of local environment, oxidation chemistry, and cooling time scales) that may perturb the radionuclide inventory captured in post-detonation debris. Work with historic debris confirmed that entrained environmental material affect debris composition, structure, and radionuclide incorporation. Complementary work utilizing a readily controllable and tunable benchtop setup (a plasma flow reactor) simulated the late cooling of a nuclear fireball (e.g., T < 6000 K) and bounded the sensitivity of actinide speciation and particle size distribution to variations in oxygen concentration and cooling rates. Concurrent laser ablation and laser heating experiments were used to investigate the chemistry and physics of processes occurring in vaporized and/or rapidly heated actinides and other elements in the presence of oxygen. A more computationally efficient microphysical model was developed for predicting and evolving size distributions of particles forming from mixed vapor terms and simulating particle formation processes under a variety of extreme conditions. Continued study of historic nuclear event film confirmed that shockwave data and physics codes agree to within the uncertainty of the data. Good agreement was achieved for thermal emission from an airburst, however the paucity of low-temperature molecular opacity data for mixtures of air, bomb debris, entrained dirt, and water vapor complicate agreement for more elaborate scenarios. A multiphysics code (ALE3D) was modified to bring the necessary physics and chemistry, including these new data and insights, onto a single platform. Code development included improved initialization of large physical systems, modernization of chemistry capabilities, and modifications to enable inclusion of particle transport.

07 ISOTOPE AND RADIATION SOURCES↗