Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modelling and simulation”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Advanced Supersonic Parachute Inflation Research and Experiment-2 (ASPIRE2) Flight Mechanics Modeling and Simulation

Introduction: The Advanced Supersonic Para-chute Inflation Research Experiment-2 (ASPIRE2) program is a sounding rocket flight test to be conducted at Wallops Island, VA in early 2025. This program is a risk mitigation exercise for the Mars Sample Retrieval Lander (MSRL) program and builds upon the success of the 2017 & 2018 ASPIRE program, [1,2]. While the ASPIRE program certified the strengthened 21.5 m diameter disk-gap-band (DGB) deployed at Mach 1.7, ASPIRE2 will certify a 24 m diameter DGB deployed at Mach 2.1. The need to certify this increased parachute performance is driven the increased lander mass for MSRL; over 50% increase in comparison to M2020. Modeling: Like its predecessor, ASPIRE2 is developing a multi-body flight dynamics simulation to predict parachute dynamics and aide in designing the flight test that will target Mars-relevant flight conditions, as shown in Figure 1. This work de-scribes the parachute modeling, flight mechanics simulations (from payload separation to splash-down), and design trades used to prepare for the 2025 ASPIRE2 flight. Discussed herein are comparisons between ASPIRE and ASPIRE2, noting key differences in the parachute modeling and vehicle configuration. A study on the attitude control system performance impacts with respect to payload section design will also be presented. This work will quantify the pre-flight parachute performance in the presence of uncertainties, such as those associated with the separation from the sounding rocket, atmosphere, the parachute system, and vehicle mass. The pre-flight predictions will include Monte Carlo analyses, powered by the flight mechanics simulations to show the ASPIRE2 vehicle performance in meeting program requirements on parachute deployment conditions (Mach & dynamic pressure), parachute loads, vehicle attitude at key milestones, and the vehicle splashdown conditions.

Entry Descent and Landing↗
DOE OSTI2018

Neutronics Phenomena Important in Modeling and Simulation of Liquid-Fuel Molten Salt Reactors

This paper discusses liquid-fuel molten salt reactors, how they will operate under normal, transient, and accident conditions, and the results of an expert elicitation to determine the corresponding neutronic phenomena important to understanding their behavior. Identifying these phenomena will enable the U.S. Nuclear Regulatory Commission (NRC) to develop or identify modeling functionalities and tools required to carry out confirmatory analyses that examine the validity and accuracy of applicants’ calculations and help determine the margin of safety in plant design. NRC frequently does an expert elicitation using a Phenomena Identification and Ranking Table (PIRT) to identify and evaluate the state of knowledge of important modeling phenomena. However, few details about the design of these reactors and the sequence of events during accidents are known, so the process used was considered a preliminary PIRT. A panel met to define phenomena that would need to be modeled and considered the impact/importance of each phenomenon with respect to specific figures-of-merit (FoMs) (e.g., power distribution, fluence, kinetics parameters and reactivity). Each FoM reflected a potential impact on radionuclide release or loss of a barrier to release. The panel considered what the path forward might be with respect to being able to model the phenomenon in a simulation code. Results are explained for both thermal and fast spectrum designs.

42 ENGINEERING↗
DOE OSTI2016

Simulink models for performance analysis of high speed DQPSK modulated optical link

This paper attempts to present the design approach for development of simulation models to study and analyze the transmission of 10 Gbps DQPSK signal over a single channel Peer to Peer link using Matlab Simulink. The simulation model considers the different optical components used in link design with their behavior represented initially by theoretical interpretation, including the transmitter topology, Mach Zehnder Modulator(MZM) module and, the propagation model for optical fibers etc. thus allowing scope for direct realization in experimental configurations. It provides the flexibility to incorporate the various photonic components as either user-defined or fixed and, can also be enhanced or removed from the model as per the design requirements. We describe the detailed operation and need of every component model and its representation in Simulink blocksets. Moreover the developed model can be extended in future to support Dense Wavelength Division Multiplexing (DWDM) system, thereby allowing high speed transmission with N × 40 Gbps systems. The various compensation techniques and their influence on system performance can be easily investigated by using such models.

42 ENGINEERING↗
NASA NTRS2001

Integrating LMINET with TAAM and SIMMOD: A Feasibility Study

LMINET is a queuing network air traffic simulation model implemented at 64 large airports and the entire National Airspace System in the United States. TAAM and SIMMOD are two widely used air traffic event-driven simulation models mostly for airports. Based on our proposed Progressive Augmented window approach, TAAM and SIMMOD are integrated with LMINET though flight schedules. In the integration, the flight schedules are modified through the flight delays reported by the other models. The benefit to the local simulation study is to let TAAM or SIMMOD take the modified schedule from LMINET, which takes into account of the air traffic congestion and flight delays at the national network level. We demonstrate the value of the integrated models by the case studies at Chicago O'Hare International Airport and Washington Dulles International Airport. Details of the integration are reported and future work for a full-blown integration is identified.

Long, Dou↗
DOE OSTI2018

Thermal-Hydraulics Phenomena Important in Modeling and Simulation of Liquid-Fuel Molten Salt Reactors

This paper discusses liquid-fuel molten salt reactors, how they will operate under normal, transient, and accident conditions, and the results of an expert elicitation to determine the corresponding thermalhydraulic phenomena important to understanding their behavior. Identifying these phenomena will enable the U.S. Nuclear Regulatory Commission (NRC) to develop or identify modeling functionalities and tools required to carry out confirmatory analyses that examine the validity and accuracy of an applicant’s calculations and help determine the margin of safety in plant design. NRC frequently does an expert elicitation using a Phenomena Identification and Ranking Table (PIRT) to identify and evaluate the state of knowledge of important modeling phenomena. However, few details about the design of these reactors and the sequence of events during accidents are known, so the process used was considered a preliminary PIRT. A panel met to define phenomena that would need to be modeled and considered the impact/importance of each phenomenon with respect to specific figures-of-merit (FoMs) (e.g., salt temperature, velocity, and composition). Each FoM reflected a potential impact on radionuclide release or loss of a barrier to release. The panel considered what the path forward might be with respect to being able to model the phenomenon in a simulation code. Results are explained for both thermal and fast spectrum designs.

42 ENGINEERING↗
DOE OSTI2016

Helium reionization simulations. I. Modeling quasars as radiation sources

We introduce a new project to understand helium reionization using fully coupled N-body, hydrodynamics, and radiative transfer simulations. This project aims to capture correctly the thermal history of the intergalactic medium as a result of reionization and make predictions about the Lyα forest and baryon temperature–density relation. The dominant sources of radiation for this transition are quasars, so modeling the source population accurately is very important for making reliable predictions. In this first paper, we present a new method for populating dark matter halos with quasars. Our set of quasar models includes two different light curves, a lightbulb (simple on/off) and symmetric exponential model, and luminosity-dependent quasar lifetimes. Our method self-consistently reproduces an input quasar luminosity function given a halo catalog from an N-body simulation, and propagates quasars through the merger history of halo hosts. After calibrating quasar clustering using measurements from the Baryon Oscillation Spectroscopic Survey, we find that the characteristic mass of quasar hosts is M{sub h}∼2.5×10{sup 12} h{sup −1} M{sub ⊙} for the lightbulb model, and M{sub h}∼2.3×10{sup 12} h{sup −1} M{sub ⊙} for the exponential model. In the latter model, the peak quasar luminosity for a given halo mass is larger than that in the former, typically by a factor of 1.5–2. The effective lifetime for quasars in the lightbulb model is 59 Myr, and in the exponential case, the effective time constant is about 15 Myr. We include semi-analytic calculations of helium reionization, and discuss how to include these quasars as sources of ionizing radiation for full hydrodynamics with radiative transfer simulations in order to study helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Dynamical simulations of the interacting galaxies in the NGC 520/UGC 957 system

Numerical simulations of the interacting galaxies in the NGC 520/UGC 957 system are presented. Two sets of models were produced to investigate the postulated three-galaxy system of two colliding disk galaxies within NGC 520 and the dwarf galaxy UGC 957. The first set of models simulated a dwarf perturbing one-disk galaxy, which tested the possibility that NGC 520 contains only one galaxy disturbed by the passage of UGC 957. The resulting morphology of the perturbed single disk in the simulation fails to reproduce the observed tidal tails and northwest mass condensation of NGC 520. A second set of models simulated two colliding disks, which tested the hypothesis that NGC 520 itself contains two galaxies in a strong collision and UGC 957 is unimportant to the interaction. These disk-disk models produced a good match to the morphology of the present NGC 520. It is concluded that (1) NGC 520 contains two colliding disk galaxies which have produced the brighter southern half of the long tidal tail and (2) UGC 957, which may originally have been a satellite of one of the disk galaxies, formed the diffuse northern tail as it orbited NGC 520.

Stanford, S. A.↗
DOE OSTI2022

Analytical modeling and simulation of electrical contact resistance for elastic rough electrode surface contact including frictional temperature rise

An improved electrical contact resistance (ECR) model for elastic rough electrode contact is proposed, incorporating the effects of asperity interactions and temperature rise by frictional and joule heating. The analytical simulation results show that the ECR decreases steeply at the beginning of the contact between Al and Cu. However, it becomes stabilized after reaching a specific contact force. It is also found that the longer elapsed sliding contact time, the higher ECR due to the increase in electrical resistivity of electrode materials by the frictional temperature rise at the interface. The effects of surface roughness parameters on ECR are studied through the 3 2 full-factorial design-of-experiment analysis. Based on the two representative roughness parameters, i.e., root-mean-square (rms) roughness and asperity radius, their individual and coupled effects on the saturated ECR are examined. The saturated ECR increases with the rms roughness for a rough machined surface condition, but it is hardly affected by the asperity radius. On the other hand, the saturated ECR increases with both the rms roughness and the asperity radius under a smooth thin film surface condition.

42 ENGINEERING↗
DOE OSTI2018

Experimental characterization and modeling of optical tweezer particle handling dynamics [Modeling and Simulation of Optical Tweezer Particle Handling Dynamics]

Here, we report a new framework for a quantitative understanding of optical trapping (OT) particle handling dynamics. We present a novel three-dimensional particle-based model that includes optical, hydrodynamic, and inter-particle forces. This semi-empirical colloid model is based on an open-source simulation code known as LAMMPS (large-scale atomic/molecular massively parallel simulator) and properly recapitulates the full OT force profile beyond the typical linear approximations valid near the trap center. Simulations are carried out with typical system parameters relevant for our experimental holographic optical trapping (HOT) system, including varied particle sizes, trap movement speeds, and beam powers. Furthermore, we present a new experimental method for measuring both the stable and metastable boundaries of the optical force profile to inform or validate the model’s underlying force profile. We show that our framework is a powerful tool for accurately predicting particle behavior in a practical experimental OT setup and can be used to characterize and predict particle handling dynamics within any arbitrary OT force profile.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Probabilistic modeling for simulation of aerodynamic uncertainties in propulsion systems

The numerical simulation of the probabilistic aerothermodynamic response of propulsion system components to randomness in their environment was explored. The reusable rocket engine turbopumps were selected as an example because of the severe cryogenic environment in which they operate. The thermal and combustion instabilities, coupled with the engine thrust requirements from start up to shut down, lead to randomness in the flow variables and uncertainties in the aerodynamic loading. The probabilistic modeling of the turbopumps aerodynamic response was accomplished using the panel method coupled with fast probability integration methods. The aerodynamic response in the form of probabilistic rotor blades and splitter loading were predicted and the results presented for specified flow coefficient and rotor preswirl variance. Possible future applications of the aerothermodynamic probabilistic modeling in engine transient simulation, condition monitoring and engine life prediction are briefly discussed.

Hamed, Awatef↗
NASA NTRS1990

Probabilistic modeling for simulation of aerodynamic uncertainties in propulsion systems

The numerical simulation of the probabilistic aerothermodynamic response of propulsion system components to randomness in their environment was explored. The reusable rocket engine turbopumps were selected as an example because of the severe cryogenic environment in which they operate. The thermal and combustion instabilities, coupled with the engine thrust requirements from start up to shut down, lead to randomness in the flow variables and uncertainties in the aerosynamic loading. The probabilistic modeling of the turbopumps' aerodynamic response was accomplished using the panel method coupled with fast probability integration methods. The aerodynamic response in the form of probabilistic rotor blades and splitter loading were predicted and the results presented for specified flow coefficient and rotor preswirl variance. Possible future applications of the aerothermodynamic probabilistic modeling in engine transient simulation, condition monitoring, and engine life prediction are briefly discussed.

Hamed, A.↗
DOE OSTI2013

Westinghouse waste simulation and optimization software tool

Applications for dynamic simulation can be found in virtually all areas of process engineering. The tangible benefits of using dynamic simulation can be seen in tighter design, smoother start-ups and optimized operation. Thus, proper implementation of dynamic simulation can deliver substantial benefits. These benefits are typically derived from improved process understanding. Simulation gives confidence in evidence based decisions and enables users to try out lots of 'what if' scenarios until one is sure that a decision is the right one. In radioactive waste treatment tasks different kinds of waste with different volumes and properties have to be treated, e.g. from NPP operation or D and D activities. Finding a commercially and technically optimized waste treatment concept is a time consuming and difficult task. The Westinghouse Waste Simulation and Optimization Software Tool will enable the user to quickly generate reliable simulation models of various process applications based on equipment modules. These modules can be built with ease and be integrated into the simulation model. This capability ensures that this tool is applicable to typical waste treatment tasks. The identified waste streams and the selected treatment methods are the basis of the simulation and optimization software. After implementing suitable equipment data into the model, process requirements and waste treatment data are fed into the simulation to finally generate primary simulation results. A sensitivity analysis of automated optimization features of the software generates the lowest possible lifecycle cost for the simulated waste stream. In combination with proven waste management equipments and integrated waste management solutions, this tool provides reliable qualitative results that lead to an effective planning and minimizes the total project planning risk of any waste management activity. It is thus the ideal tool for designing a waste treatment facility in an optimum manner, taking account of the detailed waste stream and specific requirements. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Validation of cloud forcing simulated by the National Center for Atmospheric Research Community Climate Model using observations from the Earth Radiation Budget Experiment

Satellite measurements of the effect of clouds on the top of atmosphere radiative energy budget are used to validate model simulations from the National Center for Atmospheric Research Community Climate Model (NCAR CCM). The ability of the NCAR CCM to reproduce the monthly mean global distribution and temporal variability on both daily and seasonal time scales is assessed. The comparison reveals several deficiencies in the CCM cloud representation. Most notable are the difficulties in properly simulating the effect of clouds on the planetary albedo. This problem arises from discrepancies in the model's portrayal of low-level cloudiness and leads to significant errors in the absorbed solar radiation simulated by the model. The CCM performs much better in simulating the effect of clouds on the longwave radiation emitted to space, indicating its relative success in capturing the vertical distribution of cloudiness. The daily variability of the radiative effects of clouds in both the shortwave and longwave spectral regions is systematically overestimated. Analysis of the seasonal variations illustrates a distinct lack of coupling in the seasonal changes in the radiative effects of cloudiness between the tropics and mid-latitudes and between the Northern and Southern Hemisphere. Much of this problem also arises from difficulties in simulating low-level cloudiness, placing further emphasis on the need for better model parameterizations of boundary layer clouds.

Soden, B. J.↗
DOE OSTI2019

Recrystallization and Grain Growth Simulations for Multiple-Pass Rolling and Annealing of U-10Mo

A simulation model of recrystallization and grain growth has been developed to investigate grain structure evolution during deformation and heat treatment in polycrystalline U-10wt% Mo (U-10Mo) fuel. Experimentally obtained U-10Mo microstructures were used as input for closed loop simulations of multiple rolling, intermediate heating and final annealing process. FEM calculations of deformation and Potts Model simulations of recrystallization and grain growth were used to iteratively inform each subsequent stage of simulation. Furthermore, the model was then applied to predict the grain structure evolution during multiple-pass hot rolling of U-10Mo, and benchmarked against experimentally observed U-10Mo recrystallization behavior. The results showed that our model was able to capture the coupling between deformation and recrystallization as a function of microstructure, including particle stimulated nucleation and recrystallization nucleation on grain boundaries. Additionally, we have achieved reasonable quantitative agreement with U-10Mo recrystallization and grain growth behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Development of a Stochastic Infiltration Model for PPA Implementation - 19293

Performance Assessment models often require inputs representing specific hydrologic processes and properties such as infiltration, water content, groundwater flow rate, streamflow, and sediment flux. This paper details the development of a process-level infiltration model and its use in a Probabilistic Performance Assessment (PPA) model for the West Valley Site in western New York. Additionally, this paper demonstrates a method for obtaining distributions of the long-term average hydrologic processes and properties for a PA model. For this PPA, an infiltration model was developed using HELP (Hydrologic Evaluation of Landfill Performance) to evaluate the rates of infiltration, percolation, and drainage through several site soil and proposed decommissioning cover systems, as well as the volumetric water contents of each of the unsaturated layers in a modeled soil column. A particular challenge of probabilistic modeling over long periods of time is scaling parameters appropriately so that outcomes are not grossly over- or underestimated. For these models, distributions of climatic parameters such as average precipitation rate, temperature, wind speed, humidity, and solar radiation were scaled to 1000-year averages and were used to force steady-state models that predict 1000-year average infiltration rates, drainage rates, and water contents. For the HELP model, distributions were developed for the physical and hydrologic properties of each material in a given soil column or cover system, and these properties were allowed to vary stochastically along with climatic inputs. The HELP model was run 5000 times for each of seven different soil and cover systems representing various existing and proposed conditions (i.e. native soil, backfill, existing covers, and a proposed close-in-place cover). The results of the HELP model provide ranges of 1000-year average infiltration rates, drainage rates, and water contents through and within the various layers of each of the different soil and cover conditions. The model outcomes also demonstrate the effectiveness of a proposed close-in-place cover, which essentially eliminates infiltration during its operational lifetime. Following model simulation, a sensitivity analysis was performed to establish the relative importance of each stochastic input with respect to each output. The outcomes of the model simulations and the sensitivity analysis were implemented into a system-level GoldSim model via a novel sensitivity-informed nearest neighbor approach. These infiltration rates, drainage rates, and water contents are used to parameterize mixing cells within the GoldSim Contaminant Transport module, which is the primary system-level model used in the PPA. The infiltration rates are also used as boundary conditions for a stochastic groundwater flow model that is also utilized in the PPA, but not discussed in this paper. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1998

Computational Nonlinear of Vertical-Cavity Surface-Emitting Lasers

Vertical-Cavity Surface-Emitting Lasers (VCSELs) have become increasingly important for many applications, such as optical interconnects and optical data storages. In many applications, critical requirements are very often high beam quality and/or high output powers. To help the design and optimization of VCSELs for such applications, it is essential to have a simulation model that includes the time evolution of the transverse space dependence of the beam. Such models allow arbitrary transverse profiles of laser intensity or carrier density to develop in time without any a priori assumption about the number and types of transverse modes. In addition, extensive research in the past has shown that many-body effects are important in influencing the optical gain and the refractive index of semiconductor quantum wells. It is thus important to include the key many-body effects in simulation models for semiconductor lasers. A computational model, called the Maxwell-Effective Bloch Equations (MEBE) was developed recently that simulates the above features and effects. This model combines the space-time resolution with microscopic many-body effects in a way that is both accurate enough and computationally manageable. This model has been used to simulate edge-emitting high power lasers recently. A new finite-difference algorithm has been developed for solving these equations and will be shown in the presentation. That algorithm will be used to calculate the detailed time evolution of the transverse mode profiles of VCSELs. The VCSELs are based on InGaAs/GaAs quantum well structures. Effects of pumping current and refractive index profiles will be investigated in detail. The simulation results will be compared with other results the use the rate-equation model.

Goorjian, Peter M.↗
NASA NTRS2019

TPSAS-NF1676L-34010-DND

We invert global black carbon (BC), organic carbon (OC) and desert dust (DD) aerosol emissions from POLDER/PARASOL spectral aerosol optical depth (AOD) and aerosol absorption optical depth (AAOD) using the GEOSChem inverse modelling framework. Our inverse modeling framework uses standard a priori emissions to provide a posteriori emissions that are constrained by POLDER/PARASOL AODs and AAODs. Our satellite-derived global aerosol emission database indicates an increase of BC and OC, while a decrease of DD with respect to the a priori emission inventories used in standard GEOS-Chem runs. The model simulations using a posteriori emissions (i.e. retrieved emissions) provide new estimate of global aerosol extinction and absorption constrained by observations. Additionally, the a posteriori model simulation of AOD, AAOD, single scattering albedo, Angstrom exponent, and absorption Angstrom exponent show better agreement with independent AERONET, MODIS, and OMI measurements than the a priori simulation. Thus, this study suggests that using satellite-constrained global aerosol emissions in aerosol transport models can improve the accuracy of simulated global aerosol properties and have potential to better understanding global aerosol effects on climate and human health.

Cheng Chen↗