Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initial conditions problem”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Solving time-dependent two-dimensional eddy current problems

Transient eddy current calculations are presented for an EM wave-scattering and field-penetrating case in which a two-dimensional transverse magnetic field is incident on a good (i.e., not perfect) and infinitely long conductor. The problem thus posed is of initial boundary-value interface type, where the boundary of the conductor constitutes the interface. A potential function is used for time-domain modeling of the situation, and finite difference-time domain techniques are used to march the potential function explicitly in time. Attention is given to the case of LF radiation conditions.

Lee, Min Eig↗
NASA NTRS1999

Reinforcement Learning in a Nonstationary Environment: The El Farol Problem

This paper examines the performance of simple learning rules in a complex adaptive system based on a coordination problem modeled on the El Farol problem. The key features of the El Farol problem are that it typically involves a medium number of agents and that agents' pay-off functions have a discontinuous response to increased congestion. First we consider a single adaptive agent facing a stationary environment. We demonstrate that the simple learning rules proposed by Roth and Er'ev can be extremely sensitive to small changes in the initial conditions and that events early in a simulation can affect the performance of the rule over a relatively long time horizon. In contrast, a reinforcement learning rule based on standard practice in the computer science literature converges rapidly and robustly. The situation is reversed when multiple adaptive agents interact: the RE algorithms often converge rapidly to a stable average aggregate attendance despite the slow and erratic behavior of individual learners, while the CS based learners frequently over-attend in the early and intermediate terms. The symmetric mixed strategy equilibria is unstable: all three learning rules ultimately tend towards pure strategies or stabilize in the medium term at non-equilibrium probabilities of attendance. The brittleness of the algorithms in different contexts emphasize the importance of thorough and thoughtful examination of simulation-based results.

Bell, Ann Maria↗
DOE OSTI2024

Implicit–explicit multirate infinitesimal stage-restart methods

Implicit–Explicit (IMEX) methods are flexible numerical time integration methods which solve an initial-value problem (IVP) that is split into stiff and nonstiff processes with the goal of lower computational costs than a purely implicit or explicit approach. A complementary form of flexible IVP solvers are multirate infinitesimal methods for problems split into fast- and slow-changing dynamics, that solve a multirate IVP by evolving a sequence of “fast” IVPs using any suitably accurate algorithm. This article introduces a new class of high-order implicit–explicit multirate methods that are designed for multirate IVPs in which the slow-changing dynamics are further split in an IMEX fashion. This new class, which we call implicit–explicit multirate infinitesimal stage-restart (IMEX-MRI-SR), both improves upon the previous implicit–explicit multirate infinitesimal generalized-structure additive Runge Kutta (IMEX-MRI-GARK) methods by allowing for far easier creation of new embedded methods, and extends multirate exponential Runge Kutta (MERK) methods by allowing the fast-changing dynamics to be nonlinear and the methods to be implicit. We leverage GARK theory to derive conditions for orders of accuracy up to four, and we provide second- and third-order accurate example methods, which are the first known embedded MRI methods with IMEX structure. We then perform numerical simulations demonstrating convergence rates and computational performance in both fixed-step and adaptive-step settings.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
NASA NTRS1994

Autogenic-feedback training improves pilot performance during emergency flying conditions

Studies have shown that autonomous mode behavior is one cause of aircraft fatalities due to pilot error. In such cases, the pilot is in a high state of psychological and physiological arousal and tends to focus on one problem, while ignoring more critical information. This study examined the effect of training in physiological self-recognition and regulation, as a means of improving crew cockpit performance. Seventeen pilots were assigned to the treatment and control groups matched for accumulated flight hours. The treatment group comprised three pilots of HC-130 Hercules aircraft and four HH-65 Dolphin helicopter pilots; the control group comprised three pilots of HC-130's and six Dolphin helicopter pilots. During an initial flight, physiological data were recorded for each crew member and individual crew performance was rated by an instructor pilot. Eight crewmembers were then taught to regulate their own physiological response levels using Autogenic-Feedback Training (AFT). The remaining subjects received no training. During a second flight, treatment subjects showed significant improvement in performance, while controls did not improve. The results indicate that AFT management of high states of physiological arousal may improve pilot performance during emergency flying conditions.

Kellar, Michael A.↗
NASA NTRS2004

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗
NASA NTRS2001

Autogenic Feedback Training Exercise and pilot performance: enhanced functioning under search-and-rescue flying conditions

Studies have shown that autonomous mode behavior is one cause of aircraft fatalities due to pilot error. In such cases, the pilot is in a high state of psychological and physiological arousal and tends to focus on one problem, while ignoring more critical information. This study examined the effect of training in physiological self-recognition and regulation, as a means of improving crew cockpit performance. Seventeen pilots were assigned to the treatment and control groups matched for accumulated flight hours. The treatment group contained 4 pilots from HC-130 Hercules aircraft and 4 HH-65 Dolphin helicopter pilots; the control group contained 3 pilots of HC-130s and 6 helicopter pilots. During an initial flight, physiological data were recorded on each crewmember and an instructor pilot rated individual crew performance. Eight crewmembers were then taught to regulate their own physiological response levels using Autogenic-Feedback Training Exercise (AFTE). The remaining participants received no training. During a second flight, treatment participants showed significant improvement in performance (rated by the same instructor pilot as in pretests) while controls did not improve. The results indicate that AFTE management of high states of physiological arousal may improve pilot performance during emergency flying conditions.

NASA Discipline Neuroscience↗
DOE OSTI2026

WellPINN: Accurate Well Representation for Transient Fluid Pressure Diffusion in Subsurface Reservoirs With Physics‐Informed Neural Networks

Accurate representation of pumping wells is essential for reliable reservoir characterization and simulation of operational scenarios in subsurface flow models. Physics-informed neural networks (PINNs) are emerging as a promising alternative to numerical models for reservoir modeling, offering seamless integration of monitoring data and governing physical equations. However, existing PINN-based studies face major challenges in capturing fluid pressure near wells when using a source/sink term, particularly during the early stages after pumping begins. We address this problem by introducing WellPINN, a workflow in which an initially trained PINN infers fluid pressure across the entire reservoir domain using a large equivalent well radius. This initial PINN solution is then locally refined around the well by a set of subdomain PINNs that are trained for smaller equivalent well radii. Continuity across these subdomain interfaces as well as at the initial condition is ensured by hard-constraining each PINN on its subdomain boundary. Our results demonstrate WellPINN as the first workflow of its kind to focus on accurate inference of fluid pressure from pumping rates throughout the entire injection period, significantly advancing the potential of PINNs for inverse modeling and operational scenario simulations. All data and code for this paper are openly available at https://doi.org/10.20350/DIGITALCSIC/17260.

58 GEOSCIENCES↗
DOE OSTI2021

Learning Forecasts of Rare Stratospheric Transitions from Short Simulations

Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology.

Finkel, Justin↗
DOE OSTI2021

Benchmarking of Probabilistic Fracture Mechanics Models in Grizzly

In 2020, as part of contract/task order 31310019N0006/31310020F0060, the Nuclear Regulatory Commission tasked Idaho National Laboratory (INL) with benchmarking the MOOSE/GRIZZLY code against the NRC’s Fracture Analysis of Vessels - Oak Ridge (FAVOR), version 16.1 code. Both codes can predict large light-water reactor vessel integrity deterministically and probabilistically. In fact, such predictions are the primary purpose of the FAVOR code, and are a main focus of the GRIZZLY code, although GRIZZLY is capable of modeling a wider range of geometries and problems. The vessel integrity modeling capabilities of the GRIZZLY code are still under development and being gradually expanded, and the capability to model vessel fracture by crack growth has not yet been integrated into GRIZZLY (although it does exist in pre-release versions). Thus, the present benchmark study focused on comparing predictions for the conditional probability of crack growth initiation (CPI) between GRIZZLY and FAVOR v16.1. The present benchmarking study expands on a previous study completed in 2020, which focused on single plate regions, various flaw types, and full vessel studies. The previous study was limited in scope and revealed some relatively large discrepancies, particularly in the case of weld regions. In the present study, detailed investigations into the sources of discrepancies between GRIZZLY and FAVOR were performed, with a focus on weld regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1981

The stability of solar coronal loops with realistic photospheric boundary conditions

Finite-length effects are the primary influence for the MHD stability of coronal loops. This paper gives a new physical description of the photospheric boundary conditions, and the first complete and uniformly convergent initial-perturbation set. A general energy-principle (or variational) stability analysis is then developed. The results of this calculation show that large classes of cylindrical-symmetry perturbations, those with separable contributions to the energy integral, are completely stable. Comparisons are also made with other more restrictive formulations of this problem.

Van Hoven, G.↗
NASA NTRS1982

The P-POD Project

The programmable pilot oriented display (P-POD), a multiprocessor based flight information processing and display unit, is discussed. Communication protocols between the three Z-80 processors inside P-POD are reviewed. The interface between the video processor, P2, and the main processor, P1, was tested. Hardware and software were revised. Data flow between P1 and P2 is illustrated. While running diagnostic routines to exercise the P1/P2 interface at maximum speed an intermittent failure was observed. A check of possible noise sources or race conditions in the hardware did not reveal the source of this failure. Possible software timing problems were also fully investigated. To eliminate the effects of this type of failure a performance monitor routine was used to detect these occasional failures and correct them. The data transmission cycle between P1/P2 is either allowed to continue normal operation or is re-initialized and restarted when an error is detected.

Littlefield, J. A.↗
NASA NTRS1989

A finite element formulation for the large deflection random response of thermally buckled beams

The effects of temperature and acoustic loading are included in a theoretical finite element large deflection formulation for thin, isotropic beams. Thermal loads are applied as steady-state temperature distributions, and acoustic loads are taken to be ergodic and Gaussian with zero mean and uniform magnitude and phase along the length of the beam. Material properties are considered presently to be independent of temperature. Also, inplane and rotary inertia terms are assumed to be negligible, and all inplane edge conditions are taken to be immovable. For the random response analysis, both auto- and cross-correlation terms are included. The nature of the loads leads to the solution of two separate problems. First, the problem of thermal postbuckling is solved to determine the deflections and stresses due to the thermal load only. These deflections and stresses are then used as initial deflections and stresses for the random vibration analysis. Root-mean-square (RMS) maximum deflections and strains are obtained and compared with previous classical equivalent linearization results.

Locke, James↗
DOE OSTI2022

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Identifiability and characterization of transmon qutrits through Bayesian experimental design

Robust control of a quantum system is essential to utilize the current noisy quantum hardware to its full potential, such as quantum algorithms. To achieve such a goal, a systematic search for an optimal control for any given experiment is essential. The design of optimal control pulses requires accurate numerical models and, therefore, accurate characterization of the system parameters. We present an online Bayesian approach for quantum characterization of qutrit systems, which automatically and systematically identifies optimal experiments that provide maximum information on the system parameters, thereby greatly reducing the number of experiments that need to be performed on the quantum testbed. Unlike most characterization protocols that provide point-estimates of the parameters, the proposed approach is able to estimate their probability distribution. The applicability of the Bayesian experimental design technique was demonstrated on test problems, where each experiment was defined by a parameterized control pulse. In addition to this, we also present an approach for iterative pulse extension, which is robust under uncertainties in transition frequencies and coherence times, and shot noise, despite being initialized with wide uninformative priors. Furthermore, we provide a mathematical proof of the theoretical identifiability of the model parameters and present conditions on the quantum state under which the parameters are identifiable. The proof and conditions for identifiability are presented for both closed and open quantum systems using the Schrödinger equation and the Lindblad master equation, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Learning the boundary-to-domain mapping using Lifting Product Fourier Neural Operators for partial differential equations

Neural operators such as the Fourier Neural Operator (FNO) have been shown to provide resolution-independent deep learning models that can learn mappings between function spaces. For example, an initial condition can be mapped to the solution of a partial differential equation (PDE) at a future time-step using a neural operator. Despite the popularity of neural operators, their use to predict solution functions over a domain given only data over the boundary (such as a spatially varying Dirichlet boundary condition) remains unexplored. In this paper, we refer to such problems as boundary-to-domain problems; they have a wide range of applications in areas such as fluid mechanics, solid mechanics, heat transfer etc. We present a novel FNO-based architecture, named Lifting Product FNO (or LP-FNO) which can map arbitrary boundary functions defined on the lower-dimensional boundary to a solution in the entire domain. Specifically, two FNOs defined on the lower-dimensional boundary are lifted into the higher dimensional domain using our proposed lifting product layer. We demonstrate the efficacy and resolution independence of the proposed LP-FNO for the 2D Poisson equation.

Kashi, Aditya↗
NASA NTRS1988

Rocketdyne LOX bearing tester program

The cause, or causes, for the Space Shuttle Main Engine ball wear were unknown, however, several mechanisms were suspected. Two testers were designed and built for operation in liquid oxygen to empirically gain insight into the problems and iterate solutions in a timely and cost efficient manner independent of engine testing. Schedules and test plans were developed that defined a test matrix consisting of parametric variations of loading, cooling or vapor margin, cage lubrication, material, and geometry studies. Initial test results indicated that the low pressure pump thrust bearing surface distress is a function of high axial load. Initial high pressure turbopump bearing tests produced the wear phenomenon observed in the turbopump and identified an inadequate vapor margin problem and a coolant flowrate sensitivity issue. These tests provided calibration data of analytical model predictions to give high confidence in the positive impact of future turbopump design modification for flight. Various modifications will be evaluated in these testers, since similar turbopump conditions can be produced and the benefit of the modification will be quantified in measured wear life comparisons.

Keba, J. E.↗
NASA NTRS1992

Computational aspects of helicopter trim analysis and damping levels from Floquet theory

Helicopter trim settings of periodic initial state and control inputs are investigated for convergence of Newton iteration in computing the settings sequentially and in parallel. The trim analysis uses a shooting method and a weak version of two temporal finite element methods with displacement formulation and with mixed formulation of displacements and momenta. These three methods broadly represent two main approaches of trim analysis: adaptation of initial-value and finite element boundary-value codes to periodic boundary conditions, particularly for unstable and marginally stable systems. In each method, both the sequential and in-parallel schemes are used and the resulting nonlinear algebraic equations are solved by damped Newton iteration with an optimally selected damping parameter. The impact of damped Newton iteration, including earlier-observed divergence problems in trim analysis, is demonstrated by the maximum condition number of the Jacobian matrices of the iterative scheme and by virtual elimination of divergence. The advantages of the in-parallel scheme over the conventional sequential scheme are also demonstrated.

Gaonkar, Gopal H.↗