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 253 records · Page 14

Spherically symmetric cosmological spacetimes with dust and radiation — numerical implementation

We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.

79 ASTRONOMY AND ASTROPHYSICS↗

Active Flutter Suppression Controllers Derived from Linear and Nonlinear Aerodynamics: Application to a Transport Aircraft Model

Active flutter suppression has been demonstrated in simulation by many researchers, generally using methods based on linear aerodynamics and often with simplistic geometries. In this paper, active flutter suppression is demonstrated in a simulation using a Navier-Stokes aerodynamics code, FUN3D (Fully Unstructured Navier-Stokes Three-Dimensional), and a realistic transport aircraft configuration. This is accomplished using simple observer-feedback controllers derived from linear aeroelastic models, including reduced order models built via FUN3D data. The development of these reduced order models is described here. It is shown that controllers derived from reduced order models of the nonlinear aerodynamics outperform controllers based on linear aerodynamics.

Waite, Josiah M.↗

Generating parity relations for detecting and identifying control system component failures

The monitoring of control system sensors and actuators for failures is presently undertaken by means of an exceptionally simple form of generalized parity relations based on a discrete-time model of the dynamics of linear, time-invariant systems. These generalized parity relations are constructed by recourse to a transfer matrix-description of the system that is additionally useful in the interpretation of their properties. Attention is given to a novel method for constructing the parity relation of minimum length that depends on the output of only a single sensor.

Vander Velder, Wallace E.↗

Toward the large-eddy simulations of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model (in terms of Favre-filtered fields) is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The three dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 to the third power grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. Future applications of these compressible subgrid-scale models to the large-eddy simulation of supersonic aerodynamic flows are discussed briefly.

Erlebacher, G.↗

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96 (exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine-grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly.

Erlebacher, G.↗

Monte Carlo goodness-of-fit tests for degree corrected and related stochastic blockmodels

Abstract We construct Bayesian and frequentist finite-sample goodness-of-fit tests for three different variants of the stochastic blockmodel for network data. Since all of the stochastic blockmodel variants are log-linear in form when block assignments are known, the tests for the latent block model versions combine a block membership estimator with the algebraic statistics machinery for testing goodness-of-fit in log-linear models. We describe Markov bases and marginal polytopes of the variants of the stochastic blockmodel and discuss how both facilitate the development of goodness-of-fit tests and understanding of model behaviour. The general testing methodology developed here extends to any finite mixture of log-linear models on discrete data, and as such is the first application of the algebraic statistics machinery for latent-variable models.

Karwa, Vishesh↗

Storage Viability and Optimization Web Service

Non-residential sectors offer many promising applications for electrical storage (batteries) and photovoltaics (PVs). However, choosing and operating storage under complex tariff structures poses a daunting technical and economic problem that may discourage potential customers and result in lost carbon and economic savings. Equipment vendors are unlikely to provide adequate environmental analysis or unbiased economic results to potential clients, and are even less likely to completely describe the robustness of choices in the face of changing fuel prices and tariffs. Given these considerations, researchers at Lawrence Berkeley National Laboratory (LBNL) have designed the Storage Viability and Optimization Web Service (SVOW): a tool that helps building owners, operators and managers to decide if storage technologies and PVs merit deeper analysis. SVOW is an open access, web-based energy storage and PV analysis calculator, accessible by secure remote login. Upon first login, the user sees an overview of the parameters: load profile, tariff, technologies, and solar radiation location. Each parameter has a pull-down list of possible predefined inputs and users may upload their own as necessary. Since the non-residential sectors encompass a broad range of facilities with fundamentally different characteristics, the tool starts by asking the users to select a load profile from a limited cohort group of example facilities. The example facilities are categorized according to their North American Industry Classification System (NAICS) code. After the load profile selection, users select a predefined tariff or use the widget to create their own. The technologies and solar radiation menus operate in a similar fashion. After these four parameters have been inputted, the users have to select an optimization setting as well as an optimization objective. The analytic engine of SVOW is LBNL?s Distributed Energy Resources Customer Adoption Model (DER-CAM), which is a mixed-integer linear program (MILP) written and executed in the General Algebraic Modeling System (GAMS) optimization software. LBNL has released version 1.2.0.11 of SVOW. Information can be found at http://der.lbl.gov/microgrids-lbnl/current-project-storage-viability-website.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Magnetic properties of metastable honeycomb KCoAsO 4

We present comprehensive neutron scattering data on polycrystalline samples of a new metastable honeycomb material KCoAsO 4 . Below 𝑇 𝑁 = 14 K, the system orders into a zigzag antiferromagnetic state, with spins ordered into alternating ferromagnetic chains similar to the isostructural sister compound KNiAsO 4 . In the case of KCoAsO 4 , we find the moments are robustly canted out of plane closer to the crystallographic 𝑐 axis. A combination of weak interlayer coupling, lattice strain, and inhomogeneities lead to the coexistence of two magnetic ordering wave vectors 𝑘 1 = (1.5 0 0) and 𝑘 2 = (0.5 0 0.5), where the two structures differ only in their layer stacking. Inelastic data show the presence of a spin orbital mode at 24 meV, supporting a pseudospin $\tilde{𝑆}$ = 1/2 Kramer's doublet ground state of the Co 2+ ions. We model the low energy excitations using both a conventional XXZ Hamiltonian and generalized Kitaev Heisenberg Hamiltonian within a linear spin wave limit. While either model can qualitatively reproduce the observed spectra, the lack of fine features and observation of disorder prevent a clear-cut determination of the low energy Hamiltonian. In the case of an XXZ-type model, a large easy-axis anisotropy is necessary to reproduce the gapped spectra and canting of the magnetic moments. For the generalized Kitaev model, despite the large canting of the moments away from the honeycomb layers, we find a noticeable if nondominant Kitaev term persists. In conclusion, the contrast of the magnetic properties of KCoAsO 4 to other cobalt honeycombs highlights the sensitivity of the low energy magnetic properties of Co 2+ to fine details of its crystalline environment.

Honeycomb lattice↗

A comprehensive conductivity model for drift and micro-tearing modes

The parallel electrical conductivity is a crucial parameter in the study of the linear stability of drift-modes like the micro-tearing mode (MTM). The conductivity enters by closing the electromagnetic tearing layer equations. Recent progress in the understanding of the pedestal suggests that the MTM could play an important role in its structure and evolution. For this reason, we revisit and improve previous model conductivities. This parameter is generally derived from the linearized drift kinetic equation. In the past literature, it has been computed using either simplified collision operators or neglecting the spatial dependence away from the rational surface. Here, a fully consistent expression for the conductivity that would accurately model the pedestal has not been available. By applying a novel variational procedure and with the full Fokker Plank collision operator, including electron–electron collisions, we compute a closed expression for the parallel electrical conductivity in the form of a rational function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simultaneous Evolution of Gyre and Atlantic Meridional Overturning Circulation Anomalies as an Eigenmode of the North Atlantic System

Here, the authors observe an interdecadal oscillatory mode of the North Atlantic atmosphere–ocean system in a general circulation model (GFDL CM2.1) via a linear inverse model (LIM). The oscillation mechanism is mostly embedded in the subpolar gyre: anomalous advection generates density anomalies in the eastern subpolar gyre, which propagate along the mean gyre circulation and reach the subpolar gyre center around 10 years later, when associated anomalous advection of the opposite sign starts the other half cycle. While density anomalies reach the Labrador Sea deep convection region, Atlantic meridional overturning circulation (AMOC) anomalies are also induced. Both the gyre and AMOC anomalies then propagate equatorward slowly, following the advection of density anomalies. The oscillation is further demonstrated to be more likely an ocean-only mode while excited by the atmospheric forcing; in particular, it can be approximated as a linearly driven damped oscillator that is partly excited by the North Atlantic Oscillation (NAO). The slowly evolving interdecadal oscillation significantly improves and prolongs the LIM’s prediction skill of sea surface temperature (SST) evolution.

54 ENVIRONMENTAL SCIENCES↗

Discovering the desirable alleles contributing to the lignocellulosic biomass traits in saccharum germplasm collections for energy cane improvement

Phenotyping Methods: The accessions (which includes 21 taxa and 1,177 accessions) in the World Collection of Sugarcane and Related Grasses (WCSRG) was evaluated for the following traits: arenchyma, internode length and diameter, pubescence, pith, Brix, stalk number and fiber. A core of 300 accessions that included each species in the World Collection was selected by using the Maximization Strategy in MStrat software. Results: The core had a higher diversity rating than random selections of 300 accessions. The Shannon–Weaver Diversity Index scores of the core and whole collection were similar indicating that the majority of the diversity was captured by the core collection. The ranges and medians between the core and WCSRG were similar; only two of the trait medians were not significant at P = 0.05 using the non-parametric Wilcoxon method and the coincidence rate (CR % = 96.2) was high (>80) indicating that extreme values were retained. Thus, the phenotypic diversity of these traits in the WCSRG was well represented by the core collection. Associations Methods: Genotypic and phenotypic data were collected for 1002 accessions of the WCSRG including 209 SSR markers. Association analysis was performed using both General Linear (GLM) and Maximum Likelihood (MLM) models. Different core collections with 300 accessions each were selected based on genotypic, phenotypic and combined data based on the Maximization Strategy in MStrat software. Results: A major portion of the genotyping involving SNPs is being conducted by Dr. Jianping Wang of the University of Florida under the DOE award DE-FG 02-11ER 65214 and the genotypic and phenotypic associations will be reported separately next year. In the current, study forty one and seventeen markers were found to be associated with traits using the GLM and MLM methods respectively including associations with arenchyma, internode length and diameter, pubescence, pith, and Sugar Cane Yellow Leaf Virus. The data indicates that each of the cores and the World Collection are similar to each other genotypically and phenotypically, but the core that was selected using only genotypic data was significantly different phenotypically. This indicates that there is not enough association between the genotypic and phenotypic diversity as to select using only genotypic diversity and get the full phenotypic diversity. Core Collection: Creation and Phenotyping Methods: To evaluate this germplasm for breeding purposes, a representative diversity panel selected from the WCSRG of approximately 300 accessions was planted at Canal Point, FL in three replications. These accessions were measured for stalk height and stalk number multiple times throughout the growing season and Brix and fresh biomass during harvest in 2013 and, stalk height, stalk number, stalk diameter, internode length, Brix and fresh and dry biomass was determined in the ratoon crop harvest in 2014. Results: In correlations of multiple measurements, there were higher correlations for early measurements of stalk number and stalk height with harvest traits like Brix and fresh weight. Hybrids had higher fresh mass and Brix while Saccharum spontaneum had higher stalk number and dry mass. The heritability of hybrid mass traits was lower in the ratoon crop. According to the principal component analysis, the diversity panel was divided into two groups. One group had accessions with high stalk number and high dry biomass like S. spontaneum and the other groups contained accessions with higher Brix and fresh biomass like S. officinarum. Mass traits correlated with each other as expected but hybrids had lower correlations between fresh and dry mass. Stalk number and the mass traits correlated with each other except in S. spontaneum and hybrids in the first ratoon. There were 110 accessions not significantly different in Brix from the commercial sugarcane checks including 10 S. spontaneum accessions. There were 27 dry and 6 fresh mass accessions significantly higher than the commercial sugarcane checks. Core Collection: Fiber analysis Methods: A biomass sample was taken from each accession then shredded and dried. Fiber analysis was then performed on each sample. The acetyl groups, acid insoluble lignin, acid soluble lignin, arabinan, glucan, holocellulose, total lignin, structural ash, and xylan were quantified on a % fiber basis and nonstructural ash on a % total basis. Results: There were significant, but not large differences between species for holocellulose, lignin, acetyl, acid soluble lignin, nonstructural ash, and glucan. For each trait, S. spontaneum had significantly more holocellulose, glucan, lignin, and nonstructural ash and less acetyl and acid soluble lignin than the other species. In all populations, glucan and was positively correlated with holocellulose were positively correlated and glucan and and holocellulose were negatively correlated with lignin. In hybrids, internode length correlated positively with holocellulose and nonstructural ash and negatively with lignin. The heritability estimates for each of the fiber component traits is low indicating that environment is an important factor in fiber composition. Principal component analysis indicated that a large amount of diversity exists within each of the species.

09 BIOMASS FUELS↗

Techniques in Doppler gravity inversion

The types of Doppler gravity data available for local as opposed to planetwide geophysical modeling are reviewed. Those gravity fields that are determined dynamically in orbit determination programs yield a smoothed representation of the local gravity field that may be used for quantitative modeling. An estimate of the difference between smoothed and true fields can be considered as a noise limitation in generating local gravity models. A nonlinear inversion for the geometry, depth, and density of the Mare Serenitatis mascon using an ellipsoidal model yielded a global least squares minimum in horizontal dimensions, depth, and thickness-density contrast product. It was subsequently found, by using a linear model, that there were an infinite number of solutions corresponding to various combinations of depth and lateral inhomogeneity. Linear modeling was performed by means of generalized inverse theory.

Phillips, R. J.↗

Prediction of blade vortex interaction noise from measured blade pressure

The impulsive nature of noise due to the interaction of a rotor blade with a tip vortex is studied. The time signature of this noise is calculated theoretically based on the measured blade surface pressure fluctuation of an operational load survey rotor in slow descending flight and is compared with the simultaneous microphone measurement. Particularly, the physical understanding of the characteristic features of a waveform is extensively studied in order to understand the generating mechanism and to identify the important parameters. The interaction trajectory of a tip vortex on an acoustic planform is shown to be a very important parameter for the impulsive shape of the noise. The unsteady nature of the pressure distribution at the very leading edge is also important to the pulse shape. The theoretical model using noncompact linear acoustics predicts the general shape of interaction impulse pretty well except for peak amplitude which requires more continuous pressure information along the span at the leading edge.

Nakamura, Y.↗

Theoretical Modeling for Numerical Weather Prediction

The goal is to utilize predictability theory and numerical experimentation to identify and understand some of the dynamical processes which must be modeled more realistically if large-scale numerical weather predictions are to be improved. The emphasis is on the use of relatively simple models to exlore the properties of physically comprehensive general circulation models (GCM's). A global linear quasi-geostrophic model and the Goddard Laboratory for Atmospheric Sciences (GLAS) GCM were used to investigate several mechanisms which are responsible for the decay of large-scale forecast skill in mid-latitude numerical weather predictions. Five-day forecasts for an ensemble of cases were made using First GARP Global Experiment data. It was found that forecast skill depends crucially on the specification of the stationary forcing. A lack of stationary forcing leads to spurious westwad propagation of the ultralong waves. Forecasts made with stationary forcings derived from climatological data are superior to those using forcings inferred from observations immediately preceding the forecast period. Interhemispheric forecast differences were analyzed, and the model errors were compared to errors of a simple persistence-damped-to-climatology scheme and to errors of the GLAS GCM.

Somerville, R. C. J.↗

The design and performance of a real-time self excited vocoder

This paper is concerned with a generic class of predictive speech coders that includes the newly proposed Self Excited Vocoder (SEV) and the well known Code-Excited Linear Predictive Coder (CELPC). All members of this class form an excitation sequence for a linear predictive model filter using the same general model for the excitation signal. The general excitation model is based on a block coding technique where each sequence is drawn from an ensemble of sequences. This paper reports on two developments related to this general model. The first development is a new type of excitation ensemble that can in general be populated by many different types of sequences. The second development is a means of populating this new type of ensemble based on a vector quantizer design procedure using a new distortion measure.

Rose, Richard C.↗

Linear analysis of a force reflective teleoperator

Complex force reflective teleoperation systems are often very difficult to analyze due to the large number of components and control loops involved. One mode of a force reflective teleoperator is described. An analysis of the performance of the system based on a linear analysis of the general full order model is presented. Reduced order models are derived and correlated with the full order models. Basic effects of force feedback and position feedback are examined and the effects of time delays between the master and slave are studied. The results show that with symmetrical position-position control of teleoperators, a basic trade off must be made between the intersystem stiffness of the teleoperator, and the impedance felt by the operator in free space.

Biggers, Klaus B.↗

Aspects of the Kelvin wave response to episodic wind forcing

Episodes of westerly wind are an important aspect of surface stress variability in the western Pacific. During ENSO periods, the presence of such wind episodes comprises much of the LF relaxation of the trades over the central and western Pacific. This paper describes the oceanic Kelvin pulse response to a single idealized episode of westerly wind stress, using results from linear theory as well as from a 27-level general circulation model. When stratification typical of the western and eastern Pacific is used, the conservation of energy flux predicts a reduction of surface currents associated with the first baroclinic mode and an enhancement of surface currents associated with the second baroclinic mode. The idealized wind anomaly is also used to drive an ocean general circulation model. When the wind anomaly is weak, the model Kelvin response agrees with predictions of linear theory. For more realistic strong forcing there are three important deviations from linear theory: the amplitude of low baroclinic modes increases; the amplitude of higher baroclinic modes decreases; and the phase speed increases.

Giese, Benjamin S.↗