Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling 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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The influence of dynamics on two-dimensional model results - Simulations of C-14 and stratospheric aircraft NO(x) injections

The influence of dynamics on model results is investigated via a 2D photochemical model. Three representations of atmospheric transport are tested in order to simulate total ozone and C-14 amounts after nuclear tests in the early 1960s. Three scenarios of NO(x) injections from a proposed fleet of stratospheric aircraft and their effects on ozone are simulated. The three dynamical formulations used are Dynamics A, a base dynamics used in previous work with this model; Dynamics B, a strong circulation dynamics discussed by Jackman et al. (1989); and Dynamics C, the dynamics used by Shia et al. (1989). The advective component of the stratosphere to the troposphere mass exchange rate is largest for Dynamics B (5.8 x 10 exp 17 kg/yr) and smallest for Dynamics C (1.4 x 10 exp 17 kg/yr), with the advective strat/trop exchange rate for Dynamics A (2.4 x 10 exp 17 kg/yr) being these two extremes. Ozone depletion from NO(x) injections of stratospheric aircraft showed a strong sensitivity to dynamics.

Jackman, Charles H.↗
DOE OSTI2016

High-fidelity simulation-driven model development for coarse-grained computational fluid dynamics

Nuclear reactor safety analysis requires analysis of a broad range of accident scenarios and determining their consequences. For a nuclear power plant behavior, it is impossible to obtain experimental data at sufficiently large scales to support calibration and validation of various system models. In single-phase flow convective problems, high-resolution methods of Computational Fluid Dynamics (CFD) such as Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) can provide high fidelity results under conditions where physical measurements are unavailable. However, such high-resolution simulations are computationally expensive and not suitable for simulation of long transient scenarios in nuclear reactor accidents. In this work, we investigate the use of a high fidelity simulation-driven approach to model sub-grid scale (SGS) effect in Coarse Grained Computational Fluid Dynamics CG-CFD. Notably, fine-mesh simulations are used to construct physics-informed statistical surrogate of coarse-mesh SGS model. For an initial analysis, we consider a case of turbulent natural convection in a volumetrically heated fluid layer with a thermally insulated lower boundary and isothermal upper boundary. This scenario of unstable stratification is relevant to turbulent natural convection in a severe nuclear reactor accident, as well as in containment mixing and passive cooling. For this type of flow, the SGS effect is modeled by an added turbulent diffusivity in the energy equation. It is shown that a global correction for the energy equation is sufficient to achieve a significant improvement to the CG-CFD prediction of thermal mixing in the fluid layer. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Evaluating Clouds, Aerosols, and their Interactions in Three Global Climate Models using COSP and Satellite Observations

Accurately representing aerosol-cloud interactions in global climate models is challenging. As parameterizations evolve, it is important to evaluate their performance with appropriate use of observations. In this work we compare aerosols, clouds, and their interactions in three climate models (AM3, CAM5, ModelE) to MODIS satellite observations. Modeled cloud properties were diagnosed using the CFMIP Observations Simulator Package (COSP). Cloud droplet number concentrations (N) were derived using the same algorithm for both satellite-simulated model values and observations. We find that aerosol optical depth tau simulated by models is similar to observations. For N, AM3 and CAM5 capture the observed spatial pattern of higher values in near-coast versus remote ocean regions, though modeled values in general are higher than observed. In contrast, ModelE simulates lower N in most near-coast versus remote regions. Aerosol- cloud interactions were computed as the sensitivity of N to tau for marine liquid clouds off the coasts of South Africa and Eastern Asia where aerosol pollution varies in time. AM3 and CAM5 are in most cases more sensitive than observations, while the sensitivity for ModelE is statistically insignificant. This widely used sensitivity could be subject to misinterpretation due to the confounding influence of meteorology on both aerosols and clouds. A simple framework for assessing the N – tau sensitivity at constant meteorology illustrates that observed sensitivity can change from positive to statistically insignificant when including the confounding influence of relative humidity. Satellite simulated values of N were compared to standard model output and found to be higher with a bias of 83 cm-3.

aerosols↗
DOE OSTI2013

A model to simulate the oxygen distribution in hypoxic tumors for different vascular architectures

Purpose: As hypoxic cells are more resistant to photon radiation, it is desirable to obtain information about the oxygen distribution in tumors prior to the radiation treatment. Noninvasive techniques are currently not able to provide reliable oxygenation maps with sufficient spatial resolution; therefore mathematical models may help to simulate microvascular architectures and the resulting oxygen distributions in the surrounding tissue. Here, the authors present a new computer model, which uses the vascular fraction of tumor voxels, in principle measurable noninvasively in vivo, as input parameter for simulating realistic PO2 histograms in tumors, assuming certain 3D vascular architectures.Methods: Oxygen distributions were calculated by solving a reaction-diffusion equation in a reference volume using the particle strength exchange method. Different types of vessel architectures as well as different degrees of vascular heterogeneities are considered. Two types of acute hypoxia (ischemic and hypoxemic) occurring additionally to diffusion-limited (chronic) hypoxia were implemented as well.Results: No statistically significant differences were observed when comparing 2D- and 3D-vessel architectures (p > 0.79 in all cases) and highly heterogeneously distributed linear vessels show good agreement, when comparing with published experimental intervessel distance distributions and PO2 histograms. It could be shown that, if information about additional acute hypoxia is available, its contribution to the hypoxic fraction (HF) can be simulated as well. Increases of 128% and 168% in the HF were obtained when representative cases of ischemic and hypoxemic acute hypoxia, respectively, were considered in the simulations.Conclusions: The presented model is able to simulate realistic microscopic oxygen distributions in tumors assuming reasonable vessel architectures and using the vascular fraction as macroscopic input parameter. The model may be used to generate PO2 histograms, which are needed as input in models predicting the radiation response of hypoxic tumors.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2023

Field validation of isotropic analytical models for simulating fabric shades

Fabric roller shades are common shading materials used in commercial and residential buildings. Accurately characterizing and modeling shades helps practitioners select the appropriate product and its control strategy based on climate and occupants' priorities, such as visual comfort and view to outdoors. Previous studies established a generalized method for modeling complex fenestration systems using data-driven tabulated bidirectional scattering distribution functions. However, deploying such a method at scale to all fabric shading products on the market is too costly and time-consuming. Analytical models that are based on a limited set of measurements (e.g., normal-normal and normal-hemispherical visible transmittance and reflectance, and directional cut-off angles) can be used to model the wide variety of shading products on the market. This study evaluates the performance of two isotropic analytical models, Roos-Wienold and Modified-Kotey, for modeling fabric roller shades, with a focus on the model's ability to predict occupant visual comfort. The performance evaluation was conducted through laboratory and field measurements and simulations. The results showed that both models are sufficient for predicting vertical illuminance at seated eye-level. Roos-Wienold model was able to predict binary visual comfort classification (glare/no-glare) under a wide range of luminance conditions, while Modified-Kotey model did not perform as well under high-contrast low-adaptation conditions. Both models are insufficient in predicting visual comfort at a four-point scale (e.g., imperceptible, perceptible, disturbing, intolerable). The two isotropic models become less accurate when the fabric exhibits high anisotropy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1986

A lumped parameter mathematical model for simulation of subsonic wind tunnels

Equations for a lumped parameter mathematical model of a subsonic wind tunnel circuit are presented. The equation state variables are internal energy, density, and mass flow rate. The circuit model is structured to allow for integration and analysis of tunnel subsystem models which provide functions such as control of altitude pressure and temperature. Thus the model provides a useful tool for investigating the transient behavior of the tunnel and control requirements. The model was applied to the proposed NASA Lewis Altitude Wind Tunnel (AWT) circuit and included transfer function representations of the tunnel supply/exhaust air and refrigeration subsystems. Both steady state and frequency response data are presented for the circuit model indicating the type of results and accuracy that can be expected from the model. Transient data for closed loop control of the tunnel and its subsystems are also presented, demonstrating the model's use as a control analysis tool.

Krosel, S. M.↗
NASA NTRS1992

Application of algebraic-RNG eddy viscosity model to simulation of transitional boundary layer flow

An algebraic eddy-viscosity model is derived from the renormalization group (RNG) theory of turbulence. A new length scale, based on boundary layer characteristics (displacement thickness, shape factor), is proposed. The model was applied to transitional boundary layer flow over a flat plate. Integral characteristics, such as the total skin friction coefficient, and mean velocity profile across the boundary layer, are found to be in good agreement with experimental data.

Yakhot, Alexander↗
DOE OSTI2020

Spatiotemporal dynamics of ecosystem fires and biomass burning-induced carbon emissions in China over the past two decades

Fire is a major type of disturbance that has important influences on ecosystem dynamics and carbon cycles. Yet our understanding of ecosystem fires and their carbon cycle consequences is still limited, largely due to the difficulty of large-scale fire monitoring and the complex interactions between fire, vegetation, climate, and anthropogenic factors. Here, using data from satellite-derived fire observations and ecosystem model simulations, we performed a comprehensive investigation of the spatial and temporal dynamics of China's ecosystem fire disturbances and their carbon emissions over the past two decades (1997–2016). Satellite-derived results showed that on average about 3.47 - 4.53 x 10 4 km 2 of the land was burned annually during the past two decades, among which annual burned forest area was about 0.81 - 1.25 × 10 4 km 2 , accounting for 0.33-0.51% of the forest area in China. Biomass burning emitted about 23.02 TgC per year. Compared to satellite products, simulations from the Energy Exascale Earth System Land Model (ELM) strongly overestimated China's burned area and fire-induced carbon emissions. Annual burned area and fire-induced carbon emissions were high for boreal forest in Northeast China's Daxing'anling region and subtropical dry forest in South Yunnan, as revealed by both the satellite product and the model simulations. Our results suggest that climate and anthropogenic factors play critical roles in controlling the spatial and seasonal distribution of China's ecosystem fire disturbances. Our findings highlight the importance of multiple complementary approaches in assessing ecosystem fire disturbance and its carbon consequences. Further studies are required to improve the methods of observing and modelling China's ecosystem fire disturbances, which will provide valuable information for fire management and ecosystem sustainability in an era when both human activities and the natural environment are rapidly changing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stochastic simulations of temporally and spatially variable Fe cycling within a floodplain aquifer

This data package contains model input and output files for local and global sensitivity analysis (SA) of a reactive transport model simulating redox cycling within a floodplain aquifer. The reaction network is implemented in CrunchFlow and focuses on spatially and temporally heterogeneous Fe cycling, though it contains several other species and redox pathways (37 aqueous species and 7 minerals in total). The package contains data from 600 Monte Carlo simulations (for global SA) and 104 one-at-a-time simulations (for local SA). It includes all input files required to run the CrunchFlow model, as well as limited model output – .csv files containing time series of primary species concentrations at 21 grid cells within the domain. These time series files contain the following species: H+, CO32-, SO42-, Cl-, Ca2+, Mg2+, Na+, K+, Fe(II), Fe(III), HS-, CH3COO-, O2(aq), S(aq), and Br-.These data were used to demonstrate the application of a novel form of global SA, distance-based generalized sensitivity analysis (DGSA) to a reactive transport model. The 2D model simulates the export of reduced species from pockets of fine-grained, organic-rich sediments embedded in a coarse sand aquifer. Stochastic simulations jointly varied 17 key model input parameters across several orders of magnitude, including the spatial variability of each parameter. DGSA on the model results reveal that the amplitude and variability of exported Fe(II) are most sensitive to the interaction between sand permeability and individual reaction rates. By contrast, the propagation of reducing conditions downgradient from the fine-grained lenses depends most heavily on the rate of dissolved organic carbon and sulfate release from the fines.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Rendezvous algorithms for large-scale modeling and simulation

Rendezvous algorithms encode a communication pattern that is useful when processors sending data do not know who the receiving processors should be, or vice versa. The idea is to define an intermediate decomposition where datums from different sending processors can ”rendezvous” to perform a computation, in a manner that both the senders and eventual receivers of the results can identify the appropriate rendezvous processor. Though they were originally designed for interpolating between overlaid grids with independent parallel decompositions (Plimpton et al., 2004), we have recently found rendezvous algorithms useful for a variety of operations in particle- or grid-based simulation codes when running large problems on large numbers of processors. In particular, we show they can perform well when a load-balanced intermediate decomposition is randomized and not spatial, requiring all-to-all communication to move data between processors. In this case rendezvous algorithms leverage the large bisection communication bandwidths which parallel machines provide. We describe how rendezvous algorithms work in a scientific computing context and give specific examples for molecular dynamics and Direct Simulation Monte Carlo codes which result in dramatic performance improvements versus simpler algorithms which do not scale as well. We explain how a generic rendezvous algorithm can be implemented, and also point out similarities with the MapReduce paradigm popularized by Google and Hadoop.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2014

A Computational Model for Simulating Spaceflight Induced Bone Remodeling

An overview of an initial development of a model of bone loss due to skeletal unloading in weight bearing sites is presented. The skeletal site chosen for the initial application of the model is the femoral neck region because hip fractures can be debilitating to the overall performance health of astronauts. The paper begins with the motivation for developing such a model of the time course of change in bone in order to understand the mechanism of bone demineralization experienced by astronauts in microgravity, to quantify the health risk, and to establish countermeasures. Following this, a general description of a mathematical formulation of the process of bone remodeling is discussed. Equations governing the rate of change of mineralized bone volume fraction and active osteoclast and osteoblast are illustrated. Some of the physiology of bone remodeling, the theory of how imbalance in remodeling can cause bone loss, and how the model attempts to capture this is discussed. The results of a preliminary validation analysis that was carried out are presented. The analysis compares a set of simulation results against bone loss data from control subjects who participated in two different bed rest studies. Finally, the paper concludes with outlining the current limitations and caveats of the model, and planned future work to enhance the state of the model.

Pennline, James A.↗
DOE OSTI2022

Efficient Multi-Scale, Thermo-Hydro-Mechanical Numerical Model for Simulation of Long-Term Stability of Rooms for Nuclear Waste Disposal

This project demonstrated the feasibility of developing a numerical modeling tool that provides an efficient and unique approach to analyze the geomechanical performance of a deep geological repository (DGR) for high-level nuclear waste. The numerical modeling tool is physics-based and models the thermo-hydro-mechanical (THM) processes during the required time (up to one million years) for a DGR. The unique capability of the developed code is explicit representation of fracturing and a discrete fracture network (DFN) in the model that includes the entire repository layout for simulated time of one million years. Therefore, it can potentially address one of the challenges of analysis of a DGR, which is efficient numerical simulation of a multi-scale, physics-based models.

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

Macquarie Island Cloud and Radiation Experiment (MICRE) Science Plan

Clouds over the Southern Ocean are poorly represented in present day reanalysis products and global climate model simulations. Errors in top-of-atmosphere (TOA) broadband radiative fluxes in this region are among the largest globally, with large implications for modeling both regional and global scale climate responses (e.g., Trenberth and Fasullo 2010, Ceppi et al. 2012). Recent analyses of model simulations suggest that model radiative errors in the Southern Ocean are due to a lack of low-level postfrontal clouds (including clouds well behind the front) and perhaps a lack of supercooled liquid water that contribute most to the model biases (Bodas-Salcedo et al. 2013, Huang et al. 2014). These assessments of model performance, as well as our knowledge of cloud and aerosol properties over the Southern Ocean, rely heavily on satellite data sets. Satellite data sets are incomplete in that the observations are not continuous (i.e., they are acquired only when the satellite passes nearby), generally do not sample the diurnal cycle, and view primarily the tops of cloud systems (especially for the passive instruments). This is especially problematic for retrievals of aerosol, low-cloud properties, and layers of supercooled water embedded within (rather than at the top of) clouds, as well as estimates of surface shortwave and longwave fluxes based on these properties.

Marchand, RT↗
DOE OSTI2013

Off-gas adsorption model and simulation - OSPREY

A capability of accurately simulating the dynamic behavior of advanced fuel cycle separation processes is expected to provide substantial cost savings and many technical benefits. To support this capability, a modeling effort focused on the off-gas treatment system of a used nuclear fuel recycling facility is in progress. The off-gas separation consists of a series of scrubbers and adsorption beds to capture constituents of interest. Dynamic models are being developed to simulate each unit operation involved so each unit operation can be used as a stand-alone model and in series with multiple others. Currently, an adsorption model has been developed within Multi-physics Object Oriented Simulation Environment (MOOSE) developed at the Idaho National Laboratory (INL). Off-gas Separation and Recovery (OSPREY) models the adsorption of offgas constituents for dispersed plug flow in a packed bed under non-isothermal and non-isobaric conditions. Inputs to the model include gas composition, sorbent and column properties, equilibrium and kinetic data, and inlet conditions. The simulation outputs component concentrations along the column length as a function of time from which breakthrough data can be obtained. The breakthrough data can be used to determine bed capacity, which in turn can be used to size columns. In addition to concentration data, the model predicts temperature along the column length as a function of time and pressure drop along the column length. A description of the OSPREY model, results from krypton adsorption modeling and plans for modeling the behavior of iodine, xenon, and tritium will be discussed. (author)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. Lastly, HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

58 GEOSCIENCES↗
DOE OSTI2021

Performance Testing and Assessment of Protection Scheme Using Real-Time Hardware-in-the-Loop and IEC 61850 Standard

The main challenge of the microgrid is to design a suitable protection scheme due to the complexity of the architecture of the microgrid. The importance of the proposed protection technique is threefold. First, it presents a co-simulation platform to integrate between a simulated model on power system computer aided design (PSCAD)/real time digital simulator computer aided design (RSCAD) software’s and physical devices schweitzer engineering laboratories (SEL) 421-7 relays to protect the microgrid that includes different resources connected based on inverter interface. Second, it presents a comprehensive hardware/software setup to test the protective relays in a closed loop system and shows how to configure the protective relay’s International Electrotechnical Commission 61850 communications. Third, IEEE 1588 standard is used to provide sub nanoseconds latency between the simulated model that emulated on real time digital simulator (RTDS) and the external devices. Also, the measurement signals are synchronized between RTDS and the external devices using giga-transceiver synchronization card (GTSYNC) interface card and SEL-2488 satellite-synchronized network clock. The results showed that the co-simulation infrastructure introduces a highly dependable design, analysis, and testing environment for cyber and physical data flow in the system. Besides that, the voltages at ac/dc sides and frequency at fault condition were maintained due to the energy storage device contributions at different modes of operation.

42 ENGINEERING↗
DOE OSTI2015

Theoretical investigation of the design and performance of a dual energy (kV and MV) radiotherapy imager

Purpose: In modern radiotherapy treatment rooms, megavoltage (MV) portal imaging and kilovoltage (kV) cone-beam CT (CBCT) imaging are performed using various active matrix flat-panel imager (AMFPI) designs. To expand the clinical utility of MV and kV imaging, MV AMFPIs incorporating thick, segmented scintillators and, separately, kV imaging using a beam’s eye view geometry have been investigated by a number of groups. Motivated by these previous studies, it is of interest to explore to what extent it is possible to preserve the benefits of kV and MV imaging using a single AMFPI design, given the considerably different x ray energy spectra used for kV and MV imaging. In this paper, considerations for the design of such a dual energy imager are explored through examination of the performance of a variety of hypothetical AMFPIs based on x ray converters employing segmented scintillators. Methods: Contrast, noise, and contrast-to-noise ratio performances were characterized through simulation modeling of CBCT imaging, while modulation transfer function, Swank factor, and signal performance were characterized through simulation modeling of planar imaging. The simulations were based on a previously reported hybrid modeling technique (accounting for both radiation and optical effects), augmented through modeling of electronic additive noise. All designs employed BGO scintillator material with thicknesses ranging from 0.25 to 4 cm and element-to-element pitches ranging from 0.508 to 1.016 mm. A series of studies were performed under both kV and MV imaging conditions to determine the most advantageous imager configuration (involving front or rear x ray illumination and use of a mirror or black reflector), converter design (pitch and thickness), and operating mode (pitch-binning combination). Results: Under the assumptions of the present study, the most advantageous imager design was found to employ rear illumination of the converter in combination with a black reflector, incorporate a BGO converter with a 0.508 mm pitch and a 2 cm thickness, and operate at full resolution for kV imaging and 2 × 2 binning mode for MV imaging. Such a dual energy imager design should provide soft tissue visualization at low, clinically practical doses under MV conditions, while helping to preserve the high spatial resolution and high contrast offered by kV imaging. Conclusions: The authors’ theoretical investigation suggests that a dual energy imager capable of largely preserving the desirable characteristics of both kV and MV imaging is feasible. Such an imager, when coupled to a dual energy radiation source, could facilitate simplification of current treatment room imaging systems (as well as their associated quality assurance), and facilitate more precise integration of kV and MV imaging information by virtue of reduced geometric uncertainties.

07 ISOTOPE AND RADIATION SOURCES↗