SEARCH · Engineering Papers

Results for “iterative methods”

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,747 records · Page 97

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Coupled Radiation Effects in Thermochemical Nonequilibrium Shock-Capturing Flowfield Calculations

Lunar and Mars return conditions are examined using the LAURA flow field code and the LORAN radiation code to assess the effect of radiative coupling on axisymmetric thermochemical nonequilibrium flows. Coupling of the two codes is achieved iteratively. Special treatment required to couple radiation in a shock-capturing method is discussed. Results indicate that while coupling effects are generally the same as occur in equilibrium flows, under certain conditions radiation can modify the chemical kinetics of a nonequilibrium flow and thus alter relaxation processes. Coupling effects are found to be small for all cases considered, except for a five meter diameter aerobrake returning from Mars at 13.6 kilometers per second.

Hartung, Lin C.↗

NASA NTRS2006

Thermal Model of a Current-Carrying Wire in a Vacuum

A computer program implements a thermal model of an insulated wire carrying electric current and surrounded by a vacuum. The model includes the effects of Joule heating, conduction of heat along the wire, and radiation of heat from the outer surface of the insulation on the wire. The model takes account of the temperature dependences of the thermal and electrical properties of the wire, the emissivity of the insulation, and the possibility that not only can temperature vary along the wire but, in addition, the ends of the wire can be thermally grounded at different temperatures. The resulting second-order differential equation for the steady-state temperature as a function of position along the wire is highly nonlinear. The wire is discretized along its length, and the equation is solved numerically by use of an iterative algorithm that utilizes a multidimensional version of the Newton-Raphson method.

Border, James↗

NASA NTRS2016

ISS Double-Gimbaled CMG Subsystem Simulation Using the Agile Development Method

This paper presents an evolutionary approach in simulating a cluster of 4 Control Moment Gyros (CMG) on the International Space Station (ISS) using a common sense approach (the agile development method) for concurrent mathematical modeling and simulation of the CMG subsystem. This simulation is part of Training systems for the 21st Century simulator which will provide training for crew members, instructors, and flight controllers. The basic idea of how the CMGs on the space station are used for its non-propulsive attitude control is briefly explained to set up the context for simulating a CMG subsystem. Next different reference frames and the detailed equations of motion (EOM) for multiple double-gimbal variable-speed control moment gyroscopes (DGVs) are presented. Fixing some of the terms in the EOM becomes the special case EOM for ISS's double-gimbaled fixed speed CMGs. CMG simulation development using the agile development method is presented in which customer's requirements and solutions evolve through iterative analysis, design, coding, unit testing and acceptance testing. At the end of the iteration a set of features implemented in that iteration are demonstrated to the flight controllers thus creating a short feedback loop and helping in creating adaptive development cycles. The unified modeling language (UML) tool is used in illustrating the user stories, class designs and sequence diagrams. This incremental development approach of mathematical modeling and simulating the CMG subsystem involved the development team and the customer early on, thus improving the quality of the working CMG system in each iteration and helping the team to accurately predict the cost, schedule and delivery of the software.

Inampudi, Ravi↗

NASA NTRSDate not supplied

A preliminary review and evaluation of historical unrest and eruptive activity at Misti volcano, Arequipa, Peru

Since the Spanish arrived in Arequipa in the mid-16th century, 11 potential small-volume eruptions or periods of unrest are reported at Misti. Yet, documentation of these events is scant, often based upon secondary or tertiary sources, and has never been assessed using modern volcanological methods. As a result, Misti?s historical record contains significant uncertainty and iterated errors. We assess these events with modern volcanological knowledge and when possible return to primary published and archival sources. Additionally, we scrutinized historical newspapers, photos, and contemporary travelers? accounts to better understand the volcano?s activity since the early-19th century. Misti?s earliest reported activity in 1542CE illustrates the problems associated with the volcano?s historical record. The 1542CE event first appears in the literature in the late-19th century. Despite some uncertainty of its validity in these reports, it remains listed as a potential eruption. Numerous contemporaneous accounts document a major earthquake in Arequipa on 22 January 1582CE, causing major damage and widespread fear of an eruption. Upon investigation of the volcano, such fear was quickly dispelled. We suggest that Misti?s 1542CE event did not occur, rather its origin is the fear induced by the 1582CE tectonic earthquake and an iterated typographic error in the date. Historical documents also sometimes record when Misti was explicitly not active, such as in a letter from a resident of Arequipa published in 1872CE. Events are also missing from Misti?s reported history, such as a period of enhanced degassing in 1971CE that is recorded in newspaper articles and eye-witness accounts. In summary, we consider three of Misti?s reported events to be false, four to be possible but unlikely, three are likely, and the 1985CE degassing crisis certainly occurred. We further uncovered preliminary evidence of up to four additional degassing crises at the volcano during the 19th and 20th centuries.

Christopher Harpel↗

NASA NTRS2022

Drought Monitoring with Drones: A Hundred Fields at a Time

Drought frequency and severity are likely to increase due to global warming. Droughts already have a substantial negative influence on agriculture and the economy and finding ways to reduce their effects could have a monumental impact. Although NASA already has satellites deployed to collect drought data, these satellites are more for global drought indexing than local. To alleviate droughts on a local level, this paper proposes the use of drones to map soil moisture, plant health, and other drought indicators. The proposed drone design is a fixed-wing UAV equipped with a hyperspectral camera, a LiDAR sensor, and an array of weather sensors. These tools will permit it to reliably capture the necessary data to enhance suggestions on improving drought management practices. The collected data from the drones will be deployable in many ways, including for agricultural and non-agricultural applications. The hyperspectral camera has applications in monitoring the health of crops within a field to direct relief measures to the crops most in need. Thermal and LiDAR imaging can be deployed for locating leaks, predicting shortages of water bodies, and determining a field’s water needs. To implement a drought-monitoring drone, the recommended steps include building a prototype drone design that is equipped with the outlined instruments. The prototype drone could then be deployed to collect training data to guide a neural network that would provide interpretations and predictions from the data for users. After the prototype design is iterated upon, it will be ready for deployment and inform water management methods to serve the world in our battle against drought.

climate change↗

NASA NTRSDate not supplied

Overlapping Conditions in IMPACT

Introduction The IMPACT 1.0 model assumes independence of 119 different medical conditions. That is, it is assumed that one medical event maps to a single medical condition. However, it is often the case that a single medical event may lead to several different concurrent conditions (e.g., a single accident resulting in multiple fractures, chest/abdominal trauma, and sepsis). There is a significant overlap between many of the conditions in ICL 1.0. For example, depression and anxiety are two separate conditions but often co-exist. This overlap may lead to an overestimation in the incidence calculation. In future iterations of the model, we plan to address progression of medical conditions. Methods There are two possible methods to approach progression of a medical condition. The first is to take any medical condition in ICL, calculate its incidence and determine its outcomes no matter what sequelae occur. For example, calculate the incidence of prostatitis and determine its metrics (e.g., loss of crew life, risk of medical evacuation or task time affected) no matter what the sequelae might be (sepsis, hydronephrosis, renal failure, etc.). The second approach is to take any medical condition in the ICL and calculate the outcomes based on transitions to any other condition. For example, calculate the incidence of nephrolithiasis and determine the incidence of each possible transition (e.g., hydronephrosis, renal failure, UTI, pyelonephritis, sepsis) and then determine its metrics. With perfectly informative evidence, the answers should be the same using either method. The question is which approach best reduces the opportunity for overlapping. Results The team identified 31 conditions that were the most likely to transition to a secondary condition such as sepsis and respiratory failure. It was felt that the first approach, as described above, would be the easiest to implement and would have the greatest reduction in overlapping conditions. The identified conditions were primarily infectious conditions, like UTI and pneumonia (which could lead to both respiratory failure and sepsis), toxic inhalations, cardiac arrest, seizures, and trauma. Once these conditions were identified, the conditions were divided between the four clinicians on the team to identify those that were not consistent with the others in way they were approached. There were 13 conditions that were found to be problematic. Some condition definitions needed to be changed to eliminate references to transitions or to update them in cases where the definitions were out of date. Some conditions required new evidence to correct the loss of crew life, and one condition was felt to have so much overlap that it was recommended for elimination. A second effort required modification of the sepsis incidence calculation to prevent double counting of the conditions, such as prostatitis, that could progress to sepsis. Conclusion The team made significant progress in clarifying the method for dealing with transitions from one medical condition to another. In addition, errors were found and corrected in definitions, incidence, and loss of crew life as well as a condition in which reference was made to a CLiFF that was no longer included. These modifications will provide greater fidelity in the IMPACT model and reduce overlap.

Arian Anderson↗

DOE OSTI2021

Prototype tests of the electromagnetic particle injector-2 for fast time response disruption mitigation in tokamaks

Predicting and controlling disruptions is an important and urgent issue for ITER. Some disruptions with a short warning time may be unavoidable. For these cases, a fast time response disruption mitigation method is essential. Experimental tests on a prototype system of a novel, rapid time-response disruption mitigation system being developed for tokamak-based reactors and as a backup option for ITER, referred to as the electromagnetic particle injector (EPI), have been able to verify the primary advantages of the concept. These are its ability to meet short warning time scales of <10 ms while attaining the projected high velocities for deep radiative payload penetration in reactor-scale plasmas. The EPI relies on an electromagnetic propulsion system. A metallic sabot is accelerated electromagnetically to the required velocities (>1 km s -1 ) within 2 ms, at which point it releases a radiative payload consisting of a shell pellet or well-defined microspheres. Initial experimental tests from the prototype system show attainment of over 600 m s -1 in about 1 ms. Finally, essential aspects of payload separation from the sabot and sabot capture have also been demonstrated at 200 m s -1 , and the method can be extended to over 2 km s -1 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1985

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, J. H.↗

NASA NTRS1988

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Saltz, Joel H.↗

NASA NTRS1973

Program for the solution of multipoint boundary value problems of quasilinear differential equations

Linear equations are solved by a method of superposition of solutions of a sequence of initial value problems. For nonlinear equations and/or boundary conditions, the solution is iterative and in each iteration a problem like the linear case is solved. A simple Taylor series expansion is used for the linearization of both nonlinear equations and nonlinear boundary conditions. The perturbation method of solution is used in preference to quasilinearization because of programming ease, and smaller storage requirements; and experiments indicate that the desired convergence properties exist although no proof or convergence is given.

Source record↗

NASA NTRS1987

A simultaneous viscous-inviscid interaction calculation procedure for transonic turbulent flows

A new simultaneous viscous-inviscid interaction scheme has been developed for the analysis of steady turbulent transonic separated flows. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity. The local solutions of the full potential and boundary-layer equations are treated simultaneously using the finite-difference method. The displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure of the space marching scheme. The Cebeci-Smith and Johnson-King models are used to simulate the turbulence. The computational examples showed that the simultaneous method is more efficient and robust than the semi-inverse method for transonic flows with a strong interaction.

Lee, D.↗

NASA NTRS1992

A nonlinear relaxation/quasi-Newton algorithm for the compressible Navier-Stokes equations

A highly efficient implicit method for the computation of steady, two-dimensional compressible Navier-Stokes flowfields is presented. The discretization of the governing equations is hybrid in nature, with flux-vector splitting utilized in the streamwise direction and central differences with flux-limited artificial dissipation used for the transverse fluxes. Line Jacobi relaxation is used to provide a suitable initial guess for a new nonlinear iteration strategy based on line Gauss-Seidel sweeps. The applicability of quasi-Newton methods as convergence accelerators for this and other line relaxation algorithms is discussed, and efficient implementations of such techniques are presented. Convergence histories and comparisons with experimental data are presented for supersonic flow over a flat plate and for several high-speed compression corner interactions. Results indicate a marked improvement in computational efficiency over more conventional upwind relaxation strategies, particularly for flowfields containing large pockets of streamwise subsonic flow.

Edwards, Jack R.↗

DOE OSTI2020

Novel strategies for modal-based structural material identification

Here, we present modal-based methods for model calibration in structural dynamics, and address several key challenges in the solution of gradient-based optimization problems with eigenvalues and eigenvectors, including the solution of singular Helmholtz problems encountered in sensitivity calculations, non-differentiable objective functions caused by mode swapping during optimization, and cases with repeated eigenvalues. Unlike previous literature that relied on direct solution of the eigenvector adjoint equations, we present a parallel iterative domain decomposition strategy (Adjoint Computation via Modal Superposition with Truncation Augmentation) for the solution of the singular Helmholtz problems. For problems with repeated eigenvalues we present a novel Mode Separation via Projection algorithm, and in order to address mode swapping between inverse iterations we present a novel Injective mode ordering metric. We present the implementation of these methods in a massively parallel finite element framework with the ability to use measured modal data to extract unknown structural model parameters from large complex problems. A series of increasingly complex numerical examples are presented that demonstrate the implementation and performance of the methods in a massively parallel finite element framework [7], [5], using gradient-based optimization techniques in the Rapid Optimization Library (ROL) [21].

36 MATERIALS SCIENCE↗

DOE OSTI2022

Asynchronous domain decomposition methods for nonlinear PDEs

One- and two-level parallel asynchronous methods for the numerical solution of nonlinear systems of equations, especially those arising from (nonlinear) partial differential equations, are studied. The proposed methods are based on domain decomposition techniques. Local convergence theorems are presented in several cases, with appropriate hypotheses. Computational results on a shared memory multiprocessor machine for various problems exhibiting nonlinearities are reported, illustrating the potential of these asynchronous methods, especially for heterogeneous clusters.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗

NASA NTRS1984

An algorithm for maximum likelihood estimation using an efficient method for approximating sensitivities

An algorithm for maximum likelihood (ML) estimation is developed primarily for multivariable dynamic systems. The algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). The method determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort compared with integrating the analytically determined sensitivity equations or using a finite-difference method. Different surface-fitting methods are discussed and demonstrated. Aircraft estimation problems are solved by using both simulated and real-flight data to compare MNRES with commonly used methods; in these solutions MNRES is found to be equally accurate and substantially faster. MNRES eliminates the need to derive sensitivity equations, thus producing a more generally applicable algorithm.

Murphy, P. C.↗

NASA NTRS2016

Deductive Evaluation: Formal Code Analysis With Low User Burden

We describe a framework for symbolically evaluating iterative C code using a deductive approach that automatically discovers and proves program properties. Although verification is not performed, the method can infer detailed program behavior. Software engineering work flows could be enhanced by this type of analysis. Floyd-Hoare verification principles are applied to synthesize loop invariants, using a library of iteration-specific deductive knowledge. When needed, theorem proving is interleaved with evaluation and performed on the fly. Evaluation results take the form of inferred expressions and type constraints for values of program variables. An implementation using PVS (Prototype Verification System) is presented along with results for sample C functions.

Di Vito, Ben. L↗

DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗