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 577 records · Page 32

Constructor selection system

Future space construction missions will involve both human and machine constructors. Selection of the optimum constructor mix requires a model of constructor capabilities and requirements. The database for that model is developed via extrapolation from current literature. Optimization is done via minimization of total mission cost using a linear programming approach. This prototype is the first cut at producing a general tool for choosing a near-optimum constructor mix for any space construction mission. It illuminates some significant representational and data-gathering problems with the modelling approach.

Johnson, Richard↗

Impact of ethanol blending into gasoline on aromatic compound evaporation and particle emissions from a gasoline direct injection engine

This study examined the interactions of ethanol with aromatic compounds on aromatic species evolution during gasoline evaporation and the consequent impacts on particle emissions from a single-cylinder, gasoline direct injection engine. From a chemical kinetic standpoint, the dilution of the aromatic species in gasoline by ethanol blending should reduce particle emissions because aromatics are the primary species that form particles during gasoline combustion. However, increased evaporative cooling from the ethanol and ethanol's effect on when aromatics in the fuel droplet (or surface film or pool) evaporate can drive an increase in particle emissions. The objective of this study was to quantify these competing effects. The results show that a combination of ethanol's increased evaporative cooling and impact on aromatic compound vapor-liquid equilibrium can extend droplet lifetime and cause aromatics to evaporate later in the droplet evaporation process than would be the case without ethanol. When fuels were burned in a direct-injection single-cylinder engine at high speed (2500 rpm) and relatively high load, ethanol blending caused an increase in particle emissions for fuels containing a low vapor pressure aromatic. At a lower speed (1500 rpm), no statistically significant increase in particle emissions was observed, likely because more time was available for evaporation and mixing to occur. In contrast, particle emissions from a fuel blended with a high vapor pressure aromatic were either insensitive to or reduced by ethanol blending, likewise to engine operating conditions. Current models for predicting particle emissions from gasoline engines, such as particulate matter index (PMI), generally lack variables for non-linear interactions between ethanol and gasoline hydrocarbons. To identify nonlinear interactions between predictor variables that might better explain the observed particle emissions, the least absolute shrinkage and selection operator (LASSO) regularized regression approach was employed. Among the predictor variables evaluated was the Yield Sooting Index (YSI) which is as an improved measure to quantify the aromatic's chemical tendency to generate particle emissions. Following this analysis, two variables were selected based on their frequent appearance. The first was ethanol mole percent multiplied by aromatic mole percent divided by aromatic vapor pressure [(EtOH * Aro%)/Aro VP@443 K]. The second was aromatic mole percent multiplied by aromatic YSI divided by aromatic molecular weight [(Aro% * Aro YSI)/Aro MW]. The results of linear regression including only these two combined explanatory variables achieved an r 2 = 0.959, substantially improved over the PMI model (r 2 = 0.688) and provides insight into how the PMI formalism might be modified to account for non-linear oxygenate-hydrocarbon coupling.

09 BIOMASS FUELS↗

Revising the dynamic energy budget theory with a new reserve mobilization rule and three example applications to bacterial growth

Dynamic energy budget (DEB) theory has been applied to model a wide range of organisms, including microbes. In the standard DEB model, biomass is partitioned into reserve and structural compartments, where reserve biomass is mobilized in a pseudolinear manner (while the reserve biomass density, defined as the ratio between reserve and structural biomass, decays linearly) to drive maintenance and the growth of structural biomass (and extracellular enzyme production if it is considered). However, the linear dynamics of the reserve biomass density makes the standard DEB model incapable of explaining the slowdown of microbial growth at high reserve density that is caused by macromolecular crowding effect which reduces biochemical reaction rates (a typical situation occurs when microbes are experiencing severe moisture stress) and is inconsistent with the observation that intracellular enzymatic reactions generally follow non-linear kinetics. By partitioning biomass into reserve, kinetic, and structural compartments, we show here that the Equilibrium Chemistry Approximation (ECA) kinetics can be used to represent enzymatically catalyzed reserve biomass mobilization that can then drive the kinetic and structural biomass synthesis. This revised DEB model better represents the tradeoff in ribosome allocation for structural growth and internal enzyme production, is structurally compatible with metabolic models of cell individuals, and includes the standard DEB model and the popular compromise model as special cases for representing population growth. We then applied the revised DEB model to interpret components of bacterial respiration, their dependence on substrate availability, and emergent microbial carbon use efficiency dynamics for an exponentially growing population. We found that the revised DEB model enables a better understanding of bacterial substrates use (carbon in our examples) than that can be derived from a few other models in the literature. In particular, the revised DEB model explains why carbon use efficiency may first increase, then plateau, and finally decrease with growth rate (and substrate uptake rate), as a function of proteomics. Additionally, the revised DEB model explains why the kinetic biomass compartment needs to be divided to reasonably incorporate proteomic control of microbial growth.

59 BASIC BIOLOGICAL SCIENCES↗

Integrated controls-structures design methodology development for a class of flexible spacecraft

Future utilization of space will require large space structures in low-Earth and geostationary orbits. Example missions include: Earth observation systems, personal communication systems, space science missions, space processing facilities, etc., requiring large antennas, platforms, and solar arrays. The dimensions of such structures will range from a few meters to possibly hundreds of meters. For reducing the cost of construction, launching, and operating (e.g., energy required for reboosting and control), it will be necessary to make the structure as light as possible. However, reducing structural mass tends to increase the flexibility which would make it more difficult to control with the specified precision in attitude and shape. Therefore, there is a need to develop a methodology for designing space structures which are optimal with respect to both structural design and control design. In the current spacecraft design practice, it is customary to first perform the structural design and then the controller design. However, the structural design and the control design problems are substantially coupled and must be considered concurrently in order to obtain a truly optimal spacecraft design. For example, let C denote the set of the 'control' design variables (e.g., controller gains), and L the set of the 'structural' design variables (e.g., member sizes). If a structural member thickness is changed, the dynamics would change which would then change the control law and the actuator mass. That would, in turn, change the structural model. Thus, the sets C and L depend on each other. Future space structures can be roughly divided into four mission classes. Class 1 missions include flexible spacecraft with no articulated appendages which require fine attitude pointing and vibration suppression (e.g., large space antennas). Class 2 missions consist of flexible spacecraft with articulated multiple payloads, where the requirement is to fine-point the spacecraft and each individual payload while suppressing the elastic motion. Class 3 missions include rapid slewing of spacecraft without appendages, while Class 4 missions include general nonlinear motion of a flexible spacecraft with articulated appendages and robot arms. Class 1 and 2 missions represent linear mathematical modeling and control system design problems (except for actuator and sensor nonlinearities), while Class 3 and 4 missions represent nonlinear problems. The development of an integrated controls/structures design approach for Class 1 missions is addressed. The performance for these missions is usually specified in terms of (1) root mean square (RMS) pointing errors at different locations on the structure, and (2) the rate of decay of the transient response. Both of these performance measures include the contributions of rigid as well as elastic motion.

Maghami, P. G.↗

COSMOLOGY OF CHAMELEONS WITH POWER-LAW COUPLINGS

In chameleon field theories, a scalar field can couple to matter with gravitational strength and still evade local gravity constraints due to a combination of self-interactions and the couplings to matter. Originally, these theories were proposed with a constant coupling to matter; however, the chameleon mechanism also extends to the case where the coupling becomes field dependent. We study the cosmology of chameleon models with power-law couplings and power-law potentials. It is found that these generalized chameleons, when viable, have a background expansion very close to {Lambda}CDM, but can in some special cases enhance the growth of the linear perturbations at low redshifts. For the models we consider, it is found that this region of the parameter space is ruled out by local gravity constraints. Imposing a coupling to dark matter only, the local constraints are avoided, and it is possible to have observable signatures on the linear matter perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of curvature-related DEM contact model on the macro- and micro-mechanical behaviours of granular soils

The presented paper conducted a comprehensive comparison between the Hertz-Mindlin model and the linear spring model in true triaxial shear simulations of granular soils using the discrete element method (DEM). The no-slip Hertz-Mindlin model for general elastic non-spherical particles with smooth surfaces was revisited, and implemented for superellipsoidal particles in an in-house DEM code. Three groups of specimens with a grain size distribution of Ottawa 20-30 sands, consisting of spheres, ellipsoids, and superellipsoids, respectively, were subjected to triaxial shear DEM simulations with the Hertz-Mindlin model and the linear spring model. The corresponding mechanical behaviours were examined in terms of a series of macro-and micro-parameters. It is found that the linear spring model is able to resemble the Hertz-Mindlin model on aspects of both microscopic and macroscopic mechanical behaviours of granular media with spherical and/or non-spherical particles. This finding suggests that the linear spring model can be used to investigate micro-mechanical behaviours of granular soils even with complex particle shapes.

42 ENGINEERING↗

EKAT v.1.0

E3SM Kokkos Application Toolkit (EKAT) is a collection of C++, Fortran, and CMake utilities for providing a single implementation of common kernels based on the Kokkos programming model. The library contains utilities for vectorization, tridiagonal linear system solvers, and linear interpolation as well as some general-purpose utilities such as testing utilities, parameter lists, representation of physical units, and additional interfaces. The goal is to provide a centralized implementation for high-performance computing structures and common utilities that reduce code duplication and streamline maintenance efforts. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525. SAND2022-1327 O

Bertagna, Luca↗

Numerical studies of identification in nonlinear distributed parameter systems

An abstract approximation framework and convergence theory for the identification of first and second order nonlinear distributed parameter systems developed previously by the authors and reported on in detail elsewhere are summarized and discussed. The theory is based upon results for systems whose dynamics can be described by monotone operators in Hilbert space and an abstract approximation theorem for the resulting nonlinear evolution system. The application of the theory together with numerical evidence demonstrating the feasibility of the general approach are discussed in the context of the identification of a first order quasi-linear parabolic model for one dimensional heat conduction/mass transport and the identification of a nonlinear dissipation mechanism (i.e., damping) in a second order one dimensional wave equation. Computational and implementational considerations, in particular, with regard to supercomputing, are addressed.

Banks, H. T.↗

Dynamic modeling, property investigation, and adaptive controller design of serial robotic manipulators modeled with structural compliance

Research results on general serial robotic manipulators modeled with structural compliances are presented. Two compliant manipulator modeling approaches, distributed and lumped parameter models, are used in this study. System dynamic equations for both compliant models are derived by using the first and second order influence coefficients. Also, the properties of compliant manipulator system dynamics are investigated. One of the properties, which is defined as inaccessibility of vibratory modes, is shown to display a distinct character associated with compliant manipulators. This property indicates the impact of robot geometry on the control of structural oscillations. Example studies are provided to illustrate the physical interpretation of inaccessibility of vibratory modes. Two types of controllers are designed for compliant manipulators modeled by either lumped or distributed parameter techniques. In order to maintain the generality of the results, neither linearization is introduced. Example simulations are given to demonstrate the controller performance. The second type controller is also built for general serial robot arms and is adaptive in nature which can estimate uncertain payload parameters on-line and simultaneously maintain trajectory tracking properties. The relation between manipulator motion tracking capability and convergence of parameter estimation properties is discussed through example case studies. The effect of control input update delays on adaptive controller performance is also studied.

Tesar, Delbert↗

Transonic Aerodynamic Characteristics of a 45 deg Swept Wing Fuselage Model with a Finned and Unfinned Body Pylon Mounted Beneath the Fuselage or Wing, Including Measurements of Body Loads

An investigation of a model of a standard size body in combination with a representative 45 deg swept-wing-fuselage model has been conducted in the Langley 8-foot transonic pressure tunnel over a Mach number range from 0.80 to 1.43. The body, with a fineness ratio of 8.5, was tested with and without fins, and was pylon-mounted beneath the fuselage or wing. Force measurements were obtained on the wing-fuselage model with and without the body, for an angle-of-attack range from -2 deg to approximately 12 deg and an angle-of-sideslip range from -8 deg to 8 deg. In addition, body loads were measured over the same angle-of-attack and angle-of-sideslip range. The Reynolds number for the investigation, based on the wing mean aerodynamic chord, varied from 1.85 x 10(exp 6) to 2.85 x 10(exp 6). The addition of the body beneath the fuselage or the wing increased the drag coefficient of the complete model over the Mach number range tested. On the basis of the drag increase per body, the under-fuselage position was the more favorable. Furthermore, the bodies tended to increase the lateral stability of the complete model. The variation of body loads with angle of attack for the unfinned bodies was generally small and linear over the Mach number range tested with the addition of fins causing large increases in the rates of change of normal-force coefficient and nose-down pitching-moment coefficient. The variation of body side-force coefficient with sideslip for the unfinned body beneath the fuselage was at least twice as large as the variation of this load for the unfinned body beneath the wing. The addition of fins to the body beneath either the fuselage or the wing approximately doubled the rate of change of body side-force coefficient with sideslip. Furthermore, the variation of body side-force coefficient with sideslip for the body beneath the wing was at least twice as large as the variation of this load with angle of attack.

Wornom, Dewey E.↗

Cyber-Physical Correlations for Infrastructure Resilience: A Game-Theoretic Approach

In several critical infrastructures, the cyber and physical parts are correlated so that disruptions to one affect the other and hence the whole system. These correlations may be exploited to strategically launch components attacks, and hence must be accounted for ensuring the infrastructure resilience, specified by its survival probability. We characterize the cyber-physical interactions at two levels: (i) the failure correlation function specifies the conditional survival probability of cyber sub-infrastructure given the physical sub-infrastructure as a function of their marginal probabilities, and (ii) the individual survival probabilities of both sub-infrastructures are characterized by first-order differential conditions. We formulate a resilience problem for infrastructures composed of discrete components as a game between the provider and attacker, wherein their utility functions consist of an infrastructure survival probability term and a cost term expressed in terms of the number of components attacked and reinforced. We derive Nash Equilibrium conditions and sensitivity functions that highlight the dependence of infrastructure resilience on the cost term, correlation function and sub-infrastructure survival probabilities. These results generalize earlier ones based on linear failure correlation functions and independent component failures. We apply the results to models of cloud computing infrastructures and energy grids.

cyber-physical infrastructures↗

Classification and stability of plasma motion in periodic linearly polarized relativistic waves

Based on a relativistic fluid-Maxwell model, laser-induced plasma dynamics is investigated for relativistic periodic waves. Within a one-dimensional (1D) description, the Akhiezer-Polovin model is applied to the existence of periodic, nonlinearly coupled electromagnetic and electrostatic waves, and the corresponding particle motion. Known existence criteria for periodic solutions are generalized. The corresponding stability behaviors are investigated by 1D integrators of the relativistic fluid-Maxwell model. It is shown that in contrast to the vacuum solution, linearly polarized coupled electromagnetic-electrostatic waves are unstable in plasmas. The magnitudes of the growth rates are investigated in terms of the maximum amplitudes and normalized phase velocities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic models for simulation of the 70-M antenna axis servos

Dynamic models for the various functional elements of the 70 m antenna axis servos are described. The model representing the digital position controller, the linear and nonlinear properties of the physical hardware, and the dynamics of the flexible antenna structure are encoded in six major function blocks. The general modular structure of the function blocks facilitates their adaptation to a variety of dynamic simulation studies. Model parameter values were calculated from component specifications and design data. A simulation using the models to predict limit cycle behavior produced results in excellent agreement with field test data from the DSS 14 70-m antenna.

Hill, R. E.↗

Enhanced rotor modeling tailored for rub dynamic stability analysis and simulation

New methods are presented that allow straightforward application of complex nonlinearities to finite element based rotor dynamic analyses. The key features are: (1) the methods can be implemented with existing finite element or dynamic simulation programs, (2) formulation is general for simple application to a wide range of problems, and (3) implementation is simplified because nonlinear aspects are separated from the linear part of the model. The new techniques are illustrated with examples of inertial nonlinearity and torquewhirl which can be important in rubbing turbomachinery. The sample analyses provide new understanding of these nonlinearities which are discussed.

Davis, R. R.↗

OPTICON: Pro-Matlab software for large order controlled structure design

A software package for large order controlled structure design is described and demonstrated. The primary program, called OPTICAN, uses both Pro-Matlab M-file routines and selected compiled FORTRAN routines linked into the Pro-Matlab structure. The program accepts structural model information in the form of state-space matrices and performs three basic design functions on the model: (1) open loop analyses; (2) closed loop reduced order controller synthesis; and (3) closed loop stability and performance assessment. The current controller synthesis methods which were implemented in this software are based on the Generalized Linear Quadratic Gaussian theory of Bernstein. In particular, a reduced order Optimal Projection synthesis algorithm based on a homotopy solution method was successfully applied to an experimental truss structure using a 58-state dynamic model. These results are presented and discussed. Current plans to expand the practical size of the design model to several hundred states and the intention to interface Pro-Matlab to a supercomputing environment are discussed.

Peterson, Lee D.↗

A Communication Channel Density Estimating Generative Adversarial Network

Autoencoder-based communication systems use neural network channel models to backwardly propagate message reconstruction error gradients across an approximation of the physical communication channel. In this work, we develop and test a new generative adversarial network (GAN) architecture for the purpose of training a stochastic channel approximating neural network. In previous research, investigators have focused on additive white Gaussian noise (AWGN) channels and/or simplified Rayleigh fading channels, both of which are linear and have well defined analytic solutions. Given that training a neural network is computationally expensive, channel approximation networks— and more generally the autoencoder systems—should be evaluated in communication environments that are traditionally difficult. To that end, our investigation focuses on channels that contain a combination of non-linear amplifier distortion, pulse shape filtering, intersymbol interference, frequency-dependent group delay, multipath, and non-Gaussian statistics. Each of our models are trained without any prior knowledge of the channel. We show that the trained models have learned to generalize over an arbitrary amplifier drive level and constellation alphabet. We demonstrate the versatility of our GAN architecture by comparing the marginal probability density function of several channel simulations with that of their corresponding neural network approximations

Smith, Aaron↗

First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon

Event Horizon Telescope (EHT) observations at 230 GHz have now imaged polarized emission around the supermassive black hole in M87 on event-horizon scales. This polarized synchrotron radiation probes the structure of magnetic fields and the plasma properties near the black hole. Here we compare the resolved polarization structure observed by the EHT, along with simultaneous unresolved observations with the Atacama Large Millimeter/submillimeter Array, to expectations from theoretical models. The low fractional linear polarization in the resolved image suggests that the polarization is scrambled on scales smaller than the EHT beam, which we attribute to Faraday rotation internal to the emission region. We estimate the average density n_e of order 10^4-7 cm-3, magnetic field strength B of order 1-30 G, and electron temperature Te of order (1-12) x 10^10 K of the radiating plasma in a simple one-zone emission model. We show that the net azimuthal linear polarization pattern may result from organized, poloidal magnetic fields in the emission region. In a quantitative comparison with a large library of simulated polarimetric images from general relativistic magnetohydrodynamic (GRMHD) simulations, we identify a subset of physical models that can explain critical features of the polarimetric EHT observations while producing a relativistic jet of sufficient power. The consistent GRMHD models are all of magnetically arrested accretion disks, where near-horizon magnetic fields are dynamically important. We use the models to infer a mass accretion rate onto the black hole in M87 of (3-20) x 10^-4 Msun yr-1.

79 ASTRONOMY AND ASTROPHYSICS↗