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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

All sources NASA NTRS DOE OSTI

NASA NTRS1996

An Evaluation of the Measurement Requirements for an In-Situ Wake Vortex Detection System

Results of a numerical simulation are presented to determine the feasibility of estimating the location and strength of a wake vortex from imperfect in-situ measurements. These estimates could be used to provide information to a pilot on how to avoid a hazardous wake vortex encounter. An iterative algorithm based on the method of secants was used to solve the four simultaneous equations describing the two-dimensional flow field around a pair of parallel counter-rotating vortices of equal and constant strength. The flow field information used by the algorithm could be derived from measurements from flow angle sensors mounted on the wing-tip of the detecting aircraft and an inertial navigation system. The study determined the propagated errors in the estimated location and strength of the vortex which resulted from random errors added to theoretically perfect measurements. The results are summarized in a series of charts and a table which make it possible to estimate these propagated errors for many practical situations. The situations include several generator-detector airplane combinations, different distances between the vortex and the detector airplane, as well as different levels of total measurement error.

Fuhrmann, Henri D.↗

NASA NTRS1994

A Chebyshev Collocation Method for Moving Boundaries, Heat Transfer, and Convection During Directional Solidification

Free and moving boundary problems require the simultaneous solution of unknown field variables and the boundaries of the domains on which these variables are defined. There are many technologically important processes that lead to moving boundary problems associated with fluid surfaces and solid-fluid boundaries. These include crystal growth, metal alloy and glass solidification, melting and name propagation. The directional solidification of semi-conductor crystals by the Bridgman-Stockbarger method is a typical example of such a complex process. A numerical model of this growth method must solve the appropriate heat, mass and momentum transfer equations and determine the location of the melt-solid interface. In this work, a Chebyshev pseudospectra collocation method is adapted to the problem of directional solidification. Implementation involves a solution algorithm that combines domain decomposition, finite-difference preconditioned conjugate minimum residual method and a Picard type iterative scheme.

Zhang, Yiqiang↗

NASA NTRS2004

Efficient Viscous Design of Realistic Aircraft Configurations

This paper addresses the use of the Constrained Direct Iterative Surface Curvature (CDISC) design method in the aircraft design process. A discussion of some of the requirements for practical use of CFD in the design process is followed by a description of different CFD design methods, along with their relative strengths and weaknesses. A detailed description of the CDISC design method highlights some of the aspects of the method that provide computational efficiency and portability, as well as the flow and geometry constraint capabilities. In addition, an efficient approach to multipoint design, the Weighted Averaging of Geometries (WAG) method, is described and illustrated using a couple of simple examples. The CDISC and WAG methods are then applied to a complex generic business jet geometry using an unstructured grid flow solver to demonstrate the multipoint and multicomponent design capabilities of these methods. Introduction

Campbell, Richard L.↗

NASA NTRS2020

INVESTIGATING THE FEASIBILITY AND STABILITY FOR MODELING ACOUSTIC WAVE SCATTERING USING A TIME-DOMAIN BOUNDARY INTEGRAL EQUATION WITH IMPEDANCE BOUNDARY CONDITION

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. When designing next generation quiet and environmentally friendly aircraft, it is important to be able to accurately and efficiently predict the acoustic scattering by an aircraft body from a given noise source. Acoustic liners are an effective tool for aircraft noise reduction and are characterized by a frequency-dependent impedance. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. This work considers using either an impedance or an admittance (inverse of impedance) boundary condition to allow for acoustic scattering problems to be modeled with geometries consisting of both un-lined and lined surfaces. Three acoustic liner models are discussed: the Extended Helmholtz Resonator Model, the Three-Parameter Impedance Model, and the Broadband Impedance Model. In both the Helmholtz and Three-Parameter models, liner impedance is specified at a given frequency, whereas the Broadband model allows for the investigation of multiple frequencies simultaneously. The impedance and admittance boundary conditions for acoustic liners are derived for each model and coupled with a time-domain boundary integral equation. The scattering solution is obtained iteratively using a boundary element method with constant spatial and third-order temporal basis functions. Time-domain boundary integral equations are unfortunately prone to numerical instabilities due to resonant frequencies resulting from non-trivial solutions in the interior domain. When reformulated with the Burton-Miller method, the instabilities are eliminated. Using a Burton-Miller reformulation, the stability of the boundary element method assuming a liner boundary condition is assessed using eigenvalue analysis. The stability of each liner model is discussed, and it is shown that the Three-Parameter and Broadband models are sufficient for modeling an acoustic liner on the surface of scattering bodies. The Helmholtz model demonstrates strict limitations for stability, whereas the Three-Parameter and Broadband models are stable for most cases. Also included in this work is an assessment of the spatial accuracy of the time-domain boundary element method with respect to the surface element basis functions, as well as a performance study of the numerical algorithm.

computational mathematics↗

NASA NTRS2021

A Practical Guide to Writing a Radiative Transfer Code

Using our decades-long experience in radiative transfer (RT) code development for Earth science, we endeavor to reduce the knowledge gap of bringing RT from theory to code quickly. Despite numerous classic and recent literature, it is still hard to develop anRT code from scratch within a few weeks. It is equally hard to understand, not to mention modify, an existing “monster” RT code, for which the developer is either located remotely or has retired. Following the format of “Numerical Recipes” by Press et al., we collocate in this paper small pieces of necessary theory with corresponding small pieces of RT code. These are arranged in an order that is natural for code development, which is often opposite of the natural order for laying out the theoretical basis. We focus on the transfer of unpolarized monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface. Both the surface and the atmosphere are homogeneous (uniform) at all directions. The multiple scattering is numerically solved using the deterministic method of Gauss-Seidel iterations. Except for the presented Python-Numba open-source RT code gsit, the paper does not report any new scientific results, but rather serves as an academic demonstration. If development time is an issue or the reader is familiar with basic concepts of RT theory, we recommend proceeding directly to Sec.3 “RT code development.

multiple light scattering↗

NASA NTRSDate not supplied

Efficient Parametric Uncertainty Analysis of an Earth Entry Vehicle Concept Using Least Angle Regression

The objective of this work was to outline and apply an efficient and accurate parametric un-certainty propagation approach to the analysis of convective heating on an Earth entry vehicle concept. The described approach was based on Least Angle Regression used to solve a sparse and underdetermined linear system in the point-collocation non-intrusive polynomial chaos surrogate method. This approach involved an iterative process to computing the non-zero terms of the underlying polynomial chaos model using only enough samples to converge uncertainty interval predictions and Sobol index values based global nonlinear sensitivity estimates. The Earth entry vehicle was analyzed at three points along a representative trajectory for a Mars return mission. 329 sources of uncertainty were identified in the computational fluid dynamics model used to predict the forebody convective heating. These included uncertainty in flow field chemical rates, collision integrals, heats of formation, surface finite rate char model reaction rates, wall roughness height, and the turbulent Schmidt number. Results from this study showed that convective heating uncertainty as high as 50% of the nominal was predicted with only about 50 evaluations of the computational model. This was far fewer than would be required for a sampling-based approach or a full basis polynomial chaos model, which would have required over 50,000 samples. Additionally, results showed that over 90% of the total convective heating uncertainty was due to uncertainty in the N2catalytic rate on the surface, while the remainder of the uncertainty was attributed to the turbulent Schmidt number and the wall roughness uncertainties.

Thomas K West IV↗

NASA NTRS2024

Airfoil Table Modification for Improved Calculation of Ingenuity Rotor Performance

To better characterize the behavior of the Ingenuity helicopter, multiple test campaigns were conducted. These tests included the Engineering Design Model 1 (EDM1) tests and the Transonic Rotor Test (TRT), both of which assessed Ingenuity rotor performance at various conditions. Using a comprehensive rotorcraft analysis tool, CAMRAD II, test setups of the two different Ingenuity test campaigns were simulated. The EDM1 test used coaxial test stand which focused on varying collective at various feasible Martian atmospheric densities and the TRT test used a single rotor test stand which focused on varying RPM to assess performance at different blade tip Mach numbers. Initial CAMRAD II analyses yielded performance results which differed significantly from experimental results. To reconcile the experimental results and the CAMRAD II analyses, modifications were made to the input airfoil decks for CAMRAD II. After several iterations of modifications, a specific method of modification was found that more accurately predicted experimental behavior at all tested densities for the coaxial rotor. This same modification was less successful for the single rotor test results. There are a number of differences between coaxial and single rotor aerodynamics and between the two test stands (hub and control system); also, there is a small influence of RPM (mainly tip Mach number) on dimensionless performance. These differences may be causing the coaxial rotor modifications not being successful for the single rotor.

Airfoil↗

NASA NTRS1992

Enhancements to viscous-shock-layer technique

A solution procedure is presented which considerably improves the computational efficiency of the viscous-shock-layer technique, especially for long slender bodies. The 'predictor-corrector' procedure suggested for obtaining the shock shape beyond the nose region requires only a single global pass. The accuracy of the present method is demonstrated by comparison with globally iterated results over the entire body and with ground- and flight-test data. A good comparison of the results computed with the two methods is shown for different flowfield chemistry models and axisymmetric body shapes. The new procedure results in computer run times 1/3 to 1/2 of the times required for the full-body global iteration procedure. Further, the algebraic expressions used to specify the initial shock shape eliminate the need for a shock shape generated by external means and permit immediate introduction of the full viscous-shock-layer equations. Finally, the present method of solution for the VSL equations provides to the aerothermal designer a very efficient and accurate tool for detailed flowfield as well as future technology studies.

Gupta, Roop N.↗

NASA NTRS1984

Tomography by iterative convolution - Empirical study and application to interferometry

An algorithm for computer tomography has been developed that is applicable to reconstruction from data having incomplete projections because an opaque object blocks some of the probing radiation as it passes through the object field. The algorithm is based on iteration between the object domain and the projection (Radon transform) domain. Reconstructions are computed during each iteration by the well-known convolution method. Although it is demonstrated that this algorithm does not converge, an empirically justified criterion for terminating the iteration when the most accurate estimate has been computed is presented. The algorithm has been studied by using it to reconstruct several different object fields with several different opaque regions. It also has been used to reconstruct aerodynamic density fields from interferometric data recorded in wind tunnel tests.

Vest, C. M.↗

DOE OSTI2020

Iterative X-ray spectroscopic ptychography

Spectroscopic ptychography is a powerful technique to determine the chemical composition of a sample with high spatial resolution. In spectro-ptychography, a sample is rastered through a focused X-ray beam with varying photon energy so that a series of phaseless diffraction data are recorded. Each chemical component in the material under investigation has a characteristic absorption and phase contrast as a function of photon energy. Using a dictionary formed by the set of contrast functions of each energy for each chemical component, it is possible to obtain the chemical composition of the material from high-resolution multi-spectral images. This paper presents SPA (spectroscopic ptychography with alternating direction method of multipliers), a novel algorithm to iteratively solve the spectroscopic blind ptychography problem. First, a nonlinear spectro-ptychography model based on Poisson maximum likelihood is designed, and then the proposed method is constructed on the basis of fast iterative splitting operators. SPA can be used to retrieve spectral contrast when considering either a known or an incomplete (partially known) dictionary of reference spectra. By coupling the redundancy across different spectral measurements, the proposed algorithm can achieve higher reconstruction quality when compared with standard state-of-the-art two-step methods. It is demonstrated how SPA can recover accurate chemical maps from Poisson-noised measurements, and its enhanced robustness when reconstructing reduced-redundancy ptychography data using large scanning step sizes is shown.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Evaluation of gamma-ray transmission through rectangular collimator slits for application in nuclear fuel spectrometry

Gamma-ray spectrometry is widely applied in several science fields, and in particular in non-destructive gamma scanning and gamma emission tomography of irradiated nuclear fuel. Usually, a collimator is used in the experimental setup, to selectively interrogate a region of interest in the fuel. For the optimization of instrument design, as well as for planning measurement campaigns, predictive models for the transmitted gamma-ray intensity through the collimator are needed. Commonly, Monte Carlo Radiation Transport tools are used for accurate prediction of gamma-ray transport, however, the long computation time requirements when used in low-efficiency experimental setups present challenges. Here, the full-energy peak intensity transmitted through a rectangular collimator slit was examined. A uniform planar surface source emitting isotropically was considered, and the rate of photons reaching an ideal counter plane on the opposite side of the collimator was evaluated by analytical integration. To find a closed-form primitive function, some idealizations were required, and thereby parametric models were obtained for the optical field of view, dependent on slit dimensions (length, height and width) and source-to-collimator distance. For contributions from outside the optical field of view, where a closed-form expression cannot be found, fast numerical integral methods were instead used. The results were validated using the Monte Carlo code MCNP6 and show an agreement within three percent for the numerical method. For the analytical method, deviations up to tens of percent were obtained, which is deemed to still be sufficient for instrument design and measurement planning, where often the order of magnitude of the count rate is not a priori known. The method is planned for use in iterative optimization routines in the design of Gamma Emission Tomography devices, as well as for the prediction of gamma spectra obtained in the planning of fuel inspections. An application of the proposed method was demonstrated in spectrum prediction for a short cooling-time fuel rod test from the Halden reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Cross Section Generation Capability in Griffin

The Griffin code is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code reads multigroup cross sections in the ISOXML format generated from external deterministic or Monte Carlo cross section generation codes. The implementation of the cross section generation capability in Griffin was initiated last year by plugging in the cross section application programming interface (CSAPI) and reviewing the methodologies for treating particulate fuels. The focus this year was on improving the CSAPI integration and implementing advanced self-shielding methods for applications to advanced reactor problems with TRISO fuels. First, the process for cross section library generation was updated to accurately and rigorously produce isotopic cross section data. Second, the equivalent Dancoff factor cell method performing slowing down calculations on the fly for the resonance treatment was implemented in CSAPI to improve the accuracy of effective multigroup cross sections in the resonance energy range. Third, the iterative local spatial self-shielding method was implemented under the calculation framework of the equivalent Dancoff factor cell method to accurately deal with the double heterogeneity effect of particulate fuel. The updated CSAPI with the advanced self-shielding methods, together with the cross section libraries generated based on the improved process, were tested for pin-cell, unit-cell, and fuel assembly problems with various resonance self-shielding conditions based on very high temperature reactor, high temperature test reactor, and Empire benchmark cores, indicating that the updated CSAPI in Griffin is able to produce multigroup cross sections accurately and efficiently. We also showed that the methodology worked well for pebble bed fuel from HTR-10, but the capability still needs to be fully integrated into CSAPI. In the future, further benchmark tests will be performed for various thermal reactor core problems, including particulate fuel-based pebble bed reactors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

NASA NTRS1991

Domain decomposition methods in computational fluid dynamics

The divide-and-conquer paradigm of iterative domain decomposition, or substructuring, has become a practical tool in computational fluid dynamic applications because of its flexibility in accommodating adaptive refinement through locally uniform (or quasi-uniform) grids, its ability to exploit multiple discretizations of the operator equations, and the modular pathway it provides towards parallelism. These features are illustrated on the classic model problem of flow over a backstep using Newton's method as the nonlinear iteration. Multiple discretizations (second-order in the operator and first-order in the preconditioner) and locally uniform mesh refinement pay dividends separately, and they can be combined synergistically. Sample performance results are included from an Intel iPSC/860 hypercube implementation.

Gropp, William D.↗

NASA NTRS1992

Domain decomposition methods in computational fluid dynamics

Gropp, William D.↗

NASA NTRS1991

Finite element implementation of state variable-based viscoplasticity models

The implementation of state variable-based viscoplasticity models is made in a general purpose finite element code for structural applications of metals deformed at elevated temperatures. Two constitutive models, Walker's and Robinson's models, are studied in conjunction with two implicit integration methods: the trapezoidal rule with Newton-Raphson iterations and an asymptotic integration algorithm. A comparison is made between the two integration methods, and the latter method appears to be computationally more appealing in terms of numerical accuracy and CPU time. However, in order to make the asymptotic algorithm robust, it is necessary to include a self adaptive scheme with subincremental step control and error checking of the Jacobian matrix at the integration points. Three examples are given to illustrate the numerical aspects of the integration methods tested.

Iskovitz, I.↗

NASA NTRS1990

Comparison of coupled-cluster methods which include the effects of connected triple excitations

The 'coupled cluster single, double, and triple' (CCSDT) excitation model has been used to ascertain electron correlation energies for 14 different molecules representing a variety of chemical bonds, in conjunction with several methods of this type 'CCSDT-x', which include only an approximate treatment of connected triple excitations; these methods encompass CCSDT-1a, -1b, -2, -3, and -4, as well as the novel CCSD(T). While all methods treat the effects of connected triple excitations iteratively, CCSD(T) approaches then perturbationally. For the 14 molecules considered, the CCSD(T) method's average error relative to CCSDT is substantially lower than any of the CCSDT-x methods.

Scuseria, Gustavo E.↗