Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deterministic method”

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 793 records · Page 44

Designing an Optimal LSST Deep Drilling Program for Cosmology with Type Ia Supernovae

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) is forecast to collect a large sample of Type Ia supernovae (SNe Ia) expected to be instrumental in unveiling the nature of dark energy. The feat, however, requires accurately measuring the two components of the Hubble diagram, distance modulus and redshift. Distance is estimated from SN Ia parameters extracted from light-curve fits, where the average quality of light curves is primarily driven by survey parameters. An optimal observing strategy is thus critical for measuring cosmological parameters with high accuracy. We present in this paper a three-stage analysis to assess the impact of the deep drilling (DD) strategy parameters on three critical aspects of the survey: redshift completeness, the number of well-measured SNe Ia, and cosmological measurements. We demonstrate that the current DD survey plans (internal LSST simulations) are characterized by a low completeness (z ~ 0.55–0.65), and irregular and low cadences (several days), which dramatically decrease the size of the well-measured SN Ia sample. We propose a method providing the number of visits required to reach higher redshifts. We use the results to design a set of optimized DD surveys for SN Ia cosmology taking full advantage of spectroscopic resources for host galaxy redshift measurements. The most accurate cosmological measurements are achieved with deep rolling surveys characterized by a high cadence (1 day), a rolling strategy (at least two seasons of observation per field), and ultradeep (z ≳ 0.8) and deep (z ≳ 0.6) fields. A deterministic scheduler including a gap recovery mechanism is critical to achieving a high-quality DD survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-stabilizing byzantine-fault-tolerant clock synchronization system and method

Systems and methods for rapid Byzantine-fault-tolerant self-stabilizing clock synchronization are provided. The systems and methods are based on a protocol comprising a state machine and a set of monitors that execute once every local oscillator tick. The protocol is independent of specific application specific requirements. The faults are assumed to be arbitrary and/or malicious. All timing measures of variables are based on the node's local clock and thus no central clock or externally generated pulse is used. Instances of the protocol are shown to tolerate bursts of transient failures and deterministically converge with a linear convergence time with respect to the synchronization period as predicted.

Malekpour, Mahyar R.↗

NESC GN&C TDT Workshop on 2D Image Motion Optical Transfer Functions, Pointing Performance Analysis, and Requirements

What You Will Learn: The focus is on payload imaging performance due to pointing motion. Some historical background on pointing performance analysis is given. The Optical Transfer Function (OTF) and Modulation Transfer Function (MTF) are defined. The imaging performance due to pointing motion is measured by image motion optical transfer functions (IM OTF). IM OTFs are defined for displacement, smear, and jitter motions, which are all rigorously defined. Deterministic and Statistical IM OTFs are briefly derived and graphically illustrated and compared. The IM OTFs are parameterized by pointing error metrics(PEM), which are means and covariances of displacement, smear, and jitter. Emphasis is on procedures and algorithms to evaluate the image motion optical transfer functions and pointing error metrics. Three procedures are covered, which depend on whether the pointing error data is from time-domain simulation, frequency-domain analysis, or stochastic modeling. A method to evaluate the relative contribution of disturbance sources and to identify the most significant contributors is presented. The presentation includes pertinent discussion of flexible structures and control-structure interaction. No single book can adequately cover this subject, so a book is not required for the course. A list of selected articles, reports, documents, and books is provided for reference and further study. Mathis kept to the minimum necessary to convey principles; lengthy derivations are left to the reference material. Graphics are used to illustrate concepts. As with any such learning endeavor, the knowledge gained will be retained and strengthened through actual practice.c©2019–2021 Mark E. Pittelkau— 5

NASA Engineering and Safety Center (NESC)↗

A mixed $\ell_1$ regularization approach for sparse simultaneous approximation of parameterized PDEs

We introduce and assess a novel sparse polynomial technique for the simultaneous approximation of parameterized partial differential equations (PDEs) with deterministic and stochastic inputs. Our approach treats the numerical solution as a jointly sparse reconstruction problem through the reformulation of the standard basis pursuit denoising, where the set of jointly sparse vectors is infinite. To achieve global reconstruction of sparse solutions to parameterized elliptic PDEs over both physical and parametric domains, we combine the standard measurement scheme developed for compressed sensing in the context of bounded orthonormal systems with a novel mixed-norm based $\ell_1$ regularization method that exploits both energy and sparsity. Moreover, we are able to prove that, with minimal sample complexity, error estimates comparable to the best $s$-term and quasi-optimal approximations are achievable, while requiring only {\em a priori} bounds on polynomial truncation error with respect to the energy norm.Finally, we perform extensive numerical experiments on several high-dimensional parameterized elliptic PDE models to demonstrate the superior recovery properties of the proposed approach.

97 MATHEMATICS AND COMPUTING↗

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Path Toward a Unifid Geometry for Radiation Transport

The Direct Accelerated Geometry for Radiation Analysis and Design (DAGRAD) element of the RadWorks Project under Advanced Exploration Systems (AES) within the Space Technology Mission Directorate (STMD) of NASA will enable new designs and concepts of operation for radiation risk assessment, mitigation and protection. This element is designed to produce a solution that will allow NASA to calculate the transport of space radiation through complex computer-aided design (CAD) models using the state-of-the-art analytic and Monte Carlo radiation transport codes. Due to the inherent hazard of astronaut and spacecraft exposure to ionizing radiation in low-Earth orbit (LEO) or in deep space, risk analyses must be performed for all crew vehicles and habitats. Incorporating these analyses into the design process can minimize the mass needed solely for radiation protection. Transport of the radiation fields as they pass through shielding and body materials can be simulated using Monte Carlo techniques or described by the Boltzmann equation, which is obtained by balancing changes in particle fluxes as they traverse a small volume of material with the gains and losses caused by atomic and nuclear collisions. Deterministic codes that solve the Boltzmann transport equation, such as HZETRN [high charge and energy transport code developed by NASA Langley Research Center (LaRC)], are generally computationally faster than Monte Carlo codes such as FLUKA, GEANT4, MCNP(X) or PHITS; however, they are currently limited to transport in one dimension, which poorly represents the secondary light ion and neutron radiation fields. NASA currently uses HZETRN space radiation transport software, both because it is computationally efficient and because proven methods have been developed for using this software to analyze complex geometries. Although Monte Carlo codes describe the relevant physics in a fully three-dimensional manner, their computational costs have thus far prevented their widespread use for analysis of complex CAD models, leading to the creation and maintenance of toolkit-specific simplistic geometry models. The work presented here builds on the Direct Accelerated Geometry Monte Carlo (DAGMC) toolkit developed for use with the Monte Carlo N-Particle (MCNP) transport code. The workflow for achieving radiation transport on CAD models using MCNP and FLUKA has been demonstrated and the results of analyses on realistic spacecraft/habitats will be presented. Future work is planned that will further automate this process and enable the use of multiple radiation transport codes on identical geometry models imported from CAD. This effort will enhance the modeling tools used by NASA to accurately evaluate the astronaut space radiation risk and accurately determine the protection provided by as-designed exploration mission vehicles and habitats

Lee, Kerry↗

Using interfaces to: create strongly-coupled magnetic-ferroelectrics

Starting at the level of electrons and atoms our long-term goal is to rationally design complex oxide heterostructures and interface-materials with targeted emergent behaviors. This is not a matter of simply optimizing material parameters, but rather begins with understanding a mechanism to control the interplay between the diverse microscopic degrees of freedom prevalent in complex oxides in order to create targeted macroscopic phenomena and ends with the design of new material realizations. These realizations are in turn created with atomic-layer precision, structurally assessed to see that they are the intended realization, and finally their relevant properties are measured. During this program we have developed the scientific ideas necessary to apply this materials-by-design paradigm to the creation of materials offering electrical control of magnetism. At the heart of our approach lies a partnership between theory, synthesis, and characterization. These are the areas of expertise of the three co-PIs, who have worked closely together since their arrival at Cornell. During this project bold ideas on where to put the atoms to provide electrical control of magnetism were met with the latest methods of synthesis science and electron microscopy in an attempt create targeted new multiferroics. Our work applying this materials-by-design methodology to multiferroics has enjoyed multiple successes. The Co-PIs have invented new mechanisms for multiferroics, which provide strong coupling between polarization and magnetism, and reduced them to practice. Specifically, we pioneered multiferroics based on spin-phonon coupling, rotation-driven multiferroicity, and novel geometric multiferroics. Importantly, the latter two of these mechanisms enable the deterministic switching of magnetism by an applied electric field, which is key to the application of multiferroics to future device technologies. Most recently, we performed detailed theoretical analysis on a system involving the combination of magnetoelectric and colossal magnetoresistance behaviors and predict that electrical control of a metal-insulator transition is possible in appropriately strained SmBaMn2O6. The research enabled by this grant led to 63 publications in leading refereed journals including Nature, Nature Materials, Physical Review Letters, and Advanced Materials demonstrating a materials-by-design paradigm to the creation of materials offering electrical control of magnetism.

36 MATERIALS SCIENCE↗

Investigation of HZETRN 2010 as a Tool for Single Event Effect Qualification of Avionics Systems - Part II

An accurate prediction of spacecraft avionics single event effect (SEE) radiation susceptibility is key to ensuring a safe and reliable vehicle. This is particularly important for long-duration deep space missions for human exploration where there is little or no chance for a quick emergency return to Earth. Monte Carlo nuclear reaction and transport codes such as FLUKA can be used to generate very accurate models of the expected in-flight radiation environment for SEE analyses. A major downside to using a Monte Carlo-based code is that the run times can be very long (on the order of days). A more popular choice for SEE calculations is the CREME96 deterministic code, which offers significantly shorter run times (on the order of seconds). However, CREME96, though fast and easy to use, has not been updated in several years and underestimates secondary particle shower effects in spacecraft structural shielding mass. Another modeling option to consider is the deterministic code HZETRN 20104, which includes updates to address secondary particle shower effects more accurately. This paper builds on previous work by Rojdev, et al. to compare the use of HZETRN 2010 against CREME96 as a tool to verify spacecraft avionics system reliability in a space flight SEE environment. This paper will discuss modifications made to HZETRN 2010 to improve its performance for calculating SEE rates and compare results with both in-flight SEE rates and other calculation methods.

Rojdev, Kristina↗

Application of Machine Learning Techniques to Delay Tolerant Network Routing

This dissertation discusses several machine learning techniques to improve routing in delay tolerant networks (DTNs). These are networks in which there may be long one-way trip times, asymmetric links, high error rates, and deterministic as well as non-deterministic loss of contact between network nodes, such as interplanetary satellite networks, mobile ad hoc networks and wireless sensor networks. This work uses historical network statistics to train a multi-label classifier to predict reliable paths through the network. In addition, a clustering technique is used to predict future mobile node locations. Both of these techniques are used to reduce the consumption of resources such as network bandwidth, memory and data storage that is required by replication routing methods often used in opportunistic DTN environments. Thesis contributions include: an emulation tool chain developed to create a DTN test bed for machine learning, the network and software architecture for a machine learning based routing method, the development and implementation of classification and clustering techniques and performance evaluation in terms of machine learning and routing metrics.

Dudukovich, Rachel M.↗

APPLICATION OF THE FAVOR-OCI FRACTURE MECHANICS COMPUTER PROGRAM TO ASME CODE SECTION XI, IWB-3610 FLAW ACCEPTANCE CRITERIA EVALUATIONS

This paper presents an overview of added features in a new version of the FAVOR (Fracture Analysis of Vessels Oak Ridge) computer code called FAVOR-OCI. The original FAVOR code was developed at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) under the sponsorship of the US Nuclear Regulatory Commission (NRC). FAVOR is applied by US and international nuclear power industries to perform deterministic and probabilistic fracture mechanics analyses of commercial nuclear reactor pressure vessels (RPVs). Appli-cations of FAVOR are focused on insuring that the structural integrity of aging, and increasingly embrittled, RPVs is main-tained throughout their licensed service life. Based on the final ORNL release of FAVOR, v16.1, FAVOR-OCI extends existing deterministic features of FAVOR while preserving all prev-iously-existing probabilistic capabilities of FAVOR.The objective of this paper is to describe new deterministic features in FAVOR-OCI that can be applied to analytical evaluations of planar flaws. These evaluations are consistent with the acceptance criteria of ASME Code, Section XI, Article IWB-3610, including Subarticles IWB-3611 (acceptance based on flaw size) and IWB-3612 (acceptance based on applied stress intensity factor). The linear elastic fracture mechanics (LEFM) capabilities of FAVOR-OCI also incorporate the analytical procedures presented in the Nonmandatory Appendix A, Analysis of Flaws, Article A-3000, Method of KI Determination, for both surface and subsurface (embedded) flaws.The paper describes a computational methodology for determining critical values of fracture-related parameters that satisfy ASME Code Section XI acceptance criteria for flaws exposed to multiple thermal-hydraulic transients. These com-pute-intensive analyses can be carried out with a single execution of FAVOR-OCI. The new methodology determines critical values by solving for either the point of tangency or point of intersection between applied KI versus time histories and a user-selected fracture toughness material property (e.g., ASME KIc, FAVOR Weibull KIc, or Master Curve Weibull KJc) for surface or subsurface flaws. Situations where warm prestress conditions apply can also be addressed. The paper highlights a need for this new capability via applications to a recent independent review of safety cases for RPVs in two Belgian nuclear power plants (NPPs). That review required ASME Section XI assessments of several thousand embedded, quasi-laminar flaws in the wall of each RPV. Analysis results provided by the new capability contributed to the technical bases compiled from several sources by the Belgian nuclear regulatory agency (FANC) and eventually used by FANC to justify the restart of these NPPs.

Dickson, Terry L.↗

Modeling State Transition and Head-Dependent Efficiency Curve for Pumped Storage Hydro in Look-Ahead Dispatch

As one of the most widely installed utility-scale storage facilities, pumped storage hydro (PSH) plays an essential role in providing flexibility for power systems worldwide. Thus, to accurately quantify the flexibility of PSH units in operational optimization problems is important. The conventional PSH models in the literature rarely consider detailed state transitions due to their hourly-based settings. However, it becomes imperative for operational optimization on short-term intervals, especially with increasing shares of renewable energy in power systems. To this end, this paper presents a novel deterministic PSH model that considers the transition time and trajectory between three states of PSH units in look-ahead dispatch. Moreover, to better characterize the varying efficiency of PSH units with water head and flow rate, this paper proposes to model detailed headdependent efficiency curves (hereafter called input-output curves) in look-ahead dispatch, which encounters heavy computational burdens when short time intervals are applied. In this work, a zig-zag piece-wise linear approximation method is used for input-output curve modeling. This can enable an accurate quantification for variable efficiency and head dependence in a computationally effective manner. Here, numerical results are presented to show performances of the proposed PSH model in both flexibility quantification and computation time.

13 HYDRO ENERGY↗

Quantum Foundations of Classical Reversible Computing

The reversible computation paradigm aims to provide a new foundation for general classical digital computing that is capable of circumventing the thermodynamic limits to the energy efficiency of the conventional, non-reversible digital paradigm. However, to date, the essential rationale for, and analysis of, classical reversible computing (RC) has not yet been expressed in terms that leverage the modern formal methods of non-equilibrium quantum thermodynamics (NEQT). In this paper, we begin developing an NEQT-based foundation for the physics of reversible computing. We use the framework of Gorini-Kossakowski-Sudarshan-Lindblad dynamics (a.k.a. Lindbladians) with multiple asymptotic states, incorporating recent results from resource theory, full counting statistics and stochastic thermodynamics. Important conclusions include that, as expected: (1) Landauer’s Principle indeed sets a strict lower bound on entropy generation in traditional non-reversible architectures for deterministic computing machines when we account for the loss of correlations; and (2) implementations of the alternative reversible computation paradigm can potentially avoid such losses, and thereby circumvent the Landauer limit, potentially allowing the efficiency of future digital computing technologies to continue improving indefinitely. We also outline a research plan for identifying the fundamental minimum energy dissipation of reversible computing machines as a function of speed.

97 MATHEMATICS AND COMPUTING↗

Particle-laden swirling free jets: Measurements and predictions

A theoretical and experimental investigation of single-phase and particle-laden weakly swirling jets was conducted. The jets were injected vertically downward from a 19 mm diameter tube with swirl numbers ranging from 0 to 0.33. The particle-laden jets had a single loading ratio (0.2) with particles having a SMD of 39 microns. Mean and fluctuating properties of both phases were measured using nonintrusive laser based methods while particle mass flux was measured using an isokinetic sampling probe. The continuous phase was analyzed using both a baseline kappa-epsilon turbulence model and an extended version with modifications based on the flux Richardson number to account for effects of streamline curvature. To highlight effects of interphase transport rates and particle/turbulence interactions, effects of the particles were analyzed as follows: (1) locally homogeneous flow (LHF) analysis, where interphase transport rates are assumed to be infinitely fast; (2) deterministic separated flow (DSF) analysis, where finite interphase transport rates are considered but particle/turbulence interactions are ignored; and (3) stochastic separated flow (SSF) analysis, where both effects are considered using random-walk computations.

Bulzan, D. L.↗

Particle-laden swirling free jets - Measurements and predictions

A theoretical and experimental investigation of single-phase and particle-laden weakly swirling jets was conducted. The jets were injected vertically downward from a 19 mm diameter tube with swirl numbers ranging from 0 to 0.33. The particle-laden jets had a single loading ratio (0.2) with particles having an SMD of 39 micrometers. Mean and fluctuating properties of both phases were measured using nonintrusive laser based methods while particle mass flux was measured using an isokinetic sampling probe. The continuous phase was analyzed using both a baseline k-epsilon turbulence model and an extended version with modifications based on the flux Richardson number to account for effects of streamline curvature. To highlight effects of interphase transport rates and particle/turbulence interactions, effects of the particles were analyzed as follows: (1) locally homogeneous flow (LHF) analysis, where interphase transport rates are assumed to be infinitely fast; (2) deterministic separated flow (DSF) analysis, where finite interphase transport rates are considered but particle/turbulence interactions are ignored; and (3) stochastic separated flow (SSF) analysis, where both effects are considered using random-walk computations.

Bulzan, D. L.↗

Model for radial dependence of frequency distributions for energy imparted in nanometer volumes from HZE particles

This paper develops a deterministic model of frequency distributions for energy imparted (total energy deposition) in small volumes similar to DNA molecules from high-energy ions of interest for space radiation protection and cancer therapy. Frequency distributions for energy imparted are useful for considering radiation quality and for modeling biological damage produced by ionizing radiation. For high-energy ions, secondary electron (delta-ray) tracks originating from a primary ion track make dominant contributions to energy deposition events in small volumes. Our method uses the distribution of electrons produced about an ion's path and incorporates results from Monte Carlo simulation of electron tracks to predict frequency distributions for ions, including their dependence on radial distance. The contribution from primary ion events is treated using an impact parameter formalism of spatially restricted linear energy transfer (LET) and energy-transfer straggling. We validate our model by comparing it directly to results from Monte Carlo simulations for proton and alpha-particle tracks. We show for the first time frequency distributions of energy imparted in DNA structures by several high-energy ions such as cosmic-ray iron ions. Our comparison with results from Monte Carlo simulations at low energies indicates the accuracy of the method.

Non-NASA Center↗

Aeroelastic Uncertainty Quantification of a Low-Boom Aircraft Configuration

As the state of the art in uncertainty quantification for low-boom aircraft advances, the underlying assumption of a rigid airframe must be revisited. The goal of this research is to investigate the impact of uncertainties in aeroelastic deformation of a low-boom aircraft on ground noise. Variations in structural properties and uncertainties in loading, derived from flight conditions, both factor into the overall aeroelastic deformation and subsequently the ground noise. Incorporation of these aeroelastic uncertainties in the prediction of ground noise during the design phase can lead to improved robustness. In this paper, a review of methodologies and techniques employed in low-boom uncertainty quantification will be given. In addition, methods for aeroelastic uncertainty quantification are integrated into the previous work and a generalized set of procedures is established. In a case study implementing the analysis procedures, ground noise generated from a static aeroelastic deformed low-boom aircraft increased slightly over that from the undeformed geometry for both undertrack and offtrack angles. Ground noise sensitivities to uncertainties in near field conditions and structural parameters varied significantly with atmospheric profiles. Shifts in confidence interval width in addition to shifts in deterministic values of ground noise were observed while varying atmospheric conditions.

Phillips, Benjamin D.↗

Formulations and analysis of Compton scattering for deterministic thermal radiation transport

In high-energy-density physics at low densities and high temperatures, Compton scattering can be the dominant interaction mechanism between thermal radiation and the material. Boson statistics leads to induced scattering, where the rate of scattering is proportional to the square of the photon intensity. To explore the importance of this nonlinear effect in the context of thermal radiative transfer, in this work we formulate several linear approximations to the scattering rate and compare them to the proper nonlinear rate in an infinite-medium setting. With the goal of implementing Compton scattering in a radiation transport code, we study heating-rate accuracy and iterative convergence properties for the linear methods as we vary the fidelity of the discretization in photon frequency. We find the equilibrium out-scattering formulation to be the most accurate linear method, with heating rate errors around 15% when the radiation is heating the material and up to 50% when the radiation is cooling the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Sensitivity Calculations for Systems with Polyethylene Reflector Materials Using CLUTCH

The SCALE 6.2.4 code package contains four sequences for calculating $k_{eff}$ sensitivity coefficients. Two of these sequences use deterministic transport solvers: a one-dimensional (1D) capability based on XSDRN, and a two-dimensional (2D) capability based on NEWT. These sequences are restricted to the multigroup (MG) treatment of neutron energy. The three-dimensional (3D) sequences use the KENO V.a or KENO-VI Monte Carlo transport codes and can be used to calculate sensitivity coefficients with either MG or continuous-energy (CE) transport. The 3D sensitivities are ultimately reported in an MG structure, regardless of the method used in the transport calculations. If desired, the sensitivity coefficients can be reported with very fine energy resolution from a CE calculation, but they are calculated only in the MG library structure in the MG mode. CE TSUNAMI methods are available in SCALE starting in SCALE version 6.2. Sensitivity coefficients were generated using the 3D sequences as part of the generation of the SCALE 6.2.2 Validation Report; difficulties encountered when using the CLUTCH method for thick, fissionable-material reflectors were discussed and investigated as documented in a previous paper. This paper discusses the difficulties encountered in generating accurate sensitivity coefficients using the CLUTCH technique for polyethylene reflectors for two fast spectrum benchmarks. Direct perturbation (DP) calculations were performed to confirm the accuracy of the total sensitivity coefficient for important isotopes with large sensitivities in the system. Discrepancies were detected for CLUTCH-calculated sensitivity coefficients in the reflector of a critical experiment with a radial polyethylene reflector. A simple polyethylene-reflected plutonium sphere was then used to further investigate the discrepancy. Calculations performed using the iterated fission probability (IFP) method generated accurate sensitivity coefficients in both cases. The results of this study emphasize the need to confirm CLUTCH sensitivity results with DP calculations. IFP calculations are generally less efficient but more reliable than CLUTCH calculations. Improvements to the CLUTCH methodology that retain the greater efficiency but address identified difficulties are therefore potentially useful to analysts.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗