SEARCH · Engineering Papers

Results for “reference models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,081 records · Page 60

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗

DOE OSTI2025

H2A-Lite (Hydrogen Analysis Lite Production Model) [SWR-24-69]

Within H2A-Lite, users can provide a minimal number of inputs—such as hydrogen production technology of choice —to produce estimates about characteristic scale, capital, and operations. Price projections for energy and feedstock are based on the Energy Information Administration's Annual Energy Outlook 2022, AEO2022 Reference case. The model additionally allows users to override technology default values to adapt to specific technology scales or regional energy prices. As output, H2A-Lite provides cost breakdown from rigorous financial analysis as well as greenhouse gas and criteria pollutant emissions characteristics.

Penev, Michael [National Laboratory of the Rockies↗

DOE OSTI2021

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.

42 ENGINEERING↗

DOE OSTI2021

Incorporating Valve Switching Losses Into a Static Optimal Control Algorithm for the Hybrid Hydraulic-Electric Architecture (HHEA)

The Hybrid Hydraulic-Electric Architecture (HHEA) has, in recent years, been proposed as an energy efficient alternative to conventional load-sensing architectures in mobile machines such as excavators and wheel-loaders. HHEA leverages the advantages of hydraulic power and electric power to eliminate throttling valves while also improving the energy and control performance of the system. The architecture utilizes a set of common pressure rails to provide a majority of power and and a small electric motor driven pump to modulate this power to meet the exact demand. Previous work has developed a computationally efficient Lagrange Multiplier approach for determining the optimal pressure rail selections that minimizes the energy losses in the system. The static model used considers only the energy use for each pressure rail selection but not the losses associated with the valves during the transition. This paper presents an approach to include the switching losses in the model and in the optimization procedure. To capture the switching losses, switching events between different rails and at various input and output flow rates were simulated with consideration of valve spool dynamics. A parameterized model that summarizes the losses is then obtained, allowing switching losses to be added to the previous energy analysis. The performance of the switching loss model was compared with reference data obtained from a high-fidelity simulation model. To incorporate the switching losses into optimal control algorithm, an efficient dynamic programming approach that prevents frequent switching is adopted in place of the Lagrange multiplier approach. The overall effect of switching losses on the energy consumption and optimal control decisions is presented. In general, switching losses contribute to about 9–10% of input energy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1976

An evaluation of recent quantitative magnetospheric magnetic field models

Magnetospheric field models involving dipole tilt effects are discussed, with particular reference to defined magnetopause models and boundary surface models. The models are compared with observations and with each other whenever possible. It is shown that models containing only contributions from magnetopause and tail current systems are capable of reproducing the observed quiet time field just in a qualitative way. The best quantitative agreement between models and observations take place when currents distributed in the inner magnetosphere are added to the magnetopause and tail current systems. One region in which all the models fall short is the region around the polar cusp. Obtaining physically reasonable gradients should have high priority in the development of future models.

Walker, R. J.↗

NASA NTRS1976

The NASTRAN User's Manual Level 16.0 and Supplement

The user's manual is restricted to those items related to the use of NASTRAN that are independent of the computing system being used. The features of NASTRAN described include: (1) procedures for defining and loading a structural model and a functional reference for every card that is used for structural modeling; (2) the NASTRAN data deck, including the details for each of the data cards; (3) the NASTRAN control cards that are associated with the use of the program; (4) rigid format procedures, along with specific instructions for the use of each rigid format: (5) procedures for using instructions for the use of each rigid format; (5) procedures for using the NASTRAN plotting capability; (6) procedures governing the creation of DMAP programs; and (7) the NASTRAN diagnostic messages. The NASTRAN dictionary of mnemonics, acronyms, phrases, and other commonly used NASTRAN terms is included along with a limited number of sample problems.

Source record↗

DOE OSTI2020

Influence of sun zenith angle on canopy clumping and the resulting impacts on photosynthesis

Addressing the impact of vegetation architecture on shortwave radiation transfer in land surface models is important for accurate weather forecasting, carbon budget estimates, and climate predictions. This paper investigates to what extent it is possible to retrieve structural parameters of two different parameterization schemes from direct transmittance derived from digital hemispherical photography and 3D radiative transfer modeling for two study sites with different vegetation canopy architectures. Neglecting the representation of 3D canopy structure in radiative transfer schemes leads to significant errors in shortwave radiation partitioning (up to 3.5 times more direct transmittance in the 3D model). Structural parameters, referred to as whole canopy ‘clumping indices’, were obtained in order to evaluate the impact of angular variation in clumping on shortwave radiation transfer. Impacts on photosynthesis were evaluated at site level with the UKESM land surface model, JULES. A comparison between flux tower derived and modeled photosynthesis indicates that considering zenith angular variations of structural parameters in the radiative transfer scheme of the UKESM land surface model significantly improves photosynthesis prediction in light limited ecosystems (from RMSE = 2.91 μmol CO 2 .m -2 .s -1 to RMSE = 1.51 μmol CO 2 .m -2 .s -1 , 48% smaller), typically with enhanced photosynthesis from bottom layers.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1980

Model error estimation for large space systems

In-flight estimation of large structure model errors may have to be carried out in order to detect inevitable deficiencies in large structure controller/estimator models. These error estimates can most efficiently be obtained by the minimization of a quadratic functional of the model errors and on the subsequent analysis of the resulting optimal model error estimates. An integral operator approach to estimation leads to a geometrical interpretation of the model error estimation process. One of the significant insights gained with this interpretation is that the actual but unknown model errors can be decomposed as the sum of two distinct components that are orthogonal in some sense. One of these components is a so-called minimal error vector that retains many of the significant dynamics of the actual errors. The basic ideas in the model error estimation approach are first set forth with a two-dimensional analogy that has most of the essential features of the general estimation problem. The generalized results are then established and their application to a reference large structure model illustrated.

Rodriguez, G.↗

NASA NTRS1982

A reference solar spectral irradiance for use in atmospheric modeling

The present state of knowledge concerning the absolute magnitude and temporal variability of the solar spectral irradiance is outlined with emphasis on wavelengths relevant to the mesosphere and stratosphere. Reference spectra for the wavelength region 175 to 850 nm are presented including estimates for solar maximum and solar minimum conditions. Values for the Lyman alpha emission are given separately.

Source record↗

DOE OSTI2021

Comparison of Radiation Models for a Turbulent Piloted Methane/Air Jet Flame: A Frozen-Field Study

Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations.

David, Chloe↗

NASA NTRS1987

Space Station/Shuttle Orbiter dynamics during docking

Mathematical models of a reference space station configuration (Power Tower) and a Space Shuttle Orbiter are developed and used to study the dynamic behavior of the Space Station/Orbiter system just prior to and subsequent to an impulsive docking of the two spacecraft. The physical model of the space station is a collection of rigid and flexible bodies. The orbiter is modeled as a rigid body. An algorithm developed for use in digitally simulating the dynamics of the system is described and results of its application are presented.

Fitz-Coy, N. G.↗

DOE OSTI2020

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2021

Evaluation of 18 Satellite- and Model-Based Soil Moisture Products Using in Situ Measurements From 826 Sensors

Information about the spatiotemporal variability of soil moisture is critical for many purposes, including monitoring of hydrologic extremes, irrigation scheduling, and prediction of agricultural yields. We evaluated the temporal dynamics of 18 state-of-the-art (quasi-)global near-surface soil moisture products, including six based on satellite retrievals, six based on models without satellite data assimilation (referred to hereafter as “open-loop” models), and six based on models that assimilate satellite soil moisture or brightness temperature data. Seven of the products are introduced for the first time in this study: one multi-sensor merged satellite product called MeMo (Merged soil Moisture) and six estimates from the HBV (Hydrologiska Byråns Vattenbalansavdelning) model with three precipitation inputs (ERA5, IMERG, and MSWEP) with and without assimilation of SMAPL3E satellite retrievals, respectively. As reference, we used in situ soil moisture measurements between 2015 and 2019 at 5 cm depth from 826 sensors, located primarily in the USA and Europe. The 3-hourly Pearson correlation (R) was chosen as the primary performance metric. We found that application of the Soil Wetness Index (SWI) smoothing filter resulted in improved performance for all satellite products. The best-to-worst performance ranking of the four single-sensor satellite products was SMAPL3E SWI , SMOS SWI , AMSR2 SWI , and ASCAT SWI , with the L-band-based SMAPL3E SWI (median R of 0.72) outperforming the others at 50 % of the sites. Among the two multi-sensor satellite products (MeMo and ESA-CCI SWI ), MeMo performed better on average (median R of 0.72 versus 0.67), probably due to the inclusion of SMAPL3E SWI . The best-to-worst performance ranking of the six open-loop models was HBV-MSWEP, HBV-ERA5, ERA5-Land, HBV-IMERG, VIC-PGF, and GLDAS-Noah. This ranking largely reflects the quality of the precipitation forcing. HBV-MSWEP (median R of 0.78) performed best not just among the open-loop models but among all products. The calibration of HBV improved the median R by +0.12 on average compared to random parameters, highlighting the importance of model calibration. The best-to-worst performance ranking of the six models with satellite data assimilation was HBV-MSWEP+SMAPL3E, HBV-ERA5+SMAPL3E, GLEAM, SMAPL4, HBV-IMERG+SMAPL3E, and ERA5. The assimilation of SMAPL3E retrievals into HBV-IMERG improved the median R by +0.06, suggesting that data assimilation yields significant benefits at the global scale.

Hylke E. Beck↗

NASA NTRS1977

User manual and programmer reference manual for the ATS-6 navigation model AOIPS and McIDAS versions, part 2

Development of a navigation system for a given satellite is reported. An algorithm for converting a satellite picture element location to earth location and vice versa was defined as well as a procedure for measuring the set of constants needed by the algorithm. A user manual briefly describing the current version of the navigation model and how to use the computer programs developed for it is presented.

Chatters, G. C.↗

NASA NTRS1980

The development of a self-streamlining flexible walled transonic test section

This design eliminates the uncertainties in data from conventional transonic test sections. Sidewalls are rigid, and the flexible floor and ceiling are positioned by motorized jacks controlled by on-line computer to minimize run times. The tunnel-computer combination is self-streamlining without reference to the model. Data is taken from the model only when the walls are good streamlines, and is corrected for the small, known but inevitable residual interferences. Two-dimensional validation testing in the Mach range up to about 0.85 where the walls are just supercritical shows good agreement with reference data using a height:chord ratio of 1.5. Techniques are under development to extend Mach number above 1. This work has demonstrated the feasibility of almost eliminating wall interferences, improving flow quality, and reducing power requirements or increasing Reynolds number. Extensions to three-dimensional testing are outlined.

Goodyer, M. J.↗

NASA NTRS2001

The Physics of Accretion in X-Ray Binaries

This project consists of several related investigations directed to the study of mass transfer processes in X-ray binaries. Models developed over several years incorporating highly detailed physics will be tested on a balanced mix of existing data and planned observations with both ground and space-based observatories. The extended time coverage of the observations and the existence of simultaneous X-ray, ultraviolet, and optical observations will be particularly beneficial for studying the accretion flows. These investigations, which take as detailed a look at the accretion process in X-ray binaries as is now possible, test current models to their limits, and force us to extend them. We now have the ability to do simultaneous ultraviolet/X-ray/optical spectroscopy with HST, Chandra, XMM, and ground-based observatories. The rich spectroscopy that these observations give us must be interpreted principally by reference to detailed models, the development of which is already well underway; tests of these essential interpretive tools are an important product of the proposed investigations.

Vrtilek, S.↗

NASA NTRS2003

The Physics of Accretion in X-Ray Binaries

This project consists of several related investigations directed to the study of mass transfer processes in X-ray binaries. Models developed over several years be tested on a balanced mix of existing data and planned observations with both ground and space-based observatories. The extended time coverage of the observations and the existence of simultaneous X-ray, ultraviolet, and optical observations will be particularly beneficial for studying the accretion flows. These investigations, which take as detailed a look at the accretion process in X-ray binaries as is now possible, test current models to their limits, and force us to extend them. We now have the ability to do simultaneous ultraviolet/X-ray/optical spectroscopy with HST, Chandra, XMM, and ground-based observatories. The rich spectroscopy that these observations give us must be interpreted principally by reference to detailed models, the development of which is already well underway; tests of these essential interpretive tools are an important product of the proposed investigations.

Oliversen, Ronald↗