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 145 records · Page 8

Component testing for dynamic model verification

Dynamic model verification is the process whereby an analytical model of a dynamic system is compared with experimental data, adjusted if necessary to bring it into agreement with the data, and then qualified for future use in predicting system response in a different dynamic environment. These are various ways to conduct model verification. The approach taken here employs Bayesian statistical parameter estimation. Unlike curve fitting, whose objective is to minimize the difference between some analytical function and a given quantity of test data (or curve), Bayesian estimation attempts also to minimize the difference between the parameter values of that funciton (the model) and their initial estimates, in a least squares sense. The objectives of dynamic model verification, therefore, are to produce a model which: (1) is in agreement with test data; (2) will assist in the interpretation of test data; (3) can be used to help verify a design; (4) will reliably predict performance; and (5) in the case of space structures, will facilitate dynamic control.

Hasselman, T. K.↗

A computer program for model verification of dynamic systems

Dynamic model verification is the process whereby an analytical model of a dynamic system is compared with experimental data, and then qualified for future use in predicting system response in a different dynamic environment. There are various ways to conduct model verification. The approach adopted in MOVER II employs Bayesian statistical parameter estimation. Unlike curve fitting whose objective is to minimize the difference between some analytical function and a given quantity of test data (or curve), Bayesian estimation attempts also to minimize the difference between the parameter values of that function (the model) and their initial estimates, in a least squares sense. The objectives of dynamic model verification, therefore, are to produce a model which: (1) is in agreement with test data, (2) will assist in the interpretation of test data, (3) can be used to help verify a design, (4) will reliably predict performance, and (5) in the case of space structures, facilitate dynamic control.

Chrostowski, J. D.↗

The First Circumbinary Planet Found by Microlensing: OGLE-2007-BLG-349L(AB)c

We present the analysis of the first circumbinary planet microlensing event, OGLE-2007-BLG-349. This event has a strong planetary signal that is best fit with a mass ratio of q approx. = 3.4×10(exp −4), but there is an additional signal due to an additional lens mass, either another planet or another star. We find acceptable light-curve fits with two classes of models: two-planet models (with a single host star) and circumbinary planet models. The light curve also reveals a significant microlensing parallax effect, which constrains the mass of the lens system to be M(sub L) approx. = 0.7 Stellar Mass. Hubble Space Telescope (HST) images resolve the lens and source stars from their neighbors and indicate excess flux due to the star(s) in the lens system. This is consistent with the predicted flux from the circumbinary models, where the lens mass is shared between two stars, but there is not enough flux to be consistent with the two-planet, one-star models. So, only the circumbinary models are consistent with the HST data. They indicate a planet of mass m(sub c) = 80 +/- 13 Stellar Mass, orbiting a pair of M dwarfs with masses of M(sub A) = 0.41+/- 0.07 and M(sub B) = 0.30 +/- 0.07, which makes this the lowest-mass circumbinary planet system known. The ratio of the separation between the planet and the center of mass to the separation of the two stars is approx.40, so unlike most of the circumbinary planets found by Kepler, the planet does not orbit near the stability limit.

OGLE-2007-BLG-349↗

Multi-Criterion Preliminary Design of a Tetrahedral Truss Platform

An efficient method is presented for multi-criterion preliminary design and demonstrated for a tetrahedral truss platform. The present method requires minimal analysis effort and permits rapid estimation of optimized truss behavior for preliminary design. A 14-m-diameter, 3-ring truss platform represents a candidate reflector support structure for space-based science spacecraft. The truss members are divided into 9 groups by truss ring and position. Design variables are the cross-sectional area of all members in a group, and are either 1, 3 or 5 times the minimum member area. Non-structural mass represents the node and joint hardware used to assemble the truss structure. Taguchi methods are used to efficiently identify key points in the set of Pareto-optimal truss designs. Key points identified using Taguchi methods are the maximum frequency, minimum mass, and maximum frequency-to-mass ratio truss designs. Low-order polynomial curve fits through these points are used to approximate the behavior of the full set of Pareto-optimal designs. The resulting Pareto-optimal design curve is used to predict frequency and mass for optimized trusses. Performance improvements are plotted in frequency-mass (criterion) space and compared to results for uniform trusses. Application of constraints to frequency and mass and sensitivity to constraint variation are demonstrated.

Wu, K. Chauncey↗

Multicriterion Design Methods for Aerospace Systems

An efficient method is presented for multi-criterion preliminary design and demonstrated for a tetrahedral truss platform. The present method requires minimal analysis effort and permits rapid estimation of optimized truss behavior for preliminary design. A 14-m-diameter, 3-ring truss platform represents a candidate reflector support structure for space-based science spacecraft. The truss members are divided into 9 groups by truss ring and position. Design variables are the cross-sectional area of all members in a group, and are either 1, 3 or 5 times the minimum member area. Non-structural mass represents the node and joint hardware used to assemble the truss structure. Taguchi methods are used to efficiently identify key points in the set of Pareto-optimal truss designs. Key points identified using Taguchi methods are the maximum frequency, minimum mass, and maximum frequency-to-mass ratio truss designs. Low-order polynomial curve fits through these points are used to approximate the behavior of the full set of Pareto-optimal designs. The resulting Pareto-optimal design curve is used to predict frequency and mass for optimized trusses. Performance improvements are plotted in frequency-mass (criterion) space and compared to results for uniform trusses. Application of constraints to frequency and mass and sensitivity to constraint variation are demonstrated.

taguchi methods↗

Energy from Redox Disproportionation of Sugar Carbon Drives Biotic and Abiotic Synthesis

To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1 %) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucteotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.

Weber, Arthur L.↗

Energy from Redox Disproportionation of Sugar Carbon Drives Biotic and Abiotic Synthesis

To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1%) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucleotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.

Weber, Arthur L.↗

Long-term predictive capability of erosion models

A brief overview of long-term cavitation and liquid impingement erosion and modeling methods proposed by different investigators, including the curve-fit approach is presented. A table was prepared to highlight the number of variables necessary for each model in order to compute the erosion-versus-time curves. A power law relation based on the average erosion rate is suggested which may solve several modeling problems.

Veerabhadra, P.↗

Predictive capability of long-term cavitation and liquid impingement erosion models

A brief overview of long-term cavitation and liquid impingement erosion and modeling methods proposed by different investigators, including the curve-fit approach is presented. A table was prepared to highlight the number of variables necessary for each model in order to compute the erosion-versus-time curves. A power law relation based on the average erosion rate is suggested which may solve several modeling problems. Previously announced in STAR as N83-22386

Rao, P. V.↗

JMOSFET: A MOSFET parameter extractor with geometry-dependent terms

The parameters from metal-oxide-silicon field-effect transistors (MOSFETs) that are included on the Combined Release and Radiation Effects Satellite (CRRES) test chips need to be extracted to have a simple but comprehensive method that can be used in wafer acceptance, and to have a method that is sufficiently accurate that it can be used in integrated circuits. A set of MOSFET parameter extraction procedures that are directly linked to the MOSFET model equations and that facilitate the use of simple, direct curve-fitting techniques are developed. In addition, the major physical effects that affect MOSFET operation in the linear and saturation regions of operation for devices fabricated in 1.2 to 3.0 mm CMOS technology are included. The fitting procedures were designed to establish single values for such parameters as threshold voltage and transconductance and to provide for slope matching between the linear and saturation regions of the MOSFET output current-voltage curves. Four different sizes of transistors that cover a rectangular-shaped region of the channel length-width plane are analyzed.

Buehler, M. G.↗

High speed viscous flow calculations about complex configurations

Applications of the NASA Ames Parabolized Navier-Stokes (PNS) code to a variety of complex generic configurations is presented. The algorithm, boundary conditions, initial conditions, and grid generators are discussed as applied to these configurations. The PNS code was used as the mainline procedure to numerically simulate the viscous supersonic flow over these generic configurations. The turbulence model that was used in this study is the Baldwin-Lomax model. The boundary conditions are the usual viscous no slip at the wall, and a characteristic procedure is used to fit the bow shock wave which is the outermost boundary. An elliptic grid generator is employed to discretize the flow domain. In addition, an equilibrium air capability has been incorporated into the code. It uses the curve fits of Tannehill, et al. The flow regimes vary from a Mach number of 2 up to 25. Both laminar and turbulent flow are considered. Varying angles of attack have also been computed. Configurations vary from simple cone-type bodies to lifting winged bodies, such as the space shuttle or the generic supersonic cruise fighter.

Chaussee, D. S.↗

SRB TPS closeout materials characterization

The recession characteristics of the materials used to closeout the Thermal Protection System (TPS) on the Space Shuttle Solid Rocket Booster (SRB) were assessed. Two candidate closeout materials namely, K5NA and MTA2 is the Marshall Trowelable Ablator. Test panels of each of these two materials were subjected to aerodynamic heat and shear in the NASA Hot Gas Facility. Pretest and Post-test measurements were taken and the rate of recession (R) thus determined was correlated with the heating rate. A least squares curve fit of the data together with another least squares 95% confidence level fit, was performed. The characteristics of the two materials when compared showed that both materials performed equally well at the higher heating rates (10 q 30 Btu/sq ft/sec) but MTA2 was slightly better, not receding as fast, at the lower q values. The comparison of the R versus q curves for K5NA and MTA2 with that of primary TPS P-50 cork on the SRB Aft Skirt was good.

Karu, Z. S.↗

Mechanical and swelling behaviour of well characterized polybutadiene networks

Endlinking of hydroxyl-terminated polybutadiene with the appropriate isocyanate has been used to prepare well characterized networks. Two networks have been studied with molecular weights of the prepolymers being 6100 and 2400 g/mole by g.p.c. Cylindrical specimens were prepared and the derivatives of the stored energy function with respect to the stretch invariants were determined by torque and normal force measurements in torsion. From these data the Valanis-Landel (1967) stored energy function derivatives w-prime(lambda) were determined for both networks. The stored energy function for the junction constraint model of Flory (1953, 1977, 1979, 1985) which is a special form of the Valanis-Landel function, has been fitted to that determined from the experiments. The contributions to the stored energy function from the phantom network and from the junction constraints respectively do not agree with predictions from the topologies of the networks. In spite of this, the form of w-prime(lambda) for the junction constraint model gives an excellent 'curve fit' to the data. Comparison is also made with equilibrium swelling.

Mckenna, Gregory B.↗

Fitting Prony Series To Data On Viscoelastic Materials

Improved method of fitting Prony series to data on viscoelastic materials involves use of least-squares optimization techniques. Based on optimization techniques yields closer correlation with data than traditional method. Involves no assumptions regarding the gamma'(sub i)s and higher-order terms, and provides for as many Prony terms as needed to represent higher-order subtleties in data. Curve-fitting problem treated as design-optimization problem and solved by use of partially-constrained-optimization techniques.

Hill, S. A.↗

Development and Characterization of a Rate-Dependent Three-Dimensional Macroscopic Plasticity Model Suitable for Use in Composite Impact Problems

Several key capabilities have been identified by the aerospace community as lacking in the material/models for composite materials currently available within commercial transient dynamic finite element codes such as LS-DYNA. Some of the specific desired features that have been identified include the incorporation of both plasticity and damage within the material model, the capability of using the material model to analyze the response of both three-dimensional solid elements and two dimensional shell elements, and the ability to simulate the response of composites composed with a variety of composite architectures, including laminates, weaves and braids. In addition, a need has been expressed to have a material model that utilizes tabulated experimentally based input to define the evolution of plasticity and damage as opposed to utilizing discrete input parameters (such as modulus and strength) and analytical functions based on curve fitting. To begin to address these needs, an orthotropic macroscopic plasticity based model suitable for implementation within LS-DYNA has been developed. Specifically, the Tsai-Wu composite failure model has been generalized and extended to a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The coefficients in the yield function are determined based on tabulated stress-strain curves in the various normal and shear directions, along with selected off-axis curves. Incorporating rate dependence into the yield function is achieved by using a series of tabluated input curves, each at a different constant strain rate. The non-associative flow-rule is used to compute the evolution of the effective plastic strain. Systematic procedures have been developed to determine the values of the various coefficients in the yield function and the flow rule based on the tabulated input data. An algorithm based on the radial return method has been developed to facilitate the numerical implementation of the material model. The presented paper will present in detail the development of the orthotropic plasticity model and the procedures used to obtain the required material parameters. Methods in which a combination of actual testing and selective numerical testing can be combined to yield the appropriate input data for the model will be described. A specific laminated polymer matrix composite will be examined to demonstrate the application of the model.

Impact↗

Measurements of the stochastic nature of atmospheric spectral amplitudes

The stochastic nature of the power spectral amplitudes of the neutral atmospheric boundary layer is examined. Probability density distributions and probability distributions of longitudinal and lateral power spectra amplitudes are computed from neutral atmospheric boundary layers. The statistical distributions are computed for frequencies of 0.006, 0.01, 0.03, 0.06, 0.1, and 0.5 Hz at each of the elevations of 18, 30, 60, 90, 120, and 150 m. When the probability density distributions are properly nondimensionalized, the data tend to collapse to a universal curve. An empirical curve fit to the universal nondimensionalized probability density distribution is also given. Probability distributions of individual frequency power spectral amplitudes are also presented for all elevations and frequencies. An interesting observation from the data is that greater than 10 percent of the time the power spectral amplitude at a given frequency will genrally be more than three times the temporal mean value computed by standard Fourier techniques. The standard power spectral density curves are also included in the report.

Cliff, W. C.↗

Optimal Stimulus Amplitude for Vestibular Stochastic Stimulation to Improve Sensorimotor Function

Sensorimotor changes such as postural and gait instabilities can affect the functional performance of astronauts when they transition across different gravity environments. We are developing a method, based on stochastic resonance (SR), to enhance information transfer by applying non-zero levels of external noise on the vestibular system (vestibular stochastic resonance, VSR). Our previous work has shown the advantageous effects of VSR in a balance task of standing on an unstable surface. This technique to improve detection of vestibular signals uses a stimulus delivery system that is wearable or portable and provides imperceptibly low levels of white noise-based binaural bipolar electrical stimulation of the vestibular system. The goal of this project is to determine optimal levels of stimulation for SR applications by using a defined vestibular threshold of motion detection. A series of experiments were carried out to determine a robust paradigm to identify a vestibular threshold that can then be used to recommend optimal stimulation levels for SR training applications customized to each crewmember. Customizing stimulus intensity can maximize treatment effects. The amplitude of stimulation to be used in the VSR application has varied across studies in the literature such as 60% of nociceptive stimulus thresholds. We compared subjects' perceptual threshold with that obtained from two measures of body sway. Each test session was 463s long and consisted of several 15s sinusoidal stimuli, at different current amplitudes (0-2 mA), interspersed with 20-20.5s periods of no stimulation. Subjects sat on a chair with their eyes closed and had to report their perception of motion through a joystick. A force plate underneath the chair recorded medio-lateral shear forces and roll moments. First we determined the percent time during stimulation periods for which perception of motion (activity above a pre-defined threshold) was reported using the joystick, and body sway (two standard deviation of the noise level in the baseline measurement) was detected by the sensors. The percentage time at each stimulation level for motion detection was normalized with respect to the largest value and a logistic regression curve fit was applied to these data. The threshold was defined at the 50% probability of motion detection. Comparison of threshold of motion detection obtained from joystick data versus body sway suggests that perceptual thresholds were significantly lower, and were not impacted by system noise. Further, in order to determine optimal stimulation amplitude to improve balance, two sets of experiments were carried out. In the first set of experiments, all subjects received the same level of stimuli and the intensity of optimal performance was projected back on subjects' vestibular threshold curve. In the second set of experiments, on different subjects, stimulation was administered from 20-400% of subjects' vestibular threshold obtained from joystick data. Preliminary results of our study show that, in general, using stimulation amplitudes at 40-60% of perceptual motion threshold improved balance performance significantly compared to control (no stimulation). The amplitude of vestibular stimulation that improved balance function was predominantly in the range of +/- 100 to +/- 400 micro A. We hypothesize that VSR stimulation will act synergistically with sensorimotor adaptability (SA) training to improve adaptability by increasing utilization of vestibular information and therefore will help us to optimize and personalize a SA countermeasure prescription. This combination will help to significantly reduce the number of days required to recover functional performance to preflight levels after long-duration spaceflight.

Goel, R.↗

Polynomial Expressions for Estimating Elastic Constants From the Resonance of Circular Plates

Two approaches were taken to make convenient spread sheet calculations of elastic constants from resonance data and the tables in ASTM C1259 and E1876: polynomials were fit to the tables; and an automated spread sheet interpolation routine was generated. To compare the approaches, the resonant frequencies of circular plates made of glass, hardened maraging steel, alpha silicon carbide, silicon nitride, tungsten carbide, tape cast NiO-YSZ, and zinc selenide were measured. The elastic constants, as calculated via the polynomials and linear interpolation of the tabular data in ASTM C1259 and E1876, were found comparable for engineering purposes, with the differences typically being less than 0.5 percent. Calculation of additional v values at t/R between 0 and 0.2 would allow better curve fits. This is not necessary for common engineering purposes, however, it might benefit the testing of emerging thin structures such as fuel cell electrolytes, gas conversion membranes, and coatings when Poisson s ratio is less than 0.15 and high precision is needed.

Salem, Jonathan A.↗