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

Evaluating a Commercial Dynamic Line Rating Software with the National PMU Dataset

To accelerate the development of data-driven applications for power systems, the Department of Energy (DOE) supported the collection and curation of a synchrophasor dataset spanning two years of observations from transmission utilities across the US. This National PMU Dataset (NPDS) was anonymized and distributed to awardees of a DOE research grant under nondisclosure agreements (NDAs) but has also been retained at PNNL to enable further research. Agreements with data contributors prevent the data from being shared outside the organization. However, establishing a blind research validation methodology is envisioned to maximize the value proposition of the NPDS. In this validation strategy, researchers may share algorithms/software (potentially as executables to protect intellectual property) with PNNL, and PNNL will share feedback about the software’s performance on subsets of the NPDS. Such a blind methodology ensures that sensitive information about critical infrastructure remains protected, but the value of the NPDS can be extended to research beyond PNNL. Through iterative feedback, the algorithms may be tweaked to address real-world artifacts. As the NPDS data is temporally and geographically diverse, it may capture features absent in smaller datasets used during the development of the algorithm under test. This report presents lessons learned from applying the blind validation methodology to LineID™, a synchrophasor-based dynamic line rating software developed by Topolonet Corporation. Improvements made to the software through iterative feedback, limitations of the validation methodology, as well as how the limitations of the NPDS affected the evaluation process are discussed. Observations indicate that the proposed validation methodology can be valuable for evaluating other tools in the future.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2010

Turbopump Seal Testing at Marshall Space Flight Center

The new ARES space flight program has presented many challenges to aerospace engineers and designers. One of the areas for consideration are the seals in the turbopumps that supply cryogenic propellants to the combustion chamber in the upper stage. Heritage face seals that worked in the past might not be sufficient in the newer turbopumps with increased speeds, pressures across the seals, and loads. New seal materials, engineering designs, and analysis techniques have been developed since the early use of these heritage seals, however, rub conditions and surface degradation at the sliding contact cannot be reliably predicted. Testing is required to determine the safe operating limits and verify seal wear life over the operating range. Rocketdyne in Canoga Park California entered into a task agreement with MSFC to design, fabricate, build, test, disassemble, and inspect hardware after tests of carbon materials and wear resistant coatings. The purpose of testing would be to determine the safe operating limits, empirically iterate the design, and select the best combination of materials for face seals and mating rings. This paper summarizes the many hours and efforts of individuals and teams to get the program operating successfully and presents the test results that were obtained.

Gibson, Howard G.↗
DOE OSTI2024

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

Predicting dynamic performance limits for servosystems with saturating nonlinearities

A generalized treatment for a system with a single saturating nonlinearity is presented and compared with frequency response plots obtained from an analog model of the system. Once the amplitude dynamics are predicted with the limit lines, an iterative technique is employed to determine the system phase response. The saturation limit line technique is used in conjunction with velocity and acceleration limits to predict the performance of an electro-hydraulic servosystem containing a single-stage servovalve. Good agreement was obtained between predicted performance and experimental data.

Webb, J. A., Jr.↗
DOE OSTI2021

Overcoming timestep limitations in boosted-frame particle-in-cell simulations of plasma-based acceleration

Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Loss of energetic particles due to resistive wall mode instability in ITER

Effects of an unstable n = 1 (n is the toroidal mode number) resistive wall mode (RWM) on the energetic particle (EP) confinement and loss are numerically investigated, for an ITER steady state scenario with 10 MA plasma current and 5.3 T toroidal field. The eigenfunction of the RWM is computed, with the associated three-dimensional magnetic field perturbation superposed with the 2D equilibrium field for tracing the EP drift orbits. Considered are mono-energetic EPs at 0.5 MeV and 1 MeV for deuterium ions, and 3.5 MeV for fusion-born alphas, with a range of distribution in the particle pitch angle. Modeling finds that less than 20% of EPs can be lost to the limiting surface in ITER assuming a source distribution uniform in minor radius, due to an unstable RWM that produces 100 Gauss poloidal field perturbation at the outboard mid-plane just inside the (effective) resistive wall surface. On top of the initial prompt drift orbit loss for counter-current EPs, the RWM induced particle loss occurs on a one second time scale, which is comparable to the RWM growth time in ITER. The 'wetted' area, due to the lost EPs striking the limiting surface, is generally found to be large due to the RWM. This is a favorable prediction for ITER. Here, the loss distribution in the poloidal angle is more uniform for co-current EPs. Counter-current EPs experience outward orbit drift when launched from the low-field side, and tend to more often hit the bottom region of the limiting surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The 1 → 3 massive splitting functions from QCD factorization and SCET

Splitting functions are universal functions describing the collinear dynamics of gauge theories, and as such are crucial ingredients for a wide variety of calculations in perturbative QCD. We present analytic results for the triple collinear splitting functions in QCD with a single massive parton. We derive the splitting functions using two distinct methods; first by expanding the squared matrix elements in the collinear limit, and secondly by using soft-collinear effective theory with massive quarks. We find agreement between these two approaches, providing a strong check of our results. Additionally, we also check all iterated and soft limits of our results, finding agreement with predictions from factorization. Our results provide an important ingredient for higher order perturbative calculations involving massive partons, and for the description of the collinear dynamics of heavy flavor jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

An efficient algorithm using matrix methods to solve wind tunnel force-balance equations

An iterative procedure applying matrix methods to accomplish an efficient algorithm for automatic computer reduction of wind-tunnel force-balance data has been developed. Balance equations are expressed in a matrix form that is convenient for storing balance sensitivities and interaction coefficient values for online or offline batch data reduction. The convergence of the iterative values to a unique solution of this system of equations is investigated, and it is shown that for balances which satisfy the criteria discussed, this type of solution does occur. Methods for making sensitivity adjustments and initial load effect considerations in wind-tunnel applications are also discussed, and the logic for determining the convergence accuracy limits for the iterative solution is given. This more efficient data reduction program is compared with the technique presently in use at the NASA Langley Research Center, and computational times on the order of one-third or less are demonstrated by use of this new program.

Smith, D. L.↗
DOE OSTI2024

Stochastic machine learning via sigma profiles to build a digital chemical space

This work establishes a different paradigm on digital molecular spaces and their efficient navigation by exploiting sigma profiles. To do so, the remarkable capability of Gaussian processes (GPs), a type of stochastic machine learning model, to correlate and predict physicochemical properties from sigma profiles is demonstrated, outperforming state-of-the-art neural networks previously published. The amount of chemical information encoded in sigma profiles eases the learning burden of machine learning models, permitting the training of GPs on small datasets which, due to their negligible computational cost and ease of implementation, are ideal models to be combined with optimization tools such as gradient search or Bayesian optimization (BO). Gradient search is used to efficiently navigate the sigma profile digital space, quickly converging to local extrema of target physicochemical properties. While this requires the availability of pretrained GP models on existing datasets, such limitations are eliminated with the implementation of BO, which can find global extrema with a limited number of iterations. A remarkable example of this is that of BO toward boiling temperature optimization. Holding no knowledge of chemistry except for the sigma profile and boiling temperature of carbon monoxide (the worst possible initial guess), BO finds the global maximum of the available boiling temperature dataset (over 1,000 molecules encompassing more than 40 families of organic and inorganic compounds) in just 15 iterations (i.e., 15 property measurements), cementing sigma profiles as a powerful digital chemical space for molecular optimization and discovery, particularly when little to no experimental data is initially available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1972

Stability of combustors with partial length acoustic liners

An analytical technique for the evaluation of combustion stability in rocket motors with partial length acoustic absorbers is presented. The combustors considered have concentrated combustion zones at the injector, finite mean flows, cylindrical cross sections, and acoustic liners of arbitrary length and impedance. Linear three dimensional oscillations in such combustion chambers are analyzed using an integral equation-iteration technique. Stability limits in terms of a combustion response factor are calculated for several values of Mach number, length to radius ratio, liner impedance, liner length, liner location, and nozzle admittance. Results indicate that increasing liner length increases combustor stability substantially at low Mach numbers but has a substantially smaller effect for larger Mach numbers. Increasing Mach numbers or length to radius ratio have destabilizing effects while liner location has only a minor effect on stability.

Mitchell, C. E.↗
NASA NTRS1972

Stability of combustors with partial length acoustic liners.

An analytical technique for the evaluation of combustion stability in rocket motors with partial length acoustic liners is presented. The combustors considered in this work have concentrated combustion zones at the injector, finite mean flows, cylindrical cross sections, and acoustic liners of arbitrary length and impedance. Linear three dimensional oscillations in such combustion chambers are analyzed using an integral equation-iteration technique. Stability limits in terms of the sensitive time lag model are calculated for several values of Mach number, length to radius ratio, liner impedance, liner length, liner location, and nozzle admittance. Results indicate that increasing liner length increases combustor stability substantially at low Mach numbers but has a considerably smaller effect for larger Mach numbers. Increasing Mach number or length to radius ratio has a destabilizing effect while liner location has only a minor effect on stability.

Mitchell, C. E.↗
NASA NTRS2000

A Safety Index and Method for Flightdeck Evaluation

If our goal is to improve safety through machine, interface, and training design, then we must define a metric of flightdeck safety that is usable in the design process. Current measures associated with our notions of "good" pilot performance and ultimate safety of flightdeck performance fail to provide an adequate index of safe flightdeck performance for design evaluation purposes. The goal of this research effort is to devise a safety index and method that allows us to evaluate flightdeck performance holistically and in a naturalistic experiment. This paper uses Reason's model of accident causation (1990) as a basis for measuring safety, and proposes a relational database system and method for 1) defining a safety index of flightdeck performance, and 2) evaluating the "safety" afforded by flightdeck performance for the purpose of design iteration. Methodological considerations, limitations, and benefits are discussed as well as extensions to this work.

Latorella, Kara A.↗
NASA NTRSDate not supplied

Rapid Development of Instrument Thermal Models: Perspectives and Guidelines from NASA Goddard’s Instrument Design Laboratory

- The design and development of robotic spaceflight instruments is a critical part of NASA’s vision to discover and expand knowledge for the benefit of humanity - For typical flight instrument projects, thermal engineers will develop initial instrument thermal models over weeks or months, then iterate them over a project’s lifespan – In each iteration, the engineer will: - Refine their thermal models and thermal designs in accordance with updates from other subsystems - Perform trade studies - Solve very detailed and complex analysis problems, including worst-cases and contingencies - Pick hardware and plan for testing and integration - However, prior to a project being established, or for proposal development at an early conceptual stage, the luxury of multiple instrument design iterations may be limited or nonexistent – Within a short timeline, how do you complete a thermal model or explore multiple possible instrument configurations? – What are the critical parameters for your model? Which details do you include or leave out?

Kan Yang↗
DOE OSTI2026

Modelling the limiter ramp-up of WEST for addressing the future challenges of ITER

This paper presents a joint experimental and numerical investigation into the physics of long limited plasma ramp-up in tokamaks with tungsten (W) first walls, a critical phase for ITER operation. The comparison between the average plasma quantities simulated using the SolEdge-HDG code and the measurements taken during three successive WEST discharges after boronisation shows how challenging it is to predict this phase. While simulations reproduce the general trends at the midplane, with reasonable match in density profiles, they consistently underestimate core temperatures possibly due to too large perpendicular heat conductivity. On the contrary, at the high field side (HFS) limiter, simulations overestimate the measured quantities, and highlights the limitation of using Bohm boundary conditions at grazing magnetic angles. Experimental measurements reveal that the boron layer is rapidly eroded, on a timescale comparable to a single ITER discharge. The subsequent transition from a boron-coated to a tungsten wall increases recycling and significantly degrades the core plasma, reducing the electron temperature by nearly half due to W contamination, despite wall parameters remaining stable. Furthermore, comparisons with Langmuir probes, bolometry, reflectometry, and spectroscopy indicate that the experimental far scrape-Off layer (SOL) is significantly wider than simulated. This wide SOL implies that boron erosion extends along the entire HFS limiter rather than being confined to the contact point. This work highlights some characteristics of the plasma during this phase of the discharge and emphasizes the current modelling issues that need to be resolved in order to obtain reliable predictions concerning the ITER ramp-up.

ITER↗
DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Aircraft Conflict Analysis and Real-Time Conflict Probing Using Probabilistic Trajectory Modeling

Methods for maintaining separation between aircraft in the current airspace system have been built from a foundation of structured routes and evolved procedures. However, as the airspace becomes more congested and the chance of failures or operational error become more problematic, automated conflict alerting systems have been proposed to help provide decision support and to serve as traffic monitoring aids. The problem of conflict detection and resolution has been tackled from a number of different ways, but in this thesis, it is recast as a problem of prediction in the presence of uncertainties. Much of the focus is concentrated on the errors and uncertainties from the working trajectory model used to estimate future aircraft positions. The more accurate the prediction, the more likely an ideal (no false alarms, no missed detections) alerting system can be designed. Additional insights into the problem were brought forth by a review of current operational and developmental approaches found in the literature. An iterative, trial and error approach to threshold design was identified. When examined from a probabilistic perspective, the threshold parameters were found to be a surrogate to probabilistic performance measures. To overcome the limitations in the current iterative design method, a new direct approach is presented where the performance measures are directly computed and used to perform the alerting decisions. The methodology is shown to handle complex encounter situations (3-D, multi-aircraft, multi-intent, with uncertainties) with relative ease. Utilizing a Monte Carlo approach, a method was devised to perform the probabilistic computations in near realtime. Not only does this greatly increase the method's potential as an analytical tool, but it also opens up the possibility for use as a real-time conflict alerting probe. A prototype alerting logic was developed and has been utilized in several NASA Ames Research Center experimental studies.

Yang, Lee C.↗
NASA NTRS1965

Human Tracking Performance in Uncoupled and Coupled Two-Axis Systems

This report presents tile results of an experimental and analytical study of human performance in uncoupled and coupled control systems. Human pilot performance in single and two-axis systems was mathematically modeled by linear second-order describing functions. Model parameters were determined using model matching techniques. Analysis of the models showed that the amplitude ratio and phase lead of the describing function increased with training indicating an increase in open loop bandwidth. The phase margin also decreased with training. Increasing the plant lag time constant resulted in an increase in the model lead time constant and a decrease in the zero frequency gain. No significant difference was found to exist in the normalized tracking error per axis between the two-axis tasks and the single-axis tasks. However tile model lead time constant was significantly greater in two-axis tracking. Manual tracking of two-axis systems with cross-coupling was studied experimentally and analytically. Approximate methods for modeling two-axis performance were developed and checked using a precise spectral analysis approach. Coupled and uncoupled, symmetrical and asymmetrical two-axis performance was compared. The results show that modeling of cross-coupled systems is feasible and that trained subjects are capable of decoupling the axes of some systems. A methodology study compared the identification performance of continuous, iterative, and extrapolation model matching techniques. An iterative technique employing sensitivity equations for the generation of influence coefficients was found to be the best technique due to its excellent identification accuracy and ease of implementation. Convergence in iterative techniques can be improved substantially by equalizing the parameter adjustment rates and limiting the maximum correction per iteration.

Todosiev, E. P.↗