Search NASA⌕ Search

SEARCH · Search NASA

Results for “Curve fitting”

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 379 records · Page 21

Solar Cycle Predictions

We currently lack a comprehensive model for the solar dynamo which might be used for predicting solar cycle activity levels. We can, nonetheless, predict these activity levels with some skill using a variety of techniques based on correlations between activity indicators at different times. A number of the techniques currently in use for predicting solar activity on a solar cycle time-scale have been tested with historical data. Some techniques, e.g., regression and curve fitting, work well as solar activity approaches maximum and provide a month by month description of future activity, while others, e.g., geomagnetic precursors, work well near solar minimum but only provide an estimate of the amplitude of the next cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. The combination of two uncorrelated geomagnetic precursor techniques provides a more accurate prediction for the amplitude of a solar activity cycle of a time well before activity minimum. A mathematical function dependent on the time of cycle initiation and the cycle amplitude can then be used to describe the level of solar activity month by month for the next cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between previous activity levels and this function. The success of the geomagnetic precursors in predicting the amplitude of the next solar cycle supports the idea of extended solar cycle in which the activity cycles overlap. While the current cycle is producing sunspots and active regions in the Sun's equatorial regions the next cycle is already actively influencing the geomagnetic field, The precise nature and source of this influence remains unclear.

Hathaway, D. H.↗

Tentative Identification of the 780/cm nu(sub 4) Band Q Branch of Chlorine Nitrate in High-Resolution Solar Absorption Spectra of the Stratosphere

Absorption by the Q branch of the nu(sub 4), band of ClONO2 at 780.2/cm has been tentatively identified in a series of 0.02/cm resolution balloon-borne solar absorption spectra of the stratosphere. The spectral data were recorded at sunset from a flot altitude of 33.5 km during a balloon flight from Holloman Air Force Base (32.8deg N, 106.0 deg W) near Alamogordo, New Mexico, on March 23 1998. A preliminary ClONO2 vertical profile has been determined from the stratospheric spectra by using the technique of nonlinear least squares spectral curve fitting and new spectroscopic parameters deduced from high-resolution laboratory spectra of ClONO2 and O3.

Rinsland, C. P.↗

A Synthesis of Solar Cycle Prediction Techniques

A number of techniques currently in use for predicting solar activity on a solar cycle timescale are tested with historical data. Some techniques, e.g., regression and curve fitting, work well as solar activity approaches maximum and provide a month-by-month description of future activity, while others, e.g., geomagnetic precursors, work well near solar minimum but only provide an estimate of the amplitude of the cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. A combination of two uncorrelated geomagnetic precursor techniques provides a more accurate prediction for the amplitude of a solar activity cycle at a time well before activity minimum. This combined precursor method gives a smoothed sunspot number maximum of 154 plus or minus 21 at the 95% level of confidence for the next cycle maximum. A mathematical function dependent on the time of cycle initiation and the cycle amplitude is used to describe the level of solar activity month by month for the next cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between previous activity levels and this function. This Combined Solar Cycle Activity Forecast gives, as of January 1999, a smoothed sunspot maximum of 146 plus or minus 20 at the 95% level of confidence for the next cycle maximum.

Hathaway, David H.↗

Status of Cycle 23 Forecasts

A number of techniques for predicting solar activity on a solar cycle time scale are identified, described, and tested with historical data. Some techniques, e.g,, regression and curve-fitting, work well as solar activity approaches maximum and provide a month- by-month description of future activity, while others, e.g., geomagnetic precursors, work well near solar minimum but provide an estimate only of the amplitude of the cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. A combination of two uncorrelated geomagnetic precursor techniques provides the most accurate prediction for the amplitude of a solar activity cycle at a time well before activity minimum. This precursor method gave a smoothed sunspot number maximum of 154+21 for cycle 23. A mathematical function dependent upon the time of cycle initiation and the cycle amplitude then describes the level of solar activity for the complete cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between recent activity levels and this function. This Combined Solar Cycle Activity Forecast now gives a smoothed sunspot maximum of 140+20 for cycle 23. The success of the geomagnetic precursors in predicting future solar activity suggests that solar magnetic phenomena at latitudes above the sunspot activity belts are linked to solar activity, which occurs many years later in the lower latitudes.

Hathaway, D. H.↗

Review of Skin Friction Measurements Including Recent High-Reynolds Number Results from NASA Langley NTF

This paper reviews flat plate skin friction data from early correlations of drag on plates in water to measurements in the cryogenic environment of The NASA Langley National Transonic Facility (NTF) in late 1996. The flat plate (zero pressure gradient with negligible surface curvature) incompressible skin friction at high Reynolds numbers is emphasized in this paper, due to its importance in assessing the accuracy of measurements, and as being important to the aerodynamics of large scale vehicles. A correlation of zero pressure gradient skin friction data minimizing extraneous effects between tests is often used as the first step in the calculation of skin friction in complex flows. Early data compiled by Schoenherr for a range of momentum thickness Reynolds numbers, R(sub Theta) from 860 to 370,000 contained large scatter, but has proved surprisingly accurate in its correlated form. Subsequent measurements in wind tunnels under more carefully controlled conditions have provided inputs to this database, usually to a maximum R(sub Theta) of about 40,000. Data on a large axisymmetric model in the NASA Langley National Transonic Facility extends the upper limit in incompressible R(sub Theta) to 619,800 using the van Driest transformation. Previous data, test techniques, and error sources ar discussed, and the NTF data will be discussed in detail. The NTF Preston tube and Clauser inferred data accuracy is estimated to be within -2 percent of a power-law curve fit, and falls above the Spalding theory by 1 percent at R(sub Theta) of about 600,000.

Watson, Ralph D.↗

Paralinear Oxidation of CVD SiC in Simulated Fuel-Rich Combustion

The oxidation kinetics of CVD SiC were measured by thermogravimetric analysis (TGA) in a 4H2 (central dot) 12H2O (central dot) 10CO (central dot) 7CO2 (central dot) 67N2 gas mixture flowing at 0.44 cm/s at temperatures between 1300 and 1450 C in fused quartz furnace tubes at I atm total pressure. The SiC was oxidized to form solid SiO2. At less than or = 1350 C, the SiO2 was in turn volatilized. Volatilization kinetics were consistent with the thermodynamic predictions based on SiO formation. These two simultaneous reactions resulted in overall paralinear kinetics. A curve fitting technique was used to determine the linear and parabolic rate constants from the paralinear kinetic data. Volatilization of the protective SiO2 scale resulted in accelerated consumption of SiC. Recession rates under conditions more representative of actual combustors were estimated from the furnace data.

Fox, Dennis S.↗

Optical Characterization of Bulk ZnSeTe Solid Solutions

Optical characterization was performed on wafers sliced from crystals of ZnSe, ZnTe, and ZnSe(1-x)Te(x)(0 less than x less than 0.4) grown by physical vapor transport. Energy band gaps at room temperature were determined from optical transmission measurements on 11 wafers. A best fit curve to the band gap versus composition x data gives a bowing parameter of 1.45. This number lies between the value of 1.23 determined previously on ZnSeTe bulk crystals and the value of 1.621 reported on ZnSeTe epilayers. Low-temperature photoluminescence (PL) spectra were measured on six samples. The spectra of ZnSe and ZnTe were dominated by near band edge emissions and no deep donor-acceptor pairs were observed. The PL spectrum exhibited a broad emission for each of the ZnSe(1-x)Te(x) samples, 0.09 less than x less than 0.39. For x=0.09, this emission energy is about 0.2 eV lower than the band gap energy measured at low temperature. As x increases the energy discrepancy gradually decreases and reduces to almost zero at x=0.4. The single broad PL emission spectra and the spectra measured as a function of temperature were interpreted as being associated with the exciton bound to Te clusters because of the high Te content in these samples.

Su, Ching-Hua↗

Dynamic Stiffness and Damping Characteristics of a High-Temperature Air Foil Journal Bearing

Using a high-temperature optically based displacement measurement system, a foil air bearing's stiffness and damping characteristics were experimentally determined. Results were obtained over a range of modified Sommerfeld Number from 1.5E6 to 1.5E7, and at temperatures from 25 to 538 C. An Experimental procedure was developed comparing the error in two curve fitting functions to reveal different modes of physical behavior throughout the operating domain. The maximum change in dimensionless stiffness was 3.0E-2 to 6.5E-2 over the Sommerfeld Number range tested. Stiffness decreased with temperature by as much as a factor of two from 25 to 538 C. Dimensionless damping was a stronger function of Sommerfeld Number ranging from 20 to 300. The temperature effect on damping being more qualitative, showed the damping mechanism shifted from viscous type damping to frictional type as temperature increased.

Howard, Samuel A.↗

Solar Activity Forecasting for use in Orbit Prediction

Orbital prediction for satellites in low Earth orbit (LEO) or low planetary orbit depends strongly on exospheric densities. Solar activity forecasting is important in orbital prediction, as the solar UV and EUV inflate the upper atmospheric layers of the Earth and planets, forming the exosphere in which satellites orbit. Geomagnetic effects also relate to solar activity. Because of the complex and ephemeral nature of solar activity, with different cycles varying in strength by more than 100%, many different forecasting techniques have been utilized. The methods range from purely numerical techniques (essentially curve fitting) to numerous oddball schemes, as well as a small subset, called 'Precursor techniques.' The situation can be puzzling, owing to the numerous methodologies involved, somewhat akin to the numerous ether theories near the turn of the last century. Nevertheless, the Precursor techniques alone have a physical basis, namely dynamo theory, which provides a physical explanation for why this subset seems to work. I discuss this solar cycle's predictions, as well as the Sun's observed activity. I also discuss the SODA (Solar Dynamo Amplitude) index, which provides the user with the ability to track the Sun's hidden, interior dynamo magnetic fields. As a result, one may then update solar activity predictions continuously, by monitoring the solar magnetic fields as they change throughout the solar cycle. This paper ends by providing a glimpse into what the next solar cycle (#24) portends.

Schatten, Kenneth↗

Case Studies of Water Vapor and Surface Liquid Water from AVIRIS Data Measured Over Denver, CO and Death Valley, CA

High spatial resolution column atmospheric water vapor amounts and equivalent liquid water thicknesses of surface targets are retrieved from spectral data collected by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The retrievals are made using a nonlinear least squares curve fitting technique. Two case studies from AVIRIS data acquired over Denver-Platteville area, Colorado and over Death Valley, California are presented. The column water vapor values derived from AVIRIS data over the Denver-Platteville area are compared with those obtained from radiosondes, ground level upward-looking microwave radiometers, and geostationary satellite measurements. The column water vapor image shows spatial variation patterns related to the passage of a weather front system. The column water vapor amounts derived from AVIRIS data over Death Valley decrease with increasing surface elevation. The derived liquid water image clearly shows surface drainage patterns.

Gao, B.-C.↗

Viscosity Relaxation in Molten HgZnTe

Rotating cup measurements of the viscosity of the pseudo-binary melt, HgZnTe have shown that the isothermal liquid with zinc mole fraction 0.16 requires tens of hours of equilibration time before a steady viscous state can be achieved. Over this relaxation period, the viscosity at 790 C increases by a factor of two, while the viscosity at 810 C increases by 40%. Noting that the Group VI elements tend to polymerize when molten, we suggest that the viscosity of the melt is enhanced by the slow formation of Te atom chains. To explain the build-up of linear Te n-mers, we propose a scheme, which contains formation reactions with second order kinetics that increase the molecular weight, and decomposition reactions with first order kinetics that inactivate the chains. The resulting rate equations can be solved for the time dependence of each molecular weight fraction. Using these molecular weight fractions, we calculate the time dependence of the average molecular weight. Using the standard semi-empirical relation between polymer average molecular weight and viscosity, we then calculate the viscosity relaxation curve. By curve fitting, we find that the data imply that the rate constant for n-mer formation is much smaller than the rate constant for n-mer deactivation, suggesting that Te atoms only weakly polymerize in molten HgZnTe. The steady state toward which the melt relaxes occurs as the rate of formation of an n-mer becomes exactly balanced by the sum of the rate for its deactivation and the rate for its polymerization to form an (n+1)-mer.

Su, Ching-Hua↗

Fatigue Failure Criteria for Single Crystal Nickel Superalloys

High Cycle Fatigue (HCF) induced failures in aircraft gas-turbine and rocket engine turbopump blades is a pervasive problem. Single crystal turbine blades are being utilized in rocket engine turbopumps and jet engines throughout industry and NASA because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Single-crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the pan geometry a significant factor in the overall analysis. The failure modes of single crystal turbine blades is complicated to predict due to the material orthotropy and variations in crystal orientations. A fatigue failure criteria based on the maximum shear stress amplitude [delta t max] on the 30 slip systems, is presented for single crystal nickel superalloys (FCC crystal). This criteria reduces the scatter in uniaxial LCF test data, for four different specimen orientations, for PWA 1484 at 1200 F in air, quite well. A power law curve fit of the failure parameter, delta t max, vs. cycles to failure is presented.

Arakere, Nagaraj K.↗

Crystal Growth of ZnSe and Related Ternary Compound Semiconductors by Vapor Transport

The objective of the project is to determine the relative contributions of gravity-driven fluid flows to the compositional distribution, incorporation of impurities and defects, and deviation from stoichiometry observed in the crystals grown by vapor transport as results of buoyance-driven convection and growth interface fluctuations caused by irregular fluid-flows. ZnSe and related ternary compounds, such as ZnSeS and ZnSeTe, were grown by vapor transport technique with real time in-situ non-invasive monitoring techniques. The grown crystals were characterized extensively to correlate the grown crystal properties with the growth conditions. The following are the research progress in the past two years. In-situ monitoring of partial pressure by optical absorption technique and visual observation of the growing crystal were performed during vapor growth of ZnSe. Low-temperature photoluminescence (PL) spectra and glow discharge mass spectroscopy (GDMS) were measured on ZnSe starting materials provided by various vendors and on bulk crystals grown from these starting materials by physical vapor transport (PVT) to study the effects of purification and contamination during crystal growth process. Optical characterization was performed on wafers sliced from the grown crystals of ZnSe, ZnTe and ZnSe(1-x),Te(x), (0<x<0.4). Energy band gaps at room temperature were determined from optical transmission measurements and a best fit curve to the band gap vs. composition, x, data gives a bowing parameter of 1.45. Low-temperature photoluminescence (PL) spectra of ZnSe and ZnTe were dominated by near band edge emissions and no deep donor-acceptor pairs were observed. The PL spectrum exhibited a broad emission for the ZnSe(1-x),Te(x), samples, 0.09<x<0.39. The single broad PL emission spectra and the spectra measured as a function of temperature were interpreted as being associated with the exciton bound to Te clusters because of the high Te content in these samples. To validate numerical codes, in-situ monitoring during the PVT of HgI2 was designed. Optical absorption spectra of the vapor phase over HgI2 were measured for wavelengths between 200 and 600nm at sample temperatures between 349 and 610K. The Beer's Law constants for 15 wavelengths between 200 and 440 nm were determined. From these constants the vapor pressure of HgI2 was established as a function of temperature for the liquid and the solid Beta-phases To characterize the growth conditions during the PVT growth of In-doped ZnSe the optical absorbance of the vapor phase over the In-Se system were measured and were used to obtain the partial pressures of Se2(g) and In2Se(g).

Su, Ching-Hua↗

An Investigation of Applications for Thermodynamic Work Potential Methods: Working Tables and Charts for Estimation of Thermodynamic Work Potential in Equilibrium Mixtures of Jet-A and Air

The objective of this report is to provide a tool to facilitate the application of thermodynamic work potential methods to aircraft and engine analysis. This starts with a discussion of the theoretical background underlying these methods, which is then used to derive various equations useful for thermodynamic analysis of aircraft engines. The work potential analysis method is implemented in the form of a set of working charts and tables that can be used to graphically evaluate work potential stored in high-enthalpy gas. The range of validity for these tables is 300 to 36,000 R, pressures between between 0.01 atm and 100 atm, and fuel-air ratios from zero to stoichiometric. The derivations and charts assume mixtures of Jet-A and air as the working fluid. The thermodynamic properties presented in these charts were calculated based upon standard thermodynamic curve fits.

Mavris, Dimitri↗

Prediction of Aerodynamic Coefficient using Genetic Algorithm Optimized Neural Network for Sparse Data

Wind tunnels use scale models to characterize aerodynamic coefficients, Wind tunnel testing can be slow and costly due to high personnel overhead and intensive power utilization. Although manual curve fitting can be done, it is highly efficient to use a neural network to define the complex relationship between variables. Numerical simulation of complex vehicles on the wide range of conditions required for flight simulation requires static and dynamic data. Static data at low Mach numbers and angles of attack may be obtained with simpler Euler codes. Static data of stalled vehicles where zones of flow separation are usually present at higher angles of attack require Navier-Stokes simulations which are costly due to the large processing time required to attain convergence. Preliminary dynamic data may be obtained with simpler methods based on correlations and vortex methods; however, accurate prediction of the dynamic coefficients requires complex and costly numerical simulations. A reliable and fast method of predicting complex aerodynamic coefficients for flight simulation I'S presented using a neural network. The training data for the neural network are derived from numerical simulations and wind-tunnel experiments. The aerodynamic coefficients are modeled as functions of the flow characteristics and the control surfaces of the vehicle. The basic coefficients of lift, drag and pitching moment are expressed as functions of angles of attack and Mach number. The modeled and training aerodynamic coefficients show good agreement. This method shows excellent potential for rapid development of aerodynamic models for flight simulation. Genetic Algorithms (GA) are used to optimize a previously built Artificial Neural Network (ANN) that reliably predicts aerodynamic coefficients. Results indicate that the GA provided an efficient method of optimizing the ANN model to predict aerodynamic coefficients. The reliability of the ANN using the GA includes prediction of aerodynamic coefficients to an accuracy of 110% . In our problem, we would like to get an optimized neural network architecture and minimum data set. This has been accomplished within 500 training cycles of a neural network. After removing training pairs (outliers), the GA has produced much better results. The neural network constructed is a feed forward neural network with a back propagation learning mechanism. The main goal has been to free the network design process from constraints of human biases, and to discover better forms of neural network architectures. The automation of the network architecture search by genetic algorithms seems to have been the best way to achieve this goal.

Rajkumar, T.↗

Xenon Sputter Yield Measurements for Ion Thruster Materials

In this paper, we describe a technique that was used to measure total and differential sputter yields of materials important to high specific impulse ion thrusters. The heart of the technique is a quartz crystal monitor that is swept at constant radial distance from a small target region where a high current density xenon ion beam is aimed. Differential sputtering yields were generally measured over a full 180 deg arc in a plane that included the beam centerline and the normal vector to the target surface. Sputter yield results are presented for a xenon ion energy range from 0.5 to 10 keV and an angle of incidence range from 0 deg to 70 deg from the target surface normal direction for targets consisting of molybdenum, titanium, solid (Poco) graphite, and flexible graphite (grafoil). Total sputter yields are calculated using a simple integration procedure and comparisons are made to sputter yields obtained from the literature. In general, the agreement between the available data is good. As expected for heavy xenon ions, the differential and total sputter yields are found to be strong functions of angle of incidence. Significant under- and over-cosine behavior is observed at low- and high-ion energies, respectively. In addition, strong differences in differential yield behavior are observed between low-Z targets (C and Ti) and high-Z targets (Mo). Curve fits to the differential sputter yield data are provided. They should prove useful to analysts interested in predicting the erosion profiles of ion thruster components and determining where the erosion products re-deposit.

Williams, John D.↗

Unbalance Response Prediction for Rotors on Ball Bearings Using Speed and Load Dependent Nonlinear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic analysis requires that bearing forces corresponding to the actual bearing deflection be utilized. For this work bearing forces were calculated by COBRA-AHS, a recently developed rolling-element bearing analysis code. Bearing stiffness was found to be a strong function of bearing deflection, with higher deflection producing markedly higher stiffness. Curves fitted to the bearing data for a range of speeds and loads were supplied to a flexible rotor unbalance response analysis. The rotordynamic analysis showed that vibration response varied nonlinearly with the amount of rotor imbalance. Moreover, the increase in stiffness as critical speeds were approached caused a large increase in rotor and bearing vibration amplitude over part of the speed range compared to the case of constant bearing stiffness. Regions of bistable operation were possible, in which the amplitude at a given speed was much larger during rotor acceleration than during deceleration. A moderate amount of damping will eliminate the bistable region, but this damping is not inherent in ball bearings.

Fleming, David P.↗

Failure Criteria for FRP Laminates in Plane Stress

A new set of six failure criteria for fiber reinforced polymer laminates is described. Derived from Dvorak's fracture mechanics analyses of cracked plies and from Puck's action plane concept, the physically-based criteria, denoted LaRC03, predict matrix and fiber failure accurately without requiring curve-fitting parameters. For matrix failure under transverse compression, the fracture plane is calculated by maximizing the Mohr-Coulomb effective stresses. A criterion for fiber kinking is obtained by calculating the fiber misalignment under load, and applying the matrix failure criterion in the coordinate frame of the misalignment. Fracture mechanics models of matrix cracks are used to develop a criterion for matrix in tension and to calculate the associated in-situ strengths. The LaRC03 criteria are applied to a few examples to predict failure load envelopes and to predict the failure mode for each region of the envelope. The analysis results are compared to the predictions using other available failure criteria and with experimental results. Predictions obtained with LaRC03 correlate well with the experimental results.

Davila, Carlos G.↗