Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative limiters”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Long-time simulations for fixed input states on quantum hardware

Publicly accessible quantum computers open the exciting possibility of experimental dynamical quantum simulations. While rapidly improving, current devices have short coherence times, restricting the viable circuit depth. Despite these limitations, we demonstrate long-time, high fidelity simulations on current hardware. Specifically, we simulate an XY-model spin chain on Rigetti and IBM quantum computers, maintaining a fidelity over 0.9 for 150 times longer than is possible using the iterated Trotter method. Our simulations use an algorithm we call fixed state Variational Fast Forwarding (fsVFF). Recent work has shown an approximate diagonalization of a short time evolution unitary allows a fixed-depth simulation. fsVFF substantially reduces the required resources by only diagonalizing the energy subspace spanned by the initial state, rather than over the total Hilbert space. We further demonstrate the viability of fsVFF through large numerical simulations, and provide an analysis of the noise resilience and scaling of simulation errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Realizing the potentials of density functional theory (DFT) and of the materials genome initiative (MGI)

Abstract From 1964 and 1965 to present, the wide spread utilization of an incomplete density functional theory (DFT) has led to mixed results: The second theorem of the theory asserts that the energy functional reaches its minimum if the calculation employs the ground state charge density—without providing a mechanism for finding this density. Calculations purporting to employ DFT have mostly assumed that results obtained with a judiciously selected basis set, following self-consistent iterations, are those of the ground state. The state obtained with a single basis set is a stationary one, among an infinite number of such states, with no proven relation to the actual ground state of the material. Most failures or limitations of the incomplete DFT can be traced to this error. We present results from calculations using the completed DFT. They are in excellent agreement with experiment and portend the realization of the Materials Genome Initiative.

Bagayoko, Diola↗
DOE OSTI2026

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Iterative design of one- and two-dimensional FIR digital filters

The paper describes a new iterative technique for designing FIR (finite duration impulse response) digital filters using a frequency weighted least squares approximation. The technique is as easy to implement (via FFT) and as effective in two dimensions as in one dimension, and there are virtually no limitations on the class of filter frequency spectra approximated. An adaptive adjustment of the frequency weight to achieve other types of design approximation such as Chebyshev type design is discussed.

Suk, M.↗
NASA NTRS1991

The analysis of a rocket tomography measurement of the N2+3914A emission and N2 ionization rates in an auroral arc

Techniques were developed for recovering two-dimensional distributions of auroral volume emission rates from rocket photometer measurements made in a tomographic spin scan mode. These tomographic inversion procedures are based upon an algebraic reconstruction technique (ART) and utilize two different iterative relaxation techniques for solving the problems associated with noise in the observational data. One of the inversion algorithms is based upon a least squares method and the other on a maximum probability approach. The performance of the inversion algorithms, and the limitations of the rocket tomography technique, were critically assessed using various factors such as (1) statistical and non-statistical noise in the observational data, (2) rocket penetration of the auroral form, (3) background sources of emission, (4) smearing due to the photometer field of view, and (5) temporal variations in the auroral form. These tests show that the inversion procedures may be successfully applied to rocket observations made in medium intensity aurora with standard rocket photometer instruments. The inversion procedures have been used to recover two-dimensional distributions of auroral emission rates and ionization rates from an existing set of N2+3914A rocket photometer measurements which were made in a tomographic spin scan mode during the ARIES auroral campaign. The two-dimensional distributions of the 3914A volume emission rates recoverd from the inversion of the rocket data compare very well with the distributions that were inferred from ground-based measurements using triangulation-tomography techniques and the N2 ionization rates derived from the rocket tomography results are in very good agreement with the in situ particle measurements that were made during the flight. Three pre-prints describing the tomographic inversion techniques and the tomographic analysis of the ARIES rocket data are included as appendices.

Mcdade, Ian C.↗
NASA NTRS1991

A geometric nonlinear degenerated shell element using a mixed formulation with independently assumed strain fields

A mixed formulation is chosen to overcome deficiencies of the standard displacement-based shell model. Element development is traced from the incremental variational principle on through to the final set of equilibrium equations. Particular attention is paid to developing specific guidelines for selecting the optimal set of strain parameters. A discussion of constraint index concepts and their predictive capability related to locking is included. Performance characteristics of the elements are assessed in a wide variety of linear and nonlinear plate/shell problems. Despite limiting the study to geometric nonlinear analysis, a substantial amount of additional insight concerning the finite element modeling of thin plate/shell structures is provided. For example, in nonlinear analysis, given the same mesh and load step size, mixed elements converge in fewer iterations than equivalent displacement-based models. It is also demonstrated that, in mixed formulations, lower order elements are preferred. Additionally, meshes used to obtain accurate linear solutions do not necessarily converge to the correct nonlinear solution. Finally, a new form of locking was identified associated with employing elements designed for biaxial bending in uniaxial bending applications.

Graf, Wiley E.↗
NASA NTRS2008

Post-Landing Orion Crew Survival in Warm Ocean Areas: A Case Study in Iterative Environmental Design

The Orion crew module (CM) is being designed to perform survivable land and water landings. There are many issues associated with post-landing crew survival. In general, the most challenging of the realistic Orion landing scenarios from an environmental control standpoint is the off-nominal water landing. Available power and other consumables will be very limited after landing, and it may not be possible to provide full environmental control within the crew cabin for very long after splashdown. Given the bulk and thermal insulation characteristics of the crew-worn pressure suits, landing in a warm tropical ocean area would pose a risk to crew survival from elevated core body temperatures, if for some reason the crewmembers were not able to remove their suits and/or exit the vehicle. This paper summarizes the analyses performed and conclusions reached regarding post-landing crew survival following a water landing, from the standpoint of the crew s core body temperatures.

Rains, George E.↗
NASA NTRSDate not supplied

Root Cause Analysis of the Data Refinement Process – Medical Conditions Capability Resource Tables

The medical system for spaceflight thus far has been designed to support missions in low earth orbit (LEO). Crew capabilities are limited and heavily dependent on the team of medical support staff at Mission Control Center (MCC) to guide diagnosis and management. However, missions to the Moon and Mars will suffer from several constraints that will make this ground support focused approach to care ineffective. In order to update and modify medical system design, NASA has relied on Probabilistic Risk Assessment (PRA) modeling to mitigate medical risk through trade space analysis. Specifically, capability resource tables (CRT’s) were developed to create a dataset of resources required to manage a list of accepted medical conditions significant in exploration spaceflight. With 120 conditions, this dataset contained hundreds of capabilities and thousands of resources with tens of thousands of cells of data. Initially these tables were built in excel for high throughput during development, but ultimately had to be transferred, managed, and modified into the Evidence Library database for modeling purposes. The process of collating and reviewing the Evidence Library revealed numerous errors in the dataset that had to be corrected through iterative changes. Several error types emerged during this process and can be broken into specific classifications defined as “input”, “transcription”, “structural”, “branching”, and “information”. In reviewing these error types through the root cause analysis (RCA) approach, we were able to identify the contributors to these errors which included single data review points, changing product end goals, limited software selection, time constraints and several others. By reviewing and evaluating the underlying causes we can provide possible system improvements that can be implemented for current and future data management in PRA model inputs.

A. Anderson↗
DOE OSTI2024

WEST full tungsten operation with an ITER grade divertor

The mission of WEST (tungsten-W Environment in Steady-state Tokamak) is to explore long pulse operation in a full tungsten (W) environment for preparing next-step fusion devices (ITER and DEMO) with a focus on testing the ITER actively cooled W divertor in tokamak conditions. Following the successful completion of phase 1 (2016-2021), phase 2 started in December 2022 with the lower divertor made entirely of actively cooled ITER-grade tungsten mono-blocks. A boronization prior the first plasma attempt allowed for a smooth startup with the new divertor. Despite the reduced operating window due to tungsten, rapid progress has been made in long pulse operation, resulting in discharges with a pulse length of 100 s and an injected energy of around 300 MJ per discharge. Plasma startup studies were carried out with equatorial boron nitride limiters to compare them with tungsten limiters, while Ion Cyclotron Resonance Heating assisted startup was attempted. High fluence operation in attached regime, which was the main thrust of the first campaigns, already showed the progressive build up of deposits and appearance of dust, impacting the plasma operation as the plasma fluence increased. In total, the cumulated injected energy during the first campaigns reached 43 GJ and the cumulated plasma time exceeded 5 h. Demonstration of controlled X-Point Radiator regime is also reported, opening a promising route for investigating plasma exhaust and plasma-wall interaction issues in more detached regime. This paper summarises the lessons learned from the manufacturing and the first operation of the ITER-grade divertor, describing the progress achieved in optimising operation in a full W environment with a focus on long pulse operation and plasma wall interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Density and temperature profiles after low-Z and high-Z shattered pellet injections on DIII-D

Abstract In this work we utilize the recently upgraded Thomson scattering diagnostic to resolve density and temperature plasma profiles after pure deuterium and mixed neon/deuterium shattered pellet injections (SPIs) on DIII-D. This allows us to study individual components of the staggered scheme proposed for disruption mitigation on ITER, consisting of a low-Z material SPI followed by a delayed high-Z SPI. Obtained spatio-temporal density profiles exhibit very different dynamics after dominantly neon and pure deuterium SPIs. The neon SPI causes a fast radiative plasma collapse in a few milliseconds and results in almost flat density profile once the impurity mixes with the plasma during and after the thermal quench (TQ). The deuterium SPI leads to a disruption delayed by ten and more milliseconds, but very limited core fueling can be observed before the disruption. Even during and after the TQ, the edge deuterium density significantly exceeds the core density. 1D transport modeling suggests that this poor core fueling can be explained by strong outward grad-B-induced drift of the injected deuterium. Preliminary simulations show that larger pellet shards and greater injected quantity can be used to improve the penetration of the low-Z material into the core. These results call for optimization and further evaluation of the staggered SPI on ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Recent DIII-D progress toward validating models of tungsten erosion, re-deposition, and migration for application to next-step fusion devices

Abstract Fundamental mechanisms governing the erosion and prompt re-deposition of tungsten impurities in tokamak divertors are identified and analyzed to inform the lifetime of tungsten plasma-facing components in ITER and other future devices. Various experiments conducted at DIII-D to benchmark predictive models are presented, leveraging the DiMES removable sample exposure probe capability and the Metal Rings Campaign, in which toroidally symmetric rows of tungsten-coated tiles were installed in the DIII-D divertor. In tokamak divertors, the width of the electric sheath is of the order of the main ion Larmor radius, and a vast majority of sputtered tungsten impurities are typically ionized within the sheath. Therefore, W prompt redeposition is mainly governed by the ratio of the characteristic ionization mean-free path of neutral tungsten to the width of the sheath. In-situ monitoring of the prompt redeposition of tungsten impurities in divertors is demonstrated via the use of WII/WI line ratios and the ionizations/photon (S/XB) method in L-mode discharges. Even with this relatively limited set of emission measurements, net erosion measurements were found to be a consistent upper bound to an analytic scaling based on the ratio of the W ionization length, λ iz , and the width of the magnetic sheath rather than the ratio of λ iz and the W + gyro-radius. In the far-scrape-off layer (SOL) of the ITER divertor, however, it is calculated that the measurement of photon emissions associated with the ionization of tungsten impurities up to W 5 + may be required. Finally, W deposition patterns on DiMES collector probes, interpreted via DIVIMP-WallDYN modelling, reveal the key roles of progressive W erosion/re-deposition staps and E × B drifts in regulating long-range high-Z material migration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

A finite element procedure for nonlinear prebuckling and initial postbuckling analysis

A procedure cast in a form appropriate to the finite element method is presented for geometrically nonlinear prebuckling and postbuckling structural analysis, including the identification of snap-through type of buckling. The principal features of this procedure are the use of direct iteration for solution of the nonlinear algebraic equations in the prebuckling range, an interpolation scheme for determination of the initial bifurcation point, a perturbation method in definition of the load-displacement behavior through the postbuckling regime, and extrapolation in determination of the limit point for snap-through buckling. Three numerical examples are presented in illustration of the procedure and in comparison with alternative approaches.

Mau, S. T.↗
NASA NTRS1972

A review of nonstatistical techniques for the estimation of vertical atmospheric structure from remote infrared measurements

Mathematical techniques used in atmospheric profile inversion are reviewed, with the discussion limited to methods that do not require extensive information on the statistical properties of the profiles. The basic problem of the calculation of vertical profiles from noise contaminated infrared spectral measurements is examined. Three specific examples of methods of solution are given: the iterative minimum estimation method, a second linear iterative method, and the algorithm of Chahine. Satellite data are used in an empirical comparison of the methods. Several factors are considered, including the stability of solutions, convergence behavior of the iterative methods, and dependence of the solutions on the first guess, a factor that is particularly important in the treatment of data from other planets.

Conrath, B. J.↗
NASA NTRS1975

An iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions, Addendum

New results and insights concerning a previously published iterative procedure for obtaining maximum-likelihood estimates of the parameters for a mixture of normal distributions were discussed. It was shown that the procedure converges locally to the consistent maximum likelihood estimate as long as a specified parameter is bounded between two limits. Bound values were given to yield optimal local convergence.

Peters, B. C., Jr.↗
NASA NTRS2018

Testing and Maturing a Mass Translating Mechanism for a Deep Space CubeSat

Near Earth Asteroid (NEA) Scout is a deep space satellite set to launch aboard NASA's Exploration Mission 1. The spacecraft fits within a CubeSat standard 6U (about 300 x 200 x 100 mm) and is designed to travel 1 AU over a 2.5 year mission to observe NEA VG 1991. The spacecraft will use an 86 sq.m solar sail to maneuver from lunar orbit to the NEA. One of the critical mechanisms aboard NEA Scout, the Active Mass Translator (AMT), has gone through rigorous design and test cycles since its conception in July of 2015. The AMT is a two-axis translation table required to balance the spacecraft's center of mass (CM) and solar sail center of pressure (CP) while also trimming disturbance torque created by off-nominal sail conditions. The AMT has very limited mass and volume requirements, but is still required to deliver a large translation range-about 160 x 68 mm-at sub mm accuracy and precision. The system is constrained to operate in complete exposure to space with limited power and data budgets for mechanical and thermal needs. The NEA Scout team developed and carried out a rigorous test suite for the prototype and engineering development unit (EDU). These tests uncovered numerous design failures and led to many failure investigations and iteration cycles. A paper was previously presented at the 43rd Aerospace Mechanisms Symposia entitled, "Development of a High Performance, Low Profile Translation Table with Wire Feedthrough for a Deep Space CubeSat". This paper will make note of specific lessons learned: manufacturing philosophy, testing ideologies for high-risk missions, thermal mitigation design for small, motor-driven mechanisms.

Few, Alex↗