Transient modeling and simulation of a nonisothermal sodium–sulfur cell
Not Available
SEARCH · Engineering Papers
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.
Not Available
Radio frequency (RF) waves can offer heating, current and flow drive, as well as instability control for steady state operations of fusion experiments. A particle simulation model has been designed in this work to provide a first-principles tool for studying the RF nonlinear interactions with plasmas. In this model, ions are considered as fully kinetic particles using the Vlasov equation and electrons are treated as guiding centers using the drift kinetic equation. This model has been implemented in a global gyrokinetic toroidal code using real electron-to-ion mass ratio. To confirm the model, linear simulations of ion plasma oscillation, ion Bernstein wave, and lower hybrid wave are carried out in cylindrical geometry and found to agree well with analytic predictions.
A piloted flight simulator called the Ice Contamination Effects Flight Training Device (ICEFTD), which uses low cost desktop components and a generic cockpit replication is being developed. The purpose of this device is to demonstrate the effectiveness of its use for training pilots to recognize and recover from aircraft handling anomalies that result from airframe ice formations. High-fidelity flight simulation models for various baseline (non-iced) and iced configurations were developed from wind tunnel tests of a subscale DeHavilland DHC-6 Twin Otter aircraft model. These simulation models were validated with flight test data from the NASA Twin Otter Icing Research Aircraft, which included the effects of ice on wing and tail stall characteristics. These simulation models are being implemented into an ICEFTD that will provide representative aircraft characteristics due to airframe icing. Scenario-based exercises are being constructed to give an operational-flavor to the simulation. Training pilots will learn to recognize iced aircraft characteristics from the baseline, and will practice and apply appropriate recovery procedures to a handling event.
A summary of project activities relative to the estimation of potato yields in the Columbia Basin is given. Oregon State University is using a two-pronged approach to yield estimation, one using simulation models and the other using purely empirical models. The simulation modeling approach has used satellite observations to determine key dates in the development of the crop for each field identified as potatoes. In particular, these include planting dates, emergence dates, and harvest dates. These critical dates are fed into simulation models of crop growth and development to derive yield forecasts. Two empirical modeling approaches are illustrated. One relates tuber yield to estimates of cumulative intercepted solar radiation; the other relates tuber yield to the integral under the GVI curve.
This paper describes comprehensive analytical and simulation models for the design and optimization of the electron-injection based detectors. The electron-injection detectors evaluated here operate in the short-wave infrared range and utilize a type-II band alignment in InP/GaAsSb/InGaAs material system. The unique geometry of detectors along with an inherent negative-feedback mechanism in the device allows for achieving high internal avalanche-free amplifications without any excess noise. Physics-based closed-form analytical models are derived for the detector rise time and dark current. Our optical gain model takes into account the drop in the optical gain at high optical power levels. Furthermore, numerical simulation studies of the electrical characteristics of the device show good agreement with our analytical models as well experimental data. Performance comparison between devices with different injector sizes shows that enhancement in the gain and speed is anticipated by reducing the injector size. Sensitivity analysis for the key detector parameters shows the relative importance of each parameter. The results of this study may provide useful information and guidelines for development of future electron-injection based detectors as well as other heterojunction photodetectors.
Carbon dioxide (CO 2 ) enhanced oil recovery (EOR) will likely be the primary means of geologic CO 2 storage during the early stages of commercial-scale carbon capture and storage (CCS) because of the inherent economic incentives as well as the abundant experience and demonstrated success in the United States, where CO 2 EOR has been employed since 1974. The work presented here estimates CO 2 storage efficiency factors in CO 2 EOR operations using a unique industry database of CO 2 EOR sites and 12 different reservoir simulation models. The simulation models encompass fluvial clastic and shallow shelf carbonate depositional environments for reservoir depths of 1219 and 2438 m (4000 and 8000 feet) and 7.6-, 20-, and 64-m (25-, 66-, and 209-foot)-thick pay zones. A novel statistical modeling technique incorporating the Michaelis–Menten function is used to generate empirical percentile estimates of CO 2 storage efficiency factors. West Texas San Andres dolomite water alternating gas (WAG) CO 2 flood performance data were used to derive P10, P50, and P90 CO 2 storage efficiency factors of 0.76, 1.28, and 1.74 Mscf/STB (stock tank barrel) of original oil in place. Median CO 2 storage efficiency factors from continuous CO 2 injection following conventional waterflood varied from 15% to 61% and 8% to 40% for fluvial clastic and shallow shelf carbonate simulation models, respectively, while those from WAG injection varied from 14% to 42% and 8% to 31%, respectively. Variation in the CO 2 storage efficiency factors was largely attributable to reservoir depth (a surrogate for reservoir pressure and temperature) and lithology (clastic versus carbonate). The results of this work provide practical information that can be used to quantify CO 2 storage resource estimates in oil reservoirs during CO 2 EOR operations (as opposed to storage following depletion) and the uncertainty associated with those estimates.
With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.
We compare results from six AOMIP model simulations with estimates of sea ice thickness obtained from ICESat, moored and submarine-based upward looking sensors, airborne electromagnetic measurements and drill holes. Our goal is to find patterns of model performance to guide model improvement. The satellite data is pan-arctic from 2004-2008, ice-draft data is from moored instruments in Fram Strait, the Greenland Sea and the Beaufort Sea from 1992-2008 and from submarines from 1975-2000. The drill hole data are from the Laptev and East Siberian marginal seas from 1982-1986 and from coastal stations from 1998-2009. While there are important caveats when comparing modeled results with measurements from different platforms and time periods such as these, the models agree well with moored ULS data. In general, the AOMIP models underestimate the thickness of measured ice thicker than about 2 m and overestimate thickness of ice thinner than 2 m. The simulated results are poor over the fast ice and marginal seas of the Siberian shelves. Averaging over all observational data sets, the better correlations and smaller differences from observed thickness are from the ECCO2 and UW models.
Ten 3D cloud-resolving model (CRM) simulations and four 3D limited area model (LAM) simulations of an intense mesoscale convective system observed on January 23-24, 2006 during the Tropical Warm Pool – International Cloud Experiment (TWP-ICE) are compared with each other and with observed radar reflectivity fields and dual-Doppler retrievals of vertical wind speeds in an attempt to explain published results showing a high bias in simulated convective radar reflectivity aloft. This high bias results from ice water content being large, which is a product of large, strong convective updrafts, although hydrometeor size distribution assumptions modulate the size of this bias. Snow reflectivity can exceed 40 dBZ in a two-moment scheme when a constant bulk density of 100 kg m-3 is used. Making snow mass more realistically proportional to area rather than volume should somewhat alleviate this problem. Graupel, unlike snow, produces high biased reflectivity in all simulations. This is associated with large amounts of liquid water above the freezing level in updraft cores. Peak vertical velocities in deep convective updrafts are greater than dual-Doppler retrieved values, especially in the upper troposphere. Freezing of large rainwater contents lofted above the freezing level in simulated updraft cores greatly contributes to these excessive upper tropospheric vertical velocities. Strong simulated updraft cores are nearly undiluted, with some showing supercell characteristics. Decreasing horizontal grid spacing from 900 meters to 100 meters weakens strong updrafts, but not enough to match observational retrievals. Therefore, overly intense simulated updrafts may partly be a product of interactions between convective dynamics, parameterized microphysics, and large-scale environmental biases that promote different convective modes and strengths than observed.
In preparation for next generation North American Land Data Assimilation System (NLDAS), 3 three advanced land surface models (CLM4.0, Noah-MP, and CLSM-F2.5) were run from 1979 4 to 2014 within the NLDAS-based framework. Monthly total water storage anomaly (TWSA) and 5 its individual water storage components were evaluated against satellite-based and in situ 6 observations, and reference reanalysis products at basin-wide and statewide scales. In general, all 7 three models are able to reasonably capture the monthly and interannual variability and 8 magnitudes for TWSA. However, contributions of the anomalies of individual water 9 components to TWSA are very dependent on the model and basin. A major contributor to the 10 TWSA is the anomaly of total column soil moisture content (SMCA) for CLM4.0 and Noah-MP 11 or groundwater storage anomaly (GWSA) for CLSM-F2.5 although other components such as 12 the anomaly of snow water equivalent (SWEA) also play some role. For each individual water 13 storage component, the models are able to capture broad features such as monthly and 14 interannual variability. However, there are large inter-model differences and quantitative 15 uncertainties in this study. Therefore, it should be thought of as a preliminary synthesis and 16 analysis.
Very high rates of CO2 absorption have been demonstrated using semipermeable microcapsules filled with environmentally benign sodium carbonate solutions (Vericella et al., Nature communications 6, 2015). However, such capsules have also exhibited several traits, including water loss/uptake, elastic swelling and buckling, and precipitation of solids, which make them more complicated to engineer as a carbon capture material. To address this, a mechanistic model for mass transfer and chemical reaction has been developed, which accounts for these behaviors, and can be used as a predictive tool to explore capsule performance under different CO2 capture scenarios. The model uses a modified film theory to describe CO2 and H2O mass transfer rates, a concentration based description of vapor-liquid-solid equilibrium inside the capsule, and an elastic description of the capsule shell. Model predictions compare favorably to literature data for precipitating carbonates, as well as recent data for microcapsule CO2 absorption. The model was used to assess practical issues of capsule water transfer during a temperature swing CO2 absorption/desorption cycle. Model predictions indicate that capsules filled with high weight percent, precipitating sodium carbonate solution could be significantly or completely dehydrated during CO2 absorption. Subsequent rehydration of lean capsules at elevated temperature showed strong sensitivity to both gas humidity, and capsule shell elasticity.
Abstract The Pi Cloud Chamber offers a unique opportunity to study aerosol‐cloud microphysics interactions in a steady‐state, turbulent environment. In this work, an atmospheric large‐eddy simulation (LES) model with spectral bin microphysics is scaled down to simulate these interactions, allowing comparison with experimental results. A simple scalar flux budget model is developed and used to explore the effect of sidewalls on the bulk mixing temperature, water vapor mixing ratio, and supersaturation. The scaled simulation and the simple scalar flux budget model produce comparable bulk mixing scalar values. The LES dynamics results are compared with particle image velocimetry measurements of turbulent kinetic energy, energy dissipation rates, and large‐scale oscillation frequencies from the cloud chamber. These simulated results match quantitatively to experimental results. Finally, with the bin microphysics included the LES is able to simulate steady‐state cloud conditions and broadening of the cloud droplet size distributions with decreasing droplet number concentration, as observed in the experiments. The results further suggest that collision‐coalescence does not contribute significantly to this broadening. This opens a path for further detailed intercomparison of laboratory and simulation results for model validation and exploration of specific physical processes.
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.
Applications, validation tests, and upgrades of the two- and three-dimensional system level thermal mathematical system simulation models (TMSSM) used for thermal protection system (TPS) analyses are described. The TMSSM were developed as an aid to predicting the performance requirements and configurations of the Shuttle wing leading edge (WLE) and nose cone (NC) TPS tiles. The WLE and its structure were subjected to acoustic, thermal/vacuum, and air loads tests to simulate launch, on-orbit, and re-entry behavior. STS-1, -2 and -5 flight data led to recalibration of on-board instruments and raised estimates of the thermal shock at the NC and WLE. Baseline heating data are now available for the design of future TPS.
Fusion whole device modeling simulations require comprehensive models that are simultaneously physically accurate, fast, robust, and predictive. In this paper we describe the development of two neural-network (NN) based models as a means to perform a snon-linear multivariate regression of theory-based models for the core turbulent transport fluxes, and the pedestal structure. Specifically, we find that a NN-based approach can be used to consistently reproduce the results of the TGLF and EPED1 theory-based models over a broad range of plasma regimes, and with a computational speedup of several orders of magnitudes. These models are then integrated into a predictive workflow that allows prediction with self-consistent core-pedestal coupling of the kinetic profiles within the last closed flux surface of the plasma. Finally, the NN paradigm is capable of breaking the speed-accuracy trade-off that is expected of traditional numerical physics models, and can provide the missing link towards self-consistent coupled core-pedestal whole device modeling simulations that are physically accurate and yet take only seconds to run.
Finite element (FE) computational models of the whole human body are being used to investigate the injury modalities and mechanisms arising from suited astronauts’ ExtraVehicular Activities (EVAs). These FE whole body models – simplified Elemance and THUMS pedestrian versions – have undergone a series of verification and validation (V&V) assessments for applications in automotive, domestic, and contact sports arenas. However, the FE models’ Modeling and Simulation (M&S) credibility assessments for simulating suited astronaut injuries due to fall from heights in altered gravity environments are limited. The current study evaluates the M&S credibility of simplified Elemance and THUMS pedestrian FE models based on eight credibility factors – data pedigree, input pedigree, code verification, solution verification, conceptual validation, referent validation, results uncertainty, and results robustness – as per NASA-STD-7009A [1]. These credibility factors were assigned ordinal scores ranging between 0 to 4, with a score of 0 in the case of insufficient evidence and 4 when all data for the model is known, used appropriately, and traceable to the Real-World System (RWS) with acceptable accuracy, precision, and uncertainty. The credibility factor sufficiency thresholds were assigned upon guidance from NASA EVA subject matter experts. The factor scores were ascertained through the evaluation of M&S credibility assessment procedures described in the published literature and FE model reports. The M&S credibility scores for the Elemance and THUMS FE models vary between 1 and 3 for data pedigree, input pedigree, and referent validation credibility factors, and between 0 and 1 for code verification, solution verification, conceptual validation, results uncertainty and results robustness. The results show that these FE models are credible for simulating fall from heights injury scenarios in specific kinematic and kinetic ranges that correlate with highway and contact sports-related external conditions. Multiple credibility elevation strategies are prescribed to improve the Elemance and THUMS models’ credibility for use in the EVA fall from height injury scenarios.
This slide presentation is composed of two topics. The first reviews the use of modeling and simulation (M&S) particularly as it relates to the Constellation program and discrete event simulation (DES). DES is defined as a process and system analysis, through time-based and resource constrained probabilistic simulation models, that provide insight into operation system performance. The DES shows that the cycles for a launch from manufacturing and assembly to launch and recovery is about 45 days and that approximately 4 launches per year are practicable. The second topic reviews a NASA Standard for Modeling and Simulation. The Columbia Accident Investigation Board made some recommendations related to models and simulations. Some of the ideas inherent in the new standard are the documentation of M&S activities, an assessment of the credibility, and reporting to decision makers, which should include the analysis of the results, a statement as to the uncertainty in the results,and the credibility of the results. There is also discussion about verification and validation (V&V) of models. There is also discussion about the different types of models and simulation.
A brief introduction to plasma simulation using computers and the difficulties on currently available computers is given. Through the use of an analyzing and measuring methodology - SARA, the control flow and data flow of a particle simulation model REM2-1/2D are exemplified. After recursive refinements the total execution time may be greatly shortened and a fully parallel data flow can be obtained. From this data flow, a matched computer architecture or organization could be configured to achieve the computation bound of an application problem. A sequential type simulation model, an array/pipeline type simulation model, and a fully parallel simulation model of a code REM2-1/2D are proposed and analyzed. This methodology can be applied to other application problems which have implicitly parallel nature.