SEARCH · Search NASA

Results for “Reduced Order 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,801 records · Page 100

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Coupled Inertial Navigation and Flush Air Data Sensing Algorithm for Atmosphere Estimation

This paper describes an algorithm for atmospheric state estimation based on a coupling between inertial navigation and flush air data-sensing pressure measurements. The navigation state is used in the atmospheric estimation algorithm along with the pressure measurements and a model of the surface pressure distribution to estimate the atmosphere using a nonlinear weighted least-squares algorithm. The approach uses a high-fidelity model of atmosphere stored in table-lookup form, along with simplified models propagated along the trajectory within the algorithm to aid the solution. Thus, the method is a reduced-order Kalman filter in which the inertial states are taken from the navigation solution and atmospheric states are estimated in the filter. The algorithm is applied to data from the Mars Science Laboratory entry, descent, and landing from August 2012. Reasonable estimates of the atmosphere are produced by the algorithm. The observability of winds along the trajectory are examined using an index based on the observability Gramian and the pressure measurement sensitivity matrix. The results indicate that bank reversals are responsible for adding information content. The algorithm is applied to the design of the pressure measurement system for the Mars 2020 mission. A linear covariance analysis is performed to assess estimator performance. The results indicate that the new estimator produces more precise estimates of atmospheric states than existing algorithms.

Karlgaard, Christopher D.↗

DOE OSTI2019

Enhancing Monte Carlo proxy applications on GPUs

In Monte Carlo neutron transport simulations, a computational routine commonly known as the ”cross-section lookup” has been identified as being the most computationally expensive part of these applications. A tool which is commonly used as a proxy application for these routines, named ”XSBench”, was created to simulate popular algorithms used in these routines on CPUs. Currently, however, as GPU-based HPC resources have become more widely available, there has been significant interest and efforts invested in moving these traditionally CPU-based simulations to GPUs. Unfortunately, the algorithms commonly used in the cross-section lookup routine were originally devised and developed for CPU-based platforms, and have seen limited study on GPUs to date. Additionally, platforms such as XSBench implement approximations which may have a negligible effect on CPUs, but may be quite impactful to performance on GPUs given the more resource-limited nature of the latter. As a result, we have created VEXS, a new tool for modeling the cross-section lookup routine which removes or at least reduces the approximations made by XSBench in order to provide a more realistic prediction of algorithm performance on GPUs. In this paper, we detail our efforts to remove and reduce these approximations, show the resulting improvement in performance prediction in comparison to a reference production code, Shift, and provide some basic profiling analysis of the resulting application.

Shriver, Forrest↗

DOE OSTI2016

On-the-Fly Treatment of Temperature Dependent Cross Sections in the Unresolved Resonance Regions in RMC Code

Monte Carlo method can provide high fidelity neutronics analysis of different nuclear reactors, owing to its advantages of the flexible geometry modeling and the use of continuous-energy nuclear cross sections. It can also be coupled with thermal-hydraulics codes to consider the feedback. For thermal reactors such as PWR and HTGR, two main temperature effects should be considered. The first is the effect of thermal motion of target nuclei in resolved resonance energy regions, which is known as Doppler effect. The second is the thermal scattering and bound effect in thermal energy regions. The traditional approach of pre-generated cross sections has difficulty in memory footprint for detailed temperature modelling in multi-physics calculations. Recently, on-the-fly (OTF) technique has been proposed in order to reduce the memory consumption for both resolved resonance energy and thermal energy. For fast reactors calculations, the temperature dependence of cross sections in the unresolved resonance energy regions is also important. However, less attention was paid to the unresolved resonance regions compared with resolved resonance energy and thermal energy. In this paper, the on-the-fly temperature treatment of cross sections in unresolved resonance regions was proposed and developed in RMC codes. The results show that the on-the-fly treatment has high efficiency and satisfactory fidelity. The on-the-fly treatment of unresolved resonance region was also combined with the target motion sampling (TMS) method for resolved resonance energy. The formats of libraries were changed, in which the probability tables were separated from the cross sections. Then the interpolation of probability table was performed on-the-fly during the transport routine. The results show that the on-the-fly treatment has high efficiency and satisfactory fidelity. The on-the-fly treatment of unresolved resonance region was also combined with the target motion sampling (TMS) method for resolved resonance energy, showing that TMS method is effective and accurate for cross sections treatment in resolved resonance energy regions. The memory saving of the on-the-fly treatment is also considerable. For the high-fidelity 'transport-burnup-thermal-hydraulics' coupling calculation, the memory requirement is 60 Mb. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1983

Helicopter rotor loads using discretized matched asymptotic expansions

The numerical practicality of a matched asymptotic expansion approach for the computation of unsteady three dimensional airloads on a helicopter rotor was improved. This effort utilizes a discretized repesentation of the doublet strength distribution and helical streamlines to decrease the computational requirements of the original analysis. The continuous variation of the doublet strength was approximated by piecewise constant or piecewise quadratic distributions, and the helical trajectory of a fluid particle was approximated by connected straight line segments. As a direct result of these simplified representations the computational time required for the execution of a typical flight condition was reduced by an order of magnitude with respect to the requirements of the original analysis. Airloads which were computed using the discretized method for a two bladed model rotor and a full scale four bladed rotor are in close agreement with measured results and airloads from the original asymptotic analysis. For conditions characterized by significant rotor/wake interaction the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable accuracy.

Pierce, G. A.↗

NASA NTRS2005

Thermal Analysis Methods for Aerobraking Heating

As NASA begins exploration of other planets, a method of non-propulsively slowing vehicles at the planet, aerobraking, may become a valuable technique for managing vehicle design mass and propellant. An example of this is Mars Reconnaissance Orbiter (MRO), which will launch in late 2005 and reach Mars in March of 2006. In order to save propellant, MRO will use aerobraking to modify the initial orbit at Mars. The spacecraft will dip into the atmosphere briefly on each orbit, and during the drag pass, the atmospheric drag on the spacecraft will slow it, thus lowering the orbit apoapsis. The largest area on the spacecraft, and that most affected by the heat generated during the aerobraking process, is the solar arrays. A thermal analysis of the solar arrays was conducted at NASA Langley, to simulate their performance throughout the entire roughly 6-month period of aerobraking. Several interesting methods were used to make this analysis more rapid and robust. Two separate models were built for this analysis, one in Thermal Desktop for radiation and orbital heating analysis, and one in MSC.Patran for thermal analysis. The results from the radiation model were mapped in an automated fashion to the Patran thermal model that was used to analyze the thermal behavior during the drag pass. A high degree of automation in file manipulation as well as other methods for reducing run time were employed, since toward the end of the aerobraking period the orbit period is short, and in order to support flight operations the runs must be computed rapidly. All heating within the Patran Thermal model was combined in one section of logic, such that data mapped from the radiation model and aeroheating model, as well as skin temperature effects on the aeroheating and surface radiation, could be incorporated easily. This approach calculates the aeroheating at any given node, based on its position and temperature as well as the density and velocity at that trajectory point. Run times on several different processors, computer hard drives, and operating systems (Windows versus Linux) were evaluated.

Amundsen, Ruth M.↗

NASA NTRS2014

Radiation Quality Effects on Transcriptome Profiles in 3-D Cultures After Charged Particle Irradiation

In this work, we evaluated the differential effects of low- and high-LET radiation on 3-D organotypic cultures in order to investigate radiation quality impacts on gene expression and cellular responses. Current risk models for assessment of space radiation-induced cancer have large uncertainties because the models for adverse health effects following radiation exposure are founded on epidemiological analyses of human populations exposed to low-LET radiation. Reducing these uncertainties requires new knowledge on the fundamental differences in biological responses (the so-called radiation quality effects) triggered by heavy ion particle radiation versus low-LET radiation associated with Earth-based exposures. In order to better quantify these radiation quality effects in biological systems, we are utilizing novel 3-D organotypic human tissue models for space radiation research. These models hold promise for risk assessment as they provide a format for study of human cells within a realistic tissue framework, thereby bridging the gap between 2-D monolayer culture and animal models for risk extrapolation to humans. To identify biological pathway signatures unique to heavy ion particle exposure, functional gene set enrichment analysis (GSEA) was used with whole transcriptome profiling. GSEA has been used extensively as a method to garner biological information in a variety of model systems but has not been commonly used to analyze radiation effects. It is a powerful approach for assessing the functional significance of radiation quality-dependent changes from datasets where the changes are subtle but broad, and where single gene based analysis using rankings of fold-change may not reveal important biological information.

Patel, Zarana S.↗

NASA NTRS1990

Concurrent processing adaptation of aeroplastic analysis of propfans

Discussed here is a study involving the adaptation of an advanced aeroelastic analysis program to run concurrently on a shared memory multiple processor computer. The program uses a three-dimensional compressible unsteady aerodynamic model and blade normal modes to calculate aeroelastic stability and response of propfan blades. The identification of the computational parallelism within the sequential code and the scheduling of the concurrent subtasks to minimize processor idle time are discussed. Processor idle time in the calculation of the unsteady aerodynamic coefficients was reduced by the simple strategy of appropriately ordering the computations. Speedup and efficiency results are presented for the calculation of the matched flutter point of an experimental propfan model. The results show that efficiencies above 70 percent can be obtained using the present implementation with 7 processors. The parallel computational strategy described here is also applicable to other aeroelastic analysis procedures based on panel methods.

Janetzke, David C.↗

NASA NTRS1990

Concurrent processing adaptation of aeroelastic analysis of propfans

This paper reports on a study involving the adaptation of an advanced aeroelastic analysis program to run concurrently on a shared memory multiple processor computer. The program uses a three-dimensional compressible unsteady aerodynamic model and blade normal modes to calculate aeroelastic stability and response of propfan blades. The identification of the computational parallelism within the sequential code and the scheduling of the concurrent subtasks to minimize processor idle time are discussed. Processor idle time in the calculation of the unsteady aerodynamic coefficients was reduced by the simple strategy of appropriately ordering the computations. Speedup and efficiency results are presented for the calculation of the matched flutter point of an experimental propfan model. The results show that efficiencies above 70 percent can be obtained using the present implementation with 7 processors. The parallel computational strategy described here is also applicable to other aeroelastic analysis procedures based on parallel methods.

Murthy, Durbha V.↗

NASA NTRS1992

Concurrent processing adaptation of aeroelastic analysis of propfans

Murthy, D. V.↗

DOE OSTI2020

Interplay of the Spin Density Wave and a Possible Fulde-Ferrell-Larkin-Ovchinnikov State in CeCoIn 5 in Rotating Magnetic Field

The d-wave superconductor CeCoIn 5 has been proposed as a strong candidate for supporting the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state near the low-temperature boundary of its upper critical field. Neutron diffraction, however, finds spin-density-wave (SDW) order in this part of the phase diagram for field in the a-b plane, and evidence for the SDW disappears as the applied field is rotated toward the tetragonal c axis. It is important to understand the interplay between the SDW and a possible FFLO state in CeCoIn 5 , as the mere existence of an SDW does not necessarily exclude an FFLO state. In this work, based on a model constructed on the basis of available experiments, we show that an FFLO state competes with an SDW phase. The SDW state in CeCoIn 5 is stabilized when the field is directed close to the a-b plane. When the field is rotated toward the c axis, the FFLO state emerges, and the SDW phase disappears. In the FFLO state, the nodal planes with extra quasiparticles (where the superconducting order parameter is zero) are perpendicular to the field, and in the SDW phase, the quasiparticle density of states is reduced. We test this model prediction by measuring heat transported by normal quasiparticles in the superconducting state. As a function of field, we observe a reduction of thermal conductivity for field close to the a-b plane and an enhancement of thermal conductivity when field is close to the c axis, consistent with theoretical expectations. Our modeling and experiments, therefore, indicate the existence of the FFLO state when field is parallel to the c axis.

36 MATERIALS SCIENCE↗

DOE OSTI2015

High Fidelity Forward Model Development for Nuclear Reactor Spent Fuel Technical Nuclear Forensics

The fidelity of the forward model within a spent fuel forensic analysis system was improved by using two unique methodologies. The first consisted of developing a system to create accurate one-group neutron cross-section libraries for any user specified reactor system. In such, a detailed model is developed using the depletion code MONTEBURNS. During MONTEBURNS execution, cross-section libraries are generated at every user specified burnup step in time. These libraries could be developed for many reactor systems, then housed in a database and used for analyzing unknown fuel samples. The forensic analysis system for spent fuel resulted in higher accuracy at predicting the initial uranium isotopic compositions and burnup from spent fuel samples. Using this method, the error in results was reduced from the order of 1-6% down to less than 1% when recovering a fuel sample's burnup and initial uranium isotopic composition. The second method consisted of implementing 2D/3D reactor depletion codes as the forward model within the system's framework. This method would allow the usage of potentially recoverable geometric information from an unknown sample. No predetermined cross-section library is required for the system using this method, therefore potentially reducing model error associated with the neutron flux spectrum. The accuracy of the recovered initial uranium isotopic compositions and burnup from spent fuel samples was also improved using this method, even more so than the first. For MTR reactors, the error using this method was significantly reduced and was driven to below 0.5%. However, additional research may be required to determine the ideal fission yield and recoverable energy per fission for cases where significant amounts of 239 PU are bred and burned throughout the life of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Full core LOCA analysis for BWR/6 - Methodology and first results

As part of collaboration between the Paul Scherrer Institute (PSI) and the Swiss regulator ENSI, a study aimed at core-wide estimates of the fuel rod behaviour during large-break loss-of-coolant-accidents (LB-LOCA) in BWR/6 plants was recently conducted on the basis of a Swiss reactor. To perform such full core LOCA analysis, a methodology consisting of a multi-physics sequence of interfaced best-estimate models and codes was established: CASMO-4E/SIMULATE-3 (C4/S3) for the core modelling, TRACE for the plant system analysis and FALCON for the fuel thermomechanical calculations during base irradiation (BI) as well as for the transient analyses (TA). As starting point, reference C4/S3 models validated for all operated cycles of the given plant were used to select as situation target, a recent cycle representative of modern BWR core loadings and fuel assembly designs. Based on the selected C4/S3 cycle model, PSI in-house computational modules were developed or updated in order to set up the TRACE full core model with regards to assembly geometries, fuel rods thermophysical properties and initial 3-D power distributions. Also, an interface to FALCON BI models was established for the transfer of the irradiation history, including primarily linear heat generation rates (LHGR), for an equivalent fuel pin. For each assembly, the power of the equivalent pin is defined with the maximum power of all pins at each axial position. However, in order to reduce the computational time of the fuel performance analyses, the whole core was subdivided into groups of assemblies, each represented by the most powerful assembly within the group. On this basis, the FALCON BI analyses were conducted with coupling to the PSI advanced model for gas release and swelling (GRWS-A) and an interface was developed to transfer the calculated fuel rod characteristics at end of irradiation to initialize the core fuel rods heat structures in the TRACE model. Subsequently, double-ended recirculation line rupture scenarios at hot-full-power (HFP) conditions were investigated with TRACE in order to determine a conservative transient case as function of cycle burnup, single failure criterion and availability of safety systems. For the most conservative case of this study, an interface was set-up to transfer the TRACE results as boundary conditions to the FALCON TA models, including thus primarily time and axial distributions of the clad outer temperatures. As final step, the FALCON TA were conducted for the representative assemblies in order to infer core-wide estimates of the fuel rod behaviour with specific focus on equivalent cladding reacted levels, hoop strains as well as margins to failure criteria. The objective of this paper is to present this methodology along with the obtained first results. Regarding the latter, it is found that principally, even for the most conservative LOCA scenario of this study, the clad temperature excursion remains mild enough at all core locations such that no significant rod ballooning and no cladding failure is predicted, something in line with corresponding analyses conducted earlier for BWR/4 plants. (authors)

42 ENGINEERING↗

DOE OSTI2020

Commensuration Effects in Layered Nanoparticle Solids.

We have developed a hierarchical nanoparticle transport simulator (HINTS) and adapted it to study commensuration effects in two classes of nanoparticle (NP) solids: (1) a bilayer NP solid (BNS) with an energy offset and (2) a BNS as part of a field effect transistor (FET). HINTS integrates the ab initio characterization of single NPs with the phonon-assisted tunneling transition model of the NP-NP transitions into a kinetic Monte Carlo based simulation of the charge transport in NP solids. First, we studied a BNS with an interlayer energy offset Delta, possibly caused by a fixed electric field. Our results include the following. First, in the independent energy-offset model, we observed the emergence of commensuration effects when scanning the electron filling factor (FF) across integer values. These commensuration effects were profound as they reduced the mobility by several orders of magnitude. We analyzed these commensuration effects in a five-dimensional parameter space, as a function of the on-site charging energy E-C, energy offset Delta, disorder D, the electron FF, and temperature k(B)T. We demonstrated the complexity of our model by showing that at integer FFs commensuration effects are present in some regions of the parameter space, while they vanish in other regions, thus defining distinct dynamical phases of the model. We determined the phase boundaries between these dynamical phases. Second, using these results as a foundation, we shifted our focus to the experimentally much-studied NP-FETs. NP-FETs are also characterized by an interlayer energy offset Delta, which, in contrast to our first model, is set by the gate voltage V-G and thereby related to the electron FF. We repeated many of our simulations and again demonstrated the emergence of commensuration effects and distinct dynamical phases in these NP-FETs. Notably, the commensuration effects in the NP-FETs showed many similarities to those in the independent energy-offset BNS.

Qu, Luman↗

DOE OSTI2018

A flux-balanced fluid model for collisional plasma edge turbulence: Model derivation and basic physical features

We propose a new reduced fluid model for the study of the drift wave–zonal flow dynamics in magnetically confined plasmas. Our model can be viewed as an extension of the classic Hasegawa-Wakatani (HW) model and is based on an improved treatment of the electron dynamics parallel to the field lines, to guarantee a balanced electron flux on the magnetic surfaces. Our flux-balanced HW (bHW) model contains the same drift-wave instability as previous HW models, but unlike these models, it converges exactly to the modified Hasegawa-Mima model in the collisionless limit. We rely on direct numerical simulations to illustrate some of the key features of the bHW model, such as the enhanced variability in the turbulent fluctuations and the existence of stronger and more turbulent zonal jets than the jets observed in other HW models, especially for high plasma resistivity. Our simulations also highlight the crucial role of the feedback of the third-order statistical moments in achieving a statistical equilibrium with strong zonal structures. Finally, we investigate the changes in the observed dynamics when more general dissipation effects are included and, in particular, when we include the reduced model for ion Landau damping originally proposed by Wakatani and Hasegawa.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Shape restricted additive hazards models: Monotone, unimodal, and U‐shape hazard functions

We consider estimation of the semiparametric additive hazards model with an unspecified baseline hazard function where the effect of a continuous covariate has a specific shape but otherwise unspecified. Such estimation is particularly useful for a unimodal hazard function, where the hazard is monotone increasing and monotone decreasing with an unknown mode. A popular approach of the proportional hazards model is limited in such setting due to the complicated structure of the partial likelihood. Our model defines a quadratic loss function, and its simple structure allows a global Hessian matrix that does not involve parameters. Thus, once the global Hessian matrix is computed, a standard quadratic programming method can be applicable by profiling all possible locations of the mode. However, the quadratic programming method may be inefficient to handle a large global Hessian matrix in the profiling algorithm due to a large dimensionality, where the dimension of the global Hessian matrix and number of hypothetical modes are the same order as the sample size. We propose the quadratic pool adjacent violators algorithm to reduce computational costs. The proposed algorithm is extended to the model with a time‐dependent covariate with monotone or U‐shape hazard function. In simulation studies, our proposed method improves computational speed compared to the quadratic programming method, with bias and mean square error reductions. We analyze data from a recent cardiovascular study.

Mathematical & Computational Biology↗

NASA NTRS1991

Calculations of particle-laden flows by second-moment closures

A differential second-moment closure is applied to calculate nonswirling turbulent particle-laden jet flow. The closure is effected by extending constant-density single-phase turbulence models to the continuous (carrier) phase of two-phase flow. The dispersed (discrete) phase is treated by stochastic-Lagrangian technique. To reduce the effects of numerical (false) diffusion on the predicted results, a higher order differencing scheme, namely the flux-spline, is employed. In addition, an experimental study is conducted to provide data of good quality, especially near the inlet, for model assessment. The experiment consists of a confined coaxial jet flow of air with glass beads. The beads' diameter is 105 microns and is injected from the inner jet with a mass loading ratio of 0.9. The computed results are in good agreement with experimental data, especially for the carrier phase. However, the predicted turbulence fluctuations associated with the discrete phase are lower than the data at the downstream region of the flow field.

Nikjooy, M.↗

NASA NTRS1973

Orbit determination for low-thrust spacecraft: Concepts and analysis

Earth-based orbit determination capability for SEP spacecraft in multistation tracking and in thrust subsystem error modeling is described. Five different tracking strategies are applied to a 15 day segment of an Encke rendezvous mission. Both optimal and suboptimal orbit determination performance are determined for a wide range of process noise parameter values. The multi-station tracking techniques are found to be extremely effective, reducing orbit determination errors by orders of magnitude over that obtained with conventional single-station tracking. Explicitly differenced multistation data (QVLBI) is found to be least sensitive to gross modeling errors, but if a reasonably good process noise model is available, explicit differencing is not required.

Mcdanell, J. P.↗

Refine your search to explore more results.