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,135 records · Page 63

Kinetic theory model for the flow of a simple gas from a three-dimensional axisymmetric nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.

Riley, B. R.↗

Model reduction in integrated controls-structures design

It is the objective of this paper to present a model reduction technique developed for the integrated controls-structures design of flexible structures. Integrated controls-structures design problems are typically posed as nonlinear mathematical programming problems, where the design variables consist of both structural and control parameters. In the solution process, both structural and control design variables are constantly changing; therefore, the dynamic characteristics of the structure are also changing. This presents a problem in obtaining a reduced-order model for active control design and analysis which will be valid for all design points within the design space. In other words, the frequency and number of the significant modes of the structure (modes that should be included) may vary considerably throughout the design process. This is also true as the locations and/or masses of the sensors and actuators change. Moreover, since the number of design evaluations in the integrated design process could easily run into thousands, any feasible order-reduction method should not require model reduction analysis at every design iteration. In this paper a novel and efficient technique for model reduction in the integrated controls-structures design process, which addresses these issues, is presented.

Maghami, Peiman G.↗

Improved Boundary Conditions for Cell-centered Difference Schemes

Cell-centered finite-volume (CCFV) schemes have certain attractive properties for the solution of the equations governing compressible fluid flow. Among others, they provide a natural vehicle for specifying flux conditions at the boundaries of the physical domain. Unfortunately, they lead to slow convergence for numerical programs utilizing them. In this report a method for investigating and improving the convergence of CCFV schemes is presented, which focuses on the effect of the numerical boundary conditions. The key to the method is the computation of the spectral radius of the iteration matrix of the entire demoralized system of equations, not just of the interior point scheme or the boundary conditions.

VanderWijngaart, Rob F.↗

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel.

Atkins, H. L.↗

Design of a Mach-3 Nozzle for TBCC Testing in the NASA LaRC 8-ft High Temperature Tunnel

A new nozzle is being constructed for the NASA Langley Research Center 8-Foot High Temperature Tunnel. The axisymmetric nozzle was designed with a Mach-3 exit flow for testing Turbine-Based Combined-Cycle engines at a Mach number in the vicinity of the transition from turbojet to ramjet operation. The nozzle contour was designed using the NASA Langley IMOCND computer program which solves the potential equation using the classical method of characteristics. To include viscous effects, the design procedure iterated the MOC contour generation with CFD Navier-Stokes calculations, adjusting MOC input parameters until target nozzle-exit conditions were achieved in the Navier-Stokes calculations. The design process was complicated by a requirement to use the final 29.5 inches of an existing 54.5-inch exit-diameter Mach-5 nozzle contour. This was accomplished by generating a Mach-3 contour that matched the radius of the Mach-5 contour at the match point and using a 3rd order polynomial to create a smooth transition between the two contours. During the final evaluation of the design it was realized that the throat diameter is more than half that of the upstream mixing chamber. This led to the concern that large vortical structures generated in the mixer would persist downstream, affecting nozzle-exit flow. This concern was addressed by analyzing the results of three-dimensional, viscous, numerical simulations of the entire flowfield, from the exit of the facility combustor to the nozzle exit. An analysis of the solution indicated that large scale structures do not pass through the throat and that both the total temperature and species (CO2) are well mixed in the mixer, providing uniform flow to the nozzle and subsequently the test cabin.

Gaffney, Richard L., Jr.↗

Joint Carrier-Phase Synchronization and LDPC Decoding

A method has been proposed to increase the degree of synchronization of a radio receiver with the phase of a suppressed carrier signal modulated with a binary- phase-shift-keying (BPSK) or quaternary- phase-shift-keying (QPSK) signal representing a low-density parity-check (LDPC) code. This method is an extended version of the method described in Using LDPC Code Constraints to Aid Recovery of Symbol Timing (NPO-43112), NASA Tech Briefs, Vol. 32, No. 10 (October 2008), page 54. Both methods and the receiver architectures in which they would be implemented belong to a class of timing- recovery methods and corresponding receiver architectures characterized as pilotless in that they do not require transmission and reception of pilot signals. The proposed method calls for the use of what is known in the art as soft decision feedback to remove the modulation from a replica of the incoming signal prior to feeding this replica to a phase-locked loop (PLL) or other carrier-tracking stage in the receiver. Soft decision feedback refers to suitably processed versions of intermediate results of iterative computations involved in the LDPC decoding process. Unlike a related prior method in which hard decision feedback (the final sequence of decoded symbols) is used to remove the modulation, the proposed method does not require estimation of the decoder error probability. In a basic digital implementation of the proposed method, the incoming signal (having carrier phase theta theta (sub c) plus noise would first be converted to inphase (I) and quadrature (Q) baseband signals by mixing it with I and Q signals at the carrier frequency [wc/(2 pi)] generated by a local oscillator. The resulting demodulated signals would be processed through one-symbol-period integrate and- dump filters, the outputs of which would be sampled and held, then multiplied by a soft-decision version of the baseband modulated signal. The resulting I and Q products consist of terms proportional to the cosine and sine of the carrier phase cc as well as correlated noise components. These products would be fed as inputs to a digital PLL that would include a number-controlled oscillator (NCO), which provides an estimate of the carrier phase, theta(sub c).

Simon, Marvin↗

Computational methods for structural load and resistance modeling

An automated capability for computing structural reliability considering uncertainties in both load and resistance variables is presented. The computations are carried out using an automated Advanced Mean Value iteration algorithm (AMV +) with performance functions involving load and resistance variables obtained by both explicit and implicit methods. A complete description of the procedures used is given as well as several illustrative examples, verified by Monte Carlo Analysis. In particular, the computational methods described in the paper are shown to be quite accurate and efficient for a material nonlinear structure considering material damage as a function of several primitive random variables. The results show clearly the effectiveness of the algorithms for computing the reliability of large-scale structural systems with a maximum number of resolutions.

Thacker, B. H.↗

Fast Multigrid Techniques in Total Variation-Based Image Reconstruction

Existing multigrid techniques are used to effect an efficient method for reconstructing an image from noisy, blurred data. Total Variation minimization yields a nonlinear integro-differential equation which, when discretized using cell-centered finite differences, yields a full matrix equation. A fixed point iteration is applied with the intermediate matrix equations solved via a preconditioned conjugate gradient method which utilizes multi-level quadrature (due to Brandt and Lubrecht) to apply the integral operator and a multigrid scheme (due to Ewing and Shen) to invert the differential operator. With effective preconditioning, the method presented seems to require Omicron(n) operations. Numerical results are given for a two-dimensional example.

Oman, Mary Ellen↗

Primary electron spectrometer, 18:63 UE: Electrostatic analyzer description and energy spectrum determination

The primary electron spectrometer used to detect auroral electrons on sounding rocket 18:63 UE is described. The spectrometer used exponentially decaying positive and negative voltages applied to spherical deflection plates for energy analysis. A method for determining the analyzer response which does not require the assumptions that the ratio of plate separation to mean radius, the entrance or the exit apertures are small is described. By comparison with experiment it is shown that the effect of neither entrance nor exit collimation can be ignored. The experimental and calculated values of the limiting orbits agree well. A non-iterative technique of unfolding the electron differential energy spectrum is described. This method does not require the usual assumption of a flat or histogram-type energy spectrum. The unfolded spectra using both this technique and one which assumes a flat spectrum are compared to actual input spectra. This technique is especially useful in analyzing peaked auroral electron energy spectra.

Pongratz, M. B.↗

Application of an Upwind High Resolution Finite-Differencing Scheme and Multigrid Method in Steady-State Incompressible Flow Simulations

The analysis and design of a submarine propulsor requires the ability to predict the characteristics of both laminar and turbulent flows to a higher degree of accuracy. This report presents results of certain benchmark computations based on an upwind, high-resolution, finite-differencing Navier-Stokes solver. The purpose of the computations is to evaluate the ability, the accuracy and the performance of the solver in the simulation of detailed features of viscous flows. Features of interest include flow separation and reattachment, surface pressure and skin friction distributions. Those features are particularly relevant to the propulsor analysis. Test cases with a wide range of Reynolds numbers are selected; therefore, the effects of the convective and the diffusive terms of the solver can be evaluated separately. Test cases include flows over bluff bodies, such as circular cylinders and spheres, at various low Reynolds numbers, flows over a flat plate with and without turbulence effects, and turbulent flows over axisymmetric bodies with and without propulsor effects. Finally, to enhance the iterative solution procedure, a full approximation scheme V-cycle multigrid method is implemented. Preliminary results indicate that the method significantly reduces the computational effort.

NAVIER STOKES EQUATIONS↗

Some Advanced Concepts in Discrete Aerodynamic Sensitivity Analysis

An efficient incremental-iterative approach for differentiating advanced flow codes is successfully demonstrated on a 2D inviscid model problem. The method employs the reverse-mode capability of the automatic- differentiation software tool ADIFOR 3.0, and is proven to yield accurate first-order aerodynamic sensitivity derivatives. A substantial reduction in CPU time and computer memory is demonstrated in comparison with results from a straight-forward, black-box reverse- mode application of ADIFOR 3.0 to the same flow code. An ADIFOR-assisted procedure for accurate second-order aerodynamic sensitivity derivatives is successfully verified on an inviscid transonic lifting airfoil example problem. The method requires that first-order derivatives are calculated first using both the forward (direct) and reverse (adjoint) procedures; then, a very efficient non-iterative calculation of all second-order derivatives can be accomplished. Accurate second derivatives (i.e., the complete Hessian matrices) of lift, wave-drag, and pitching-moment coefficients are calculated with respect to geometric- shape, angle-of-attack, and freestream Mach number

Taylor, Arthur C., III↗

Power System Optimization Codes Modified

A major modification of and addition to existing Closed Brayton Cycle (CBC) space power system optimization codes was completed. These modifications relate to the global minimum mass search driver programs containing three nested iteration loops comprising iterations on cycle temperature ratio, and three separate pressure ratio iteration loops--one loop for maximizing thermodynamic efficiency, one for minimizing radiator area, and a final loop for minimizing overall power system mass. Using the method of steepest ascent, the code sweeps through the pressure ratio space repeatedly, each time with smaller iteration step sizes, so that the three optimum pressure ratios can be obtained to any desired accuracy for each of the objective functions referred to above (i.e., maximum thermodynamic efficiency, minimum radiator area, and minimum system mass). Two separate options for the power system heat source are available: 1. A nuclear fission reactor can be used. It is provided with a radiation shield 1. (composed of a lithium hydride (LiH) neutron shield and tungsten (W) gamma shield). Suboptions can be used to select the type of reactor (i.e., fast spectrum liquid metal cooled or epithermal high-temperature gas reactor (HTGR)). 2. A solar heat source can be used. This option includes a parabolic concentrator and heat receiver for raising the temperature of the recirculating working fluid. A useful feature of the code modifications is that key cycle parameters are displayed, including the overall system specific mass in kilograms per kilowatt and the system specific power in watts per kilogram, as the results for each temperature ratio are computed. As the minimum mass temperature ratio is encountered, a message is printed out. Several levels of detailed information on cycle state points, subsystem mass results, and radiator temperature profiles are stored for this temperature ratio condition and can be displayed or printed by users.

Juhasz, Albert J.↗

Construction of surface pressure field from scatterometer wind field

An account of the construction of surface pressure fields from Seasat-A satellite scatterometer (SASS) winds as carried out by different methods, and the comparison of these pressure fields with those derived from in situ ship observations is presented. On the assumption that the pressure adjusts itself instantaneously to the motion field, it may be computed by various methods. One of these makes use of planetary boundary theory, and of the possible techniques in this category a two-layer iterative scheme admitting of the parametrization of diabatic and baroclinic effects and of secondary flow was chosen. A second method involves the assumption of zero two-dimensional divergence, leading to a Laplace's equation (the balance equation) in pressure, with the wind field serving as a forcing function. This method does not accommodate adiabatic or baroclinic effects, and requires a knowledge of the pressure at all boundary points. Two comparison fields are used for validation: the conventional operational analyses of the US National Meteorological Center (NMC), and the special analyses of the Gulf of Alaska Experiment (GOASEX), which were done by hand. The results of the computations were as follows: (1) The pressure fields, as computed from the SASS winds alone, closely approximated the NMC fields in regions where reasonable in situ coverage was available (typically, one or two mb differences over most of the chart, three to four mb in extreme cases); (2) In some cases the SASS-derived pressure fields displayed high-resolution phenomena not detected by the NMC fields, but evident in the GOASEX data; and, (3) As expected, the pressure fields derived from the balance equation were much smoother and less well resolved than the SASS-derived or NMC fields. The divergence as measured from the SASS winds is smaller than, but of the same order of magnitude as, the vorticity.

Wurtele, Morton G.↗

DYNGEN: A program for calculating steady-state and transient performance of turbojet and turbofan engines

The DYNGEN, a digital computer program for analyzing the steady state and transient performance of turbojet and turbofan engines, is described. The DYNGEN is based on earlier computer codes (SMOTE, GENENG, and GENENG 2) which are capable of calculating the steady state performance of turbojet and turbofan engines at design and off-design operating conditions. The DYNGEN has the combined capabilities of GENENG and GENENG 2 for calculating steady state performance; to these the further capability for calculating transient performance was added. The DYNGEN can be used to analyze one- and two-spool turbojet engines or two- and three-spool turbofan engines without modification to the basic program. A modified Euler method is used by DYNGEN to solve the differential equations which model the dynamics of the engine. This new method frees the programmer from having to minimize the number of equations which require iterative solution. As a result, some of the approximations normally used in transient engine simulations can be eliminated. This tends to produce better agreement when answers are compared with those from purely steady state simulations. The modified Euler method also permits the user to specify large time steps (about 0.10 sec) to be used in the solution of the differential equations. This saves computer execution time when long transients are run. Examples of the use of the program are included, and program results are compared with those from an existing hybrid-computer simulation of a two-spool turbofan.

Sellers, J. F.↗

Best Merge Region Growing with Integrated Probabilistic Classification for Hyperspectral Imagery

A new method for spectral-spatial classification of hyperspectral images is proposed. The method is based on the integration of probabilistic classification within the hierarchical best merge region growing algorithm. For this purpose, preliminary probabilistic support vector machines classification is performed. Then, hierarchical step-wise optimization algorithm is applied, by iteratively merging regions with the smallest Dissimilarity Criterion (DC). The main novelty of this method consists in defining a DC between regions as a function of region statistical and geometrical features along with classification probabilities. Experimental results are presented on a 200-band AVIRIS image of the Northwestern Indiana s vegetation area and compared with those obtained by recently proposed spectral-spatial classification techniques. The proposed method improves classification accuracies when compared to other classification approaches.

Tarabalka, Yuliya↗

Matching Inviscid/Boundary Layer Flowfields

New boundary-layer equations are developed and the solutions match all the boundary-layer properties, except for the normal velocity, exactly with the corresponding inviscid properties. The numerical procedure solves tridiagonal matrices at each marching station. As part of the solution an inviscid transpiration velocity at the surface is calculated from the boundary-layer solution. This transpiration velocity could be used as a boundary condition to calculate a new inviscid solution. If the inviscid/boundary layer solutions are iterated, then the normal velocity from the boundary layer solution will match the inviscid values exactly. Solutions are calculated for shear flows over a flat plate. Results from the present method compared well with Navier-Stokes solutions for incompressible constant shear flow and sinusoidal shear flow. Iterations of the inviscid/boundary layer solutions were not necessary. Compressible flow of an exponential shear at Mach 4 was calculated for cold, adiabatic, and hot walls. It was found that inviscid shear flows have a significant effect on skin friction and heating rates. The present method should yield accurate boundary-layer solutions at Reynolds numbers lower than the traditional boundary-layer equations.

DeJarnette, Fred R.↗

Iterative solution of transonic flows over airfoils and wings, including flows at Mach 1

A new method of calculating transonic flows based on a 'rotated' difference scheme is described. It is suitable for the calculation of both two- and three-dimensional flows without restriction on the speed at infinity and is well adapted to computer use. The Murman procedure is modified to eliminate any assumptions about the direction of flow when constructing the difference scheme. The proper directional property is obtained by rotating the difference scheme to conform with the local stream direction. In the hyperbolic region retarded difference formulas are used for all contributions to the streamwise second derivative, producing a correctly oriented positive artificial viscosity. In the absence of a simple implicit scheme in the hyperbolic and elliptic regions, the concept of iterations as steps in artificial time is introduced. Computer testing of this procedure provides numerical confirmation of the existence and uniqueness of weak solutions of the potential equation when a suitable entropy inequality is enforced.

Jameson, A.↗

Neural network training by integration of adjoint systems of equations forward in time

A method and apparatus for supervised neural learning of time dependent trajectories exploits the concepts of adjoint operators to enable computation of the gradient of an objective functional with respect to the various parameters of the network architecture in a highly efficient manner. Specifically, it combines the advantage of dramatic reductions in computational complexity inherent in adjoint methods with the ability to solve two adjoint systems of equations together forward in time. Not only is a large amount of computation and storage saved, but the handling of real-time applications becomes also possible. The invention has been applied it to two examples of representative complexity which have recently been analyzed in the open literature and demonstrated that a circular trajectory can be learned in approximately 200 iterations compared to the 12000 reported in the literature. A figure eight trajectory was achieved in under 500 iterations compared to 20000 previously required. The trajectories computed using our new method are much closer to the target trajectories than was reported in previous studies.

Toomarian, Nikzad↗