Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation software.”

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

Advanced Reactors-Intermediate Heat Exchanger (IHX) Coupling: Theoretical Modeling and Experimental Validation

The overall goal of the research project was to model the behavior of the advanced reactorintermediate heat exchange system and to develop advanced control techniques for off-normal conditions. The specific objectives defined for the project were: 1. To develop the steady-state thermal hydraulic design of the intermediate heat exchanger (IHX); 2. To develop mathematical models to describe the advanced nuclear reactor-IHX-chemical process/power generation coupling during normal and off-normal operations, and to simulate models using multiphysics software; 3. To develop control strategies using genetic algorithm or neural network techniques and couple these techniques with the multiphysics software; 4. To validate the models experimentally The project objectives were accomplished by defining and executing four different tasks corresponding to these specific objectives. The first task involved selection of IHX candidates and developing steady state designs for those. The second task involved modeling of the transient and offnormal operation of the reactor-IHX system. The subsequent task dealt with the development of control strategies and involved algorithm development and simulation. The last task involved experimental validation of the thermal hydraulic performances of the two prototype heat exchangers designed and fabricated for the project at steady state and transient conditions to simulate the coupling of the reactor- IHX-process plant system. The experimental work utilized the two test facilities at The Ohio State University (OSU) including one existing High-Temperature Helium Test Facility (HTHF) and the newly developed high-temperature molten salt facility.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Integrated Reflection Seismic Monitoring and Reservoir Modeling for Geologic CO 2 Sequestration

The US DOE/NETL CCS MVA program funded a project with Fusion Petroleum Technologies Inc. (now SIGMA) to model the proof of concept of using sparse seismic data in the monitoring of CO 2 injected into saline aquifers. The goal of the project was to develop and demonstrate an active source reflection seismic imaging strategy based on deployment of spatially sparse surface seismic arrays. The primary objective was to test the feasibility of sparse seismic array systems to monitor the CO 2 plume migration injected into deep saline aquifers. The USDOE/RMOTC Teapot Dome (Wyoming) 3D seismic and reservoir data targeting the Crow Mountain formation was used as a realistic proxy to evaluate the feasibility of the proposed methodology. Though the RMOTC field has been well studied, the Crow Mountain as a saline aquifer has not been studied previously as a CO 2 sequestration (storage) candidate reservoir. A full reprocessing of the seismic data from field tapes that included prestack time migration (PSTM) followed by prestack depth migration (PSDM) was performed. A baseline reservoir model was generated from the new imaging results that characterized the faults and horizon surfaces of the Crow Mountain reservoir. The 3D interpretation was integrated with the petrophysical data from available wells and incorporated into a geocellular model. The reservoir structure used in the geocellular model was developed using advanced inversion technologies including Fusion's ThinMAN™ broadband spectral inversion. Seal failure risk was assessed using Fusion's proprietary GEOPRESS™ pore pressure and fracture pressure prediction technology. CO 2 injection was simulated into the Crow Mountain with a commercial reservoir simulator. Approximately 1.2MM tons of CO 2 was simulated to be injected into the Crow Mountain reservoir over 30 years and subsequently let 'soak' in the reservoir for 970 years. The relatively small plume developed from this injection was observed migrating due to gravity to the apexes of the double anticline in the Crow Mountain reservoir of the Teapot dome. Four models were generated from the reservoir simulation task of the project which included three saturation models representing snapshots at different times during and after simulated CO 2 injection and a fully saturated CO 2 fluid substitution model. The saturation models were used along with a Gassmann fluid substitution model for CO 2 to perform fluid volumetric substitution in the Crow Mountain formation. The fluid substitution resulted in a velocity and density model for the 3D volume at each saturation condition that was used to generate a synthetic seismic survey. FPTI's (Fusion Petroleum Technologies Inc.) proprietary SeisModelPRO™ full acoustic wave equation software was used to simulate acquisition of a 3D seismic survey on the four models over a subset of the field area. The simulated acquisition area included the injection wells and the majority of the simulated plume area.

02 PETROLEUM↗

The C 2 H 4 O isomers in the oxidation of ethylene

A combined rotational spectroscopy, thermochemistry, and kinetic modeling study explores the mechanism of acetaldehyde (ethanal), vinyl alcohol (ethenol), and ethylene oxide (oxirane) formation in the oxidation of ethylene (ethene). Multiplexed quantitative detection of oxidation intermediates and products exiting a SiO 2 / SiC microreactor at 1700 K is demonstrated with the BrightSpec W-band chirped-pulse Fourier transform millimeter-wave spectrometer. The broadband rotational spectrum contains transitions of formaldehyde (CH 2 O), methoxy (CH 3 O), methanol (CH 3 OH), ketene (CH 2 CO), acetaldehyde (CH 3 CHO), syn-vinyl alcohol (syn- CH 2 CHOH), anti-vinyl alcohol (anti-CH 2 CHOH), oxirane (c-CH 2 OCH 2 ), propyne (CH 3 CCH), and syn-propanal (syn-CH 3 CH 2 CHO). We focus on the three C 2 H 4 O species and deduce their concentration ratio [CH 3 CHO]:[CH 2 CHOH]:[c-CH 2 OCH 2 ] = 1:0.7(2):0.06(2) by comparing the observed line intensities to those simulated with the PGOPHER software. Detailed thermochemistry of the C 2 H 4 O isomers and conformers is provided by Active Thermochemical Tables. The observed excess concentrations of vinyl alcohol and oxirane relative to the more stable acetaldehyde compared to the equilibrium ratio at 1700 K (1:0.087:0.000024), point to direct chemical pathways to these higher energy isomers. The mechanism for the formation of the three C 2 H 4 O isomers is analyzed using a 0-D homogenous reactor kinetics simulation for ethylene oxidation. The ratios of the C 2 H 4 O isomers concentrations predicted by the kinetic model are compared to the experimental values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-time Distribution Simulation and Application Development for Power Systems Education

To help bridge the gap between traditional power system engineering instruction and emerging industry needs, a new set of classroom and research tools are needed. Real-time simulation tools emulating power system control room software present an opportunity to introduce students to the array of operational considerations, technical challenges, and decision-making associated with power system operations. The GridAPPS-D platform is proposed to support coursework and academic research as it provides an open-source platform for simulation, application development, and software integration. The GridAPPS-D platform, simulation capabilities, development environment, and interactive training are discussed in the context of lessons-learned from implementation for undergraduate and graduate students in the US and India. A series of potential GridAPPS-D supported academic capabilities are introduced.

Active distribution networks, open educational res↗

Evaluation of Computational Method of High Reynolds Number Slurry Flow for Caverns Backfilling

The abandonment of salt caverns used for brining or product storage poses a significant environmental and economic risk. Risk mitigation can in part be addressed by the process of backfilling which can improve the cavern geomechanical stability and reduce the risk of fluid loss to the environment. This study evaluates a currently available computational tool, Barracuda, to simulate such processes as slurry flow at high Reynolds number with high particle loading. Using Barracuda software, a parametric sequence of simulations evaluated slurry flow at Reynolds number up to 15000 and loading up to 25%. Limitations come into the long time required to run these simulations due in particular to the mesh size requirement at the jet nozzle. This study has found that slurry-jet width and centerline velocities are functions of Reynolds number and volume fraction. The solid phase was found to spread less than the water-phase with a spreading rate smaller than 1, dependent on the volume fraction. Particle size distribution does seem to have a large influence on the jet flow development. This study constitutes a first step to understand the behavior of highly loaded slurries and their ultimate application to cavern backfilling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Stand Alone Battery Thermal Management System

The objective of this project is research, development and demonstration of innovative thermal management concepts that reduce the cell or battery weight, complexity (component count) and/or cost by at least 20%. The project addresses two issues that are common problems with current state of the art lithium ion battery packs used in vehicles; low power at cold temperatures and reduced battery life when exposed to high temperatures. Typically, battery packs are “oversized” to satisfy the two issues mentioned above. The first phase of the project was spent making a battery pack simulation model using AMEsim software. The battery pack used as a benchmark was from the Fiat 500EV. FCA and NREL provided vehicle data and cell data that allowed an accurate model to be created that matched the electrical and thermal characteristics of the actual battery pack. The second phase involved using the battery model from the first phase and evaluate different thermal management concepts. In the end, a gas injection heat pump system was chosen as the dedicated thermal system to both heat and cool the battery pack. Based on the simulation model. The heat pump system could use 50% less energy to heat the battery pack in -20°C ambient conditions, and by keeping the battery cooler at hot climates, the battery pack size could be reduced by 5% and still meet the warranty requirements. During the final phase, the actual battery pack and heat pump system were installed in a test bench at DENSO to validate the simulation results. Also during this phase, the system was moved to NREL where testing was also done to validate the results. In conclusion, the heat pump system can improve “fuel economy” (for electric vehicle) by 12% average in cold climates. Also, the battery pack size, or capacity, could be reduced 5%, or if pack size is kept constant, the pack life could be increased by two years. Finally, the total battery pack and thermal system cost could be reduced 5% only if the system is integrated with the vehicle cabin air conditioning system. The reason why we were not able to achieve the 20% reduction target is because of the natural decay of the battery cell due to the number of cycles. Perhaps newer battery chemistries that are not so sensitive to cycling would have more potential for reducing the battery size due to thermal issues.

25 ENERGY STORAGE↗

Three-Dimensional Permeability of Thick-Section Glass Fabric Reinforced Polymer Composite by Vacuum-Assisted Resin Infusion Molding

Abstract Determination of permeability of thick-section glass fabric preforms with fabric layers of different architectures is critical for manufacturing large, thick composite structures with complex geometry, such as wind turbine blades. The thick-section reinforcement permeability is inherently three-dimensional and needs to be obtained for accurate composite processing modeling and analysis. Numerical simulation of the liquid stage of vacuum-assisted resin infusion molding (VARIM) is important to advance the composite manufacturing process and reduce processing-induced defects. In this research, the 3D permeability of thick-section E-glass fabric reinforcement preforms is determined, and the results are validated by a comparison between flow front progressions from experiments and from numerical simulations using ansys fluent software. The orientation of the principal permeability axes were unknown prior to experiments. The approach used in this research differs from those in literature in that the through-thickness permeability is determined as a function of flow front positions along the principal axes and the in-plane permeabilities and is not dependent on the inlet radius. The approach was tested on reinforcements with fabric architectures which vary through-the-thickness direction, such as those in a spar cap of a wind turbine blade. The computational simulations of the flow-front progression through-the-thickness were consistent with experimental observations.

Engineering↗

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗

Statistical variability and confidence intervals for planar dose QA pass rates

Purpose: The most common metric for comparing measured to calculated dose, such as for pretreatment quality assurance of intensity-modulated photon fields, is a pass rate (%) generated using percent difference (%Diff), distance-to-agreement (DTA), or some combination of the two (e.g., gamma evaluation). For many dosimeters, the grid of analyzed points corresponds to an array with a low areal density of point detectors. In these cases, the pass rates for any given comparison criteria are not absolute but exhibit statistical variability that is a function, in part, on the detector sampling geometry. In this work, the authors analyze the statistics of various methods commonly used to calculate pass rates and propose methods for establishing confidence intervals for pass rates obtained with low-density arrays. Methods: Dose planes were acquired for 25 prostate and 79 head and neck intensity-modulated fields via diode array and electronic portal imaging device (EPID), and matching calculated dose planes were created via a commercial treatment planning system. Pass rates for each dose plane pair (both centered to the beam central axis) were calculated with several common comparison methods: %Diff/DTA composite analysis and gamma evaluation, using absolute dose comparison with both local and global normalization. Specialized software was designed to selectively sample the measured EPID response (very high data density) down to discrete points to simulate low-density measurements. The software was used to realign the simulated detector grid at many simulated positions with respect to the beam central axis, thereby altering the low-density sampled grid. Simulations were repeated with 100 positional iterations using a 1 detector/cm{sup 2} uniform grid, a 2 detector/cm{sup 2} uniform grid, and similar random detector grids. For each simulation, %/DTA composite pass rates were calculated with various %Diff/DTA criteria and for both local and global %Diff normalization techniques. Results: For the prostate and head/neck cases studied, the pass rates obtained with gamma analysis of high density dose planes were 2%-5% higher than respective %/DTA composite analysis on average (ranging as high as 11%), depending on tolerances and normalization. Meanwhile, the pass rates obtained via local normalization were 2%-12% lower than with global maximum normalization on average (ranging as high as 27%), depending on tolerances and calculation method. Repositioning of simulated low-density sampled grids leads to a distribution of possible pass rates for each measured/calculated dose plane pair. These distributions can be predicted using a binomial distribution in order to establish confidence intervals that depend largely on the sampling density and the observed pass rate (i.e., the degree of difference between measured and calculated dose). These results can be extended to apply to 3D arrays of detectors, as well. Conclusions: Dose plane QA analysis can be greatly affected by choice of calculation metric and user-defined parameters, and so all pass rates should be reported with a complete description of calculation method. Pass rates for low-density arrays are subject to statistical uncertainty (vs. the high-density pass rate), but these sampling errors can be modeled using statistical confidence intervals derived from the sampled pass rate and detector density. Thus, pass rates for low-density array measurements should be accompanied by a confidence interval indicating the uncertainty of each pass rate.

61 RADIATION PROTECTION AND DOSIMETRY↗

Polarity effects and apparent ion recombination in microionization chambers

Purpose: Microchambers demonstrate anomalous voltage-dependent polarity effects. Existing polarity and ion recombination correction factors do not account for these effects. As a result, many commercial microchamber models do not meet the specification of a reference-class ionization chamber as defined by the American Association of Physicists in Medicine. The purpose of this investigation is to determine the cause of these voltage-dependent polarity effects. Methods: A series of microchamber prototypes were produced to isolate the source of the voltage-dependent polarity effects. Parameters including ionization-chamber collecting-volume size, stem and cable irradiation, chamber assembly, contaminants, high-Z materials, and individual chamber components were investigated. Measurements were performed with electrodes coated with graphite to isolate electrode conductivity. Chamber response was measured as the potential bias of the guard electrode was altered with respect to the collecting electrode, through the integration of additional power supplies. Ionization chamber models were also simulated using COMSOL Multiphysics software to investigate the effect of a potential difference between electrodes on electric field lines and collecting volume definition. Results: Investigations with microchamber prototypes demonstrated that the significant source of the voltage-dependent polarity effects was a potential difference between the guard and collecting electrodes of the chambers. The voltage-dependent polarity effects for each prototype were primarily isolated to either the guard or collecting electrode. Polarity effects were reduced by coating the isolated electrode with a conductive layer of graphite. Polarity effects were increased by introducing a potential difference between the electrodes. COMSOL simulations further demonstrated that for a given potential difference between electrodes, the collecting volume of the chamber changed as the applied voltage was altered, producing voltage-dependent polarity effects in the chamber response. Ionization chamber measurements and COMSOL simulations demonstrated an inverse relationship between the chamber collecting volume size and the severity of voltage-dependent polarity effects on chamber response. The effect of a given potential difference on chamber polarity effects was roughly ten times greater for microchambers as compared to Farmer-type chambers. Stem and cable irradiations, chamber assembly, contaminants, and high-Z materials were not found to be a significant source of the voltage-dependent polarity effects. Conclusions: A potential difference between the guard and collecting electrodes was found to be the primary source of the voltage-dependent polarity effects demonstrated by microchambers. For a given potential difference between electrodes, the relative change in the collecting volume is smaller for larger-volume chambers, illustrating why these polarity effects are not seen in larger-volume chambers with similar guard and collecting electrode designs. Thus, for small-volume chambers, it is necessary to reduce the potential difference between the guard and collecting electrodes in order to reduce polarity effects for reference dosimetry measurements.

60 APPLIED LIFE SCIENCES↗

A Software Toolkit to Study Systematic Uncertainties of the Physics Models of the Geant4 Simulation Package

The Geant4 toolkit is used to model interactions between particles and matter. Geant4 employs a set of validated physics models that span a wide range of interaction energies. These models are tuned to cover a large variety of possible applications. This raises the critical question of what uncertainties are associated with the Geant4 physics model, or group of models, involved in a simulation project. To address the challenge, we have designed and implemented a comprehen- sive, modular, user-friendly software toolkit that allows the variation of one or more parameters of one or more Geant4 physics models involved in simulation studies. It also enables analysis of multiple variants of the resulting physics observables of interest in order to estimate the uncertain- ties associated with the simulation model choices. Key functionalities of the toolkit are presented in this paper and are illustrated with selected results.

Genser, Krzysztof↗

Designing a Framework for Solving Multiobjective Simulation Optimization Problems

Multiobjective simulation optimization (MOSO) problems are optimization problems with multiple conflicting objectives, where evaluation of at least one of the objectives depends on a black-box numerical code or real-world experiment, which we refer to as a simulation. Whereas an extensive body of research is dedicated to developing new algorithms and methods for solving these and related problems, it is challenging and time-consuming to integrate these techniques into real-world production-ready solvers. This is partly because of the diversity and complexity of modern state-of-the-art MOSO algorithms and methods and partly because of the complexity and specificity of many real-world problems and their corresponding computing environments. The complexity of this problem is only compounded when introducing potentially complex and/or domain-specific surrogate-modeling techniques, problem formulations, design spaces, and data acquisition functions. Here, this paper carefully surveys the current state of the art in MOSO algorithms, techniques, and solvers, as well as problem types and computational environments where MOSO is commonly applied. We then present several key challenges in the design of a parallel multiobjective simulation optimization framework (ParMOO) and how they have been addressed. Finally, we provide two case studies demonstrating how customized ParMOO solvers can be quickly built and deployed to solve real-world MOSO problems.

engineering design optimization↗

Open Source Platform for Live Processing of High Data Rate Electron Microscopy (SBIR Phase I Award Topic 15, Subtopic B, Final Scientific / Technical Report)

4D scanning transmission electron microscopy (STEM) is experiencing a revolution both in terms of data size, and data veracity, as are many other areas of experimental data acquisition. This presents enormous opportunities to make new scientific discoveries, but is also hamstrung by aging software that is no longer fit for purpose. Experimental data acquisition platforms must make the same investments in software infrastructure that the simulation community made decades ago to take advantage of high-performance computing (HPC) if they hope to realize the full potential of the latest advances made in detector technology. This requires rethinking the design and deployment of hardware, moving away from small, fixed commodity computing capabilities with closed interfaces to open, modular, scalable systems designed for parallelism and throughput. The technical approach will leverage best-in-class open source technologies, and develop new approaches strategically to build an open source platform for experimental data acquisition targeting high data rate electron microscopy. The project will develop a software platform that uses the power and performance of modern HPC resources to ingest the large volume of data produced by electron microscopes using direct electron detectors. A set of configurable pipelines will be used to enable the processing to be tailored to the problem being investigated with builtin visualization of data at interactive rates. This project will leverage investments made by the DOE and other agencies to provide a powerful software platform for the experimental data community working with most powerful electron microscopes, such as 4D STEM detectors. It will create a platform that serves as a reference for the community. Permissive open source licensing will make the platform available to all, with a simple deployment strategy offering the ability to keep data close to where it is generated. This advanced software platform will require customization for each deployment, creating opportunities for software services that are typically offered to extend open source software platforms.

97 MATHEMATICS AND COMPUTING↗

Simulation of optimum parameters for GaN MSM UV photodetector

In this study the optimum parameters of GaN M-S-M photodetector are discussed. The evaluation of the photodetector depends on many parameters, the most of the important parameters the quality of the GaN film and others depend on the geometry of the interdigited electrode. In this simulation work using MATLAB software with consideration of the reflection and absorption on the metal contacts, a detailed study involving various electrode spacings (S) and widths (W) reveals conclusive results in device design. The optimum interelectrode design for interdigitated MSM-PD has been specified and evaluated by effect on quantum efficiency and responsivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Software-Defined Ultrasonic Communication System Based on Time-reversal Signal Processing

Ultrasonic communication through solid channels suffers from channel dispersion, attenuation, and multipath effect. For example, in elbow shape pipes or complex structures, the elastic waves are reflected and reverberated in multipath resulting in a major dispersion of the signals. In solid channels such as flat bar, multipath propagation is the main reason for bit-wavelet spreading and intersymbol interference (ISI). The conventional modulation method like OOK (on-off-keying) or PSK (phase-shift-keying) fails to provide a useful Signal-to-Noise Ratio (SNR). Time-reversal (TR) signal processing method is a promising technique that can compress the spreading of the bit-wavelet and strengthen the peak signal energy on the receiver. In this study, Finite Element Analysis (FEA) using ABAQUS software tools is used to characterize the bit-wavelet spreading effect caused by dispersive communication channels. Also, this study presents the impact of TR ultrasonic communication as a function of excitation pulse width, channel length, and carrier frequency through a flat-bar channel. A Software-Defined System-on-Chip (SD-SoC) platform is examined as a viable solution for the low-complexity and energy-efficient ultrasonic communication system. The proposed TR based SD-SoC ultrasonic communication system can make full use of signal dispersion for peak performance by minimizing the ISI within a dispersive and reverberant metal channel.

FEM simulation↗

OpenFAST [SWR-08-12]

OpenFAST is a wind turbine simulation tool which builds on NREL's FAST v8 software. NREL's FAST.Farm software extends the capability of OpenFAST to simulate multi-turbine wind farms. They were created with the goal of being community models developed and used by research laboratories, academia, and industry. They are managed by a dedicated team at the National Renewable Energy Laboratory (NREL). Our objective is to ensure that OpenFAST and FAST.Farm are sustainable software that are well tested and well documented. If you'd like to contribute, see the Developer Documentation and any open GitHub issues with the Help Wanted tag.

Jonkman, Jason↗

Geomagnetically Induced Current Field Test on Large Grid-Connected Power Transformers: Analysis, Model Development, and Simulations

Geomagnetic-induced current (GIC) flow in power grids can cause undesirable effects such as transformer overheating, harmonics, higher reactive power demand, etc. Many simulation models have been developed to study these effects, but real-world verification on modern transformer designs is rare. Here, this paper presents the first long-duration GIC field test in the U.S. performed on high-voltage, grid-connected transformers featuring winding clamps and tie rods instead of conventional tie bars. Field measurements were taken to evaluate GIC effects. These measurements also aided in developing and validating thermal and electromagnetic transient (EMT) models of the transformers. During the test, significant current and voltage distortions were observed along with considerable transformer reactive power losses. Analysis of the field measurements showed that the transformers’ hottest spot was at the inner windings, and their k-factors were close to factory test and software default values. Thermal simulations indicated that the transformers would not violate their thermal limits even for a GIC waveform that peaks at about 200 A/phase. EMT simulations revealed that increased transformer loading may reduce GIC-induced reactive power demand and harmonics in certain scenarios. The study also highlighted potential inaccuracies in using the k-factor method to calculate transformer reactive power losses.

EMTDC↗