Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LEAST SQUARE FIT”

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 613 records · Page 34

Jig-Shape Optimization of a Low-Boom Supersonic Aircraft

A simple approach for optimizing the jig-shape is proposed in this study. This simple approach is based on an unconstrained optimization problem and applied to a low-boom supersonic aircraft. In this study, the jig-shape optimization is performed using the two-step approach. First, starting design variables are computed using the least squares surface fitting technique. Next, the jig-shape is further tuned using a numerical optimization procedure based on in-house object-oriented optimization tool.

numerical analysis↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Data Processing of Miniaturized Laser Heterodyne Radiometer (mini-LHR) Ground Instrument Retrievals

The Miniaturized Laser Heterodyne Radiometer (mini-LHR) is a passive ground instrument that observes the mole fraction of carbon dioxide (CO) and methane (CH)in the atmospheric column by measuring their absorption of sunlight at 1.6 microns. A laser heterodyne radiometer is similar in design to the super heterodyne radio receiver that is well known by ham radio enthusiasts. While not previously a commercial technique, laser heterodyne radiometers have a history of measuring atmospheric trace gases that started in the 1960s. With the commercial availability of inexpensive, low-power, thumbnail-sized lasers developed for the telecommunications industry, it was possible to miniaturize this technique and ultimately commercialize it. The mini-LHR has been under development at NASA GSFC since 2009. During that time, in addition to signicant technical improvements, processing has also evolved and been streamlined. Here we present details of the processing approach for raw mini-LHR data to produce 30- and 60-minute data products of CH and CO column mole fractions. Processing occurs in two general stages: a python-based pre-processing of raw data, followed by ingestion into a Planetary Spectrum Generator (PSG) retrieval algorithm. Raw data processing involves removal of outliers, correcting for changes in air mass throughout the day, averaging scans, and ultimately converting averaged scans into transmittance vs. wavelength. The PSG retrieval simulates a spectra for the time/day/location of the scan with meteorological inputs from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data set and then perturbs the concentrations of CO and CH to obtain a t based on an iterative least-squares curve fitting procedure.

Giancarlo Roberto Zambrano↗

An Initial Electric Motor Rotor Vibration Model

Ongoing integration of outrunner brushless electric motors into drones and Advanced Air Mobility aircraft presents a need for accurate acoustic predictions derived from prediction of the motor rotor vibrations. In outrunner motors, the external rotor vibrates and drives the acoustic field. At resonant frequencies of the rotor, the rotor displacements will be the largest. The vibrations can lead to acoustic tones that are relevant to the overall acoustic design of these types of vehicles. This study performs a finite element analysis to assess the modes in two electric motor rotors. This model is then refined into a simple parameterized geometry that can make the same predictions for motors of the same class. A parametric sweep and least-squares-curve fit to the finite element analysis data results in a series of simple curves that can predict the mode shapes and frequencies that are likely to appear in this class of motors without the need for any simulation. The simulations are validated by comparison to acoustic and experimental modal-analysis data.

Electric Motor Noise↗

An Initial Electric Motor Rotor Vibration Model

Ongoing integration of outrunner brushless electric motors into drones and Advanced Air Mobility aircraft presents a need for accurate acoustic predictions derived from prediction of the motor rotor vibrations. In outrunner motors, the external rotor vibrates and drives the acoustic field. At resonant frequencies of the rotor, the rotor displacements will be the largest. The vibrations can lead to acoustic tones that are relevant to the overall acoustic design of these types of vehicles. This study performs a finite element analysis to assess the modes in two electric motor rotors. This model is then refined into a simple parameterized geometry that can make the same predictions for motors of the same class. A parametric sweep and least-squares-curve fit to the finite element analysis data results in a series of simple curves that can predict the mode shapes and frequencies that are likely to appear in this class of motors without the need for any simulation. The simulations are validated by comparison to acoustic and experimental modal-analysis data.

Electric Motor Noise↗

Simultaneous Observations of Surface Layer Profiles of Humidity, Temperature, and Wind Using Scanning Lidar Instruments

We demonstrate the combination of three Doppler lidars (DLs), a water vapor differential absorption lidar and a temperature rotational Raman lidar for the investigation of the interactions between the land-surface and the atmospheric boundary layer. This combination of scanning lidars was operated for the first time during the Land-Atmosphere Feedback Experiment at the Atmospheric Radiation Measurement program's Southern Great Plain site, Oklahoma, USA, in August 2017, and provided simultaneous surface layer profiles of horizontal wind, humidity and temperature. The horizontal wind profiles were determined using the dual-Doppler method with two DLs. The scans were performed above four towers providing atmospheric variables and fluxes. These combined lidar data allowed for the estimation of the friction velocity as well as the surface latent and sensible heat fluxes. For this purpose, profiles calculated with Monin-Obukhov similarity theory were fitted with the Levenberg-Marquardt nonlinear least squares curve-fitting method to the measured surface-layer lidar profiles. We present case studies of three 50-min periods to illustrate the new method. The measurements resulted in friction velocities of 0.31, 0.29, and 0.38 ms –1 , sensible surface heat fluxes of 312.1, 234.1, and 183.3 Wm –2 and latent surface heat fluxes of 251.8, 227.3, and 274.4 Wm –2 . These values were compared with the in-situ measurements of the towers. Considering the sampling differences and the error analyses, the results agree within 25.1%, 42.5%, and 28.8%, respectively, demonstrating that this synergy of scanning active remote sensing systems can be used to derive surface fluxes with reasonable accuracy.

54 ENVIRONMENTAL SCIENCES↗

A Double-Wiebe Function for Reactivity Controlled Compression Ignition Combustion Using Reformate Diesel

Abstract Reactivity controlled compression ignition (RCCI) combustion has previously been proposed as a method to achieve high fuel conversion efficiency and reduce engine emissions. A single-fuel RCCI combustion strategy can have decreased fuel system complexity by using a reformate fuel for port fuel injection and the parent fuel (diesel) for direct injection. This paper presents a one-dimensional computational model of a compression ignition engine with single-fuel RCCI. A Wiebe function is used to predict the combustion process by representing the mass fraction burned (MFB) on a crank angle resolved basis. One single-Wiebe function (SWF) and two double-Wiebe functions(DWFs) were fitted to experimentally derive MFB data using the least-square method. The fitted results were compared with MFBs calculated from experimental data to verify the accuracy. The SWF did not fully capture the MFB curve with high fidelity while the detailed DWF captured the MFB curve within a root mean square error of 1.4%. The reduced double-Wiebe function (RDWF) also resulted in a predicted combustion profile with similar accuracy. Hence, the RDWF was used in a GT-power thermodynamic study to understand the effects of the low-temperature heat release (LTHR) fraction and combustion phasing on combustion characteristics. At optimum phasing of 5–10 crank angle degree after the top dead center, increasing the LTHR fraction from 20% to 60% resulted in the fuel conversion efficiency increasing from 39.5% to 41.1%, thus suggesting that the reformate fuel-based RCCI strategy is viable to unlock improved combustion performance.

Energy & Fuels↗

Spherical fitting for oscilloscope film trace reading

This report describes a spherical fitting algorithm that can be used for estimating the oscilloscope film trace coordinates. This algorithm was adopted from a least square algebraic spherical fitting algorithm developed by the author earlier for photovoltaic velocimetry applications. Preliminary assessment suggests that the algorithm may be used for optimizing the coordinates manually digitized by the user, or may be used as an uncertainty estimator for the user’s film reading uncertainties.

97 MATHEMATICS AND COMPUTING↗

Parameterization of solar cells

The aggregation (sorting) of the individual solar cells into an array is commonly based on a single operating point on the current-voltage (I-V) characteristic curve. An alternative approach for cell performance prediction and cell screening is provided by modeling the cell using an equivalent electrical circuit, in which the parameters involved are related to the physical phenomena in the device. These analytical models may be represented by a double exponential I-V characteristic with seven parameters, by a double exponential model with five parameters, or by a single exponential equation with four or five parameters. In this article we address issues concerning methodologies for the determination of solar cell parameters based on measured data points of the I-V characteristic, and introduce a procedure for screening of solar cells for arrays. We show that common curve fitting techniques, e.g., least squares, may produce many combinations of parameter values while maintaining a good fit between the fitted and measured I-V characteristics of the cell. Therefore, techniques relying on curve fitting criteria alone cannot be directly used for cell parameterization. We propose a consistent procedure which takes into account the entire set of parameter values for a batch of cells. This procedure is based on a definition of a mean cell representing the batch, and takes into account the relative contribution of each parameter to the overall goodness of fit. The procedure is demonstrated on a batch of 50 silicon cells for Space Station Freedom.

Appelbaum, J.↗

Some observations regarding the statistical determination of stress rupture regression lines

Observations concerning the statistical evaluation of creep data are presented. Methods currently employed in the determination of stress rupture regression lines can result in conflicting and necessarily invalid results. Anomalous behavior is principally associated with the selection of the dependent variable. However, it is the least squares method of curve fitting which introduces regression bias. Methods to improve the validity of least squares regressions are suggested.

Pizzo, P. P.↗

A curve-fitting technique for the determination of gain and saturation intensity in homogeneously-broadened gas lasers: Application to a waveguide CO2 laser

A least squares curve-fitting algorithm is derived which allows the simultaneous estimation of the small signal gain and the saturation intensity from an arbitrary number of data points relating power output to the incidence angle of an internal coupling plate. The method is used to study the dependence of the two parameters on tube pressure and discharge current in a waveguide CO2 laser having a 2 mm diameter capillary. It is found that, at pressures greater than 28 torr, rising CO2 temperature degrades the small signal gain at current levels as low as three milliamperes.

Degnan, J. J.↗

Fitting Polynomial Equations to Curves and Surfaces

FIT is computer program for interactively determining least-squares polynomial equations that fit user-supplied data. Finds leastsquares fits for functions of two independent variables. Interactive graphical and editing capabilities in FIT enables user to control polynomial equations to be fitted to data arising from most practical applications. FIT written in FORTRAN and COMPASS.

Arbuckle, P. D.↗

An interactive approach to surface-fitting complex geometries for flowfield applications

Numerical flowfield methods require a geometry subprogram which can calculate body coordinates, slopes, and radii of curvature for typical aircraft and spacecraft configurations. The objective of this paper is to develop a new surface-fitting technique which addresses two major problems with existing geometry packages: computer storage requirements and the time required of the user for the initial set-up of the geometry model. In the present method, coordinates of cross sections are fit in a least-squares sense using segments of general conic sections. After fitting each cross section, the next step is to blend the cross-sectional curve-fits in the longitudinal direction using general conics to fit specific meridional half-planes. For the initial setup of the geometry model, an interactive, completely menu-driven computer code has been developed to allow the user to make modifications to the initial fit for a given cross section or meridional cut. Graphic displays are provided to assist the user in the visualization of the effect of each modification. The completed model may be viewed from any angle using the code's three-dimensional graphics package. Geometry results for the modeling of the Space Shuttle and a proposed Aeroassist Flight Experiment (AFE) geometry are presented, in addition to calculated heat-transfer rates based on these models.

Dejarnette, Fred R.↗

Multiwavelength Pyrometry To Correct For Reflections

Computerized curve fitting yields more and better information on thermal radiation. Multiwavelength pyrometry involves measurement of spectrum of thermal radiation emitted by and reflected from specimen. Auxiliary source enables determination of spectral reflectance of specimen. Spectral reflectance and temperature of specimen obtained from spectral measurements by nonlinear least-squares curve-fitting routine.

Ng, Daniel↗

Dynamic Differential Image Circle Diameter Measurement Precision Assessment: Application to Burning Droplets

Dynamic measurement precision assessment has been achieved for a differential circle measurement application. Differential circle diameter measurement, in image analysis, typically requires fitting a circle model that optimizes for image distortions, defects or occlusions. The differential task occurs when precise measurements of diameter change are required given object size variation with time. An automated system was designed to provide diameter measurements and associated measurement precision of images of a fuel droplet undergoing combustion in zero gravity for the FLEX-2 dataset. An image gradient-based, least-squares boundary point fitting method to a circle or ellipse model is used for diameter measurement. The presence of soot aggregates poses significant challenges for diameter measurements when it occludes part of the droplet boundary. The precision of the diameter measurements depends upon the image quality. Using synthetic image simulations that model the soot behavior, we developed a model based on image quality measures that assesses the measurement precision for each individual diameter measurement. Thus, diameter measurements with precision assessments were made available for follow-up scientific analysis. As a result, the algorithm's success rate for measurable runs was 98%. In cases of limited occlusion, a measurement precision of ±0.2 pixels for the FLEX-2 dataset was achieved.

42 ENGINEERING↗

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗