Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,333 records · Page 74

Visualization software for a global three-dimensional upper atmosphere model

A software environment has been developed to assist in analysis and visualization of numerical outputs from NCAR's thermospheric general circulation models. The codes produce vector graphics, raster imagery, and time dependent animation of model results, with provision for comparison with empirical models and relevant data from satellite and ground based instruments. Separate programs are executed on the NCAR CRAY YMP8/864 and HAO divisional Sun workstations, providing both large volume highspeed batch and custom interactive computing environments. The visualization system assembles a suite of existing software packages including NCAR Graphics, IDL, netCDF, and the Xt intrinsics library with the Athena widget set.

Foster, Benjamin T.↗

CEDAR Electrodynamics Thermosphere Ionosphere (ETI) Challenge for Systematic Assessment of Ionosphere/Thermosphere Models: NmF2, hmF2, and Vertical Drift Using Ground-Based Observations

Objective quantification of model performance based on metrics helps us evaluate the current state of space physics modeling capability, address differences among various modeling approaches, and track model improvements over time. The Coupling, Energetics, and Dynamics of Atmospheric Regions (CEDAR) Electrodynamics Thermosphere Ionosphere (ETI) Challenge was initiated in 2009 to assess accuracy of various ionosphere/thermosphere models in reproducing ionosphere and thermosphere parameters. A total of nine events and five physical parameters were selected to compare between model outputs and observations. The nine events included two strong and one moderate geomagnetic storm events from GEM Challenge events and three moderate storms and three quiet periods from the first half of the International Polar Year (IPY) campaign, which lasted for 2 years, from March 2007 to March 2009. The five physical parameters selected were NmF2 and hmF2 from ISRs and LEO satellites such as CHAMP and COSMIC, vertical drifts at Jicamarca, and electron and neutral densities along the track of the CHAMP satellite. For this study, four different metrics and up to 10 models were used. In this paper, we focus on preliminary results of the study using ground-based measurements, which include NmF2 and hmF2 from Incoherent Scatter Radars (ISRs), and vertical drifts at Jicamarca. The results show that the model performance strongly depends on the type of metrics used, and thus no model is ranked top for all used metrics. The analysis further indicates that performance of the model also varies with latitude and geomagnetic activity level.

Shim, J. S.↗

A Realization of Bias Correction Method in the GMAO Coupled System

Over the past several decades, a tremendous effort has been made to improve model performance in the simulation of the climate system. The cold or warm sea surface temperature (SST) bias in the tropics is still a problem common to most coupled ocean atmosphere general circulation models (CGCMs). The precipitation biases in CGCMs are also accompanied by SST and surface wind biases. The deficiencies and biases over the equatorial oceans through their influence on the Walker circulation likely contribute the precipitation biases over land surfaces. In this study, we introduce an approach in the CGCM modeling to correct model biases. This approach utilizes the history of the model's short-term forecasting errors and their seasonal dependence to modify model's tendency term and to minimize its climate drift. The study shows that such an approach removes most of model climate biases. A number of other aspects of the model simulation (e.g. extratropical transient activities) are also improved considerably due to the imposed pre-processed initial 3-hour model drift corrections. Because many regional biases in the GEOS-5 CGCM are common amongst other current models, our approaches and findings are applicable to these other models as well.

Chang, Yehui↗

International Space Station Bipropellant Plume Contamination Model Update for Short Thruster Pulse Widths

The International Space Station (ISS) Bipropellant Plume Contamination Model has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS but was not intended to be used for very short thruster pulse widths. This paper presents an updated model that ensures full start-up and shut-down phases are modeled for each thruster firing, when the majority of liquid phase contaminant mass is released. The updated ISS Bipropellant Plume Contamination Model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations and provides a way to take advantage of thruster start-up and shut-down data performance data gathered during thruster test programs, if available. The International Space Station (ISS) Bipropellant Plume Contamination Model developed by Boeing Space Environments is a semi-empirical model anchored in flight experiment data and has been a vital tool for characterizing the thruster plume-induced contamination environment at the ISS.[1] The current model utilizes flight experiment data from the Plume Impingement Contamination (PIC) and Shuttle Plume Impingement Flight Experiment (SPIFEX) studies, which include Orbiter 3870 N Primary Reaction Control System (PRCS) and Russian 130 N thrusters operating in pulse mode with 80-100 ms.[2] As the next generation of crew and cargo visiting vehicles are developed and arrive at ISS, minimum pulse widths of vehicle thrusters used for ISS proximity operations have decreased significantly below 80 ms. An update to the ISS Bipropellant Contamination Model is needed to prevent potential under-prediction of thruster plume-induced contamination for these very short pulse widths. Contamination due to thruster plumes occurs in the liquid phase (i.e., unburned or partially burned propellant in the plume). Liquid phase releases primarily occur during thruster start-up and shut-down phases, with the steady state phase contributing a small amount to the total contaminant mass released. The current ISS Bipropellant Plume Contamination Model includes functional dependencies on thruster parameters (thrust, mass flow rate, specific impulse (Isp)), distance to receiver surface, and angle off plume centerline. For modeling purposes, contaminant mass released scales linearly with pulse width. This is a conservative approach for pulse widths greater than 80 ms but could under-predict contamination for very short pulse widths (i.e., with little or no steady state phase). This paper presents a model update to add a functional dependency on commanded pulse width to ensure the portion of the pulse spent in start-up and shut-down phases is appropriately modeled, when the majority of liquid phase contaminant mass is released. The updated model prevents potential under-prediction of thruster plume-induced contamination due to visiting vehicle proximity operations while providing a way to take advantage of thruster start-up and shut-down data gathered during test programs, if available. Developing the plume model for a specific thruster using the ISS Bipropellant Plume Contamination Model must be done in consideration of all available thruster performance data. Example thruster data will be used to illustrate this point and options for implementing the model update with existing ISS Bipropellant Plume Contamination Model code will be discussed.

plume contamination↗

Validation and Sensitivity Analysis of a 1-D Lake Model Across Global Lakes

Lakes have important influence on weather and climate from local to global scales. However, their prediction using numerical models is notoriously difficult because lakes are highly heterogeneous across the globe, but observations are sparse. In this study, we assessed the performance of a 1-D lake model in simulating the thermal structures of 58 lakes with diverse morphometric and geographic characteristics by following the phase 2a local lake protocol of the Inter-sectoral Impact Model Intercomparison Project (ISIMIP2a). After calibration, the root-mean-square errors (RMSE) were below 2 °C for 70% and 75% of the lakes for epilimnion and full-profile temperature simulations, with an average of 1.71 °C and 1.43 °C, respectively. The model performance mainly depended on lake shape rather than location, supporting the possibility of grouping model parameters by lake shape for global applications. Furthermore, through machine-learning based parameter sensitivity tests, we identified turbulent heat fluxes, wind-driven mixing and water transparency as the major processes controlling lake thermal and mixing regimes. Snow density was also important for modeling the ice phenology of high-latitude lakes. The relative influence of the key processes and the corresponding parameters mainly depended on lake latitude and depth. Turbulent heat fluxes showed a decreasing importance in affecting epilimnion temperature with increasing latitude. Wind-driven mixing was less influential to lake stratification for deeper lakes while the impact of light extinction, on the contrary, showed a positive correlation with depth. Our findings may guide improvements in 1-D lake model parameterizations to achieve higher fidelity in simulating global lake thermal dynamics.

54 ENVIRONMENTAL SCIENCES↗

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

A study of plasmaspheric density distributions for diffusive equilibrium conditions

The plasmaspheric density distribution has been modeled for a range of solar cycle, seasonal and diurnal conditions with a magnetic flux tube dependent diffusive equilibrium model by using experimentally determined values of ionospheric parameters at 675 km as boundary conditions. Data is presented in terms of plasmaspheric H(+) and He(+) density contours, total flux tube content, and equatorial plasma density for a range of L-values from 1.15 to 3.0. The variation of equatorial density with L-value shows good agreement with the 1/L exp 4 dependence observed experimentally. The results show that the model predicts larger solar cycle and diurnal variation in equatorial plasma density than observed using whistler techniques. However, the whistler method requires a model to deduce the equatorial density and is therefore open to interpretation. Seasonal variations are rather artificial since in this general model no attempt has been made to match equatorial densities for flux tubes emanating from the winter and summer hemispheres.

Li, W.↗

Bridging microscale to macroscale mechanical property measurements of FeCrAl alloys by crystal plasticity modeling

FeCrAl alloys are candidates for accident tolerant fuel cladding of light water reactors. In this work, a microstructure- and temperature-dependent crystal plasticity model is employed to bridge microscale to macroscale mechanical property measurements of FeCrAl alloys. With the visco-plastic self-consistent (VPSC) polycrystal plasticity framework, a mechanism-based single crystal plasticity (MSCP) model adopts the Arrhenius type rate equation to describe the dependence of the critical resolved shear stress for dislocation slips on their temperature-dependent intrinsic frictional resistance and the microstructure-dependent irradiation hardening. The intrinsic frictional resistance associated with {110}<111> and {112}<111> slip systems were measured by in-situ micromechanical testing on unirradiated/irradiated samples at 25-500 °C. The irradiation hardening is estimated by the Bacon-Kocks-Scattergood (BKS) model with density and size of radiation-induced defects measured from microstructural characterization. Several features associated with thermo-mechanical behavior of unirradiated/irradiated polycrystalline FeCrAl alloys are captured. High density of deformation-induced dislocations and radiation-induced defects results in obvious hardening at room temperature, which is weakened at high temperature, and facilitates damage evolution during deformation. Moreover, both high temperature and radiation-induced defects, which facilitate dislocation multiplication, trigger large hardening rate. Finally, the proposed method together with application of accelerator-based ion irradiation technique is a surrogate approach to simulate neutron damage, improving the efficiency associated with evaluation of mechanical properties of FeCrAl alloys exposed to temperature, stress and radiation conditions.

36 MATERIALS SCIENCE↗

Time dependent plasma injection by Io

A two parameter model of time-dependent, flux-tube interchange diffusion is fit to the Voyager 1 plasma data obtained in the Io plasma disk. The interpretation of the parameters required to achieve the fit is that plasma injection increased suddenly and substantially (by more than an order of magnitude) at some time prior to the arrival of Voyager 1 (between 1 and 100 days prior). The injection rate was about 2 x 10 to the (29 plus or minus 1) power ion/sec. At this rate, the centrifugally driven interchange instability dominated outward diffusion, causing the outward diffusion rate to be about a factor of 50 greater than the inward diffusion rate. The material diffusing inward had time to cool by radiation, possibly accounting for the observed temperature drop inside the orbit of Io.

Richardson, J. D.↗

A semikinetic model for early stage plasmasphere refilling. I - Effects of Coulomb collisions

A collisionless, time-dependent, kinetic plasma model is applied to the problem of baseline plasmasphere refilling of an initially depleted flux tube, without regard for the effects of wave-particle interactions. Refilling calculations for various flux tubes and for different ionospheric plasma fluxes and temperatures are performed. In each case considered, the same set of events occurs. Initially, two polar wind outflows develop from each hemisphere and set up counterstreaming beams. With time the vacant phase space region between these beams fills, primarily because of collision-induced particle diffusion but also because of lowering ambipolar potential drops from the increasing density in the plasmasphere. In contrast to all previous hydrodynamic approaches, no formation of shocks was found. The plasma first evolves an isotropic, nearly Maxwellian velocity distribution in a region that starts near the ionosphere and moves outward toward the equator. For reasonable topside ionospheric temperatures and fluxes, the thermal plasma all along an L shell is found to become nearly isotropic in 6 to 30 hr, consistent with the observations of Horwitz et al. (1984).

Wilson, G. R.↗

Designing a repository in domal salt: The influence of design variants in different modelling environments

To understand the long-term environmental impact of disposing radioactive waste of in a deep geological repository and to optimise its design, performance assessments are used. In this study, four teams (COVRA, GRS, Quintessa, and DOE) modified the previously developed generic repository of DECOVALEX task F2 to identify commonalities and differences between the teams for specific changes in repository design. The teams tested six design modifications: (1) Replacing concrete abutments with run-of-mine salt; (2) Replacing the salt seal with a concrete abutment and using run-of-mine salt instead for the two concrete abutments in each drift seal; (3) Halving the size of the infrastructure area; (4) Using run-of-mine salt instead of gravel for backfilling the infrastructure area; (5) Disposal of spent nuclear fuel without the POLLUX-10 containers (6); Lower initial saturation of the spent nuclear fuel and vitrified waste disposal drifts. Despite different modelling strategies used, models agreed that a smaller infrastructure area has a limited effect on radionuclide transport. Responses to the absence of the two concrete abutments in each seal, the use of single large concrete abutments (200 m each), or the use of run-of-mine salt in the infrastructure area differ between teams due to differing modelling assumptions. Based on these results, the estimated efficacy of containment depends strongly on the model assumptions of each team. More specifically, it appears to depend on the compaction model used and therefore on the backfill material used in different areas of the repository. However, the drift seal appears to be a critical design element in all models, effectively limiting radionuclide transport by hydrologically disconnecting sections of the repository. Additional beneficial design choices include the use of dry salt in disposal drifts to limit radionuclide transport and reducing the infrastructure area costs and minimizing host rock damage.

DECOVALEX↗

Data-driven model validation for neutrino-nucleus cross section measurements

Neutrino-nucleus cross section measurements are needed to improve interaction modeling to meet the precision needs of neutrino experiments in efforts to measure oscillation parameters and search for physics beyond the Standard Model. We review the difficulties associated with modeling neutrino-nucleus interactions that lead to a dependence on event generators in oscillation analyses and cross section measurements alike. We then describe data-driven model validation techniques intended to address this model dependence. The method relies on utilizing various goodness-of-fit tests and the correlations between different observables and channels to probe the model for defects in the phase space relevant for the desired analysis. These techniques shed light on relevant mismodeling, allowing it to be detected before it begins to bias the cross section results. We compare more commonly used model validation methods which directly validate the model against alternative ones to these data-driven techniques and show their efficacy with fake data studies. These studies demonstrate that employing data-driven model validation in cross section measurements represents a reliable strategy to produce robust results that will stimulate the desired improvements to interaction modeling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Data-driven model validation for neutrino-nucleus cross section measurements

Neutrino-nucleus cross section measurements are needed to improve interaction modeling to meet the precision needs of neutrino experiments in efforts to measure oscillation parameters and search for physics beyond the Standard Model. We review the difficulties associated with modeling neutrino-nucleus interactions that lead to a dependence on event generators in oscillation analyses and cross section measurements alike. We then describe data-driven model validation techniques intended to address this model dependence. The method relies on utilizing various goodness-of-fit tests and the correlations between different observables and channels to probe the model for defects in the phase space relevant for the desired analysis. These techniques shed light on relevant mismodeling, allowing it to be detected before it begins to bias the cross section results. We compare more commonly used model validation methods which directly validate the model against alternative ones to these data-driven techniques and show their efficacy with fake data studies. These studies demonstrate that employing data-driven model validation in cross section measurements represents a reliable strategy to produce robust results that will stimulate the desired improvements to interaction modeling.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A framework for feasibility-level validation of high-resolution wave hindcast models

The value of long-term wave hindcasts for investigating wave climates, wave energy resources, and extreme wave conditions has motivated research developing, calibrating and validating wave hindcast models. Past hindcast model validation studies examined the accuracy in modeling bulk wave parameters of overall sea states without considering the dependency of the model's skill within different sea states. In the present study, a framework for wave hindcast model validation is developed by examining the model accuracy for the most frequently occurring sea states, sea states contributing the most energy to total wave power, sea states associated with hurricane events, and those with the largest model error. Here, validations using bulk wave parameters and frequency-directional spectra at these key sea states and extreme wave conditions based on univariate and bivariate-contour methods provide insights to improve model accuracy, identifying the model's strong and weak points, and pathways for improvement, e.g., modeling wave-current interactions and adjusting wind data. This study adds to a growing body of research demonstrating that a carefully calibrated and verified spectral wave hindcast model can be used to estimate key wave energy parameters over a wide range of wave energy climates, as well as their spatial, temporal, frequency, directional, and probabilistic distributions.

54 ENVIRONMENTAL SCIENCES↗

Anomalous transport in high-mobility superconducting SrTiO 3 thin films

The study of subtle effects on transport in semiconductors requires high-quality epitaxial structures with low defect density. Using hybrid molecular beam epitaxy (MBE), SrTiO 3 films with a low-temperature mobility exceeding 42,000 cm 2 V -1 s -1 at a low carrier density of 3 × 10 17 cm -3 were achieved. A sudden and sharp decrease in residual resistivity accompanied by an enhancement in the superconducting transition temperature were observed across the second Lifshitz transition where the third band becomes occupied, revealing dominant intraband scattering. These films further revealed an anomalous behavior in the Hall carrier density as a consequence of the antiferrodistortive (AFD) transition and the temperature dependence of the Hall scattering factor. Using hybrid MBE growth, phenomenological modeling, temperature-dependent transport measurements, and scanning superconducting quantum interference device imaging, we provide critical insights into the important role of inter- versus intraband scattering and of AFD domain walls on normal-state and superconducting properties of SrTiO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantifying the Impact of Vertical Resolution on the Representation of Marine Boundary Layer Physics for Global-Scale Models

Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement. Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

Meteorology & Atmospheric Sciences↗

Hierarchical Conceptual Model Event and Activity Domains for Forecasting State-Sponsored Civil Nuclear Power Activities

In FY20, the Savannah River National Laboratory and the Sanghani Center for Artificial Intelligence and Data Analytics (Virginia Tech) entered a collaboration funded by the Department of Energy’s Office of Defense Nuclear Nonproliferation Research and Development to develop a machine learning based modeling pipeline to extract proliferation events of interest from open data sources. The prototype modeling pipeline that was developed relies on the use of time dependent word embedding models to identify contextual shifts in key words and phrases that act as indicators of events of interest. The FY20- 21 efforts were focused on a narrow topical domain of forecasting “fissile core fabrication” at the Savannah River Site prior to its official announcement in 2018. In FY22, the research team was funded to continue development of the modeling pipeline by applying it to the problem of forecasting new, and/or significant changes to existing, civil nuclear power reactors around the world. In this effort, the development will focus on proving applicability to a broader topical domain and in data environments that may contain more sparse information, relative to the United States. In this report, the worldwide landscape of civil nuclear reactors is outlined, the timeline of interest is defined, and a hierarchical conceptual model is established to identify the activity domains of interest and the event domains of interest. This preliminary effort will guide the curation of a glossary of key terms for data acquisition, as well as the downstream modeling efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗