Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Generalized Linear Model”

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 307 records · Page 17

A decision and control multi-axis pilot model based on an urgency for action concept

An all digital multiple axis, multiple input, multiple output pilot model flies five different tasks in aircraft of different classes while encountering turbulence represented by the Dryden spectral model. Six degree of freedom linearized aircraft equations of motion are used together with a generalized stability augmentation system. The pilot model consists of two parts: (1) a decision maker, and (2) a control action implementer. This decision maker selects the critical variable and predicts the error at action implementation. The control action implementer decides the magnitude of control to be applied and applies it. This system has been successfully compared with hybrid/analog man in the loop simulations of the F-5, A-7, T-33, and 707 aircraft, thus validating the all digital simulation and the decision and control pilot model.

Pollard, J. J.↗

Progressive Aerodynamic Model Identification From Dynamic Water Tunnel Test of the F-16XL Aircraft

Development of a general aerodynamic model that is adequate for predicting the forces and moments in the nonlinear and unsteady portions of the flight envelope has not been accomplished to a satisfactory degree. Predicting aerodynamic response during arbitrary motion of an aircraft over the complete flight envelope requires further development of the mathematical model and the associated methods for ground-based testing in order to allow identification of the model. In this study, a general nonlinear unsteady aerodynamic model is presented, followed by a summary of a linear modeling methodology that includes test and identification methods, and then a progressive series of steps suggesting a roadmap to develop a general nonlinear methodology that defines modeling, testing, and identification methods. Initial steps of the general methodology were applied to static and oscillatory test data to identify rolling-moment coefficient. Static measurements uncovered complicated dependencies of the aerodynamic coefficient on angle of attack and sideslip in the stall region making it difficult to find a simple analytical expression for the measurement data. In order to assess the effect of sideslip on the damping and unsteady terms, oscillatory tests in roll were conducted at different values of an initial offset in sideslip. Candidate runs for analyses were selected where higher order harmonics were required for the model and where in-phase and out-of-phase components varied with frequency. From these results it was found that only data in the angle-of-attack range of 35 degrees to 37.5 degrees met these requirements. From the limited results it was observed that the identified models fit the data well and both the damping-in-roll and the unsteady term gain are decreasing with increasing sideslip and motion amplitude. Limited similarity between parameter values in the nonlinear model and the linear model suggest that identifiability of parameters in both terms may be a problem. However, the proposed methodology can still be used with careful experiment design and carefully selected values of angle of attack, sideslip, amplitude, and frequency of the oscillatory data.

Murphy, Patrick C.↗

A surrogate-based uncertainty quantification with quantifiable errors

Surrogate models are often employed to reduce the computational cost required to complete uncertainty quantification, where one is interested in propagating input parameters uncertainties throughout a complex engineering model to estimate responses uncertainties. An improved surrogate construction approach is introduced here which places a premium on reducing the associated computational cost. Unlike existing methods where the surrogate is constructed first, then employed to propagate uncertainties, the new approach combines both sensitivity and uncertainty information to render further reduction in the computational cost. Mathematically, the reduction is described by a range finding algorithm that identifies a subspace in the parameters space, whereby parameters uncertainties orthogonal to the subspace contribute negligible amount to the propagated uncertainties. Moreover, the error resulting from the reduction can be upper-bounded. The new approach is demonstrated using a realistic nuclear assembly model and compared to existing methods in terms of computational cost and accuracy of uncertainties. Although we believe the algorithm is general, it will be applied here for linear-based surrogates and Gaussian parameters uncertainties. The generalization to nonlinear models will be detailed in a separate article. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Two dimensional thick center vortex model

The potential between static color source is calculated in the SU (3) gauge group by introducing a two dimensional vortex flux. To generalize the model, the length of the Wilson loop is equal to R oriented along the x axis, and the vortex flux is considered as a function of x and y. The comparison between the generalized model and the original one shows that the intermediate linear regime is increased significantly and better agreement with Casimir scaling is achieved. Furthermore, the model is applied to calculate the potential between baryons.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Polyhedral approximation in mixed-integer convex optimization

Generalizing both mixed-integer linear optimization and convex optimization, mixed-integer convex optimization possesses broad modeling power but has seen relatively few advances in general-purpose solvers in recent years. Here, we intend to provide a broadly accessible introduction to our recent work in developing algorithms and software for this problem class. Our approach is based on constructing polyhedral outer approximations of the convex constraints, resulting in a global solution by solving a finite number of mixed-integer linear and continuous convex subproblems. The key advance we present is to strengthen the polyhedral approximations by constructing them in a higher-dimensional space. In order to automate this extended formulation we rely on the algebraic modeling technique of disciplined convex programming (DCP), and for generality and ease of implementation we use conic representations of the convex constraints. Although our framework requires a manual translation of existing models into DCP form, after performing this transformation on the MINLPLIB2 benchmark library we were able to solve a number of unsolved instances and on many other instances achieve superior performance compared with state-of-the-art solvers like Bonmin, SCIP, and Artelys Knitro.

97 MATHEMATICS AND COMPUTING↗

ACCRUE—An Integral Index for Measuring Experimental Relevance in Support of Neutronic Model Validation

A key challenge for the introduction of any design changes, e.g., advanced fuel concepts, first-of-a-kind nuclear reactor designs, etc., is the cost of the associated experiments, which are required by law to validate the use of computer models for the various stages, starting from conceptual design, to deployment, licensing, operation, and safety. To achieve that, a criterion is needed to decide on whether a given experiment, past or planned, is relevant to the application of interest. This allows the analyst to select the best experiments for the given application leading to the highest measures of confidence for the computer model predictions. The state-of-the-art methods rely on the concept of similarity or representativity, which is a linear Gaussian-based inner-product metric measuring the angle—as weighted by a prior model parameters covariance matrix—between two gradients, one representing the application and the other a single validation experiment. This manuscript emphasizes the concept of experimental relevance which extends the basic similarity index to account for the value accrued from past experiments and the associated experimental uncertainties, both currently missing from the extant similarity methods. Accounting for multiple experiments is key to the overall experimental cost reduction by prescreening for redundant information from multiple equally-relevant experiments as measured by the basic similarity index. Accounting for experimental uncertainties is also important as it allows one to select between two different experimental setups, thus providing for a quantitative basis for sensor selection and optimization. The proposed metric is denoted by ACCRUE, short for Accumulative Correlation Coefficient for Relevance of Uncertainties in Experimental validation. Using a number of criticality experiments for highly enriched fast metal systems and low enriched thermal compound systems with accident tolerant fuel concept, the manuscript will compare the performance of the ACCRUE and basic similarity indices for prioritizing the relevance of a group of experiments to the given application.

97 MATHEMATICS AND COMPUTING↗

Modeling Aircraft Wing Loads from Flight Data Using Neural Networks

Neural networks were used to model wing bending-moment loads, torsion loads, and control surface hinge-moments of the Active Aeroelastic Wing (AAW) aircraft. Accurate loads models are required for the development of control laws designed to increase roll performance through wing twist while not exceeding load limits. Inputs to the model include aircraft rates, accelerations, and control surface positions. Neural networks were chosen to model aircraft loads because they can account for uncharacterized nonlinear effects while retaining the capability to generalize. The accuracy of the neural network models was improved by first developing linear loads models to use as starting points for network training. Neural networks were then trained with flight data for rolls, loaded reversals, wind-up-turns, and individual control surface doublets for load excitation. Generalization was improved by using gain weighting and early stopping. Results are presented for neural network loads models of four wing loads and four control surface hinge moments at Mach 0.90 and an altitude of 15,000 ft. An average model prediction error reduction of 18.6 percent was calculated for the neural network models when compared to the linear models. This paper documents the input data conditioning, input parameter selection, structure, training, and validation of the neural network models.

Allen, Michael J.↗

First Sagittarius A* Event Horizon Telescope Results. VIII. Physical Interpretation of the Polarized Ring

Abstract In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

79 ASTRONOMY AND ASTROPHYSICS↗

Closed Form Solution for Minimum Norm Model-Validating Uncertainty

A methodology in which structured uncertainty models are directly constructed from measurement data for use in robust control design of multivariable systems is proposed. The formulation allows a general linear fractional transformation uncertainty structure connections with respect to a given nominal model. Existence conditions are given, and under mild assumptions, a closed-form expression for the smallest norm structured uncertainty that validates the model is given. The uncertainty bound computation is simple and is formulated for both open and closed loop systems.

Lim, Kyong Been↗

Gravity Waves and Linear Inflation From Axion Monodromy

Wrapped branes in string compactifications introduce a monodromy that extends the field range of individual closed-string axions to beyond the Planck scale. Furthermore, approximate shift symmetries of the system naturally control corrections to the axion potential. This suggests a general mechanism for chaotic inflation driven by monodromy-extended closed-string axions. We systematically analyze this possibility and show that the mechanism is compatible with moduli stabilization and can be realized in many types of compactifications, including warped Calabi-Yau manifolds and more general Ricci-curved spaces. In this broad class of models, the potential is linear in the canonical inflaton field, predicting a tensor to scalar ratio r {approx} 0.07 accessible to upcoming cosmic microwave background (CMB) observations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Correlation function induced by a generalized diffusion equation with the presence of a harmonic potential

An integro-differential diffusion equation with linear force, based on the continuous time random walk model, is considered. The equation generalizes the ordinary and fractional diffusion equations, which includes short, intermediate and long-time memory effects described by the waiting time probability density function. Analytical expression for the correlation function is obtained and analyzed, which can be used to describe, for instance, internal motions of proteins. The result shows that the generalized diffusion equation has a broad application and it may be used to describe different kinds of systems. - Highlights: • Calculation of the correlation function. • The correlation function is connected to the survival probability. • The model can be applied to the internal dynamics of proteins.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy exchange in strongly coupled plasmas with electron drift

In this paper, the generalized viscoelastic collisional quantum hydrodynamic model is employed in order to investigate the linear dielectric response of a quantum plasma in the presence of strong electron-beam plasma interactions. The generalized Chandrasekhar's relativistic degeneracy pressure together with the electron-exchange and Coulomb interaction effects are taken into account in order to extend current research to a wide range of plasma number density relevant to big planetary cores and astrophysical compact objects. The previously calculated shear viscosity and the electron-ion collision frequencies are used for strongly coupled ion fluid. The effect of the electron-beam velocity on complex linear dielectric function is found to be profound. This effect is clearly interpreted in terms of the wave-particle interactions and their energy-exchange according to the sign of the imaginary dielectric function, which is closely related to the wave attenuation coefficient in plasmas. Such kinetic effect is also shown to be in close connection with the stopping power of a charged-particle beam in a quantum plasma. The effect of many independent plasma parameters, such as the ion charge-state, electron beam-velocity, and relativistic degeneracy, is shown to be significant on the growing/damping of plasma instability or energy loss/gain of the electron-beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Feedback stabilization and control of linear neutral systems

The first problem treated here is the realization and stabilization of linear neutral systems with discrete delays. It is shown that any autonomous linear neutral system with discrete delays is zero-state equivalent to an abstract linear system over a local ring of operators. Using the abstract model, the basic existence question for neutral realization is then settled. For general infinite dimensional linear systems, there is no precise analog of the finite dimensional state space isomorphism theorem. Because of this, the notion of spectral minimality must be introduced. For the case of single input-single output systems, realizations are obtained that are both minimal and spectrally minimal. Using the Cruz-Hale theory of stable D-operators, conditions are given that ensure that any poles introduced into the realization are strictly contained in the left half plane and indeed are characterized as characteristic values of the D-operator. The problem of the feedback stabilization of neutral systems is then considered using the abstract model. It is shown that, for neutral systems with commensurable delays and a stable D-operator in the sense of Cruz and Hale, Morses theorem (1976) on pole assignment over a PID implies stabilizability in the reachable case.

Tarn, T. J.↗

Electron kinetic effects on interferometry, polarimetry and Thomson scattering measurements in burning plasmas (invited)

At anticipated high electron temperatures in ITER, the effects of electron thermal motion on Thomson scattering (TS), toroidal interferometer/polarimeter (TIP), and poloidal polarimeter (PoPola) diagnostics will be significant and must be accurately treated. The precision of the previous lowest order linear in τ = T{sub e}/m{sub e}c{sup 2} model may be insufficient; we present a more precise model with τ{sup 2}-order corrections to satisfy the high accuracy required for ITER TIP and PoPola diagnostics. The linear model is extended from Maxwellian to a more general class of anisotropic electron distributions that allows us to take into account distortions caused by equilibrium current, ECRH, and RF current drive effects. The classical problem of the degree of polarization of incoherent Thomson scattered radiation is solved analytically exactly without any approximations for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. The results are discussed in the context of the possible use of the polarization properties of Thomson scattered light as a method of T{sub e} measurement relevant to ITER operational scenarios.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Electron kinetic effects on optical diagnostics in fusion plasmas

At anticipated high electron temperatures in ITER, the effects of electron thermal motion on Thomson scattering (TS), toroidal interferometer/polarimeter (TIP) and poloidal polarimeter (PoPola) diagnostics will be significant and must be accurately treated. We calculate electron thermal corrections to the interferometric phase and polarization state of an EM wave propagating along tangential and poloidal chords (Faraday and Cotton-Mouton polarimetry) and perform analysis of the degree of polarization for incoherent TS. The precision of the previous lowest order linear in τ = T{sub e}/m{sub e}c{sup 2} model may be insufficient; we present a more precise model with τ{sup 2}-order corrections to satisfy the high accuracy required for ITER TIP and PoPola diagnostics. The linear model is extended from Maxwellian to a more general class of anisotropic electron distributions that allows us to take into account distortions caused by equilibrium current, ECRH and RF current drive effects. The classical problem of degree of polarization of incoherent Thomson scattered radiation is solved analytically exactly without any approximations for the full range of incident polarizations, scattering angles, and electron thermal motion from non-relativistic to ultra-relativistic. The results are discussed in the context of the possible use of the polarization properties of Thomson scattered light as a method of T{sup e} measurement relevant to ITER operational scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ion thermal and dispersion effects in Farley-Buneman instabilities

Farley-Buneman modes are an example of the collisional instability, which is thought to be the dominant mechanism for the irregularities in low ionosphere region. Despite high collisionality due to electron-neutral and ion-neutral collisions, the kinetic effects associated with finite temperature are important for determination of the mode frequencies and growth rate. This is especially important for ion component that is largely unmagnetized due to low ion cyclotron frequency. The ion thermal effects are strongly pronounced for shorter wavelengths and are crucial for the growth rate cut-off at high wavenumbers. We develop an extended fluid model for ion dynamics to incorporate the effects of ion thermal motion. The model is based on the extended MHD model that includes the evolution equations for higher order moments such as ion viscosity and ion heat flux. We also develop the generalized Chapman-Enskog closure model that provides exact linear closures based on the linearized kinetic equation. The results of these models are compared and tested against the linear kinetic model. The dispersion of Farley-Buneman modes and growth rate behavior are investigated in the short wavelength region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spin-electron acoustic soliton and exchange interaction in separate spin evolution quantum plasmas

Separate spin evolution quantum hydrodynamics is generalized to include the Coulomb exchange interaction, which is considered as interaction between the spin-down electrons being in quantum states occupied by one electron. The generalized model is applied to study the non-linear spin-electron acoustic waves. Existence of the spin-electron acoustic soliton is demonstrated. Contributions of concentration, spin polarization, and exchange interaction to the properties of the spin electron acoustic soliton are studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Wave propagation in nonlinear disordered media

We analyze mechanisms and regimes of wave packet spreading in nonlinear disordered media. We discuss resonance probabilities, and predict a dynamical crossover from strong to weak chaos. The crossover is controlled by the ratio of nonlinear frequency shifts and the average eigenvalue spacing of eigenstates of the linear equations within one localization volume. We consider generalized models in higher lattice dimensions and obtain critical values for the nonlinearity power, the dimension, and norm density, which influence possible dynamical outcomes in a qualitative way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗