SEARCH · Engineering Papers

Results for “Reduced Order 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,639 records · Page 91

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Simulation of hardened cement degradation and estimation of uncertainty in predicted failure times with peridynamics

Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, here we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Numerical integrators for Lagrangian oceanography

Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Quantifying the benefits of improved satellite remote-sensing observations for inverse modeling of NOx and NMVOC emissions

This study aims to demonstrate the benefits of using novel high spatiotemporal retrieval products from newer satellites for top-down emission estimates of nitrogen oxides (NO x ) and non-methane volatile organic compounds (NMVOCs) for the summer of 2019 over the contiguous United States. Recent satellite retrievals have not only advanced spatiotemporal resolution but also greatly reduced error and uncertainty due to reduced noise in the retrievals compared to spaceborne sensors launched in the past. We applied inverse modeling techniques using tropospheric nitrogen dioxide (NO 2 ) and formaldehyde (H-CHO) column retrieval products from the Ozone Monitoring Instrument (OMI) and TROPOspheric Monitoring Instrument (TROPOMI) in conjunction with the Weather Research Forecast and Community Multiscale Air Quality Modeling system (WRF-CMAQ). In order to provide a better representation of background chemical composition and avoid misalignment of emission adjustment, we applied monthly scaling factors for ozone (O 3 ) and CO boundary concentrations in addition to the inclusion of lightning and aviation emissions. Satellite-constrained NO x and NMVOCs posterior emissions showed a mitigated discrepancy between observed and modeled columns. The improvement in the model performance was greater when using TROPOMI, primarily benefiting from reduced errors/biases of the satellite retrievals that enabled us to explore corresponding changes in O 3 concentrations and production sensitivity regimes using the ratio of H-CHO and NO 2 .

remote sensing↗

DOE OSTI2023

Elucidating Design Rules toward Enhanced Solid-State Charge Transport in Oligoether-Functionalized Dioxythiophene-Based Alternating Copolymers

This study investigates the solid-state charge transport properties of the oxidized forms of dioxythiophene-based alternating copolymers consisting of an oligoether-functionalized 3,4-propylenedioxythiophene (ProDOT) copolymerized with different aryl groups, dimethyl ProDOT (DMP), 3,4-ethylenedioxythiophene (EDOT), and 3,4-phenylenedioxythiophene (PheDOT), respectively, to yield copolymers P(OE3)-D, P(OE3)-E, and P(OE3)-Ph. At a dopant concentration of 5 mM FeTos 3 , the electrical conductivities of these copolymers vary significantly (ranging between 9 and 195 S cm –1 ) with the EDOT copolymer, P(OE3)-E, achieving the highest electrical conductivity. UV–vis–NIR and X-ray spectroscopies show differences in both susceptibility to oxidative doping and extent of oxidation for the P(OE3) series, with P(OE3)-E being the most doped. Wide-angle X-ray scattering measurements indicate that P(OE3)-E generally demonstrates the lowest paracrystallinity values in the series, as well as relatively small π–π stacking distances. The significant (i.e., order of magnitude) increase in electrical conductivity of doped P(OE3)-E films versus doped P(OE3)-D or P(OE3)-Ph films can therefore be attributed to P(OE3)-E exhibiting both the highest carrier ratios in the P(OE3) series, along with good π–π overlap and local ordering (low paracrystallinity values). Furthermore, these trends in the extent of doping and paracrystallinity are consistent with the reduced Fermi energy level and transport function prefactor parameters calculated using the semilocalized transport (SLoT) model. Observed differences in carrier ratios at the transport edge (c t ) and reduced Fermi energies [η(c)] suggest a broader electronic band (better overlap and more delocalization) for the EDOT-incorporating P(OE3)-E polymer relative to P(OE3)-D and P(OE3)-Ph. Ultimately, we rationalize improvements in electrical conductivity due to microstructural and doping enhancements caused by EDOT incorporation, a structure–property relationship worth considering in the future design of highly electrically conductive systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Taylor wave solution for a general equation of state

This document describes a solution procedure for calculating the Taylor wave behind an unsupported Chapman–Jouguet (CJ) detonation in planar, cylindrical, and spherical geometries given a general equation of state. The resulting semi-analytic solution can be utilized to examine new equation of state models for detonation products and during the verification of hydrodynamic codes. The governing partial differential equations are reduced to ordinary differential equations in both characteristic and self-similar forms. The first-order systems corresponding to each geometry are amenable to solution numerically using commonly available methods. A difficulty arises at the CJ point in radial coordinates where the similarity equations become singular. Two separate strategies are proposed to integrate the first-order system. The first one uses an asymptotic approximation near the CJ point that can be used to perturb the boundary conditions. The second one applies a change of variables which removes the singularity at the expense of an additional equation to be integrated. A test problem is provided for the Davis products equation of state to illustrate the qualitative features of the Taylor wave in each geometric configuration and compared with a Lagrangian hydrodynamics research code. A Python code listing gives an implementation using the SciPy library to assists users in generating the results.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

A fast, reliable algorithm for computing frequency responses of state space models

Computation of frequency responses for large order systems described by time invariant state space systems often provides a bottleneck in control system analysis. It is shown that banding the A-matrix in the state space model can effectively reduce the computation time for such systems while maintaining reliability in the results produced.

Wette, Matt↗

NASA NTRS1993

Optimal fixed-finite-dimensional compensator for Burgers' equation with unbounded input/output operators

The problem of using reduced order dynamic compensators to control a class of nonlinear parabolic distributed parameter systems was considered. Concentration was on a system with unbounded input and output operators governed by Burgers' equation. A linearized model was used to compute low-order-finite-dimensional control laws by minimizing certain energy functionals. Then these laws were applied to the nonlinear model. Standard approaches to this problem employ model/controller reduction techniques in conjunction with linear quadratic Gaussian (LQG) theory. The approach used is based on the finite dimensional Bernstein/Hyland optimal projection theory which yields a fixed-finite-order controller.

Burns, John A.↗

NASA NTRS2009

Bayes Analysis and Reliability Implications of Stress-Rupture Testing a Kevlar/Epoxy COPV Using Temperature and Pressure Acceleration

Composite Overwrapped Pressure Vessels (COPVs) that have survived a long service time under pressure generally must be recertified before service is extended. Flight certification is dependent on the reliability analysis to quantify the risk of stress rupture failure in existing flight vessels. Full certification of this reliability model would require a statistically significant number of lifetime tests to be performed and is impractical given the cost and limited flight hardware for certification testing purposes. One approach to confirm the reliability model is to perform a stress rupture test on a flight COPV. Currently, testing of such a Kevlar49 (Dupont)/epoxy COPV is nearing completion. The present paper focuses on a Bayesian statistical approach to analyze the possible failure time results of this test and to assess the implications in choosing between possible model parameter values that in the past have had significant uncertainty. The key uncertain parameters in this case are the actual fiber stress ratio at operating pressure, and the Weibull shape parameter for lifetime; the former has been uncertain due to ambiguities in interpreting the original and a duplicate burst test. The latter has been uncertain due to major differences between COPVs in the database and the actual COPVs in service. Any information obtained that clarifies and eliminates uncertainty in these parameters will have a major effect on the predicted reliability of the service COPVs going forward. The key result is that the longer the vessel survives, the more likely the more optimistic stress ratio model is correct. At the time of writing, the resulting effect on predicted future reliability is dramatic, increasing it by about one "nine," that is, reducing the predicted probability of failure by an order of magnitude. However, testing one vessel does not change the uncertainty on the Weibull shape parameter for lifetime since testing several vessels would be necessary.

Phoenix, S. Leigh↗

NASA NTRS2011

Base Flow Model Validation

A method was developed of obtaining propulsive base flow data in both hot and cold jet environments, at Mach numbers and altitude of relevance to NASA launcher designs. The base flow data was used to perform computational fluid dynamics (CFD) turbulence model assessments of base flow predictive capabilities in order to provide increased confidence in base thermal and pressure load predictions obtained from computational modeling efforts. Predictive CFD analyses were used in the design of the experiments, available propulsive models were used to reduce program costs and increase success, and a wind tunnel facility was used. The data obtained allowed assessment of CFD/turbulence models in a complex flow environment, working within a building-block procedure to validation, where cold, non-reacting test data was first used for validation, followed by more complex reacting base flow validation.

Sinha, Neeraj↗

DOE OSTI2022

Evaluation of Physical Security Risk for Potential Implementation of FLEX using Dynamic Simulation Methods

The requirements for United States nuclear power plants to maintain a large onsite physical security force contribute to their large operational costs. The cost of maintaining the current physical security posture is approximately 10% of the overall operation and maintenance budget for commercial nuclear power plants. The goal of the Light Water Reactor Sustainability Program Physical Security Pathway is to develop tools, methods, and technologies and provide the technical basis for an optimized physical security posture. This pathway will analyze and minimize the conservatisms built into current security postures in order to reduce security costs while still ensuring adequate security and operational safety. The research performed at Idaho National Laboratory within this pathway has successfully developed a dynamic force-on-force (FOF) modeling framework using various computer simulation tools and integrated them with the dynamic assessment Event Modeling Risk Assessment using Linked Diagrams (EMRALD) tool. This document provides an overview of lessons learned in applying a dynamic computational framework that links results from a commercially available FOF simulation tool, a commercially available thermal-hydraulic tool, and EMRALD to an operating commercial nuclear power plant. This process of including plant procedures and multiple analysis results is being called Modeling and Analysis for Safety Security using Dynamic EMRALD Framework. Previous reports described how a user could integrate their plant-specific FOF models with the dynamic simulation tool EMRALD, model operator actions, integrate with probabilistic risk assessment tools, such as Computer Aided Fault Tree Analysis System or Systems Analysis Programs for Hands-on Integrated Reliability Evaluations, and with thermal-hydraulic tools, such as RELAP-5. Previous reports applied various combinations of available simulations codes with EMRALD using generic plant models to demonstrate how to perform the analysis. This report documents the results of applying the dynamic computational framework to an actual nuclear facility using their security scenarios and timelines. The purpose of this study was to verify that results achieved using generic models are similar to actual plant results and to refine our guidance of the use of the framework. Such an assessment enables further analysis, such as what-if scenarios and staff-reduction evaluation, thereby optimizing physical security at plants. NOTE: The work performed in this report is based on a generic EMRALD model with actual plant data used for the analysis. However, only the generic model and general results of the analysis are in the report. No plant’s sensitive information is discussed in this report. The discussion shows examples of insights that can be obtained from the MASS-DEF methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1977

NASTRAN hydroelastic modal studies. Volume 1: Introduction, theory, and results

The implementation of a general three-dimensional hydroelastic capability in the NASTRAN computer program is described. A method for analyzing the combined mode shapes of arbitrary fluid and structure finite element models was developed. The fluid was modeled with three-dimensional solid elements with options for tetrahedron, wedge, and hexahedron shapes. The elements were then connected to fluid grid points which defined the pressure in the fluid at the specified location. Matrices were formulated to provide efficient solutions for large-order problems. The structure matrices were processed separately and could be reduced using matrix condensation procedures or through a modal formulation using the normal modes of the empty structure as solution coordinates. The fluid/structure interface and the free surface were defined by the user with special NASTRAN boundary elements. A special purpose mesh generator program was provided to generate the actual NASTRAN data cards for the fluid, the structure, and the boundary elements for typical tank-type models.

Source record↗

NASA NTRS2014

Influence of Triply-Charged Ions and Ionization Cross-Sections in a Hybrid-PIC Model of a Hall Thruster Discharge

The sensitivity of xenon ionization rates to collision cross-sections is studied within the framework of a hybrid-PIC model of a Hall thruster discharge. A revised curve fit based on the Drawin form is proposed and is shown to better reproduce the measured crosssections at high electron energies, with differences in the integrated rate coefficients being on the order of 10% for electron temperatures between 20 eV and 30 eV. The revised fit is implemented into HPHall and the updated model is used to simulate NASA's HiVHAc EDU2 Hall thruster at discharge voltages of 300, 400, and 500 V. For all three operating points, the revised cross-sections result in an increase in the predicted thrust and anode efficiency, reducing the error relative to experimental performance measurements. Electron temperature and ionization reaction rates are shown to follow the trends expected based on the integrated rate coefficients. The effects of triply-charged xenon are also assessed. The predicted thruster performance is found to have little or no dependence on the presence of triply-charged ions. The fraction of ion current carried by triply-charged ions is found to be on the order of 1% and increases slightly with increasing discharge voltage. The reaction rates for the 0→III, I→III, and II→III ionization reactions are found to be of similar order of magnitude and are about one order of magnitude smaller than the rate of 0→II ionization in the discharge channel.

Hall Thrusters↗

DOE OSTI2024

Source term analysis of FeCrAl accident tolerant fuel using MELCOR

It has been established that incremental improvements in beyond design basis accident performance can be achieved through accident tolerant fuel (ATF). However, they have the potential to recover margin with respect to conventional fuel and therefore enhance plant economics through uprate or cycle length increase. To realize this potential, it is necessary to quantify the reduction in source term due to use of ATF, and correspondingly how this is affected by increasing the cycle length and/or burnup. This requires development of a risk informed analysis methodology for ATF under high burnup conditions, which is being developed within LWRS. To this end, a MELCOR simulation of FeCrAl ATF in a recovered Large Break LOCA (LBLOCA) scenario has been developed, using a model based on the Zion Pressurized Water Reactor (PWR). The new user defined material capability, along with the inclusion of detailed neutronics- and depletion-derived parameters such as core power distribution, decay heat behavior, and fission product inventories, allows a more detailed simulation of FeCrAl clad material properties and behavior than has previously been possible using MELCOR. These detailed FeCrAl results were compared with a zircaloy clad model to investigate the differences between the two clad types and quantify the benefits of FeCrAl ATF with an 18-month cycle. Next, the fuel cycle with ATF was extended to 24 months, to determine whether any additional safety margin provided by FeCrAl ATF could be leveraged to implement high burnup FeCrAl-clad fuel while retaining the same operating and safety limits as current zircaloy-clad fuel. For the particular scenario analyzed, the delay to fuel failure from using ATF was of the same order as the LOCA recovery time, and hence significant in reducing fission product release, with the 24-month FeCrAl core performing better than the 18-month zircaloy core. It is noted that for other transients, the reduction in release due to using FeCrAl may be less significant. Furthermore, the material model developed here can be used in such further studies in support of determining the overall source term.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗

NASA NTRS1983

Magnetic drifts at Io - Depletion of 10-MeV electrons at Voyager 1 encounter due to a forbidden zone

A model for magnetic field lines is used to track the Voyager spacecraft trajectory past Io back to the equatorial plane in order to determine the cause of an observed depletion of 10 MeV electrons in Io's orbit. The field lines are in the vicinity of an Alfven wing, a noncompressive field perturbation. It is suggested that the Jovian magnetospheric drift features influence the behavior of particles around Io, which also has a zone of reduced convection velocity. An analytical model is developed that demonstrates that a forbidden zone may exist around Io due to the gradient drifts and reduced convection velocity, and projections of particles convecting to the Io surface agree well with the Voyager data for particle concentrations. Discrepancies in the measured 9 MeV fluxes with respect to Pioneer data indicate that a possible change in the Io conductance in the interval between the passages of the two spacecraft.

Goldstein, B. E.↗

NASA NTRS1984

Implicit Numerical Methods in Meteorology

The development of a fully implicit finite-difference model, whose time step is chosen solely to resolve accurately the physical flow of interest is discussed. The method is based on an operator factorization which reduces the dimensionality of the implicit approach: at each time step only (spatially) one-dimensional block-tridiagonal linear systems must be solved. The scheme uses two time levels and is second-order accurate in time. Compact implicit spatial differences are used, yielding fourth-order accuracy both vertically and horizontally. In addition, the development of a fully interactive computer code is discussed. With this code the user will have a choice of models, with various levels of accuracy and sophistication, which are imbedded, as subsets of the fully implicit 3D code.

Augenbaum, J.↗

NASA NTRS1996

A Computational and Experimental Study of Nonlinear Aspects of Induced Drag

Despite the 80-year history of classical wing theory, considerable research has recently been directed toward planform and wake effects on induced drag. Nonlinear interactions between the trailing wake and the wing offer the possibility of reducing drag. The nonlinear effect of compressibility on induced drag characteristics may also influence wing design. This thesis deals with the prediction of these nonlinear aspects of induced drag and ways to exploit them. A potential benefit of only a few percent of the drag represents a large fuel savings for the world's commercial transport fleet. Computational methods must be applied carefully to obtain accurate induced drag predictions. Trefftz-plane drag integration is far more reliable than surface pressure integration, but is very sensitive to the accuracy of the force-free wake model. The practical use of Trefftz plane drag integration was extended to transonic flow with the Tranair full-potential code. The induced drag characteristics of a typical transport wing were studied with Tranair, a full-potential method, and A502, a high-order linear panel method to investigate changes in lift distribution and span efficiency due to compressibility. Modeling the force-free wake is a nonlinear problem, even when the flow governing equation is linear. A novel method was developed for computing the force-free wake shape. This hybrid wake-relaxation scheme couples the well-behaved nature of the discrete vortex wake with viscous-core modeling and the high-accuracy velocity prediction of the high-order panel method. The hybrid scheme produced converged wake shapes that allowed accurate Trefftz-plane integration. An unusual split-tip wing concept was studied for exploiting nonlinear wake interaction to reduced induced drag. This design exhibits significant nonlinear interactions between the wing and wake that produced a 12% reduction in induced drag compared to an equivalent elliptical wing at a lift coefficient of 0.7. The performance of the split-tip wing was also investigated by wing tunnel experiments. Induced drag was determined from force measurements by subtracting the estimated viscous drag, and from an analytical drag-decomposition method using a wake survey. The experimental results confirm the computational prediction.

Smith, Stephen C.↗

DOE OSTI2021

PDG: A Composite Method Based on the Resolution of the Identity

The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller–Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles–doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the “standard” G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. Furthermore, the efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4–5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗