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,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Employing Sensitivity Derivatives to Estimate Uncertainty Propagation in CFD

Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.

Putko, Michele M.↗
NASA NTRS2012

Deconvolution Methods and Systems for the Mapping of Acoustic Sources from Phased Microphone Arrays

Mapping coherent/incoherent acoustic sources as determined from a phased microphone array. A linear configuration of equations and unknowns are formed by accounting for a reciprocal influence of one or more cross-beamforming characteristics thereof at varying grid locations among the plurality of grid locations. An equation derived from the linear configuration of equations and unknowns can then be iteratively determined. The equation can be attained by the solution requirement of a constraint equivalent to the physical assumption that the coherent sources have only in phase coherence. The size of the problem may then be reduced using zoning methods. An optimized noise source distribution is then generated over an identified aeroacoustic source region associated with a phased microphone array (microphones arranged in an optimized grid pattern including a plurality of grid locations) in order to compile an output presentation thereof, thereby removing beamforming characteristics from the resulting output presentation.

Brooks, Thomas F.↗
NASA NTRS1974

Multiple scattering from finite inhomogeneous media

Utilizing the characteristic information concerning the apparent phase constant difference between the electric and magnetic fields propagating in an inhomogeneous medium, a theoretical analysis of the multiple scattering of electromagnetic waves in finite inhomogeneous media is presented. The solution is obtained by first approximating the coefficients of a pair of exact coupled first-order differential equations and then solving the equations by first-order iteration. The present first-order approximate solution with multiple scattering considerations is shown to be more accurate than the WKB solution. Methods to improve the accuracy of the first-order solution further are discussed. Application of the solution to slowly varying finite media with periodic properties demonstrates the validity of the solution. The same approach can be extended to frequencies in the optical region by retaining additional terms in the coefficients of the coupled differential equations.

Tang, C. C. H.↗
NASA NTRS1977

A reexamination of two-fluid solar wind models

The two-fluid solar-wind equations have been solved by a method which is approximately 50 times faster than any previously developed, through the use of asymptotic expansions which are self-consistently iterated upon to find a solution that passes through the critical point. The energy assumptions in two-fluid solar-wind models are reexamined, and the conclusions are as follows: (1) proton thermal conduction may not be neglected, (2) the Coulomb logarithm must be calculated as a function of radius, and (3) the electron and proton temperatures at the base need not be equal, even when the time scale for energy exchange between the species is an order of magnitude smaller than the expansion time at the base. It is possible to reproduce reasonable quiet-time solar-wind parameters at 1 AU, but only if the proton temperature is approximately twice the electron temperature at 1 solar radius. This may indicate that extended proton heating is important in the outer solar corona. Winds with velocities at 1 AU of 450 km/s are generated without nonthermal energy deposition but require high proton temperatures as well as very low densities at the base. Higher-velocity solutions are not possible in a spherically symmetric geometry for reasonable particle fluxes at 1 AU, and it is suggested that these higher-velocity states probably require additional heating, acceleration mechanisms, or nonradial flow.

Nerney, S.↗
NASA NTRS1981

Dynamic stability of a rotor blade using finite element analysis

The aeroelastic stability of flap bending, lead-lag bending, and torsion of a helicopter rotor blade in hover is examined using a finite element formulation based on the principle of virtual work. Quasi-steady two-dimensional airfoil theory is used to obtain the aerodynamic loads. The rotor blade is discretized into beam elements, each with ten modal degrees of freedom. The resulting nonlinear equations of motion are solved for steady-state blade deflections through an iterative procedure. The flutter solution is calculated assuming blade motion to be a small perturbation about the steady solution. The normal mode method based on the coupled rotating natural modes about the steady deflections is used to reduce the number of equations in the flutter eigenanalysis. Results are presented for hingeless and articulated rotor blade configurations.

Sivaneri, N. T.↗
NASA NTRS1981

The Stokes problem for the ellipsoid using ellipsoidal kernels

A brief review of Stokes' problem for the ellipsoid as a reference surface is given. Another solution of the problem using an ellipsoidal kernel, which represents an iterative form of Stokes' integral, is suggested with a relative error of the order of the flattening. On studying of Rapp's method in detail the procedures of improving its convergence are discussed.

Zhu, Z.↗
NASA NTRS1984

Radiative-convective equilibrium models of Jupiter and Saturn

Radiative-convective equilibrium models for Jupiter and Saturn have been produced in a study concentrating on the stratospheric energy balance and the possible role of aerosol heating. These models are compared directly with the thermal structure profiles obtained from Voyager radio occultation measurements. The method is based on a straightforward flux divergence formulation derived from earlier work. The balance between absorbed and emitted energies is computed iteratively at each level in the atmosphere, assuming local thermodynamic equilibrium and employing a standard treatment of opacities. Results for Jupiter indicate that a dust-free model furnishes a good mean thermal profile for the stratosphere when compared with the Voyager 1 radio occultation measurements. Observations of the equatorial region exhibit periodic vertical structure. The Saturn models are relatively simple and in good agreement with the Voyager 2 radio occultation temperature profiles at all levels. Aerosol heating played a minor role in Saturn's midlatitude stratospheric energy balance at the time of the Voyager 2 encounter.

Appleby, J. F.↗
NASA NTRS1985

Computation of transonic vortex flows past delta wings Integral equation approach

The steady full-potential equation is written in the form of Poisson's equation, and the solution of the velocity field is expressed in terms of an integral equation. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a volume around the wing and its free-vortex sheets. The solution is obtained through successive iteration cycles. The source distribution is computed by using a mixed finite-difference scheme of the Murman-Cole type. The method is applied to delta wings. Numerical examples show that a conical shock is captured on the suction side of the wing. It is attached to the lower surface of the leading-edge vortex but does not necessarily reach to the wing surface.

Kandil, O. A.↗
NASA NTRS1990

Static aeroelastic analysis of composite wing

A static aeroelastic analysis capability that can predict aerodynamic loads for the deformed shape of the composite wing has been developed. The finite element method (FEM) was used for composite plate structural analysis, and the linear vortex lattice method (VLM) was used for steady aerodynamic analysis. The final deformed shape of the wing due to applied forces is determined by iterative manner using FEM and VLM. FEM and VLM analysis are related by a surface spline interpolation procedure. The wing with Gr/Ep composite material has been investigated to see the wing deformation effect. Aerodynamic load change due to wing flexibility has been investigated. Also, the effect of fiber orientation and sweep angle on the deformation pattern and aerodynamic coefficients are examined. For a certain fiber orientation, the deflection and aerodynamic loading of the composite wing is very much reduced. The swept forward wing has more significant effect of wing flexibility on aerodynamic coefficient than the swept back wing does.

Lee, IN↗
NASA NTRS1992

A new method for recognizing quadric surfaces from range data and its application to telerobotics and automation

Pose and orientation of an object is one of the central issues in 3-D recognition problems. Most of today's available techniques require considerable pre-processing such as detecting edges or joints, fitting curves or surfaces to segment images, and trying to extract higher order features from the input images. We present a method based on analytical geometry, whereby all the rotation parameters of any quadric surface are determined and subsequently eliminated. This procedure is iterative in nature and was found to converge to the desired results in as few as three iterations. The approach enables us to position the quadric surface in a desired coordinate system, and then to utilize the presented shape information to explicitly represent and recognize the 3-D surface. Experiments were conducted with simulated data for objects such as hyperboloid of one and two sheets, elliptic and hyperbolic paraboloid, elliptic and hyperbolic cylinders, ellipsoids, and quadric cones. Real data of quadric cones and cylinders were also utilized. Both of these sets yielded excellent results.

Alvertos, Nicolas↗
NASA NTRS1997

Application of the Two-Step Filter to Process Ranging Measurements for Relative Navigation in an Elliptical Orbit

This estimator breaks a nonlinear estimation problem into a set of over determined 'first step' states which are linear in the observations and 'second step' states which are ultimately the states of interest. Linear estimation methods are applied to filter the observations and produce the optimal first step state estimate. The 'second step' states are obtained through iterative nonlinear parameter estimation considering the first step states as observations. It has been shown that this process exactly minimizes the least squares cost function for static problems and provides a better solution than the iterated extended Kalman filter (EKF) for dynamic problems. The two step filter is applied in this paper to process range and range rate measurements between the two spacecraft. Details of the application of the two step estimator to this problem will be given, highlighting the use of a test for ill-conditioned covariance estimates that can result from the first order covariance propagation. A comparison will be made between the performance of the two step filter and the IEKF.

Garrison, James L.↗
NASA NTRS2004

Estimated spectrum adaptive postfilter and the iterative prepost filtering algirighms

The invention presents The Estimated Spectrum Adaptive Postfilter (ESAP) and the Iterative Prepost Filter (IPF) algorithms. These algorithms model a number of image-adaptive post-filtering and pre-post filtering methods. They are designed to minimize Discrete Cosine Transform (DCT) blocking distortion caused when images are highly compressed with the Joint Photographic Expert Group (JPEG) standard. The ESAP and the IPF techniques of the present invention minimize the mean square error (MSE) to improve the objective and subjective quality of low-bit-rate JPEG gray-scale images while simultaneously enhancing perceptual visual quality with respect to baseline JPEG images.

Linares, Irving↗
NASA NTRS2017

A Trilateration Scheme for Relative Positioning

We introduce a trilateration scheme that evaluates the 3-dimensional (3-D) relative position between a reference spacecraft and a target spacecraft using raw-range measurements from a distance baseline of known locations, which we call “anchors”. The anchors can be antennas of a ground-based network (e.g., Deep Space Network (DSN) or Near Earth Network (NEN) stations), or satellites of a spacebased network (e.g., global positioning system (GPS) or tracking and data relay satellite (TDRS)). We define raw-range as the range that includes all the systematic errors that occur during range measurements. A unique feature of this approach is that accurate relative position is derived from a “differencing function” of raw-range measurements of the reference spacecraft and target spacecraft, thereby eliminating most of the systematic errors, such as media effects, ephemeris errors, instrument delays, clock bias, etc. There can be an arbitrary number of target spacecraft, and relative positioning of target spacecraft with respect to the reference spacecraft can be done simultaneously. In this paper, we first assume an idealized system in which clocks on the reference and target spacecraft are synchronized, with clocks of the anchors synchronized as well. We develop a novel iterative algorithm that computes the relative position of the target spacecraft with respect to the reference spacecraft. We illustrate the relative positioning method using the scenario of a network of three ground stations (i.e., the anchors) at Goldstone, California, USA, Madrid, Spain, and Marlargue, Argentina tracking two spacecraft at geosynchronous orbit distance. We demonstrate that the algorithm converges to submeter accuracy in estimating the relative position, in the presence of random errors and systematic errors in raw-range measurements, and in the presence of angular errors in estimating the pointing vectors between the anchors and the reference spacecraft. Next, we relax the requirement of perfect time synchronization between spacecraft, and show that by using an additional anchor, one can estimate and remove the clock biases between the reference and target spacecraft. We add a ground station at Kourou to the above example of three ground stations of Goldstone, Madrid, and Marlargue, and demonstrate that the updated algorithm also converges to meter-level accuracy (submeter in some cases) in the presence of clock biases in addition to the random errors, systematic errors, and angular errors as shown in the above case. We compare this scheme with a similar trilateration scheme for relative positioning scheme first proposed by Montenbruck in 2002.

Cheung, Kar-Ming↗
NASA NTRS2023

Split-Step Simulations to Assess the Effects of Atmospheric Boundary Layer Turbulence on the Dose Variability of N-Waves and Shaped Booms

The effects of atmospheric boundary layer turbulence on the loudness variability of a sonic boom N-wave and shaped boom are examined with split-step simulations. The shaped boom is representative of a design iteration of the NASA X-59 aircraft. Inhomogeneous atmospheric boundary layer turbulence is generated in the computational domain by a Fourier synthesis method. The N-wave and shaped boom are propagated through turbulent fields representing eight different convection levels measured at the NASA Kennedy Space Center and the NASA Armstrong Flight Research Center. Probability density functions of the formation of caustic regions along the propagation direction are computed from the N-wave results, and a parameter to collapse the caustic PDFs that accounts for both fluctuation intensities and length scales is proposed. Statistical results concerning loudness metric variability are presented, and the standard deviations of several metrics are shown to collapse across different convection levels of turbulence for small nondimensional propagation distances. The loudness metric distributions are observed to be well approximated by a normal distribution for a given range of propagation distances, and become increasingly skewed as distance increases. A model function for the dose variability is proposed, and the function parameters are found to be related to the convection level of the turbulence. The model for the dose variability distribution is compared to simulation data that were not used to find the regression parameters of the model. At several nondimensional propagation distances, agreement is observed between the model and the simulation data. These results indicate that the model may be suitable for providing quick estimates of noise dose variability in the primary carpet region across a wide range of atmospheric boundary layer conditions.

Sonic boom↗
NASA NTRS2023

Requirements Engineering Scorecard and the Next-Generation Space Suit

The objective for a NASA contractor, the performing organization in this case study, is to provide engineering services to develop and deliver the next generation space suit to NASA, the customer in this case study. A case study with qualitative and quantitative analyses regarding a new process and approach to requirements engineering is described, with the intent that if utilized, these tools may have contributed to improvements across the project in terms of meeting cost, scope, budget and quality while appropriately accounting for risk management. The procedure entails a research method in which the current state of the project, current state of the art, and the identified systems engineering challenges are evaluated. Iterative models are tempered through development by continual improvements by engineering evaluation of engineers on the project. The current results have produced a prototype of a requirements engineering scorecard with implementations of FMEA and quantitative analysis to (i) identify root cause of underdeveloped requirements and (ii) project management impacts with regards to project risk. Forward work includes customer, performing organization, acceptance against applicable INCOSE community accepted practices.

INCOSE↗
NASA NTRS1980

Three-layer interactive method for computing supersonic laminar separated flows

An interactive model for numerical computation of complicated two-dimensional flowfields including regions of reversed flow is proposed. The present approach is one of dividing the flowfield into three regions, in each of which a simplified mathematical model is applied: (1) outer, supersonic flow for which the full potential equation (hyperbolic) is used; (2) viscous, laminar layer in which the compressible boundary-layer model (parabolic) is used; and (3) recirculating flow modeled by the incompressible Navier-Stokes equations (elliptic). For matching of the numerical solutions in the three layers, two interaction models are developed: one for pressure interaction, the other for interaction between the shear layer and the recirculating flow. The uniform solution for the whole flowfield is then obtained by iteration of the local solutions under the constraints imposed by matching. The three-layer interactive model is used for solution of the flowfield past an asymmetric cavity. The method is shown to be capable of dealing with backflow without encountering problems at separation, characteristic to the boundary-layer approach.

Brandeis, J.↗
NASA NTRS1975

A coupled fluid-structure dynamic analysis for water impact loads

In this paper a nonlinear transient, hydroelastic method is developed for response analysis of the Space Shuttle recoverable booster on water impact. At each time step, coupled equations of fluid flow and structure dynamics are solved with an iterative, self-correcting process. Rigid body motions are represented by large angle equations. Local deformations are represented by a NASTRAN-generated model of the three-dimensional structure. The fluid is represented by series solutions of the potential flow equation which include the effects of the local structure motions. The resulting water loads and structural stresses obtained will aid in qualification of the final structural design. Comparisons are made to experimental impact data to validate the method.

Herting, D. N.↗