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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A note on higher-order and nonlinear limiting approaches for continuously bounds-preserving discontinuous Galerkin methods

In Dzanic (2024), a limiting approach for high-order discontinuous Galerkin schemes was introduced which allowed for imposing constraints on the solution continuously (i.e., everywhere within the element). While exact for linear constraint functionals, this approach only imposed a sufficient (but not the minimum necessary) amount of limiting for nonlinear constraint functionals. This short note shows how this limiting approach can be extended to allow exactness for general nonlinear quasiconcave constraint functionals through a nonlinear limiting procedure, reducing unnecessary numerical dissipation. Finally, some examples are shown for nonlinear pressure and entropy constraints in the compressible gas dynamics equations, where both analytic and iterative approaches are used.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

2009 Space Shuttle Probabilistic Risk Assessment Overview

Loss of a Space Shuttle during flight has severe consequences, including loss of a significant national asset; loss of national confidence and pride; and, most importantly, loss of human life. The Shuttle Probabilistic Risk Assessment (SPRA) is used to identify risk contributors and their significance; thus, assisting management in determining how to reduce risk. In 2006, an overview of the SPRA Iteration 2.1 was presented at PSAM 8 [1]. Like all successful PRAs, the SPRA is a living PRA and has undergone revisions since PSAM 8. The latest revision to the SPRA is Iteration 3. 1, and it will not be the last as the Shuttle program progresses and more is learned. This paper discusses the SPRA scope, overall methodology, and results, as well as provides risk insights. The scope, assumptions, uncertainties, and limitations of this assessment provide risk-informed perspective to aid management s decision-making process. In addition, this paper compares the Iteration 3.1 analysis and results to the Iteration 2.1 analysis and results presented at PSAM 8.

Hamlin, Teri L.↗
DOE OSTI2022

A new discovery of edge localized modes suppression using ICRH

Here, the high-confinement mode (H-mode) is very important for the fusion performance of ITER and future fusion reactors. The H-mode pedestal is a region featuring strongly reduced turbulence and transport just inside the limiting flux surface with strong plasma gradients, which drive edge localized modes (ELMs) in tokamaks. The ELMs, however, are quasiperiodic relaxations of the pedestal, resulting in a series of hot plasma eruptions similar to magnetospheric substorms, solar and stellar flares, which could potentially damage the ITER divertor plates and first walls. To mitigate and suppress the ELMs in H-mode plasmas, the fusion community has tried to develop effective actuators over the past two decades. Two separate excellent inventions have been awarded to recognize their achievements. The 2018 APS excellence award in Plasma Physics Research was given to Drs. Todd E. Evans, Max E. Fenstermacher, and Richard Alan Moyer for their first experimental demonstration of the stabilization of ELMs in high-confinement diverted discharges by application of very small edge-resonant magnetic perturbations, which led to the adoption of suppression coils in the ITER design. The 2018 Fusion Technology Award from the IEEE’s Nuclear and Plasma Science Society (NPSS) was given to Dr. Larry Baylor for his work designing the fueling, pumping and disruption mitigation system for the US ITER Project and for the implementation of pellet ELM mitigation for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗
DOE OSTI2020

Identification of a network of nonlinear interactions as a mechanism triggering the onset of edge localized modes

Edge localized modes (ELMs) remain a critical issue for the lifetime of in-vessel components of fusion reactors such as ITER. We study ELMs triggered when the pedestal parameters are below the peeling-ballooning limit. These ELMs occur while the pedestal recovers after an ELM crash. We identify a mechanism by which a non-linear perturbation involving many active triads—a network of non-linear interactions – leads to the triggering of these ELMs. Each active triad leads to the non-linear transfer of energy between three waves, such that the network of non-linear interactions lead to a transfer of energy between many waves. The triggering of this non-linear perturbation systematically coincides with a pedestal perturbation induced by the neutral beam injection beams and localized near the magnetic surface of safety factor q = 5. Among the non-linear interactions detected during the ELM by a bi-coherence analysis, the most active triad involves a slow frequency that is qualitatively consistent with a geodesic acoustic mode located at the same location, i.e., near the q = 5 surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quantum Electromagnetic Transients Program

This work devises a quantum electromagnetic transients program (QEMTP) which is the first attempt to tackle the computational challenges in solving EMTP through quantum computing. The main contributions lie in: (1) A quantum-enabled EMTP formulation with Dommel's model encoded in quantum states; 2) An Harrow-Hassidim-Lloyd (HHL)-based QEMTP to solve the discrete-time nodal equations; 3) An iteration-based HHL revision for mitigating temporal errors to achieve high accuracy of QEMTP with limited quantum resources. Case studies verify the correctness and efficacy of QEMTP in simulating both fast and slow dynamics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1980

Spacecraft loads methodology - Transient vs. shock spectra method

The methodology for the establishment of spacecraft loads is strongly influenced by project constraints which include the cost, schedule and allowable weight. The most rigorous approach is the transient loads analysis which requires a composite mathematical model of the spacecraft and launch vehicle. The structural member loads for the entire composite structure are computed by applying the forcing functions, which represent various dynamic environments during the mission, to the composite model. Although this method ideally leads to a lightweight design, it is costly and time consuming due to complex interfaces involving many organizations. To reduce complexity and cost a shock spectra method has been used to design spacecraft structures. This method utilizes envelopes of shock spectra of launch vehicle accelerations obtained from analysis and/or flight measurements. Since only limited information on the launch vehicle model is involved in this process the design loads iteration cycle can be rapidly performed within the payload organization. In the present paper, these two methods will be evaluated by comparing the loads for several spacecraft. Flight measured loads will also be used in the evaluation.

Chen, J. C.↗
NASA NTRS2002

Investigation of Fundamental Modeling and Thermal Performance Issues for a Metallic Thermal Protection System Design

A study was performed to develop an understanding of the key factors that govern the performance of metallic thermal protection systems for reusable launch vehicles. A current advanced metallic thermal protection system (TPS) concept was systematically analyzed to discover the most important factors governing the thermal performance of metallic TPS. A large number of relevant factors that influence the thermal analysis and thermal performance of metallic TPS were identified and quantified. Detailed finite element models were developed for predicting the thermal performance of design variations of the advanced metallic TPS concept mounted on a simple, unstiffened structure. The computational models were also used, in an automated iterative procedure, for sizing the metallic TPS to maintain the structure below a specified temperature limit. A statistical sensitivity analysis method, based on orthogonal matrix techniques used in robust design, was used to quantify and rank the relative importance of the various modeling and design factors considered in this study. Results of the study indicate that radiation, even in small gaps between panels, can reduce significantly the thermal performance of metallic TPS, so that gaps should be eliminated by design if possible. Thermal performance was also shown to be sensitive to several analytical assumptions that should be chosen carefully. One of the factors that was found to have the greatest effect on thermal performance is the heat capacity of the underlying structure. Therefore the structure and TPS should be designed concurrently.

Blosser, Max L.↗
DOE OSTI2020

Fuel Property-Informed Process Design for the Direct Catalytic Conversion of Cellulosics

The direct catalytic conversion of cellulosics (DC3) with supercritical methanol presents a promising pathway to completely solubilize woody biomass and produce >60% mass yield of C2-C6 oxygenates for light duty vehicle applications. Although the oxygenate composition can be tuned by modifying catalyst and process conditions, limited efforts to date have examined DC3 fuel properties for light duty fuel applications to iteratively inform conversion process design. In this work, we designed and evaluated multi-component surrogate bioblendstocks that represent major DC3 light oxygenates and refined catalyst formulations and distillation parameters to improve the resulting light duty fuel quality. Surrogate design for this novel fuel production pathway was based on a preliminary product slate that was systematically adjusted to determine impacts on fuel properties of interest including heating value, oxidative stability, and heat of vaporization, among others. The resulting predictions and measurements of these surrogates were then used to inform improvement of conversion process, and specifically, reformulation of the catalyst which resulted in production of a higher energy density fuel product primarily comprised of saturated alcohols. Fuel testing outcomes also informed downstream oxygenate separations requirements. This resulted in the production of a promising DC3 light-duty fuel blendstock with >15% greater energy density relative to ethanol, while maintaining high octane number and octane sensitivity.

BIOMASS FUELS↗
NASA NTRS1984

Bearing tester data compilation, analysis and reporting and bearing math modeling, volume 1

Thermal and mechanical models of high speed angular contact ball bearings operating in LOX and LN2 were developed and verified with limited test data in an effort to further understand the parameters that determine or effect the SSME turbopump bearing operational characteristics and service life. The SHABERTH bearing analysis program which was adapted to evaluate shaft bearing systems in cryogenics is not capable of accommodating varying thermal properties and two phase flow. A bearing model with this capability was developed using the SINDA thermal analyzer. Iteration between the SHABERTH and the SINDA models enable the establishment of preliminary bounds for stable operation in LN2. These limits were established in terms of fluid flow, fluid inlet temperature, and axial load for a shaft speed of 30,000 RPM.

Marshall, D. D.↗
NASA NTRS1976

A computer program to calculate the longitudinal aerodynamic characteristics of wing-flap configurations with externally blown flaps

A vortex lattice lifting-surface method is used to model the wing and multiple flaps. Each lifting surface may be of arbitrary planform having camber and twist, and the multiple-slotted trailing-edge flap system may consist of up to ten flaps with different spans and deflection angles. The engine wakes model consists of a series of closely spaced vortex rings with circular or elliptic cross sections. The rings are normal to a wake centerline which is free to move vertically and laterally to accommodate the local flow field beneath the wing and flaps. The two potential flow models are used in an iterative fashion to calculate the wing-flap loading distribution including the influence of the waves from up to two turbofan engines on the semispan. The method is limited to the condition where the flow and geometry of the configurations are symmetric about the vertical plane containing the wing root chord. The calculation procedure starts with arbitrarily positioned wake centerlines and the iterative calculation continues until the total configuration loading converges within a prescribed tolerance. Program results include total configuration forces and moments, individual lifting-surface load distributions, including pressure distributions, individual flap hinge moments, and flow field calculation at arbitrary field points.

Mendenhall, M. R.↗
DOE OSTI2023

Accurate real space iterative reconstruction (RESIRE) algorithm for tomography

Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2009

Decentralized Control of Sound Radiation Using Iterative Loop Recovery

A decentralized model-based control strategy is designed to reduce low-frequency sound radiation from periodically stiffened panels. While decentralized control systems tend to be scalable, performance can be limited due to modeling error introduced by the unmodeled interaction between neighboring control units. Since bounds on modeling error are not known in advance, it is difficult to ensure the decentralized control system will be robust without making the controller overly conservative. Therefore an iterative approach is suggested, which utilizes frequency-shaped loop recovery. The approach accounts for modeling error introduced by neighboring control loops, requires no communication between subsystems, and is relatively simple. The control strategy is evaluated numerically using a model of a stiffened aluminum panel that is representative of the sidewall of an aircraft. Simulations demonstrate that the iterative approach can achieve significant reductions in radiated sound power from the stiffened panel without destabilizing neighboring control units.

Schiller, Noah H.↗
DOE OSTI2026

Beyond Contrast Transfer: Spectral SNR as a Finite-Dose Metric for STEM Phase Retrieval

The contrast transfer function (CTF) is widely used to evaluate phase retrieval methods in scanning transmission electron microscopy (STEM), including center-of-mass imaging, parallax imaging, direct ptychography, and iterative ptychography. However, the CTF reflects only the maximum usable signal, neglecting the effects of finite electron fluence and the Poisson-limited nature of detection. As a result, it can significantly overestimate practical performance, especially in low-dose regimes. Here, we employ the spectral signal-to-noise ratio (SSNR), as a finite-dose statistical framework to evaluate the recoverable signal as a function of spatial frequency. Using numerical reconstructions of white-noise objects, we show that center-of-mass, parallax, and direct ptychography exhibit dose-independent SSNRs, with close-form analytic expressions. In contrast, iterative ptychography exhibits a surprising dose dependence: at low fluence, its SSNR converges to that of direct ptychography; at high fluence, it saturates at a value consistent with the maximum detective quantum efficiency predicted by recent quantum Fisher information bounds. The results highlight the limitations of CTF-based evaluation and motivate SSNR as a more accurate, finite-dose metric for assessing STEM phase retrieval methods.

STEM phase retrieval↗
DOE OSTI2021

Inter-ELM pedestal evolution and the role of edge fluctuations in the C-Mod and DIII-D tokamaks

Edge localized modes (ELMs) lead to a cyclical behaviour of the H-mode pedestal density, temperature and pressure. Substantial evidence exists associating the ELM crash with violation of ideal MHD stability, via current driven kink/peeling modes and/or pressure gradient driven ballooning modes. Recovery from ELM crashes and buildup to the next ELM crash is an active and critical area of edge physics research for ITER projections. In one set of theories, the pedestal pressure profile in H-mode discharges is predicted to be limited by micro-instabilities. One such micro-instability is the kinetic ballooning mode (KBM), which is hypothesized to provide a ‘soft’ limit that regulates the edge transport, thereby restricting the local pressure gradient. In this picture, e.g. in the EPED model, the pedestal width expands slowly until the ideal MHD stability limit is reached. Experiments were recently performed on the C-Mod and DIII-D devices to search for instabilities correlated with the pedestal evolution between ELMs. The results show correlations between the onset of quasi-coherent fluctuations between type I ELMs and a critical temperature gradient. Linear gyrokinetic calculations have shown that these fluctuations have characteristics similar to those expected for KBMs. These results provide additional data toward validation of the EPED model, increasing confidence in ITER projections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Adapting CLUTCH methodology to multigroup TSUNAMI-3D for eigenvalue sensitivity calculations

The sensitivity of the eigenvalue to uncertainties in nuclear data and its evaluation are important for nuclear criticality safety. TSUNAMI-3D sequences within the SCALE code system offer several options to the user community for calculating eigenvalue sensitivity coefficients with multigroup (MG) and continuous energy (CE) 3D transport capabilities. TSUNAMI-3D sequences implement the adjoint-based perturbation theory with MG KENO code, the Contributon Linked eigenvalue sensitivity/Uncertainty estimation via Track length importance CHaracterization (CLUTCH) method with CE KENO code, and the Iterated Fission Probability (IFP) method with CE KENO and Shift codes. Each method has benefits and limitations depending on the problem that is run. The work presented here aims to adapt the CLUTCH method, which enables the Contributon method's mesh-free, memory-efficient approach for calculating adjoint-weighted tallies for sensitivity calculations, to the MG TSUNAMI-3D sequence. This application would eliminate the explicit adjoint KENO calculation, as well as the memory-consuming mesh flux moment tallies required by the conventional MG TSUNAMI-3D. Smaller memory footprints in the CLUTCH methodology and relatively shorter runtimes in MG KENO transport can make MG TSUNAMI-3D a viable method for some complex problems. Moreover, this adaptation allows MG sensitivity calculations with Shift, ORNL's next-generation high-performance Monte Carlo transport code, which currently does not offer any sensitivity capabilities with MG particle transport simulations. Initial implementation of the new MG TSUNAMI-3D sequence and its preliminary results with a selected critical benchmark experiment in the Verified, Archived Library of Inputs and Data (VALID) are presented in this study.

KENO↗
DOE OSTI2020

Development of the finite slab geometry in the NORA code

In this work we develop the complete linearization formalism for the radiation transfer equation in slab geometry. This development was motivated by an observed convergence instability arising from the NORA code’s inner boundary while modeling a suite of simple test problems in slab geometry presented at the Line Transfer Workshop held in LLNL in May of 2019. The convergence issue was seen to emerge early on in Newton-Raphson iterations. These test problems required parameters in the former NORA inner boundary to achieve the asymptotic limit of zero incoming intensity. In this work we derive the limiting case directly to gain insight into the former algorithm’s pathology.

97 MATHEMATICS AND COMPUTING↗