Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “factorization 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

Newton like: Minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

Newton-like minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

Design and Analysis Tools for Supersonic Inlets

Computational tools are being developed for the design and analysis of supersonic inlets. The objective is to update existing tools and provide design and low-order aerodynamic analysis capability for advanced inlet concepts. The Inlet Tools effort includes aspects of creating an electronic database of inlet design information, a document describing inlet design and analysis methods, a geometry model for describing the shape of inlets, and computer tools that implement the geometry model and methods. The geometry model has a set of basic inlet shapes that include pitot, two-dimensional, axisymmetric, and stream-traced inlet shapes. The inlet model divides the inlet flow field into parts that facilitate the design and analysis methods. The inlet geometry model constructs the inlet surfaces through the generation and transformation of planar entities based on key inlet design factors. Future efforts will focus on developing the inlet geometry model, the inlet design and analysis methods, a Fortran 95 code to implement the model and methods. Other computational platforms, such as Java, will also be explored.

Slater, John W.↗

On Variable Geometric Factor Systems for Top-Hat Electrostatic Space Plasma Analyzers

Even in the relatively small region of space that is the Earth's magnetosphere, ion and electron fluxes can vary by several orders of magnitude. Top-hat electrostatic analyzers currently do not possess the dynamic range required to sample plasma under all conditions. The purpose of this study was to compare, through computer simulation, three new electrostatic methods that would allow the sensitivity of a sensor to be varied through control of its geometric factor (GF) (much like an aperture on a camera). The methods studied were inner filter plates, split hemispherical analyzer (SHA) and top-cap electrode. This is the first discussion of the filter plate concept and also the first study where all three systems are studied within a common analyzer design, so that their relative merits could be fairly compared. Filter plates were found to have the important advantage that they facilitate the reduction in instrument sensitivity whilst keeping all other instrument parameters constant. However, it was discovered that filter plates have numerous disadvantages that make such a system impracticable for a top-hat electrostatic analyzer. It was found that both the top-cap electrode and SHA are promising variable geometric factor system (VGFS) concepts for implementation into a top-hat electrostatic analyzer, each with distinct advantages over the other.

Kataria, Dhiren O.↗

Mapping Mangrove Extent and Change: A Globally Applicable Approach

This study demonstrates a globally applicable method for monitoring mangrove forest extent at high spatial resolution. A 2010 mangrove baseline was classified for 16 study areas using a combination of ALOS PALSAR and Landsat composite imagery within a random forests classifier. A novel map-to-image change method was used to detect annual and decadal changes in extent using ALOS PALSAR/JERS-1 imagery. The map-to-image method presented makes fewer assumptions of the data than existing methods, is less sensitive to variation between scenes due to environmental factors (e.g., tide or soil moisture) and is able to automatically identify a change threshold. Change maps were derived from the 2010 baseline to 1996 using JERS-1 SAR and to 2007, 2008 and 2009 using ALOS PALSAR. This study demonstrated results for 16 known hotspots of mangrove change distributed globally, with a total mangrove area of 2,529,760 ha. The method was demonstrated to have accuracies consistently in excess of 90% (overall accuracy: 92.2–93.3%, kappa: 0.86) for mapping baseline extent. The accuracies of the change maps were more variable and were dependent upon the time period between images and number of change features. Total change from 1996 to 2010 was 204,850 ha (127,990 ha gain, 76,860 ha loss), with the highest gains observed in French Guiana (15,570 ha) and the highest losses observed in East Kalimantan, Indonesia (23,003 ha). Changes in mangrove extent were the consequence of both natural and anthropogenic drivers, yielding net increases or decreases in extent dependent upon the study site. These updated maps are of importance to the mangrove research community, particularly as the continual updating of the baseline with currently available and anticipated spaceborne sensors. It is recommended that mangrove baselines are updated on at least a 5-year interval to suit the requirements of policy makers.

Thomas, Nathan↗

Computing a Stability Spectrum by Use of the HHT

The Hilbert-Huang transform (HHT) is part of the mathematical basis of a method of calculating a stability spectrum. This method can be regarded as an extended and improved version of a prior HHT-based method of calculating a damping spectrum. In the prior method, information on positive damping (which leads to stability) and negative damping (which leads to instability) becomes mixed into a single squared damping loss factor. Hence, there is no way to distinguish between stability and instability by examining a damping spectrum. In contrast, in the present stabilityspectrum method, information on the mathematical sign of the damping is retained, making it possible to identify reThe Hilbert-Huang transform (HHT) is part of the mathematical basis of a method of calculating a stability spectrum. This method can be regarded as an extended and improved version of a prior HHT-based method of calculating a damping spectrum. In the prior method, information on positive damping (which leads to stability) and negative damping (which leads to instability) becomes mixed into a single squared damping loss factor. Hence, there is no way to distinguish between stability and instability by examining a damping spectrum. In contrast, in the present stabilityspectrum method, information on the mathematical sign of the damping is retained, making it possible to identify regions of instability in a spectrum. This method is expected to be especially useful for analyzing vibration-test data for the purpose of predicting vibrational instabilities in structures (e.g., flutter in airplane wings).

Huang, Norden E.↗

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗

Insulin and insulin-like growth factor-1 induce pronounced hypertrophy of skeletal myofibers in tissue culture

Skeletal myofibers differentiated from primary avian myoblasts in tissue culture can be maintained in positive nitrogen balance in a serum-free medium for at least 6 to 7 days when embedded in a three dimensional collagen gel matrix. The myofibers are metabolically sensitive to physiological concentrations of insulin but these concentrations do not stimulate cell growth. Higher insulin concentrations stimulate both cell hyperplasia and myofiber hypertrophy. Cell growth results from a long term 42 percent increase in total protein synthesis and a 38 percent increase in protein degradation. Myofiber diameters increase by 71 to 98 percent after 6 to 7 days in insulin-containing medium. Insulin-like growth factor-1 but not insulin-like growth factor-2, at 250 ng/ml, is as effective as insulin in stimulating cell hyperplasia and myofiber hypertrophy. This model system provides a new method for studying the long-term anabolic effects of insulin and insulin-like growth factors on myofiber hypertrophy under defined tissue culture conditions.

Vandenburgh, Herman H.↗

Filtering flight data prior to aerodynamic system identification

An algorithm for processing flight-test data to provide state estimates and instrument calibrations for aerodynamic or hydrodynamic system identification by the equation error estimation method is developed and demonstrated on synthesized data. The extended-Kalman-filter algorithm employs a locally level, north-pointing frame of reference, accounts for rotating ellipsoidal earth effects, and estimates sensor bias, scale factors, wind components, and process noise levels by maximum-likelihood parameters. The method is found to be most effective with navigation quality inertial input data. The algorithm is applied to data from a six-degree-of-freedom F-4 aircraft simulation and shown to produce state estimates in good agreement with the simulation values.

Trankle, T. L.↗

Probabilistic finite elements for fracture and fatigue analysis

The fusion of the probabilistic finite element method (PFEM) and reliability analysis for probabilistic fracture mechanics (PFM) is presented. A comprehensive method for determining the probability of fatigue failure for curved crack growth was developed. The criterion for failure or performance function is stated as: the fatigue life of a component must exceed the service life of the component; otherwise failure will occur. An enriched element that has the near-crack-tip singular strain field embedded in the element is used to formulate the equilibrium equation and solve for the stress intensity factors at the crack-tip. Performance and accuracy of the method is demonstrated on a classical mode 1 fatigue problem.

Liu, W. K.↗

A Rotor Tip Vortex Tracing Algorithm for Image Post-Processing

A neurite tracing algorithm, originally developed for medical image processing, was used to trace the location of the rotor tip vortex in density gradient flow visualization images. The tracing algorithm was applied to several representative test images to form case studies. The accuracy of the tracing algorithm was compared to two current methods including a manual point and click method and a cross-correlation template method. It is shown that the neurite tracing algorithm can reduce the post-processing time to trace the vortex by a factor of 10 to 15 without compromising the accuracy of the tip vortex location compared to other methods presented in literature.

Overmeyer, Austin D.↗

Resin Systems and Chemistry-Degradation Mechanisms and Durability in Long-Term Durability of Polymeric Matrix Composites

In choosing a polymer-matrix composite material for a particular application, a number of factors need to be weighed. Among these are mechanical requirements, fabrication method (e.g. press-molding, resin infusion, filament winding, tape layup), and use conditions. Primary among the environmental exposures encountered in aerospace structures are moisture and elevated temperatures, but certain applications may require resistance to other fluids and solvents, alkaline agents, thermal cycling, radiation, or rapid, localized heating (for example, lightning strike). In this chapter, the main classes of polymer resin systems found in aerospace composites will be discussed. Within each class, their responses to environmental factors and the associated degradation mechanisms will be reviewed.

Hinkley, Jeffrey A.↗

Probabilistic finite elements for fracture mechanics

The probabilistic finite element method (PFEM) is developed for probabilistic fracture mechanics (PFM). A finite element which has the near crack-tip singular strain embedded in the element is used. Probabilistic distributions, such as expectation, covariance and correlation stress intensity factors, are calculated for random load, random material and random crack length. The method is computationally quite efficient and can be expected to determine the probability of fracture or reliability.

Besterfield, Glen↗

Matrix Method of Determining the Longitudinal-Stability Coefficients and Frequency Response of an Aircraft from Transient Flight Data

A matrix method is presented for determining the longitudinal-stability coefficients and frequency response of an aircraft from arbitrary maneuvers. The method is devised so that it can be applied to time-history measurements of combinations of such simple quantities as angle of attack, pitching velocity, load factor, elevator angle, and hinge moment to obtain the over-all coefficients. Although the method has been devised primarily for the evaluation of stability coefficients which are of primary interest in most aircraft loads and stability studies, it can be used also, with a simple additional computation, to determine the frequency-response characteristics. The entire procedure can be applied or extended to other problems which can be expressed by linear differential equations.

Donegan, James J↗

A simple method for estimating minimum autorotative descent rate of single rotor helicopters

Flight test results of minimum autorotative descent rate are compared with calculations based on the minimum power required for steady level flight. Empirical correction factors are derived that account for differences in energy dissipation between these two flight conditions. A method is also presented for estimating the minimum power coefficient for level flight for any helicopter for use in the empirical estimation procedure of autorotative descent rate.

Talbot, P. D.↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

The tip of the red giant branch as a distance indicator for resolved galaxies. 2: Computer simulations

Based on both empirical data for the nearby galaxies, and on computer simulations, we show that measuring the position of the tip of the first-ascent red-giant branch provides a means of obtaining the distances to nearby galaxies with a precision and accuracy comparable to using Cepheids and/or RR Lyrae variables. We present an analysis of synthetic I vs (V-I) color magnitude diagrams of Population 2 systems to investigate the use of the observed discontinuity in the I-band luminosity function as a primary distance indicator. In the simulations we quantify the effects (1) signal to noise, (2) crowding, (3) population size, and (4) non-giant-branch-star contamination, on the method adopted for detecting the discontinuity,, measuring its luminosity, and estimating its uncertainity. We discuss sources of systematic error in the context of observable parameters, such as the signal-to-noise ratio and/or surface brightness. The simulations are then scaled to observed color-magnitude diagrams. It is concluded, that from the ground the tip of the red-giant-branch method can be sucessfully used to determine distances accurate to +/- 10% for galaxies out to 3 Mpc (mu approximately 27.5 mag); and from space a factor of four further in distance (mu approximately 30.6 mag) can be reached using HST. This method can be applied whereever a metal-poor population (-2.0 less than Z less than -0.7) of red-giant stars is detected (whose age is in the range 7-17 Gyr), whether that population resides in the halo of a spiral galaxy, the extended outer disk of a dwarf irregular, or in the outer periphery of an elliptical galaxy.

Madore, Barry F.↗

Implicit methods for computing chemically reacting flow

The backward Euler scheme was used to solve a large system of inviscid flow and chemical rate equations in three spatial coordinates. The flow equations were integrated simultaneously in time by a conventional ADI factorization technique, then the species equations were solved by either simultaneous or successive techniques. The methods were evaluated in their efficiency and robustness for a hypersonic flow problem involving an aerobrake configuration. It was found that both implicit methods can effectively reduce the stiffness associated with the chemical production term and that the successive solution for the species was as stable as the simultaneous solution. The latter method is more economical because the computation time varies linearly with the number of species.

Li, C. P.↗