Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative 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 91 records · Page 5

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗

Comparison of Electrical Output Format Options for the Analysis of Strain-Gage Balance Calibration Data

Fundamental characteristics of three gage output formats are discussed that may be used to both interpret and analyze wind tunnel strain-gage balance data. The first output format uses raw outputs, i.e., absolute voltage measurements, for the analysis. This choice requires an intercept term in the regression model of the outputs if the Iterative Method is chosen for the load prediction. Raw outputs can be used for the analysis of all known output characteristics as long as the Iterative Method is applied. However, raw outputs cannot be used to analyze data of a balance with bi-directional outputs if the Non-Iterative Method is chosen for the analysis. The second output format, i.e., difference type 1, uses the difference between raw outputs and the natural zeros of the balance gages for the analysis. In this case, the use of the intercept term becomes optional. Analysis results for difference type 1 will exactly match analysis results obtained by using raw outputs as long as identical math models are used for the regression analysis. In addition, difference type 1 may be used to analyze data of a balance with bi directional outputs if an analyst prefers to process data using the Non-Iterative Method. The third output format, i.e., difference type 2, uses the difference between raw outputs and the outputs of the zero load point of a load series for the data analysis. The application of this format is only recommended if the maximum magnitude of all tare loads of a given balance calibration data set is less than two percent of capacity. Data from the calibration of a force balance is used to illustrate the application of the three gage output formats.

Wind Tunnel Strain-Gage Balance↗

Improved Convergence and Robustness of USM3D Solutions on Mixed-Element Grids

Several improvements to the mixed-elementUSM3Ddiscretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Method, has been developed and implemented. The Hierarchical Adaptive Nonlinear Iteration Method provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier-Stokes equations and a nonlinear control of the solution update. Two variants of the Hierarchical Adaptive Nonlinear Iteration Method are assessed on four benchmark cases, namely, a zero-pressure-gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the preconditioner-alone method representing the baseline solver technology.

Pandya, Mohagna J.↗

Improved Convergence and Robustness of USM3D Solutions on Mixed-Element Grids

Several improvements to the mixed-element USM3D discretization and defect-correction schemes have been made. A new methodology for nonlinear iterations, called the Hierarchical Adaptive Nonlinear Iteration Method, has been developed and implemented. The Hierarchical Adaptive Nonlinear Iteration Method provides two additional hierarchies around a simple and approximate preconditioner of USM3D. The hierarchies are a matrix-free linear solver for the exact linearization of Reynolds-averaged Navier-Stokes equations and a nonlinear control of the solution update. Two variants of the Hierarchical Adaptive Nonlinear Iteration Method are assessed on four benchmark cases, namely, a zero-pressure-gradient flat plate, a bump-in-channel configuration, the NACA 0012 airfoil, and a NASA Common Research Model configuration. The new methodology provides a convergence acceleration factor of 1.4 to 13 over the preconditioner-alone method representing the baseline solver technology.

Pandya, Mohagna J.↗

Analytic Solution to the Problem of Aircraft Electric Field Mill Calibration

It is by no means a simple task to retrieve storm electric fields from an aircraft instrumented with electric field mill sensors. The presence of the aircraft distorts the ambient field in a complicated way. Before retrievals of the storm field can be made, the field mill measurement system must be "calibrated". In other words, a relationship between impressed (i.e., ambient) electric field and mill output must be established. If this relationship can be determined, it is mathematically inverted so that ambient field can be inferred from the mill outputs. Previous studies have primarily focused on linear theories where the relationship between ambient field and mill output is described by a "calibration matrix" M. Each element of the matrix describes how a particular component of the ambient field is enhanced by the aircraft. For example the product M(sub ix), E(sub x), is the contribution of the E(sub x) field to the i(th) mill output. Similarly, net aircraft charge (described by a "charge field component" E(sub q)) contributes an amount M(sub iq)E(sub q) to the output of the i(th) sensor. The central difficulty in obtaining M stems from the fact that the impressed field (E(sub x), E(sub y), E(sub z), E(sub q) is not known but is instead estimated. Typically, the aircraft is flown through a series of roll and pitch maneuvers in fair weather, and the values of the fair weather field and aircraft charge are estimated at each point along the aircraft trajectory. These initial estimates are often highly inadequate, but several investigators have improved the estimates by implementing various (ad hoc) iterative methods. Unfortunately, none of the iterative methods guarantee absolute convergence to correct values (i.e., absolute convergence to correct values has not been rigorously proven). In this work, the mathematical problem is solved directly by analytic means. For m mills installed on an arbitrary aircraft, it is shown that it is possible to solve for a single 2m-vector that provides all other needed variables (i.e., the unknown fair weather field, the unknown aircraft charge, and the unknown matrix M). Numerical tests of the solution, effects of measurement errors, and studies of solution non-uniqueness are ongoing as of this writing.

Koshak, William↗

A Universal Algorithm for the Detection of Bi-directional Gage Output Characteristics

A universal algorithm was developed that may be used to assess the bi-directional characteristics of the gage outputs of a wind tunnel strain-gage balance. The algorithm assumes that balance loads and gage outputs are described in the design format of the balance. It can also be applied to balance calibration data that is processed by using either the Iterative Method or the Non-Iterative Method. The algorithm uses an estimate of the bi-directional part of a gage output at load capacity as input. In addition, the statistical significance of the principle absolute value term in the regression model of either the gage output or the related primary load component is determined. A gage output is assumed to be bi-directional if two conditions are fulfilled: the bi-directional part of the output at load capacity exceeds 0.5 percent of the maximum output at load capacity; the p-value of the principle absolute value term of the regression model of the balance data is less than the threshold of 0.001. Data from the calibration of two six-component force balances and one five-component semi-span balance are used to illustrate the application of the universal detection algorithm.

wind tunnel test↗

A parallel iterative solution method for systems of nonlinear hyperbolic equations

An iterative algorithm suitable for the solution of a system of nonlinear hyperbolic partial differentiation equations in multiple dimensions is discussed. Current numerical methods for systems of nonlinear PDEs have limited parallelism due to strong coupling between the equations. This method decouples the PDEs by linearizing the convention coefficient for a space-time domain. This provides large grain parallelism. The linearization also allows the treatment of some terms in the equations as source terms, providing more freedom to choose from a wider variety of numerical methods. Smaller grain parallelism may be exploited within the solves for each equation. Thus, the method has potential for parallelism at several levels.

Scroggs, Jeffrey S.↗

Iterative PNS method for attached flows with upstream influence

A stable global-iteration procedure is developed by utilizing successive sweeps, from inflow to outflow boundaries, of a parabolized Navier-Stokes code for attached, steady supersonic flow. It is shown that the procedure converges in about ten or fewer iterations, and allows for the upstream influence within the subsonic region of a supersonic boundary layer. An implicit forward difference is employed in the subsonic region to evaluate the pressure-gradient term in the streamwise momentum equation. The pressure-gradient term is normally approximated in standard single-sweep parabolized methods, suppressing the upstream influence. The iterative results obtained demonstrate the validity of the single-sweep method for weak interactions. For hypersonic viscous interaction on a flat plate at a Mach number of 5.8, it is found that the single-sweep method has a small error that vanishes as the Reynolds number of the flow is increased. However, the skin friction has an error of about 10 percent for low Reynolds numbers and for hot wall conditions. The application of this method to laminar two-dimensional flow over weak expansion and compression corners shows that the pressure and skin friction results in the vicinity of the expansion corner compare favorably with a time-dependent Navier-Stokes numerical solution.

Rakich, J. V.↗

A Study of Morrison's Iterative Noise Removal Method

Morrison's iterative noise removal method is studied by characterizing its effect upon systems of differing noise level and response function. The nature of data acquired from a linear shift invariant instrument is discussed so as to define the relationship between the input signal, the instrument response function, and the output signal. Fourier analysis is introduced, along with several pertinent theorems, as a tool to more thorough understanding of the nature of and difficulties with deconvolution. In relation to such difficulties the necessity of a noise removal process is discussed. Morrison's iterative noise removal method and the restrictions upon its application are developed. The nature of permissible response functions is discussed, as is the choice of the response functions used.

Ioup, G. E.↗

Finite element solvers for incompressible fluid flows and heat transfer

Two different finite-element solvers for incompressible viscous flow, i.e., the mixed interpolation method and the SIMPLE-type iterative method, are compared and tested with some benchmark problems. The advantages of the SIMPLE-type iterative method are the decoupling of the governing equations and the use of equal-order interpolation functions for both velocity and pressure. Even though there is a significant difference between the two methods in terms of the pressure field, similar solutions are obtained for the velocity field.

Sohn, Jeong L.↗

Application of a Temperature-Dependent Load Prediction Method to a RUAG Six-Component Block-Type Balance

Temperature-dependent data of a RUAG six-component block-type balance was analyzed to assess the accuracy of two load prediction methods for temperature-dependent balance data. The supplied data was prepared for the analysis by splitting it into calibration and check load data subsets. The first calibration data subset was obtained at a temperature of 294 Kelvin. The second calibration data subset was obtained at a temperature of 315 Kelvin. A subset of 38 points was extracted from the second data set and used as check loads so that the accuracy of the two load prediction methods could be tested. First, the Iterative Method in combination with an extended independent and dependent variable set was used for the balance load prediction. This approach fits electrical outputs as a function of loads and the temperature and, afterwards, constructs a load iteration scheme from the regression coefficients so that loads can be predicted from outputs and the temperature during a wind tunnel test. The Non-Iterative Method was also used for the load prediction. This alternate method can more easily be implemented in a data system as loads are directly fitted as a function of electrical outputs and the temperature. Analysis results for the axial force are only discussed in the paper as similar results were obtained for the other five load components. Results for both methods clearly show that the cross-product term constructed from either a primary gage load or a primary gage output and the temperature explains the majority of the temperature-dependent part of the predicted balance load. This term models the temperature dependent nature of the gage sensitivity. Therefore, it is recommended to apply primary gage loadings at different temperatures during a balance calibration whenever temperature effects need to be described. These loadings will contain information about the temperature-dependent nature of the gage sensitivities that can be quantified by related cross-product terms in regression models of the data.

strain-gage balance↗

Wind Tunnel Strain-Gage Balance Calibration Data Analysis Using a Weighted Least Squares Approach

A new approach is presented that uses a weighted least squares fit to analyze wind tunnel strain-gage balance calibration data. The weighted least squares fit is specifically designed to increase the influence of single-component loadings during the regression analysis. The weighted least squares fit also reduces the impact of calibration load schedule asymmetries on the predicted primary sensitivities of the balance gages. A weighting factor between zero and one is assigned to each calibration data point that depends on a simple count of its intentionally loaded load components or gages. The greater the number of a data point's intentionally loaded load components or gages is, the smaller its weighting factor becomes. The proposed approach is applicable to both the Iterative and Non-Iterative Methods that are used for the analysis of strain-gage balance calibration data in the aerospace testing community. The Iterative Method uses a reasonable estimate of the tare corrected load set as input for the determination of the weighting factors. The Non-Iterative Method, on the other hand, uses gage output differences relative to the natural zeros as input for the determination of the weighting factors. Machine calibration data of a six-component force balance is used to illustrate benefits of the proposed weighted least squares fit. In addition, a detailed derivation of the PRESS residuals associated with a weighted least squares fit is given in the appendices of the paper as this information could not be found in the literature. These PRESS residuals may be needed to evaluate the predictive capabilities of the final regression models that result from a weighted least squares fit of the balance calibration data.

calibration analysis↗

Eigenvector determination by iterative optical methods

Three power algorithms are considered by which an iterative optical processor can be used to compute the eigenvalues and eigenvectors of a matrix. The algorithms are appropriate for three applications: (1) calculating the largest eigenvalue; (2) calculating the eigenvalues in order of decreasingly dominant eigenvalues; and (3) calculating the eigenvector whose eigenvalue is closest to a given value.

Vijaya Kumar, B. V. K.↗

Radiation pattern synthesis of planar antennas using the iterative sampling method

A synthesis method is presented for determining an excitation of an arbitrary (but fixed) planar source configuration. The desired radiation pattern is specified over all or part of the visible region. It may have multiple and/or shaped main beams with low sidelobes. The iterative sampling method is used to find an excitation of the source which yields a radiation pattern that approximates the desired pattern to within a specified tolerance. In this paper the method is used to calculate excitations for line sources, linear arrays (equally and unequally spaced), rectangular apertures, rectangular arrays (arbitrary spacing grid), and circular apertures. Examples using these sources to form patterns with shaped main beams, multiple main beams, shaped sidelobe levels, and combinations thereof are given.

Stutzman, W. L.↗

Assessment of Temperature-Dependent Regression Model Terms of a RUAG Six-Component Block-Type Balance

A metric called the percent contribution was applied to regression models of temperature-dependent calibration data of a RUAG six-component block-type balance in order to assess the influence of temperature-dependent regression model terms on the balance load prediction. Regression models were examined that are needed if either the Iterative or the Non-Iterative Method is used for the load prediction. Computed values of the percent contribution confirmed that the cross-product term defined by a primary load and the temperature difference is the most influential temperature-dependent term of the regression model of a primary output that the Iterative Method needs. Similarly, the analysis showed that the cross-product term defined by a primary output and the temperature difference is the most influential temperature-dependent term of the regression model of a primary load that the Non-Iterative Method needs. Computed results support conclusions that were reported in an earlier theoretical study. This study asserted that the cross-product term defined by a primary load or output and the temperature difference models the temperature-dependent shift of the gage sensitivity. The influence of other temperature-dependent terms used in the regression models of the calibration data of RUAG's balance was negligible. This observation may be explained by the fact that RUAG's block-type balances have highly linear characteristics. Overall, the percent contribution has proven itself to be a reliable and easy-to-implement metric that may also be used for the assessment of the influence of temperature-dependent regression model terms on the load prediction of a six-component strain-gage balance.

Ulbrich, N.↗

Upwind relaxation methods for the Navier-Stokes equations using inner iterations

A supersonic and subsonic test problem is solved using an upwind line relaxation algorithm for the Navier-Stokes equations which employs inner iterations. In the supersonic test case, the inner iteration method was found to require only about 65 percent of the CPU time which was required by the most efficient line relaxation method without inner iterations. In the subsonic test case, poor matrix conditioning resulted in less efficiency than line relaxation methods employing a more conventional CPU-saving strategy.

Taylor, Arthur C., III↗

Discrete fourier transform (DFT) analysis for applications using iterative transform methods

According to various embodiments, a method is provided for determining aberration data for an optical system. The method comprises collecting a data signal, and generating a pre-transformation algorithm. The data is pre-transformed by multiplying the data with the pre-transformation algorithm. A discrete Fourier transform of the pre-transformed data is performed in an iterative loop. The method further comprises back-transforming the data to generate aberration data.

Dean, Bruce H.↗

Development of parallel algorithms for electrical power management in space applications

The application of parallel techniques for electrical power system analysis is discussed. The Newton-Raphson method of load flow analysis was used along with the decomposition-coordination technique to perform load flow analysis. The decomposition-coordination technique enables tasks to be performed in parallel by partitioning the electrical power system into independent local problems. Each independent local problem represents a portion of the total electrical power system on which a loan flow analysis can be performed. The load flow analysis is performed on these partitioned elements by using the Newton-Raphson load flow method. These independent local problems will produce results for voltage and power which can then be passed to the coordinator portion of the solution procedure. The coordinator problem uses the results of the local problems to determine if any correction is needed on the local problems. The coordinator problem is also solved by an iterative method much like the local problem. The iterative method for the coordination problem will also be the Newton-Raphson method. Therefore, each iteration at the coordination level will result in new values for the local problems. The local problems will have to be solved again along with the coordinator problem until some convergence conditions are met.

Berry, Frederick C.↗