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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

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↗
NASA NTRS1972

Numerical method for the solution of large systems of differential equations of the boundary layer type

A numerical method for the solution of large systems of nonlinear differential equations of the boundary-layer type is described. The method is a modification of the technique for satisfying asymptotic boundary conditions. The present method employs inverse interpolation instead of the Newton method to adjust the initial conditions of the related initial-value problem. This eliminates the so-called perturbation equations. The elimination of the perturbation equations not only reduces the user's preliminary work in the application of the method, but also reduces the number of time-consuming initial-value problems to be numerically solved at each iteration. For further ease of application, the solution of the overdetermined system for the unknown initial conditions is obtained automatically by applying Golub's linear least-squares algorithm. The relative ease of application of the proposed numerical method increases directly as the order of the differential-equation system increases. Hence, the method is especially attractive for the solution of large-order systems. After the method is described, it is applied to a fifth-order problem from boundary-layer theory.

Green, M. J.↗
NASA NTRS1993

On eigensystem assignment with dissipativity constraints

Stability of the closed-loop system is of vital importance when controller synthesis is being performed using a reduced order model of the open-loop plant. For passive systems such as flexible space structures, closed-loop stability is often ensured by low-authority dissipative controllers using collocated and compatible sensors and actuators. This paper considers the problem of eigensystem assignment for flexible structures using constant-gain output feedback, under the dissipativity constraints requiring that the gain matrices be nonnegative definite. A novel method is developed in which a sequential procedure is used to assign one complex-conjugate pair of closed-loop eigenvalues in each step. At each step, an iterative technique is employed for generating a nonnegative definite gain matrix that assigns a new pair while keeping the previously assigned pairs unchanged. This method can assign up to m closed-loop system eigenvalues, where m is the number of control inputs. A numerical example of damping enhancement for a large flexible structure is presented to demonstrate the approach.

Maghami, Peiman G.↗
NASA NTRS1981

An analytical approach to airfoil icing

An analytical procedure has been developed to predict rime ice growth on unprotected airfoil sections and to evaluate the aerodynamic performance. A time stepping method is used in which: (1) water droplet trajectories are calculated, (2) a rime ice shape determined, (3) the flowfield around the iced airfoil is recalculated, and (4) the build-up process iterated upon until the desired icing time is reached. The performance of the iced airfoil shapes are then determined from existing analytic methods. Rime ice shapes determined in the NASA Lewis Icing Research Tunnel on a modified NACA 64 series airfoil agree well with the shapes predicted by the analytical method. Measured and predicted increases in drag due to the rime ice also agree favorably. A simplified scaling analysis is also presented and verified which provides the duplication of full scale results of rime ice accretions in small scale model tests.

Bragg, M. B.↗
NASA NTRS2012

Agile Development Methods for Space Operations

Main stream industry software development practice has gone from a traditional waterfall process to agile iterative development that allows for fast response to customer inputs and produces higher quality software at lower cost. How can we, the space ops community, adopt state of the art software development practice, achieve greater productivity at lower cost, and maintain safe and effective space flight operations? At NASA Ames, we are developing Mission Control Technologies Software, in collaboration with Johnson Space Center (JSC) and, more recently, the Jet Propulsion Laboratory (JPL).

Trimble, Jay↗
NASA NTRS2003

An Iterative Decambering Approach for Post-Stall Prediction of Wing Characteristics using known Section Data

An iterative decambering approach for the post stall prediction of wings using known section data as inputs is presented. The method can currently be used for incompressible .ow and can be extended to compressible subsonic .ow using Mach number correction schemes. A detailed discussion of past work on this topic is presented first. Next, an overview of the decambering approach is presented and is illustrated by applying the approach to the prediction of the two-dimensional C(sub l) and C(sub m) curves for an airfoil. The implementation of the approach for iterative decambering of wing sections is then discussed. A novel feature of the current e.ort is the use of a multidimensional Newton iteration for taking into consideration the coupling between the di.erent sections of the wing. The approach lends itself to implementation in a variety of finite-wing analysis methods such as lifting-line theory, discrete-vortex Weissinger's method, and vortex lattice codes. Results are presented for a rectangular wing for a from 0 to 25 deg. The results are compared for both increasing and decreasing directions of a, and they show that a hysteresis loop can be predicted for post-stall angles of attack.

Mukherjee, Rinku↗
NASA NTRS2014

Selection of Thermal Worst-Case Orbits via Modified Efficient Global Optimization

Efficient Global Optimization (EGO) was used to select orbits with worst-case hot and cold thermal environments for the Stratospheric Aerosol and Gas Experiment (SAGE) III. The SAGE III system thermal model changed substantially since the previous selection of worst-case orbits (which did not use the EGO method), so the selections were revised to ensure the worst cases are being captured. The EGO method consists of first conducting an initial set of parametric runs, generated with a space-filling Design of Experiments (DoE) method, then fitting a surrogate model to the data and searching for points of maximum Expected Improvement (EI) to conduct additional runs. The general EGO method was modified by using a multi-start optimizer to identify multiple new test points at each iteration. This modification facilitates parallel computing and decreases the burden of user interaction when the optimizer code is not integrated with the model. Thermal worst-case orbits for SAGE III were successfully identified and shown by direct comparison to be more severe than those identified in the previous selection. The EGO method is a useful tool for this application and can result in computational savings if the initial Design of Experiments (DoE) is selected appropriately.

Moeller, Timothy M.↗
NASA NTRS2009

Hybrid diversity method utilizing adaptive diversity function for recovering unknown aberrations in an optical system

A method of recovering unknown aberrations in an optical system includes collecting intensity data produced by the optical system, generating an initial estimate of a phase of the optical system, iteratively performing a phase retrieval on the intensity data to generate a phase estimate using an initial diversity function corresponding to the intensity data, generating a phase map from the phase retrieval phase estimate, decomposing the phase map to generate a decomposition vector, generating an updated diversity function by combining the initial diversity function with the decomposition vector, generating an updated estimate of the phase of the optical system by removing the initial diversity function from the phase map. The method may further include repeating the process beginning with iteratively performing a phase retrieval on the intensity data using the updated estimate of the phase of the optical system in place of the initial estimate of the phase of the optical system, and using the updated diversity function in place of the initial diversity function, until a predetermined convergence is achieved.

Dean, Bruce H.↗
NASA NTRS1992

A decomposition method for optimization of large-system reliability

This paper considers a decomposition approach for optimization of the reliability of a large system with a general network structure. A 3-level methodology is developed for optimal allocation of available resources among subsystems in order to ensure maximization of system-reliability. The decentralized nature of this methodology greatly reduces the complexity of the large problem and facilitates seeking the optimal solution. Two examples show that the complexity of a large system can be greatly reduced by solving several smaller-dimensional sub-problems iteratively. Sub-problems whose dimensions are small can be efficiently solved by any existing nonlinear programming method. Another important feature of this approach is the possible simplification of the objective function during the solution. This leads in some cases to an analytic solution for the lower-level optimization problems in a 3-level decomposition solution.

Li, Duan↗
NASA NTRS1982

An implicit formulation for the one-dimensional two-phase multi-interface Stefan Problem

In recent years there has been an intensified effort in the analysis of the Stefan Problem. Particular interest is directed towards determining the location of the interface between the solid and liquid phases of the substance being analyzed. A method has been developed using an implicit finite difference scheme which provides an accurate calculation of the phase change interface locations for a problem with simultaneous multiple phases. This is achieved through an iterative technique which requires the phase change energy to be correctly distributed within the nodes which are undergoing a phase transition. The new method predicts the interface locations without requiring time step size control or node size modification at the interfaces.

Williams, S. D.↗
NASA NTRS2006

On the Development of an Efficient Parallel Hybrid Solver with Application to Acoustically Treated Aero-Engine Nacelles

A finite element solution to the convected Helmholtz equation in a nonuniform flow is used to model the noise field within 3-D acoustically treated aero-engine nacelles. Options to select linear or cubic Hermite polynomial basis functions and isoparametric elements are included. However, the key feature of the method is a domain decomposition procedure that is based upon the inter-mixing of an iterative and a direct solve strategy for solving the discrete finite element equations. This procedure is optimized to take full advantage of sparsity and exploit the increased memory and parallel processing capability of modern computer architectures. Example computations are presented for the Langley Flow Impedance Test facility and a rectangular mapping of a full scale, generic aero-engine nacelle. The accuracy and parallel performance of this new solver are tested on both model problems using a supercomputer that contains hundreds of central processing units. Results show that the method gives extremely accurate attenuation predictions, achieves super-linear speedup over hundreds of CPUs, and solves upward of 25 million complex equations in a quarter of an hour.

Watson, Willie R.↗
NASA NTRS1983

Gravity, topography, and crustal evolution of Venus

A method for estimating radial gravity anomalies from line-of-sight acceleration data, their interpolation, and the use of iteration for improved radial anomaly estimates is outlined. Results based on spherical harmonic coefficients for Venus' gravitational field are first discussed, and the procedure is then presented for determining radial acceleration estimates from the line-of-sight vectors and for interpolating the results at preselected locations using the method of objective mapping. A preliminary gravity anomaly map of much of the surface of Venus based on the first estimation step of this technique is presented, and the results are discussed.

Bowin, C.↗
NASA NTRS1972

A finite element procedure for nonlinear prebuckling and initial postbuckling analysis

A procedure cast in a form appropriate to the finite element method is presented for geometrically nonlinear prebuckling and postbuckling structural analysis, including the identification of snap-through type of buckling. The principal features of this procedure are the use of direct iteration for solution of the nonlinear algebraic equations in the prebuckling range, an interpolation scheme for determination of the initial bifurcation point, a perturbation method in definition of the load-displacement behavior through the postbuckling regime, and extrapolation in determination of the limit point for snap-through buckling. Three numerical examples are presented in illustration of the procedure and in comparison with alternative approaches.

Mau, S. T.↗
NASA NTRS1974

Ambiguity resolution in systems using Omega for position location

The lane ambiguity problem prevents the utilization of the Omega system for many applications such as locating buoys and balloons. The method of multiple lines of position introduced herein uses signals from four or more Omega stations for ambiguity resolution. The coordinates of the candidate points are determined first through the use of the Newton iterative procedure. Subsequently, a likelihood function is generated for each point, and the ambiguity is resolved by selecting the most likely point. The method was tested through simulation.

Frenkel, G.↗
NASA NTRS1986

Force-free magnetic fields - The magneto-frictional method

The problem under discussion is that of calculating magnetic field configurations in which the Lorentz force j x B is everywhere zero, subject to specified boundary conditions. We choose to represent the magnetic field in terms of Clebsch variables in the form B = grad alpha x grad beta. These variables are constant on any field line so that each field line is labeled by the corresponding values of alpha and beta. When the field is described in this way, the most appropriate choice of boundary conditions is to specify the values of alpha and beta on the bounding surface. We show that such field configurations may be calculated by a magneto-frictional method. We imagine that the field lines move through a stationary medium, and that each element of magnetic field is subject to a frictional force parallel to and opposing the velocity of the field line. This concept leads to an iteration procedure for modifying the variables alpha and beta, that tends asymptotically towards the force-free state. We apply the method first to a simple problem in two rectangular dimensions, and then to a problem of cylindrical symmetry that was previously discussed by Barnes and Sturrock (1972). In one important respect, our new results differ from the earlier results of Barnes and Sturrock, and we conclude that the earlier article was in error.

Yang, W. H.↗
NASA NTRS1989

The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes

The present numerical method for the solution of unsteady, incompressible three-dimensional flow Navier-Stokes equations using velocity-vorticity variables and irregular Cartesian grids proceeds by solving: (1) equations of Cauchy-Riemann type for the velocity; and (2) transport-diffusion equations for the vorticity, whose solenoidal vorticity components are generated by a Poisson equation for an appropriate scalar potential. Iterations are used to solve the finite difference equations, facilitating the use of vector and parallel-computing methods; numerical experiments have verified the method's second-order spatial and temporal accuracy.

Gatski, T. B.↗
NASA NTRSDate not supplied

Performance Optimization Methods for a Memory-Bound, Unstructured-Grid CFD Application on Massively Parallel GPU Platforms

Computational performance of the FUN3D unstructured-grid computational fluid dynamics (CFD) application on massively parallel GPU environments is memory-bound and highly dependent upon efficient reads from and atomic updates to the irregular cell-, edge-, and node-based data structures. In this talk, we present recent efforts into optimizing select performance-critical kernels on NVIDIA Tesla V100 and A100 GPUs and AMD CDNA MI100 GPUs. A novel use of L2 cache residency controls and asynchronous loads into on-chip shared memory are explored on the A100 GPU for the sparse iterative solver, which is dominated by mixed-precision, sparse matrix vector multiplication. Demonstrations show that these methods improve global memory bandwidth utilization by 13.5% on the A100 GPU. Several techniques are also presented that use registers and/or shared memory to facilitate array transposition and aggregation which combine to reduce the frequency and increase the cache efficiency of floating-point atomic updates to the irregular data structures. These methods are demonstrated to improve the kernel throughput by nearly 500% on select kernels on the AMD MI100 over atomic updates directly to global memory. Overall, both V100 and A100 GPUs outperformed the MI100 GPU on kernels dominated by double-precision atomic updates; however, the techniques demonstrated here reduced the performance gap and improved the MI100 performance.

GPU CPU unstructured CFD memory↗
NASA NTRS1976

Free vibration analysis of spinning flexible space structures

Efficient computation of natural frequencies and associated free vibration modes of spinning flexible structures is required to accurately determine the nature of interaction between the flexible structure and the attitude control system, which is vital in relating control torques to attitude angles. While structural discretization is effected by the finite element method, the resulting eigenvalue problem is solved by a combined Sturm sequence and inverse iteration procedure that yields a few specified roots and associated vectors. An eigenvalue procedure, based on a simultaneous iteration technique, provides efficient computation of the first few roots and vectors; a modal synthesis procedure proves to be useful for eigenproblem solutions of unusual structures such as spacecraft.

Gupta, K. K.↗