Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear least squares”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

High Resolution Spectroscopy to Support Atmospheric Measurements

The major research activities performed during the cooperative agreement enhanced our spectroscopic knowledge of molecules of atmospheric interest such as H2O (water vapor), O3 (ozone), HCN (hydrogen cyanide), CH4 (methane), NO2 (nitrogen dioxide) and CO (carbon monoxide). The data required for the analyses were obtained from two different Fourier Transform Spectrometers (FTS); one of which is located at the National Solar Observatory (NSO) on Kitt Peak, Arizona and the other instrument is located at the Pacific Northwest National Laboratories (PNNL) at Richland, Washington. The data were analyzed using a modified multispectrum nonlinear least squares fitting algorithm developed by Dr. D. Chris Benner of the College of William and Mary. The results from these studies made significant improvements in the line positons and intensities for these molecules. The measurements of pressure broadening and pressure induced line shift coefficients and the temperature dependence of pressure broadening and pressure induced shift coefficients for hundreds of infrared transitions of HCN, CO3 CH4 and H2O were also performed during this period. Results from these studies have been used for retrievals of stratospheric gas concentration profiles from data collected by several Upper Atmospheric Research Satellite (UARS) infrared instruments as well as in the analysis of high resolution atmospheric spectra such as those acquired by space-based, ground-based, and various balloon- and aircraft-borne experiments. Our results made significant contributions in several updates of the HITRAN (HIgh resolution TRANsmission) spectral line parameters database. This database enjoys worldwide recognition in research involving diversified scientific fields. The research conducted during the period 2003-2006 has resulted in publications given in this paper. In addition to Journal publications, several oral and poster presentations were given at various Scientific conferences within the United States and abroad.

Venkataraman, Malathy Devi↗
NASA NTRS2006

Temperature-Dependence of Air-Broadened Line Widths and Shifts in the nu3 Band of Ozone

The 9.6-micron bands of O3 are used by many remote-sensing experiments for retrievals of terrestrial atmospheric ozone concentration profiles. Line parameter errors can contribute significantly to the total errors in these retrievals, particularly for nadir-viewing. The McMath-Pierce Fourier transform spectrometer at the National Solar Observatory on Kitt Peak was used to record numerous high-resolution infrared absorption spectra of O3 broadened by various gases at temperatures between 160 and 300 K. Over 30 spectra were analyzed simultaneously using a multispectrum nonlinear least squares fitting technique to determine Lorentz air-broadening and pressure-induced shift coefficients along with their temperature dependences for selected transitions in the 3 fundamental band of (16)O3. We compare the present results with other measurements reported in the literature and with the ozone parameters on the 2000 and 2004 editions of the HITRAN database.

Smith, Mary A. H.↗
NASA NTRS2006

Processing TES Level-2 Data

TES Level 2 Subsystem is a set of computer programs that performs functions complementary to those of the program summarized in the immediately preceding article. TES Level-2 data pertain to retrieved species (or temperature) profiles, and errors thereof. Geolocation, quality, and other data (e.g., surface characteristics for nadir observations) are also included. The subsystem processes gridded meteorological information and extracts parameters that can be interpolated to the appropriate latitude, longitude, and pressure level based on the date and time. Radiances are simulated using the aforementioned meteorological information for initial guesses, and spectroscopic-parameter tables are generated. At each step of the retrieval, a nonlinear-least-squares- solving routine is run over multiple iterations, retrieving a subset of atmospheric constituents, and error analysis is performed. Scientific TES Level-2 data products are written in a format known as Hierarchical Data Format Earth Observing System 5 (HDF-EOS 5) for public distribution.

Poosti, Sassaneh↗
NASA NTRS2004

Time-Average Measurement of Velocity, Density, Temperature, and Turbulence Using Molecular Rayleigh Scattering

Measurement of time-averaged velocity, density, temperature, and turbulence in gas flows using a nonintrusive, point-wise measurement technique based on molecular Rayleigh scattering is discussed. Subsonic and supersonic flows in a 25.4-mm diameter free jet facility were studied. The developed instrumentation utilizes a Fabry-Perot interferometer to spectrally resolve molecularly scattered light from a laser beam passed through a gas flow. The spectrum of the scattered light contains information about velocity, density, and temperature of the gas. The technique uses a slow scan, low noise 16-bit depth CCD camera to record images of the fringes formed by Rayleigh scattered light passing through the interferometer. A kinetic theory model of the Rayleigh scattered light is used in a nonlinear least squares fitting routine to estimate the unknown parameters from the fringe images. The ability to extract turbulence information from the fringe image data proved to be a challenge since the fringe is broadened by not only turbulence, but also thermal fluctuations and aperture effects from collecting light over a range of scattering angles. Figure 1 illustrates broadening of a Rayleigh spectrum typical of flow conditions observed in this work due to aperture effects and turbulence for a scattering angle, chi(sub s), of 90 degrees, f/3.67 collection optics, mean flow velocity, u(sub k), of 300 m/s, and turbulent velocity fluctuations, sigma (sub uk), of 55 m/s. The greatest difficulty in processing the image data was decoupling the thermal and turbulence broadening in the spectrum. To aid in this endeavor, it was necessary to seed the ambient air with smoke and dust particulates; taking advantage of the turbulence broadening in the Mie scattering component of the spectrum of the collected light (not shown in the figure). The primary jet flow was not seeded due to the difficulty of the task. For measurement points lacking particles, velocity, density, and temperature information could reliably be recovered, however the turbulence estimates contained significant uncertainty. Resulting flow parameter estimates are presented for surveys of Mach 0.6, 0.95, and 1.4 jet flows. Velocity, density, and temperature were determined with accuracies of 5 m/s, 1.5%, and 1%, respectively, in flows with no particles present, and with accuracies of 5 m/s, 1-4%, and 2% in flows with particles. Comparison with hotwire data for the Mach 0.6 condition demonstrated turbulence estimates with accuracies of about 5 m/s outside the jet core where Mie scattering from dust/smoke particulates aided in the estimation of turbulence. Turbulence estimates could not be recovered with any significant accuracy for measurement points where no particles were present.

Mielke, Amy F.↗
NASA NTRS2004

Time-Average Molecular Rayleigh Scattering Technique for Measurement of Velocity, Denisty, Temperature, and Turbulence Intensity in High Speed Nozzle Flows

A molecular Rayleigh scattering based flow diagnostic is developed to measure time average velocity, density, temperature, and turbulence intensity in a 25.4-mm diameter nozzle free jet facility. The spectrum of the Rayleigh scattered light is analyzed using a Fabry-Perot interferometer operated in the static imaging mode. The resulting fringe pattern containing spectral information of the scattered light is recorded using a low noise CCD camera. Nonlinear least squares analysis of the fringe pattern using a kinetic theory model of the Rayleigh scattered light provides estimates of density, velocity, temperature, and turbulence intensity of the gas flow. Resulting flow parameter estimates are presented for an axial scan of subsonic flow at Mach 0.95 for comparison with previously acquired pitot tube data, and axial scans of supersonic flow in an underexpanded screeching jet. The issues related to obtaining accurate turbulence intensity measurements using this technique are discussed.

Mielke, Amy F.↗
NASA NTRS2012

Spectral Line Parameters Including Temperature Dependences of Self- and Air-Broadening in the 2 (left arrow) 0 Band of CO at 2.3 micrometers

Temperature dependences of pressure-broadened half-width and pressure-induced shift coefficients along with accurate positions and intensities have been determined for transitions in the 2<--0 band of C-12 O-16 from analyzing high-resolution and high signal-to-noise spectra recorded with two different Fourier transform spectrometers. A total of 28 spectra, 16 self-broadened and 12 air-broadened, recorded using high- purity (greater than or equal to 99.5% C-12-enriched) CO samples and CO diluted with dry air(research grade) at different temperatures and pressures, were analyzed simultaneously to maximize the accuracy of the retrieved parameters. The sample temperatures ranged from 150 to 298K and the total pressures varied between 5 and 700 Torr. A multispectrum nonlinear least squares spectrum fitting technique was used to adjust the rovibrational constants (G, B, D, etc.) and intensity parameters (including Herman-Wallis coefficients), rather than determining individual line positions and intensities. Self-and air-broadened Lorentz half-width coefficients, their temperature dependence exponents, self- and air-pressure-induced shift coefficients, their temperature dependences, self- and air-line mixing coefficients, their temperature dependences and speed dependence have been retrieved from the analysis. Speed-dependent line shapes with line mixing employing off-diagonal relaxation matrix element formalism were needed to minimize the fit residuals. This study presents a precise and complete set of spectral line parameters that consistently reproduce the spectrum of carbon monoxide over terrestrial atmospheric conditions.

Devi, V. Malathy↗
NASA NTRS2013

Air- and Self-Broadened Half Widths, Pressure-Induced Shifts, and Line Mixing in the Nu(sub 2) Band of (12)CH4

Lorentz self- and air-broadened half width and pressure-induced shift coefficients and their dependences on temperature have been measured from laboratory absorption spectra for nearly 130 transitions in the nu(sub 2) band of (12)CH4. In addition line mixing coefficients (using the relaxation matrix element formalism) for both self- and airbroadening were experimentally determined for the first time for a small number of transitions in this band. Accurate line positions and absolute line intensities were also determined. These parameters were obtained by analyzing high-resolution (approx. 0.003 to 0.01 per cm) laboratory spectra of high-purity natural CH4 and air-broadened CH4 recorded at temperatures between 226 and 297 K using the McMath-Pierce Fourier transform spectrometer (FTS) located at the National Solar Observatory on Kitt Peak, Arizona. A multispectrum nonlinear least squares technique was used to fit short (5-15 per cm) spectral intervals in 24-29 spectra simultaneously. Parameters were determined for nu(sub 2) transitions up to J" = 16. The variations of the measured broadening and shift parameters with the rotational quantum number index and tetrahedral symmetry species are examined. The present results are also compared with previous measurements available in the literature.

Smith, M. A. H.↗
NASA NTRS1979

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗
DOE OSTI2023

Estimating and Calibrating DER Model Parameters Using Levenberg–Marquardt Algorithm in Renewable Rich Power Grid

The proliferation of inverter-based distributed energy resources (IBDERs) has increased the number of control variables and dynamic interactions, leading to new grid control challenges. For stability analysis and designing appropriate protection controls, it is important that IBDER models are accurate. This paper focuses on the accurate estimation and parameter calibration of DER_A, a recently proposed aggregated IBDER model. In particular, we focus on the parameters of the reactive power–voltage regulation module. We formulate the problem of parameter tuning as a non-linear least square minimization problem and solve it using the Levenberg–Marquardt (LM) method. The LM method is primarily chosen due to its flexibility in adaptively selecting between the steepest descent and Gauss–Newton methods through a damping parameter. The LM approach is used to minimize the error between the actual measurements and the estimated response of the model. Further, the computational challenges posed by the numerical calculation of the Jacobian are tackled using a quasi-Newton root-finding approach. The proposed method is validated on a real feeder model in the northeastern part of the United States. The feeder is modeled in OpenDSS and the measurements thus obtained are fed to the DER_A model for calibration. The simulation results indicate that our approach is able to successfully calibrate the relevant model parameters quickly and with high accuracy, with a total sum of square error of 3.57 × 10 –7 .

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS1999

Optimized System Identification

In system identification, one usually cares most about finding a model whose outputs are as close as possible to the true system outputs when the same input is applied to both. However, most system identification algorithms do not minimize this output error. Often they minimize model equation error instead, as in typical least-squares fits using a finite-difference model, and it is seen here that this distinction is significant. Here, we develop a set of system identification algorithms that minimize output error for multi-input/multi-output and multi-input/single-output systems. This is done with sequential quadratic programming iterations on the nonlinear least-squares problems, with an eigendecomposition to handle indefinite second partials. This optimization minimizes a nonlinear function of many variables, and hence can converge to local minima. To handle this problem, we start the iterations from the OKID (Observer/Kalman Identification) algorithm result. Not only has OKID proved very effective in practice, it minimizes an output error of an observer which has the property that as the data set gets large, it converges to minimizing the criterion of interest here. Hence, it is a particularly good starting point for the nonlinear iterations here. Examples show that the methods developed here eliminate the bias that is often observed using any system identification methods of either over-estimating or under-estimating the damping of vibration modes in lightly damped structures.

Juang, Jer-Nan↗
NASA NTRS1985

An adaptive learning control system for large flexible structures

The objective of the research has been to study the design of adaptive/learning control systems for the control of large flexible structures. In the first activity an adaptive/learning control methodology for flexible space structures was investigated. The approach was based on using a modal model of the flexible structure dynamics and an output-error identification scheme to identify modal parameters. In the second activity, a least-squares identification scheme was proposed for estimating both modal parameters and modal-to-actuator and modal-to-sensor shape functions. The technique was applied to experimental data obtained from the NASA Langley beam experiment. In the third activity, a separable nonlinear least-squares approach was developed for estimating the number of excited modes, shape functions, modal parameters, and modal amplitude and velocity time functions for a flexible structure. In the final research activity, a dual-adaptive control strategy was developed for regulating the modal dynamics and identifying modal parameters of a flexible structure. A min-max approach was used for finding an input to provide modal parameter identification while not exceeding reasonable bounds on modal displacement.

Thau, F. E.↗
NASA NTRS2007

Numerical Technique for Analyzing Rotating Rake Mode Measurements in a Duct With Passive Treatment and Shear Flow

A technique is presented for the analysis of measured data obtained from a rotating microphone rake system. The system is designed to measure the interaction modes of ducted fans. A Fourier analysis of the data from the rotating system results in a set of circumferential mode levels at each radial location of a microphone inside the duct. Radial basis functions are then least-squares fit to this data to obtain the radial mode amplitudes. For ducts with soft walls and mean flow, the radial basis functions must be numerically computed. The linear companion matrix method is used to obtain both the eigenvalues of interest, without an initial guess, and the radial basis functions. The governing equations allow for the mean flow to have a boundary layer at the wall. In addition, a nonlinear least-squares method is used to adjust the wall impedance to best fit the data in an attempt to use the rotating system as an in-duct wall impedance measurement tool. Simulated and measured data are used to show the effects of wall impedance and mean flow on the computed results.

Dahl, Milo D.↗
NASA NTRS2000

Modeling Step-Strain Relaxation and Cyclic Deformations of Elastomers

Data for step-strain relaxation and cyclic compressive deformations of highly viscous short elastomer cylinders are modeled using a large strain rubber viscoelastic constitutive theory with a rate-independent friction stress term added. In the tests, both small and large amplitude cyclic compressive strains, in the range of 1% to 10%, were superimposed on steady state compressed strains, in the range of 5% to 20%, for frequencies of 1 and 10 Hz. The elastomer cylinders were conditioned prior to each test to soften them. The constants in the viscoclastic-friction constitutive theory are determined by employing a nonlinear least-squares method to fit the analytical stresses for a Maxwell model, which includes friction, to measured relaxation stresses obtained from a 20% step-strain compression test. The simulation of the relaxation data with the nonlinear model is successful at compressive strains of 5%, 10%, 15%, and 20%. Simulations of hysteresis stresses for enforced cyclic compressive strains of 20% +/- 5% are made with the model calibrated by the relaxation data. The predicted hysteresis stresses are lower than the measured stresses.

Johnson, A.R.↗
NASA NTRS1991

Least-squares/parabolized Navier-Stokes procedure for optimizing hypersonic wind tunnel nozzles

A new procedure is demonstrated for optimizing hypersonic wind-tunnel-nozzle contours. The procedure couples a CFD computer code to an optimization algorithm, and is applied to both conical and contoured hypersonic nozzles for the purpose of determining an optimal set of parameters to describe the surface geometry. A design-objective function is specified based on the deviation from the desired test-section flow-field conditions. The objective function is minimized by optimizing the parameters used to describe the nozzle contour based on the solution to a nonlinear least-squares problem. The effect of the changes in the nozzle wall parameters are evaluated by computing the nozzle flow using the parabolized Navier-Stokes equations. The advantage of the new procedure is that it directly takes into account the displacement effect of the boundary layer on the wall contour. The new procedure provides a method for optimizing hypersonic nozzles of high Mach numbers which have been designed by classical procedures, but are shown to produce poor flow quality due to the large boundary layers present in the test section. The procedure is demonstrated by finding the optimum design parameters for a Mach 10 conical nozzle and a Mach 6 and a Mach 15 contoured nozzle.

Korte, John J.↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Subspace Methods in Multi-Parameter Seismic Full Waveform Inversion

In full waveform inversion (FWI) high-resolution subsurface model parameters are sought. FWI is normally treated as a nonlinear least-squares inverse problem, in which the minimum of the corresponding misfit function is found by updating the model parameters. When multiple elastic or acoustic properties are solved for, simple gradient methods tend to confuse parameter classes. This is referred to as parameter cross-talk; it leads to incorrect model solutions, poor convergence and strong dependence on the scaling of the different parameter types. Determining step lengths in a subspace domain, rather than directly in terms of gradients of different parameters, is a potentially valuable approach to address this problem. The particular subspace used can be defined over a span of different sets of data or different parameter classes, provided it involves a small number of vectors compared to those contained in the whole model space. Additionally, in a subspace method, the basis vectors are defined first, and a local minimum is found in the space spanned by these. We examine the application of the subspace method within acoustic FWI in determining simultaneously updates for velocity and density. We first discuss the choice of basis vectors to construct the spanned space, from linear updates by distinguishing only the contributions of different parameter classes towards nonlinear updates by adding the contributions of higher-order perturbations of each parameter class. The numerical character of FWI solutions generated via subspace methods involving different basis vectors is then analyzed and compared with traditional FWI methods. The subspace methods can provide better reconstructions of the model, especially for the velocity, as well as improved convergence rates, while the computational costs are still comparable with the traditional FWI methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗