Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Fast, Two-dimensional Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization

Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of data points has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions where, typically, the second dimension is small. This method scales linearly with the total number of data points when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2006

An Empirical Comparison between Two Recursive Filters for Attitude and Rate Estimation of Spinning Spacecraft

The advantages of inducing a constant spin rate on a spacecraft are well known. A variety of science missions have used this technique as a relatively low cost method for conducting science. Starting in the late 1970s, NASA focused on building spacecraft using 3-axis control as opposed to the single-axis control mentioned above. Considerable effort was expended toward sensor and control system development, as well as the development of ground systems to independently process the data. As a result, spinning spacecraft development and their resulting ground system development stagnated. In the 1990s, shrinking budgets made spinning spacecraft an attractive option for science. The attitude requirements for recent spinning spacecraft are more stringent and the ground systems must be enhanced in order to provide the necessary attitude estimation accuracy. Since spinning spacecraft (SC) typically have no gyroscopes for measuring attitude rate, any new estimator would need to rely on the spacecraft dynamics equations. One estimation technique that utilized the SC dynamics and has been used successfully in 3-axis gyro-less spacecraft ground systems is the pseudo-linear Kalman filter algorithm. Consequently, a pseudo-linear Kalman filter has been developed which directly estimates the spacecraft attitude quaternion and rate for a spinning SC. Recently, a filter using Markley variables was developed specifically for spinning spacecraft. The pseudo-linear Kalman filter has the advantage of being easier to implement but estimates the quaternion which, due to the relatively high spinning rate, changes rapidly for a spinning spacecraft. The Markley variable filter is more complicated to implement but, being based on the SC angular momentum, estimates parameters which vary slowly. This paper presents a comparison of the performance of these two filters. Monte-Carlo simulation runs will be presented which demonstrate the advantages and disadvantages of both filters.

Harman, Richard R.↗
DOE OSTI2023

Network-Wide Traffic Signal Control Using Bilinear System Modeling and Adaptive Optimization

This study proposes a new multi-input multi-output optimal bilinear signal control method in which a bilinear dynamic model approximation is used to capture the nonlinear dynamics of the urban traffic networks. With signal green time splits as the control input and traffic delay changes as the output for each intersections in the network, a bilinear system model was developed, which, on the basis of linear system modeling, takes interactions among traffic delays and signal timing splits into consideration. Based on the bilinear system modeling framework, we conducted two steps in each time interval to derive traffic control strategies: (1) we used the normalized least-squared algorithm to estimate system parameters; and (2) we solved an online optimization problem to obtain the updated traffic control inputs for the signal timing that minimizes future traffic delays. We evaluated the proposed method in a microscopic traffic simulation environment (VISSIM) with a 35-intersection network of Bellevue city in Washington. Two different traffic demand patterns: (1) normal traffic demands; and (2) time-varying traffic demands were simulated to compare the performance of different control strategies. Experimental results show that (1) the proposed bilinear system model can better describe traffic system dynamics than linear-model based methods, such as our previously developed linear-quadratic regulator control; and (2) the proposed method outperforms the state-of-the-art signal control strategies, namely the max-pressure and the self-organizing traffic light control methods. We have also shown that the proposed method is applicable to all other possible network layouts and signal controller phasing structures.

42 ENGINEERING↗
NASA NTRS2003

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗
NASA NTRS1992

Preliminary subsonic aerodynamic model for simulation studies of the HL-20 lifting body

A nonlinear, six-degree-of-freedom aerodynamic model for an early version of the HL-20 lifting body is described and compared with wind tunnel data upon which it is based. Polynomial functions describing most of the aerodynamic parameters are given and tables of these functions are presented. Techniques used to arrive at these functions are described. Basic aerodynamic coefficients were modeled as functions of angles of attack and sideslip. Vehicle lateral symmetry was assumed. Compressibility (Mach) effects were ignored. Control-surface effectiveness was assumed to vary linearly with angle of deflection and was assumed to be invariant with the angle of sideslip. Dynamic derivatives were obtained from predictive aerodynamic codes. Landing-gear and ground effects were scaled from Space Shuttle data. The model described is provided to support pilot-in-the-loop simulation studies of the HL-20. By providing the data in tabular format, the model is suitable for the data interpolation architecture of many existing engineering simulation facilities. Because of the preliminary nature of the data, however, this model is not recommended for study of the absolute performance of the HL-20.

Jackson, E. Bruce↗
NASA NTRS1995

Performance analysis of a GPS Interferometric attitude determination system for a gravity gradient stabilized spacecraft

The performance of an unaided attitude determination system based on GPS interferometry is examined using linear covariance analysis. The modelled system includes four GPS antennae onboard a gravity gradient stabilized spacecraft, specifically the Air Force's RADCAL satellite. The principal error sources are identified and modelled. The optimal system's sensitivities to these error sources are examined through an error budget and by varying system parameters. The effects of two satellite selection algorithms, Geometric and Attitude Dilution of Precision (GDOP and ADOP, respectively) are examined. The attitude performance of two optimal-suboptimal filters is also presented. Based on this analysis, the limiting factors in attitude accuracy are the knowledge of the relative antenna locations, the electrical path lengths from the antennae to the receiver, and the multipath environment. The performance of the system is found to be fairly insensitive to torque errors, orbital inclination, and the two satellite geometry figures-of-merit tested.

Stoll, John C.↗
NASA NTRS1973

Concepts for a theoretical and experimental study of lifting rotor random loads and vibrations (identification of lifting rotor system parameters from transient response data), Phase 7-B

System identification methods have been applied to rotorcraft to estimate stability derivatives from transient flight control response data. While these applications assumed a linear constant coefficient representation of the rotorcraft, the computer experiments used transient responses in flap-bending and torsion of a rotor blade at high advance ratio which is a rapidly time varying periodic system. It was found that a simple system identification method applying a linear sequential estimator also called least square estimator or equation of motion estimator, is suitable for this periodic system and can be used directly if only the acceleration data are noise polluted. In the case of noise being present also in the state variable data the direct application of the estimator gave poor results.

Hohenemser, K. H.↗
DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Thermal Characterization of Functionally Graded Materials: Design of Optimum Experiments

This paper is a study of optimal experiment design applied to the measure of thermal properties in functionally graded materials. As a first step, a material with linearly-varying thermal properties is analyzed, and several different tran- sient experimental designs are discussed. An optimality criterion, based on sen- sitivity coefficients, is used to identify the best experimental design. Simulated experimental results are analyzed to verify that the identified best experiment design has the smallest errors in the estimated parameters. This procedure is general and can be applied to design of experiments for a variety of materials.

Cole, Kevin D.↗
DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Models of bars − II. Exponential profiles

ABSTRACT We present a new model for galactic bars with exponentially falling major axis luminosity profiles and Gaussian cross-sections. This is based on the linear superposition of Gaussian potential–density pairs with an exponential weight function, using an extension of the method originally introduced by Long & Murali. We compute the density, potential, and forces, using Gaussian quadrature. These quantities are given as explicit functions of position. There are three independent scaled bar parameters that can be varied continuously to produce bespoke bars of a given mass and shape. We categorize the effective potential by splitting a reduced parameter space into six regions. Unusually, we find bars with three stable Lagrange points on the major axis are possible. Our model reveals a variety of unexpected orbital structure, including a bifurcating x1 orbit coexisting with a stable x4 orbit. Propeller orbits are found to play a dominant role in the orbital structure, and we find striking similarities between our bar configuration and the model of Kaufmann & Contopoulos. We find a candidate orbital family, sired from the propeller orbits, that may be responsible for the observed high-velocity peaks in the Milky Way’s bar. As a cross-check, we inspect, for the first time, the proper motions of stars in the high-velocity peaks, which also match our suggested orbital family well. This work adds to the increasing body of evidence that real galactic bars may be supported at least partly by propeller orbits rather than solely by elliptical-like orbits of the x1 family.

McGough, D. P.↗
NASA NTRS1961

A Threshold Theory for Phase-Locked Loops

A model of a phase-locked loop has been developed which is valid for all signal-to-noise ratios. The model is in the form of a nonlinear feedback system with randomly time-varying parameters. The analysis considers two operating regions. In low signal-to-noise ratio regions, the important consideration is stability. We want to study the asymptotic stability in the mean of a nonlinear system. It follows directly that a necessary condition for asymptotic stability of any nonlinear system is that a linearized model about some equilibrium point be asymptotically stable. By considering all possible equilibrium points, we can find an upper bound on the value of noise density which makes the system unstable. This upper bound represents a threshold value for system operation. In high signal-to-noise ratio regions, our results provide an exact statistical. description of system behavior. Therefore, knowledge of the spectrum of the signal and noise may be used to optimize the system configuration.

Van Trees, H. L.↗
NASA NTRS1985

Crystallinity determination by curvefit procedure for a semi-crystalline polymer

Wide angle X-ray scattering (WAXS) data from poly(etheretherketone) (PEEK) was resolved into a crystalline contribution represented as four reflections and an amorphous contribution represented as a broad, smoothly varying curve. The crystalline scatter is described as a linear combination of Cauchy and Gaussian functions while that of the amorphous halo is expressed as a cubic polynomial. Statistical analysis of the measured scattered intensity from an amorphous specimen with that calculated from the cubic polynomial, as a function of the combination parameter (fraction of Cauchy and Gaussian functions), suggests that the crystalline fraction of the polymer specimen studied is about 0.39.

Wakelyn, N. T.↗
NASA NTRS2005

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗
DOE OSTI2024

Bayesian reduced-order deep learning surrogate model for dynamic systems described by partial differential equations

We propose a reduced-order deep-learning surrogate model for dynamic systems described by time-dependent partial differential equations. This method employs space–time Karhunen–Loève expansions (KLEs) of the state variables and space-dependent KLEs of space-varying parameters to identify the reduced (latent) dimensions. Subsequently, a deep neural network (DNN) is used to map the parameter latent space to the state variable latent space. An approximate Bayesian method is developed for uncertainty quantification (UQ) in the proposed KL-DNN surrogate model. The KL-DNN method is tested for the linear advection–diffusion and nonlinear diffusion equations, and the Bayesian approach for UQ is compared with the deep ensembling (DE) approach, commonly used for quantifying uncertainty in DNN models. It was found that the approximate Bayesian method provides a more informative distribution of the PDE solutions in terms of the coverage of the reference PDE solutions (the percentage of nodes where the reference solution is within the confidence interval predicted by the UQ methods) and log predictive probability. The DE method is found to underestimate uncertainty and introduce bias. For the nonlinear diffusion equation, we compare the KL-DNN method with the Fourier Neural Operator (FNO) method and find that KL-DNN is 10% more accurate and needs less training time than the FNO method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Process parameters and system responses in friction extrusion

We report friction extrusion is an emerging manufacturing process that imparts improved properties and unique microstructures to materials. In this study, the effects of die features (i.e., the absence or presence of die face scrolling) and process parameters (i.e., die rotational speeds and feed rates) on response variables (e.g., extrusion force, torque, power, and temperature) were investigated. During this work, wires with good integrity and similar geometries were extruded over a range of experimental trials with different rotational speeds and feed rates at a constant extrusion ratio. The extrusion force linearly varies with die advance per revolution, which is the ratio of feed rate and rotational speed. The extrusion energy per kilogram mass of extrudate varies non-linearly and asymptotically such that energy consumption is highest at the lowest feed rate and lowest at the highest feed rate. Additionally relationships among other variables (i.e., die advance per revolution, torque, power heat input, and temperature) also were evaluated to elucidate their interdependences.

42 ENGINEERING↗
NASA NTRS1991

Herman-Wallis factors in the C2H2nu5 infrared fundamental near 14 microns

The presence of acetylene has been confirmed for some time in the atmospheres of the outer planets Jupiter, Saturn, Neptune, and Saturn's satellite Titan. For these atmospheres, the determination of C2H2 abundances using its strong nu5 fundamental requires laboratory line position and intensity measurements. The 1-m Fourier transform spectrometer at McMath solar telescope of Kitt Peak National Observatory was used to measure C2H2 at an unapodized spectral resolution of 0.0025/cm. Synthetic spectra are generated by convolving a Voigt line shape with an instrument function and varying intensity parameters by means of a nonlinear least squares technique. Intensities of 37 nu5 lines spanning P18 through R20 were measured using 0.123 torr of gas in a 1-cm cell. A Herman-Wallis intensity correction parameter of 1.3(4) x 10 to the -3rd has been derived using a least squares linear fit.

Maguire, William C.↗
DOE OSTI2024

Induced Polarization of Clayey Rocks and Soils: Non‐Linear Complex Conductivity Models

Abstract The past decades have witnessed the increased applications of induced polarization (IP) method in the critical zone studies with ubiquitous clay minerals. Although IP outperforms traditional electrical and electromagnetic methods through its unique ability to measure quadrature conductivity, the nonlinearity that quadrature conductivity behaves with salinities and frequencies greatly tortures IP practitioners, as (a) salinity‐dependency makes the quadrature conductivity a varyingly unstable parameter to quantitatively estimate hydraulic properties and clay content; (b) frequency‐dependent Cole‐Cole and Debye/Warburg decomposition models, although mathematically sound, physically mingle the properties of pore water and clay minerals and are empirical in nature. From basic principles, we demonstrate that quadrature conductivity remains a hybrid property involving both clay and water, and develop relevant models to distinguish them. Our models are validated by theories, experiments, simulations, and comparisons, all of which proclaim considerable advantages over previous models and offer the prospect of quantitative applications.

58 GEOSCIENCES↗