Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration method”

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 1,765 records · Page 98

Dancoff-based Wigner-Seitz approximation for the subgroup resonance self-shielding in the VERA neutronic simulator MPACT

The MPACT neutronics module of the Virtual Environment for Reactor Analysis (VERA) has used the subgroup method for resonance self-shielding calculation, for which two-dimensional (2D) fixed-source transport calculations are performed using the method of characteristics for resonance energy groups. When considering thermal feedbacks, the subgroup calculation must be performed at each outer iteration. Therefore, the computing time for cross section processing is a significant burden for computational efficiency. The Dancoff-based Wigner-Seitz approximation (DWA) capability has been implemented into MPACT in conjunction with the subgroup method, which has been used in SCALE/XSProc since SCALE version 6.0 and has recently been called an equivalent Dancoff-factor cell (EDC) method. The issue of relatively large reactivity bias in DWA for the gadolinia rods was resolved by introducing multiple Dancoff factors. Benchmark results for the VERA pressurized and boiling water reactor benchmark suites show that the DWA capability would significantly enhance computational efficiency with comparable accuracy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Experiments with a self-correcting wind tunnel

The feasibility of controlling the flow actively through the walls of a transonic, porous wall wind tunnel in order to minimize wall interference effects on a test model is demonstrated. The method is based upon measuring the components of the disturbance velocity at discrete points along an imaginary surface in the flow field within the tunnel. A mathematical formulation of the flow field exterior to the surface including the boundary condition for unconfined flow, i.e., that all disturbance vanish at infinite, is used to determine if these measured velocity components are consistent with that boundary condition. If they are not, the theory provides a better approximation to the velocity component for unconfined flow, and the flow through the tunnel walls is readjusted iteratively until the measured quantities are consistent with unconfined flow. A brief review of theoretical methods is followed by a description of the Calspan self correcting wind tunnel design and operation, calibration with and without active wall control. Typical results obtained by approximating a conventional porous wall wind tunnel for an 0012 airfoil show that active wall control largely reproduces the correct shock wave position, eliminates wall interference of lift and drag, and reduces the interference effects on pitching moment to 10%.

Vidal, R. J.↗

Numerical treatment of nonlinear latent heat boundary conditions at moving interfaces in genuine two dimensional solidification problems

The proposed method for the treatment of two dimensional solidification problems is based on quasilinearization of the transformed heat conduction equation and latent heat condition at the interface and an iterative sequence in which these are solved simultaneously. Modern algorithms for solving such sparse systems mean that most of the storage advantage of other methods are reduced and the speed of solution can be improved.

Beckett, P. M.↗

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a fine numerical resolution. Here, the main challenge of such computations is the high memory effort of classical solvers, which limits the accuracy of chosen discretizations. In this work, we derive a method for the computation of effective eigenvalues when the underlying solution has a low-rank structure. This is accomplished by utilizing dynamical low-rank approximation (DLRA), which is an efficient strategy to derive time evolution equations for low-rank solution representations. The main idea is to interpret the iterates of the classical inverse power iteration as pseudo-time steps and apply the DLRA concepts in this framework. In our numerical experiment, we demonstrate that our method significantly reduces memory requirements while achieving the desired accuracy. Analytic investigations show that the proposed iteration scheme inherits the convergence speed of the inverse power iteration, at least for a simplified setting.

97 MATHEMATICS AND COMPUTING↗

Bayesian Optimization of Catalysis with In-Context Learning

Large language models (LLMs) can perform accurate classification with zero or few examples through in-context learning (ICL), allowing the model to observe query-relevant examples at inference time and eliminating the need for additional weight updates to generalize beyond its original training data. We extend this capability to regression with uncertainty estimation using frozen LLMs (e.g., GPT-4o, Gemini), enabling Bayesian optimization (BO) in natural language without explicit model training or feature engineering. We apply this to materials discovery by representing materials as synthesis and testing procedures for use in natural language prompts. This Bayesian, design-first approach prioritizes optimization toward target material properties before detailed characterization, in contrast to conventional experimental workflows that often emphasize characterization of suboptimal materials. On benchmarks like aqueous solubility and oxidative coupling of methane (OCM), BO-ICL matches or outperforms Gaussian processes. In live experiments on the reverse water–gas shift (RWGS) reaction, BO-ICL identifies multimetallic catalysts that approach equilibrium CO yield within 6 and 10 iterations from a pool of 3,700 and 360,000 candidates, respectively. Our method redefines materials representation and accelerates discovery, with broad applications across catalysis, materials science, and AI.

Calibration↗

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

Matilda v1.0: An R package for probabilistic climate projections using a reduced complexity climate model

A primary advantage to using reduced complexity climate models (RCMs) has been their ability to quickly conduct probabilistic climate projections, a key component of uncertainty quantification in many impact studies and multisector systems. Providing frameworks for such analyses has been a target of several RCMs used in studies of the future co-evolution of the human and Earth systems. In this paper, we present Matilda, an open-science R software package that facilitates probabilistic climate projection analysis, implemented here using the Hector simple climate model in a seamless and easily applied framework. The primary goal of Matilda is to provide the user with a turn-key method to build parameter sets from literature-based prior distributions, run Hector iteratively to produce perturbed parameter ensembles (PPEs), weight ensembles for realism against observed historical climate data, and compute probabilistic projections for different climate variables. This workflow gives the user the ability to explore viable parameter space and propagate uncertainty to model ensembles with just a few lines of code. The package provides significant freedom to select different scoring criteria and algorithms to weight ensemble members, as well as the flexibility to implement custom criteria. Additionally, the architecture of the package simplifies the process of building and analyzing PPEs without requiring significant programming expertise, to accommodate diverse use cases. We present a case study that provides illustrative results of a probabilistic analysis of mean global surface temperature as an example of the software application.

54 ENVIRONMENTAL SCIENCES↗

Integration of sensors through additive manufacturing leading to increased efficiencies of gas turbines for power generation and propulsion

To realize the full capability of additively manufactured components in complex energy systems, it is imperative to minimize early component failures during development phases and during operation. Traditional field feedback timelines and offline inspection protocols significantly reduce the design-manufacturing iteration times. To address this specific question, the project developed and demonstrated a method for the integration of sensors into complex components through additive manufacturing. The team used gas turbine engines as a platform, which meets the need of both power generation and propulsion and offer opportunities for cost reductions and efficiency increases. The innovation of this intelligent integration of sensors into complex components uniquely customized to address questions of integrity and durability for additively manufactured components. With real-time sensing data from additively manufactured components, turbine manufacturers will realize higher efficiencies, reduced component failures, and a 30-50% acceleration in product deployment of high efficiency gas turbine components due to a faster reduction in component risk assessment under actual operating conditions. This is a transformative shift towards a data-driven design and qualification of additively manufactured gas turbine components. To directly integrate sensors into additively manufactured components with all the complexities of actual hardware, powder bed fusion (direct metal laser sintering) and laser metal deposition technologies was developed. Validation took take place in two university laboratories both of which contain actual engine hardware and closely simulate a gas turbine prior to demonstrating the technology in a turbine development test. Indeed, two major technologies from this research cold impact turbine systems in the near future: (1) higher efficiency materials and designs enabled by additive manufacturing with 50% faster design to manufacturing cycle time, to enable faster time-to-market targets; and (2) integration of sensors into additively manufactured components enabling broad health and condition based prognostics for faster component and engine risk reduction.

33 ADVANCED PROPULSION SYSTEMS↗

Reachable sets for linear dynamical systems.

The properties of reachable sets for linear dynamical systems for specified control sets are discussed. Iterative procedures for determining numerical approximations of the reachable set are suggested and methods of obtaining an admissible control function which transfers an initial state to as near a prescribed target as possible is described. The problem of reachability with multiple control constraints is discussed and certain aspects of reachability for time-invariant systems with adjustable parameters is considered.

Pecsvaradi, T.↗

Inversion of infrared limb emission measurements for temperature and trace gas concentrations

Limb emission measurements are characterized by sharp weighting functions at high altitudes, and for temperature determinations, strongly nonlinear dependence of the weighting function on the temperature. Several methods for inverting this type of measurement have been described and used, including iterative, statistical, nonlinear and approximate direct approaches. These approaches are described and advantages and disadvantages of each are outlined.

Gille, J. C.↗

High-Resolution Numerical Simulation Of Shock Waves

NASA technical memorandum compares results of upwind and symmetric shock-capturing methods in numerical simulation of gas-dynamic flows. Both methods find shocks as sharp variations in fluid properties over few grid points. Methods differ in type of artificial viscosity introduced to stabilize computations. Symmetric method shown as accurate as upwind method, but with fewer and simpler time steps in transient case or fewer iteration steps in steady-state case.

Yee, H. C.↗

Finite element large deflection random response of thermally buckled plates

The governing nonlinear equations of motion are derived for a rectangular isotropic plate subjected to a combined thermal-acoustic loading. The thermal loading is taken to be a steady state temperature distribution, and the acoustic loading is considered to be a stationary Gaussian pressure with zero mean and uniform magnitude and phase over the plate surface. With these types of loads, the thermal-acoustic analysis can be separated into two problems: (1) thermal post-buckling; and (2) random vibration of the thermally buckled plate. Since both analyses are nonlinear, iterative techniques are used to solve each. For the random vibration analysis, the method of equivalent linearization is used; both auto- and cross-correlation terms are included. Results obtained using the present formulation are compared with previous classical equivalent linearization results.

Locke, J. E.↗

Registration of video sequences from multiple sensors

In this paper, we describe an approach for registration of video sequences from a suite of multiple sensors including television, infrared and radar. Video sequences generated by these sensors may contain abrupt changes in local contrast and inconsistent image features, which pose additional difficulties for registration. Our approach to registration addresses the difficulties caused by using multiple sensors. We use a representation for registration that is invariant to local contrast changes, followed by smoothing of the resulting error measure used for registration, for robust estimation of registration parameters. We use an iterative procedure to reduce the effect of inconsistent features. Finally, we describe a method that uses same-sensor registration to aide in performing registration of sequences of video frames across multiple sensors.

Sharma, Ravi K.↗

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

Simplified Abrasion Test Methodology for Candidate EVA Glove Lay-Ups

During the Apollo Program, space suit outer‐layer fabrics were badly abraded after performing just a few extravehicular activities (EVAs). For example, the Apollo 12 commander reported abrasive wear on the boots that penetrated the outer‐layer fabric into the thermal protection layers after less than 8 hrs of surface operations. Current plans for the exploration planetary space suits require the space suits to support hundreds of hours of EVA on a lunar or Martian surface, creating a challenge for space suit designers to utilize materials advances made over the last 40 years and improve on the space suit fabrics used in the Apollo Program. Over the past 25 years the NASA Johnson Space Center Crew and Thermal Systems Division has focused on tumble testing as means of simulating wear on the outer layer of the space suit fabric. Most recently, in 2009, testing was performed on 4 different candidate outer layers to gather baseline data for future use in design of planetary space suit outer layers. In support of the High Performance EVA Glove Element of the Next Generation Life Support Project, testing a new configuration was recently attempted in which require 10% of the fabric per replicate of that need in 2009. The smaller fabric samples allowed for reduced per sample cost and flexibility to test small samples from manufacturers without the overhead to have a production run completed. Data collected from this iteration was compared to that taken in 2009 to validate the new test method. In addition the method also evaluated the fabrics and fabric layups used in a prototype thermal micrometeoroid garment (TMG) developed for EVA gloves under the NASA High Performance EVA Glove Project. This paper provides a review of previous abrasion studies on space suit fabrics, details methodologies used for abrasion testing in this particular study, results of the validation study, and results of the TMG testing.

Rabel, Emily↗

Study of Geometric Porosity on Static Stability and Drag Using Computational Fluid Dynamics for Rigid Parachute Shapes

This paper explores use of computational fluid dynamics to study the e↵ect of geometric porosity on static stability and drag for NASA's Multi-Purpose Crew Vehicle main parachute. Both of these aerodynamic characteristics are of interest to in parachute design, and computational methods promise designers the ability to perform detailed parametric studies and other design iterations with a level of control previously unobtainable using ground or flight testing. The approach presented here uses a canopy structural analysis code to define the inflated parachute shapes on which structured computational grids are generated. These grids are used by the computational fluid dynamics code OVERFLOW and are modeled as rigid, impermeable bodies for this analysis. Comparisons to Apollo drop test data is shown as preliminary validation of the technique. Results include several parametric sweeps through design variables in order to better understand the trade between static stability and drag. Finally, designs that maximize static stability with a minimal loss in drag are suggested for further study in subscale ground and flight testing.

Greathouse, James S.↗

Efficacy of FTIR Analysis in Determining CO2 Loading on Diglycolamine

In support of advanced air revitalization technologies to enable human spaceflight beyond low earth orbit, performance studies have been conducted using a liquid amine, Diglycolamine (DGA) between teams at NASA's Johnson Spaceflight Center (JSC) and Ames Research Center (ARC). Liquid amines have been used in regenerable earth-based systems to remove CO2 from industrial systems as well as for closed-environment air revitalization because they can be regenerated at lower temperatures than solid sorbent systems. As an additional advantage to solid sorbent-based systems, liquid sorbents can be cycled between an adsorbing contactor and degassing chamber, thereby reducing system complexity by operation in a continuous loop. In an effort to inform a regeneration system design for micro-gravity applications, ARC has performed a number of tests to characterize the degas mechanics of DGA. In order to accurately measure the amount of CO2 captured or released by the amine, methods such as gravimetric weighing and chemical desorption are reasonable, however the first iteration test setup for a scaled down degas system required analysis on small sample sizes. Fourier-transform infrared spectroscopy (FTIR) analysis was experimentally evaluated to analyze CO2 concentration because it can produce measurements with sample sizes on the order of 100's of μL. Calibration against chemical desorption showed relatively good correlation and test data showed reasonable adherence to expected trends, however more extensive testing should be conducted to fully validate the usage of FTIR to determine CO2 loading on DGA.

Diglycolamine↗

Development of a Maxwell X-57 High Lift Motor Reference Design

NASA's all-electric X-57 airplane will utilize 14 electric motors, of which 12 are exclusively for lift augmentation during takeoff and landing. This report covers the design and development process taken to create an open reference model representative of the 12 lift augmenting motors. A combined worst case scenario was used as the design point, which represents the simultaneously occurring worst case aspects of thermal, static stress, electromagnetic, and rotor dynamic conditions. This work also highlights the tightly coupled nature of aerospace electric motor design, requiring constant iteration between all disciplines involved. Further adding to the uniqueness is the cooling method, which is limited to nacelle skin forced convection cooling only, no internal air flow is permitted. The stator outer diameter limit of 156.45 mm greatly impacts the degree of coupling between the electromagnetic design with the thermal analysis. The permanent magnet synchronous motor developed here operates between 385 V and 538 V, at a peak current of 50 A. Detailed electromagnetic, thermal, static load, and rotordynamic analysis was completed for this electric motor; all of which are required for a full design. The rotordynamic analysis took into consideration the motor housing which is designed specifically for this motor. The final electric motor has a mass of 2.34 kg, produces 24.1 Nm of torque with a specific power of 5.56 kW/kg, and has an efficiency of 96.61% at the combined worst case design point.

Electric Aircraft↗