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 127 records · Page 7

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗

Unsteady Aerodynamic Force Sensing from Measured Strain

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm. A cantilevered rectangular wing built and tested at the NASA Langley Research Center (Hampton, Virginia, USA) in 1959 is used to validate the simple approach. Unsteady aerodynamic forces as well as wing deflections, velocities, accelerations, and strains are computed using the CFL3D computational fluid dynamics (CFD) code and an MSC/NASTRAN code (MSC Software Corporation, Newport Beach, California, USA), and these CFL3D-based results are assumed as measured quantities. Based on the measured strains, wing deflections, velocities, accelerations, and aerodynamic forces are computed using the proposed approach. These computed deflections, velocities, accelerations, and unsteady aerodynamic forces are compared with the CFL3D/NASTRAN-based results. In general, computed aerodynamic forces based on the lifting surface theory in subsonic speeds are in good agreement with the target aerodynamic forces generated using CFL3D code with the Euler equation. Excellent aeroelastic responses are obtained even with unsteady strain data under the signal to noise ratio of -9.8dB. The deflections, velocities, and accelerations at each sensor location are independent of structural and aerodynamic models. Therefore, the distributed strain data together with the current proposed approaches can be used as distributed deflection, velocity, and acceleration sensors. This research demonstrates the feasibility of obtaining induced drag and lift forces through the use of distributed sensor technology with measured strain data. An active induced drag control system thus can be designed using the two computed aerodynamic forces, induced drag and lift, to improve the fuel efficiency of an aircraft. Interpolation elements between structural finite element grids and the CFD grids and centroids are successfully incorporated with the unsteady aeroelastic computation scheme. The most critical technology for the success of the proposed approach is the robust on-line parameter estimator, since the least-squares curve fitting method depends heavily on aeroelastic system frequencies and damping factors.

shape sensing↗

Aerosol Particle Deposition on a Spent Nuclear Fuel Assembly Spacer Grid

The flow and particle deposition patterns on surfaces in an idealized spacer grid for a 17x17 pressurized water reactor (PWR) assembly in a spent fuel canister are modeled using computational fluid dynamics (CFD) with laminar flow. The effects of gravitational settling, non-Stokesian flow, and particle slip are first rigorously analyzed. From the analysis, non-Stokesian effects and slip may be neglected for the particle sizes and conditions expected in a canister. For particles that do not settle out, a swirling flow pattern at the corners of a spacer grid channel directs particles to the leeward side of the flow vanes where much of the deposition occurs. Particle deposition increases with increasing particle diameter. Deposition also increases with decreasing flow velocity as this provides more time for particles to settle and deposit on the leeward side of the flow vanes. The fraction of particles that are transmitted through a spacer grid is determined as a function of inlet gas velocity and particle diameter by running the CFD calculation for each set of conditions and for each particle diameter. Curve fits of the transmission curve as a function of particle diameter for a specified spacer grid and flow velocity are applied to a lognormal particle mass density function for the inlet particles. The resulting mass density function and aerosol mass fraction that passes through the spacer grid can be determined analytically without resorting to numerical iteration. A sample calculation of the analytical solution is demonstrated for a lognormal particle mass density function.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Flexible honeycomb structure can bend to fit compound curves

For flexibility in forming a curved surface, a honeycomb configuration using multiple pleats has proved superior to the usual core structures. The partial pleats formed in individual cell walls permit movements to and from the central axis without tearing.

Carmody, R. J.↗

Approximate Solutions Of Equations Of Steady Diffusion

Rigorous analysis yields reliable criteria for "best-fit" functions. Improved "curve-fitting" method yields approximate solutions to differential equations of steady-state diffusion. Method applies to problems in which rates of diffusion depend linearly or nonlinearly on concentrations of diffusants, approximate solutions analytic or numerical, and boundary conditions of Dirichlet type, of Neumann type, or mixture of both types. Applied to equations for diffusion of charge carriers in semiconductors in which mobilities and lifetimes of charge carriers depend on concentrations.

Edmonds, Larry D.↗

An improved method to determine vibration damping of materials and structural components

The present, improved frequency-domain technique for vibration-damping estimation involves the fitting of a 'best curve' for the measured frequency-response data near a resonance by means of a least-squares error criterion. The damping ratio and the undamped natural frequency are then computed using the coefficients of the rational fraction polynomial that had been determined on the basis of the curve-fitting technique. Attention is given to the results of this method's application to damping measurements for graphite-reinforced epoxy composites. The results obtained are superior to those of the 'half-power points' method.

Rao, Mohan D.↗

Experimental validation of solid oxide fuel cell polarization modeling: An LSM-YSZ/YSZ/Ni-YSZ case study

Polarization modeling based on the measured I-V curves is a useful tool for analyzing SOFC performance and diagnosing performance limitations. Using the polarization model, the overpotential of a cell is usually separated into different polarization contributions, including area specific ohmic polarization, activation polarization, and anodic and cathodic concentration polarizations. In this study, we aim to experimentally validate the accuracy of the polarization model. As a case study, anode-supported cells with LSM-YSZ/YSZ/Ni-YSZ configuration are used. By curve-fitting the I-V curves of a cell tested under different cathode oxygen partial pressures into the polarization model, the contributions of different polarization losses are quantified. To validate the polarization model, independent experiments/studies including electrochemical impedance spectroscopy, cathode kinetics study, ex-situ diffusivity measurement, and tortuosity measurement using 3D reconstructed anode are conducted. The results from polarization modeling and from validation studies are in good agreement, thereby validating the polarization model for SOFC performance analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Signal Processing in BOTDA: Utilizing PCA and PCA-Based Neural Networks for Temperature Monitoring

This work presents a comparative analysis of the various signal processing techniques used in the Brillouin gain spectrum (BGS) peak estimation. Traditional fitting methods such as Lorentzian curve fitting (LCF) are slow and less effective in noisy data. PCA-based methods were tested on the experimental data: A Euclidian distance-based approach, and a probabilistic deep neural network (PDNN) based approach, both using 5 principal components to represent a single BGS. Both methods significantly reduce computational time with respect to LCF, whereas PDNN offers uncertainty insights along with the parameter value. Measuring a range of temperatures, analyzing accuracy, and speed, it can be concluded that PCA trained PDNN outperforms other methods, and appears to be helpful in scenario where large datasets are generated.

Brillouin optical time domain analysis↗

Meteoroid and debris special investigation group; status of 3-D crater analysis from binocular imagery

During the 3 month deintegration of the LDEF, the M&D SIG generated approximately 5000 digital color stereo image pairs of impact related features from all space exposed surfaces. Currently, these images are being processed at JSC to yield more accurate feature information. Work is currently underway to determine the minimum number of data points necessary to parametrically define impact crater morphologies in order to minimize the man-hour intensive task of tie point selection. Initial attempts at deriving accurate crater depth and diameter measurements from binocular imagery were based on the assumption that the crater geometries were best defined by paraboloid. We made no assumptions regarding the crater depth/diameter ratios but instead allowed each crater to define its own coefficients by performing a least-squares fit based on user-selected tiepoints. Initial test cases resulted in larger errors than desired, so it was decided to test our basic assumptions that the crater geometries could be parametrically defined as paraboloids. The method for testing this assumption was to carefully slice test craters (experimentally produced in an appropriate aluminum alloy) vertically through the center resulting in a readily visible cross-section of the crater geometry. Initially, five separate craters were cross-sectioned in this fashion. A digital image of each cross-section was then created, and the 2-D crater geometry was then hand-digitized to create a table of XY position for each crater. A 2nd order polynomial (parabolic) was fitted to the data using a least-squares approach. The differences between the fit equation and the actual data were fairly significant, and easily large enough to account for the errors found in the 3-D fits. The differences between the curve fit and the actual data were consistent between the caters. This consistency suggested that the differences were due to the fact that a parabola did not sufficiently define the generic crater geometry. Fourth and 6th order equations were then fitted to each crater cross-section, and significantly better estimates of the crater geometry were obtained with each fit. Work is presently underway to determine the best way to make use of this new parametric crater definition.

Sapp, Clyde A.↗

Scaling laws for light weight optics, studies of light weight mirrors mounting and dynamic mirror stress, and light weight mirror and mount designs

Scaling laws for light-weight optical systems are examined. A cubic relationship between mirror diameter and weight has been suggested and used by many designers of optical systems as the best description for all light-weight mirrors. A survey of existing light-weight systems in the open literature was made to clarify this issue. Fifty existing optical systems were surveyed with all varieties of light-weight mirrors including glass and beryllium structured mirrors, contoured mirrors, and very thin solid mirrors. These mirrors were then categorized and weight to diameter ratio was plotted to find a best curve for each case. A best fitting curve program tests nineteen different equations and ranks a goodness-to-fit for each of these equations. The resulting relationship found for each light-weight mirror category helps to quantify light-weight optical systems and methods of fabrication and provides comparisons between mirror types.

Vukobratovich, Daniel↗

Lithospheric flexure beneath the Freyja Montes foredeep, Venus - Constraints on lithospheric thermal gradient and heat flow

Analysis of Venera 15 and 16 radar images and topographic data from the Freyja Montes region on Venus suggests that this mountain belt formed as a result of a sequence of underthrusts of the lithosphere of the North Polar Plains beneath the highlands of Ishtar Terra. The lithospheric flexure hypothesis is tested, and the effective thickness Te of the elastic lithosphere of the underthrusting portion of the North Polar Plains is estimated by fitting individual topographic profiles to deflection curves for a broken elastic plate. The theoretical curves fit the observed topographic profiles to within measurement error for values of flexural rigidity D in the range (0.8-3) x 10 to the 22nd N m, equivalent to Te in the range 11-18 km.

Solomon, Sean C.↗

Wing Shape Sensing from Measured Strain

A new two step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to extrapolate the deflection and slope of the entire structure through the use of System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular wing. It is then applied to test data from a cantilevered swept wing model.

system equivalent reduction and expansion process↗

Designing an Optimal LSST Deep Drilling Program for Cosmology with Type Ia Supernovae

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) is forecast to collect a large sample of Type Ia supernovae (SNe Ia) expected to be instrumental in unveiling the nature of dark energy. The feat, however, requires accurately measuring the two components of the Hubble diagram, distance modulus and redshift. Distance is estimated from SN Ia parameters extracted from light-curve fits, where the average quality of light curves is primarily driven by survey parameters. An optimal observing strategy is thus critical for measuring cosmological parameters with high accuracy. We present in this paper a three-stage analysis to assess the impact of the deep drilling (DD) strategy parameters on three critical aspects of the survey: redshift completeness, the number of well-measured SNe Ia, and cosmological measurements. We demonstrate that the current DD survey plans (internal LSST simulations) are characterized by a low completeness (z ~ 0.55–0.65), and irregular and low cadences (several days), which dramatically decrease the size of the well-measured SN Ia sample. We propose a method providing the number of visits required to reach higher redshifts. We use the results to design a set of optimized DD surveys for SN Ia cosmology taking full advantage of spectroscopic resources for host galaxy redshift measurements. The most accurate cosmological measurements are achieved with deep rolling surveys characterized by a high cadence (1 day), a rolling strategy (at least two seasons of observation per field), and ultradeep (z ≳ 0.8) and deep (z ≳ 0.6) fields. A deterministic scheduler including a gap recovery mechanism is critical to achieving a high-quality DD survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolution of wide-angle X-ray scattering from a thermoplastic composite

Wide-angle X-ray scattering from a polyetheretherketone/carbonaceous-fiber composite has been resolved using a curve-fit procedure. This resolution yielded under-curve areas for crystalline, amorphous, and fiber scatter which could be used as the basis for an estimation of crystallinity.

Wakelyn, N. T.↗

Wing Shape Sensing from Measured Strain

A new two-step theory is investigated for predicting the deflection and slope of an entire structure using strain measurements at discrete locations. In the first step, a measured strain is fitted using a piecewise least-squares curve fitting method together with the cubic spline technique. These fitted strains are integrated twice to obtain deflection data along the fibers. In the second step, computed deflection along the fibers are combined with a finite element model of the structure in order to interpolate and extrapolate the deflection and slope of the entire structure through the use of the System Equivalent Reduction and Expansion Process. The theory is first validated on a computational model, a cantilevered rectangular plate wing. The theory is then applied to test data from a cantilevered swept-plate wing model. Computed results are compared with finite element results, results using another strain-based method, and photogrammetry data. For the computational model under an aeroelastic load, maximum deflection errors in the fore and aft, lateral, and vertical directions are -3.2 percent, 0.28 percent, and 0.09 percent, respectively; and maximum slope errors in roll and pitch directions are 0.28 percent and -3.2 percent, respectively. For the experimental model, deflection results at the tip are shown to be accurate to within 3.8 percent of the photogrammetry data and are accurate to within 2.2 percent in most cases. In general, excellent matching between target and computed values are accomplished in this study. Future refinement of this theory will allow it to monitor the deflection and health of an entire aircraft in real time, allowing for aerodynamic load computation, active flexible motion control, and active induced drag reduction..

system equivalent reduction and expansion process↗

Development and Validation of an Initial Electric Motor Rotor Vibration Model

An overview of the development of a low fidelity prediction tool for electric motor rotor resonant modes and frequencies. These predictions are useful for the prediction of UAM motor noise during the design phase. Finite Element Analysis of motor rotors has led to the development of a parameterized model to predict these modes shapes. This has in turn allowed to a curve-fit based model from parametric sweeps. These curves can predict resonant frequencies and shapes for 1-5kW class motors without running a simulation.

acoustics↗

Modal vector estimation for closely spaced frequency modes

Techniques for obtaining improved modal vector estimates for systems with closely spaced frequency modes are discussed. In describing the dynamical behavior of a complex structure modal parameters are often analyzed: undamped natural frequency, mode shape, modal mass, modal stiffness and modal damping. From both an analytical standpoint and an experimental standpoint, identification of modal parameters is more difficult if the system has repeated frequencies or even closely spaced frequencies. The more complex the structure, the more likely it is to have closely spaced frequencies. This makes it difficult to determine valid mode shapes using single shaker test methods. By employing band selectable analysis (zoom) techniques and by employing Kennedy-Pancu circle fitting or some multiple degree of freedom (MDOF) curve fit procedure, the usefulness of the single shaker approach can be extended.

Craig, R. R., Jr.↗