Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Incremental Computing”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The DAO Retrospective Data Assimilation System

The fixed-lag Kalman smoother of Cohn et al. (1994), or a version of a more computationally feasible approximation of it developed by Todling et al. (1998), is under implementation at the Data Assimilation Office (DAO) to build a retrospective data assimilation system (RDAS). The initial version of the RDAS uses the physical-space statistical analysis system and a modification of the incremental analysis update (IAU) procedure of Bloom et al. (1996); a follow up version of the RDAS Nill require the use of the adjoint of the DAO general circulation model (developed by Y. Yang and M Navon). The retrospective procedure is designed to produce improved analyses as well as improved assimilated fields consequently providing an improved climate representation through data assimilation.

Todling, Richardo↗
NASA NTRS2006

Using Quasiparticle Poisoning To Detect Photons

According to a proposal, a phenomenon associated with excitation of quasiparticles in certain superconducting quantum devices would be exploited as a means of detecting photons with exquisite sensitivity. The phenomenon could also be exploited to perform medium-resolution spectroscopy. The proposal was inspired by the observation that Coulomb blockade devices upon which some quantum logic gates are based are extremely sensitive to quasiparticles excited above the superconducting gaps in their leads. The presence of quasiparticles in the leads can be easily detected via the charge states. If quasiparticles could be generated in the leads by absorption of photons, then the devices could be used as very sensitive detectors of electromagnetic radiation over the spectral range from x-rays to submillimeter waves. The devices in question are single-Cooper-pair boxes (SCBs), which are mesoscopic superconducting devices developed for quantum computing. An SCB consists of a small superconducting island connected to a reservoir via a small tunnel junction and connected to a voltage source through a gate capacitor. An SCB is an artificial two-level quantum system, the Hamiltonian of which can be controlled by the gate voltage. One measures the expected value of the charge of the eigenvectors of this quantum system by use of a radio-frequency single-electron transistor. A plot of this expected value of charge as a function of gate voltage resembles a staircase that, in the ideal case, consists of steps of height 2 e (where e is the charge of one electron). Experiments have shown that depending on the parameters of the device, quasiparticles in the form of "broken" Cooper pairs present in the reservoir can tunnel to the island, giving rise to steps of 1 e. This effect is sometimes called "poisoning." Simulations have shown that an extremely small average number of quasiparticles can generate a 1-e periodic signal. In a device according to the proposal, this poisoning would be turned to advantage. Depending on the wavelength, an antenna or other component would be used to couple radiation into the reservoir, wherein the absorption of photons would break Cooper pairs, thereby creating quasiparticles that, in turn, would tunnel to the island, creating a 1-e signal. On the basis of conservative estimates of device parameters derived from experimental data and computational simulations that fit the data, it has been estimated that the noise equivalent power of a device according to the proposal could be as low as 6 10(exp -22) W/Hz(exp 1/2). It has also been estimated that the spectroscopic resolution (photon energy divided by increment of photon energy) of such a device in visible light would exceed 100.

Echternach, Pierre↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Faster Algorithm For Computation Of Incompressible Flow

Improved algorithm yields faster numerical solutions of Navier-Stokes equations of steady or unsteady three-dimensional flow of incompressible fluid. In artificial-compressibility method, unsteady flow treated as incompressible in advancing from one time step to next, but at each time step (or in steady state), fluid treated as having variable compressibility enabling propagation of flow field, and subiterations performed in increments of pseudotime until effects of compressibility subside. Directly couples pressure and velocity fields at same time step and converts elliptic incompressible Navier-Stokes equations to hyperbolic form more amenable to numerical integration.

Rogers, S. E.↗
NASA NTRS1988

Numerical solution methods for viscoelastic orthotropic materials

Numerical solution methods for viscoelastic orthotropic materials, specifically fiber reinforced composite materials, are examined. The methods include classical lamination theory using time increments, direction solution of the Volterra Integral, Zienkiewicz's linear Prony series method, and a new method called Nonlinear Differential Equation Method (NDEM) which uses a nonlinear Prony series. The criteria used for comparison of the various methods include the stability of the solution technique, time step size stability, computer solution time length, and computer memory storage. The Volterra Integral allowed the implementation of higher order solution techniques but had difficulties solving singular and weakly singular compliance function. The Zienkiewicz solution technique, which requires the viscoelastic response to be modeled by a Prony series, works well for linear viscoelastic isotropic materials and small time steps. The new method, NDEM, uses a modified Prony series which allows nonlinear stress effects to be included and can be used with orthotropic nonlinear viscoelastic materials. The NDEM technique is shown to be accurate and stable for both linear and nonlinear conditions with minimal computer time.

Gramoll, K. C.↗
NASA NTRS1988

BRYNTRN: A baryon transport computer code, computation procedures and data base

The development is described of an interaction data base and a numerical solution to the transport of baryons through the arbitrary shield material based on a straight ahead approximation of the Boltzmann equation. The code is most accurate for continuous energy boundary values but gives reasonable results for discrete spectra at the boundary with even a relatively coarse energy grid (30 points) and large spatial increments (1 cm in H2O).

Wilson, John W.↗
NASA NTRS1988

Dynamic analysis of nonlinear rotor-housing systems

Nonlinear analysis methods are developed which will enable the reliable prediction of the dynamic behavior of the space shuttle main engine (SSME) turbopumps in the presence of bearing clearances and other local nonlinearities. A computationally efficient convolution method, based on discretized Duhamel and transition matrix integral formulations, is developed for the transient analysis. In the formulation, the coupling forces due to the nonlinearities are treated as external forces acting on the coupled subsystems. Iteration is utilized to determine their magnitudes at each time increment. The method is applied to a nonlinear generic model of the high pressure oxygen turbopump (HPOTP). As compared to the fourth order Runge-Kutta numerical integration methods, the convolution approach proved to be more accurate and more highly efficient. For determining the nonlinear, steady-state periodic responses, an incremental harmonic balance method was also developed. The method was successfully used to determine dominantly harmonic and subharmonic responses fo the HPOTP generic model with bearing clearances. A reduction method similar to the impedance formulation utilized with linear systems is used to reduce the housing-rotor models to their coordinates at the bearing clearances. Recommendations are included for further development of the method, for extending the analysis to aperiodic and chaotic regimes and for conducting critical parameteric studies of the nonlinear response of the current SSME turbopumps.

Noah, Sherif T.↗
NASA NTRS2001

Geometrically Nonlinear Finite Element Analysis of a Composite Space Reflector

Lightweight aerospace structures, such as low areal density composite space reflectors, are highly flexible and may undergo large deflection under applied loading, especially during the launch phase. Accordingly, geometrically nonlinear analysis that takes into account the effect of finite rotation may be needed to determine the deformed shape for a clearance check and the stress and strain state to ensure structural integrity. In this study, deformation of the space reflector is determined under static conditions using a geometrically nonlinear solid shell finite element model. For the solid shell element formulation, the kinematics of deformation is described by six variables that are purely vector components. Because rotational angles are not used, this approach is free of the limitations of small angle increments. This also allows easy connections between substructures and large load increments with respect to the conventional shell formulation using rotational parameters. Geometrically nonlinear analyses were carried out for three cases of static point loads applied at selected points. A chart shows results for a case when the load is applied at the center point of the reflector dish. The computed results capture the nonlinear behavior of the composite reflector as the applied load increases. Also, they are in good agreement with the data obtained by experiments.

Lee, Kee-Joo↗
NASA NTRS2002

Contingency Planning for Planetary Rovers

There has been considerable work in AI on planning under uncertainty. But this work generally assumes an extremely simple model of action that does not consider continuous time and resources. These assumptions are not reasonable for a Mars rover, which must cope with uncertainty about the duration of tasks, the power required, the data storage necessary, along with its position and orientation. In this paper, we outline an approach to generating contingency plans when the sources of uncertainty involve continuous quantities such as time and resources. The approach involves first constructing a "seed" plan, and then incrementally adding contingent branches to this plan in order to improve utility. The challenge is to figure out the best places to insert contingency branches. This requires an estimate of how much utility could be gained by building a contingent branch at any given place in the seed plan. Computing this utility exactly is intractable, but we outline an approximation method that back propagates utility distributions through a graph structure similar to that of a plan graph.

Dearden, Richard↗
NASA NTRS2003

Incremental Contingency Planning

There has been considerable work in AI on planning under uncertainty. However, this work generally assumes an extremely simple model of action that does not consider continuous time and resources. These assumptions are not reasonable for a Mars rover, which must cope with uncertainty about the duration of tasks, the energy required, the data storage necessary, and its current position and orientation. In this paper, we outline an approach to generating contingency plans when the sources of uncertainty involve continuous quantities such as time and resources. The approach involves first constructing a "seed" plan, and then incrementally adding contingent branches to this plan in order to improve utility. The challenge is to figure out the best places to insert contingency branches. This requires an estimate of how much utility could be gained by building a contingent branch at any given place in the seed plan. Computing this utility exactly is intractable, but we outline an approximation method that back propagates utility distributions through a graph structure similar to that of a plan graph.

Dearden, Richard↗
NASA NTRS2011

Sleep Disruption Medical Intervention Forecasting (SDMIF) Module for the Integrated Medical Model

The NASA Integrated Medical Model (IMM) assesses the risk, including likelihood and impact of occurrence, of all credible in-flight medical conditions. Fatigue due to sleep disruption is a condition that could lead to operational errors, potentially resulting in loss of mission or crew. Pharmacological consumables are mitigation strategies used to manage the risks associated with sleep deficits. The likelihood of medical intervention due to sleep disruption was estimated with a well validated sleep model and a Monte Carlo computer simulation in an effort to optimize the quantity of consumables. METHODS: The key components of the model are the mission parameter program, the calculation of sleep intensity and the diagnosis and decision module. The mission parameter program was used to create simulated daily sleep/wake schedules for an ISS increment. The hypothetical schedules included critical events such as dockings and extravehicular activities and included actual sleep time and sleep quality. The schedules were used as inputs to the Sleep, Activity, Fatigue and Task Effectiveness (SAFTE) Model (IBR Inc., Baltimore MD), which calculated sleep intensity. Sleep data from an ISS study was used to relate calculated sleep intensity to the probability of sleep medication use, using a generalized linear model for binomial regression. A human yes/no decision process using a binomial random number was also factored into sleep medication use probability. RESULTS: These probability calculations were repeated 5000 times resulting in an estimate of the most likely amount of sleep aids used during an ISS mission and a 95% confidence interval. CONCLUSIONS: These results were transferred to the parent IMM for further weighting and integration with other medical conditions, to help inform operational decisions. This model is a potential planning tool for ensuring adequate sleep during sleep disrupted periods of a mission.

Lewandowski, Beth↗
DOE OSTI2021

Streaming Generalized Canonical Polyadic Tensor Decompositions

In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

First simultaneous observation of co- and counter-current fast-ion losses in the ASDEX Upgrade tokamak

In ITER and future fusion power plants, the source of the fusion born alpha particles is almost isotropic in pitch angle, thus having co- and counter-current populations. For trapped ions, the co-current side of the orbit corresponds to its outer leg, while the counter-current side corresponds to the inner leg. Understanding the mechanisms responsible for the fast-ion losses (FILs) is critical for future magnetically confined fusion power plants. To further study the interplay of fast ions with plasma instabilities, a double pinhole collimator has been developed for a Fast-Ion Loss Detector (FILD) in the ASDEX Upgrade tokamak (AUG). This new FILD opens the operational window to simultaneous measurements of the co- and counter-current ion velocity-space. In this paper, the first results for the AUG double collimator FILD detector are shown. The commissioning of this new probe is carried out in H-mode plasmas with an on-axis magnetic field $B_0 = -2.5$ T, and a plasma current $I_{\textrm{p}} = 0.7$ MA. Simultaneous co- and counter-current FILs have been measured. Both have shown a similar dependence on Ion Cyclotron Resonance Heating (ICRH) power, where the main difference is the intensity of the losses, with the co-losses being an order of magnitude larger. Toroidal Alfvén eigenmode-coherent ICRH-only losses have been identified for the co-current ions. Additionally, the presence of Edge Localized Modes during the discharge were shown to increment Neutral Beam Injection prompt losses, while partially mitigating ICRH-driven losses on both co- and counter- sides of the velocity-space. Finally, a very trapped and high gyroradius losses, with an unclear origin, have been measured in the co- and counter-current velocity-space. The computed ion trajectories show that these ions remain permanently near the vessel wall, suggesting that they are accelerated within the scrape-off layer.

ELM↗
NASA NTRS1971

Gravity flow of powder in a lunar environment. Part 2: Analysis of flow initiation

A small displacement-small strain finite element technique utilizing the constant strain triangle and incremental constitutive equations for elasticplastic (media nonhardening and obeying a Coulomb yield condition) was applied to the analysis of gravity flow initiation. This was done in a V-shaped hopper containing a powder under lunar environmental conditions. Three methods of loading were examined. Of the three, the method of computing the initial state of stress in a filled hopper prior to drawdown, by adding material to the hopper layer by layer, was the best. Results of the analysis of a typical hopper problem show that the initial state of stress, the elastic moduli, and the strength parameters have an important influence on material response subsequent to the opening of the hopper outlet.

Pariseau, W. G.↗
DOE OSTI2025

Micropolar Elastoplasticity Using a Fast Fourier Transform‐Based Solver

ABSTRACT This work presents a micromechanical spectral formulation for obtaining the full‐field and homogenized response of elastoplastic micropolar composites. A closed‐form radial‐return mapping is derived from thermodynamics‐based micropolar elastoplastic constitutive equations to determine the increment of plastic strain necessary to return the generalized stress state to the yield surface, and the algorithm implementation is verified using the method of numerically manufactured solutions. Then, size‐dependent material response and micro‐plasticity are shown as features that may be efficiently simulated in this micropolar elastoplastic framework. The computational efficiency of the formulation enables the generation of large datasets in reasonable computing times.

42 ENGINEERING↗
NASA NTRS1982

Implicit numerical integration for periodic solutions of autonomous nonlinear systems

A change of variables that stabilizes numerical computations for periodic solutions of autonomous systems is derived. Computation of the period is decoupled from the rest of the problem for conservative systems of any order and for any second-order system. Numerical results are included for a second-order conservative system under a suddenly applied constant load. Near the critical load for the system, a small increment in load amplitude results in a large increase in amplitude of the response.

Thurston, G. A.↗
NASA NTRS1984

Simulation Experiments on the Relative Accuracy of Inferred Atmospheric States from Idealized Wind and Temperature Profiling Systems

The First GARP Global Experiment (FGGE) led to the planning, design and implementation of a global observing system. Technologies as coherent CO2 LIDAR systems and highly sensitive IR detectors offer the potential of more accurate global spaceborne observation systems. The advent of super computers allows general circulation modeling and more sophisticated data analysis schemes. The consideration of advanced spaceborne systems better suited to meet the new emerging requirements are suggested. The Global Weather Experiment provides us with a new baseline to assess data accuracy and forecasting capabilities which is used for study of the incremental performance.

Halem, M.↗
NASA NTRS2018

An Algorithm for Trajectory Generation in Redundant Manipulators with Joint Transmission Accommodation

Trajectory generation for manipulators involves generating incremental updates of joint variables to achieve a desired end–effector motion. The Jacobian matrix maps incremental joint motion to incremental end–effector motion in a linear fashion and is typically used in manipulator trajectory generation algorithms. In the case of a redundant manipulator, the Jacobian matrix is not square or invertible and therefore algorithms based on pseudoinverses and their variations are commonly used for trajectory generation. These methods either are computationally not efficient or do not utilize all the joints in motion generation and therefore do not completely exploit the redundancy of the manipulator. The method presented in this paper is a simple method that maximizes transmission of all joint variables onto a desired end–effector motion trajectory. The method is based on aligning the null–space of an augmented Jacobian matrix with the path of the desired end–effector motion, from which a linear combination of joints that projects fully onto the desired end–effector trajectory is obtained. In this manner, all joints of the redundant manipulator are used to generate the end–effector trajectory accommodating the ability of each joint in terms of its motion transmission.

Ravani, Bahram↗