Search NASA⌕ Search

SEARCH · Search NASA

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 469 records · Page 26

GVE-Based Dynamics and Control for Formation Flying Spacecraft

Formation flying is an enabling technology for many future space missions. This paper presents extensions to the equations of relative motion expressed in Keplerian orbital elements, including new initialization techniques for general formation configurations. A new linear time-varying form of the equations of relative motion is developed from Gauss Variational Equations and used in a model predictive controller. The linearizing assumptions for these equations are shown to be consistent with typical formation flying scenarios. Several linear, convex initialization techniques are presented, as well as a general, decentralized method for coordinating a tetrahedral formation using differential orbital elements. Control methods are validated using a commercial numerical propagator.

Breger, Louis↗

Properties of finite difference models of non-linear conservative oscillators

Finite-difference (FD) approaches to the numerical solution of the differential equations describing the motion of a nonlinear conservative oscillator are investigated analytically. A generalized formulation of the Duffing and modified Duffing equations is derived and analyzed using several FD techniques, and it is concluded that, although it is always possible to contstruct FD models of conservative oscillators which are themselves conservative, caution is required to avoid numerical solutions which do not accurately reflect the properties of the original equation.

Mickens, R. E.↗

On optimal modeling of systems.

A procedure for modeling a linear system by an optimal finite-dimensional approximation is developed on the basis of Sard's (1967) generalized spline. Error bounds are given, and the application of the procedure is illustrated by two presented examples. The procedure may be employed, not only in the modeling of conventional lumped-parameter and distributed-parameter continuous dynamical systems, but also in mathematical operations such as those involved in pattern recognition and picture enhancement problems.

De Figueiredo, R. J. P.↗

Attribution of Snowpack Errors to Simulated Temperature and Precipitation in E3SMv1 Over the Contiguous United States

Abstract Snow water equivalent (SWE), temperature, and precipitation biases and trends are evaluated in the atmosphere‐land simulations of the Energy Exascale Earth System Model version 1 (E3SMv1) in comparison to the Community Earth System Model version 2 (CESM2) and two other models using the ground measurement‐based University of Arizona (UA) snow product. SWE, temperature, and precipitation biases are highest in magnitude in the Western contiguous United States (CONUS). SWE errors are attributed to temperature and precipitation through multiple linear regressions of normalized errors, the coefficients of which represent the sensitivities to temperature and precipitation errors. SWE errors are more sensitive to temperature errors throughout the CONUS. Model SWE and temperature trends are generally opposite from UA product trends in the Western CONUS. SWE trend errors are also attributed to temperature and precipitation trend errors using multiple linear regressions of normalized trend errors. SWE trend errors are more sensitive to those of precipitation at higher elevations (>1,500 m) in the Western CONUS in these simulations. Thus, the sensitivity to temperature and precipitation differ for SWE errors and its trend errors. Furthermore, the SWE trend errors are more sensitive to temperature and precipitation in the atmosphere‐ocean coupled simulations in which the atmosphere‐land is coupled to an active ocean model. These results suggest that both errors in simulated temperature and precipitation contribute to SWE errors.

Meteorology & Atmospheric Sciences↗

An Adaptive Sequential Design for Model Discrimination and Parameter Estimation in Non-Linear Nested Models

The DKL-optimality criterion has been recently proposed for the dual problem of model discrimination and parameter estimation, for the case of two rival models. A sequential version of the DKL-optimality criterion is herein proposed in order to discriminate and efficiently estimate more than two nested non-linear models. Our sequential method is inspired by the procedure of Biswas and Chaudhuri (2002), which is however useful only in the set up of nested linear models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global impact of the Antarctic ozone hole: Simulations with a 3-D chemical transport model

A study of the Antarctic ozone hole was made with a 3-D chemical transport model using linearized photochemistry for ozone based on observed distribution. The tracer model uses the winds and convection from the GISS general circulation model (8 deg x 10 deg x 23 layers). A 3-year control run of the ozone distribution is compared with the observed climatology. In two experiments, a hypothetical Antarctic ozone hole is induced on October 1 and on November 1; the tracer model is integrated for 1 year with the standard linearized chemistry. The initial depletion, 90 percent of the O sub 3 poleward of 70 S between 25 and 180 mbar, amounts to about 5 percent of the total O sub 3 in the Southerm Hemisphere. As the vortex breaks down and the hole is dispersed, significant depletions to column ozone, of order 10 D.U., occur as far north as 36 S during Austral summer. One year later, about 25 percent of the original depletion remains, mostly below 100 mbar and poleward of 30 S. Details of the calculations are shown, along with a budget analysis showing the fraction of the hole filled in by photochemistry versus that transported into the troposhere.

Prather, Michael J.↗

Transfer of contrast sensitivity in linear visual networks

Contrast sensitivity is a useful measure of the ability of an observer to distinguish contrast signals from noise. Although usually applied to human observers, contrast sensitivity can also be defined operationally for individual visual neurons. In a model linear neuron consisting of a filter and noise source, this operational measure is a function of filter gain, noise power spectrum, signal duration, and a performance criterion. This definition allows one to relate the sensitivities of linear neurons at different levels in the visual pathway. Mathematical formulas describing these relationships are derived, and the general model is applied to the specific problem of relating the sensitivities of parvocellular LGN neurons and cortical simple cells in the primate.

Watson, Andrew B.↗

Soil Respiration Variability and Correlation Across a Wide Range of Temporal Scales

Abstract The high temporal variability of the soil‐to‐atmosphere CO 2 flux (soil respiration, R S ) has been studied at hourly to multiannual time scales but remains less well understood than R S spatial variability. How R S fluxes vary and are autocorrelated at various time lags has practical implications for sampling and more fundamentally for our understanding of its abiotic and biotic underlying mechanisms. We examined the variability, correlation, and sampling requirements of R S over a wide range of temporal scales in a temperate deciduous forest in eastern Maryland, USA, using both automated (temporally continuous, N = 30,036 over 10 months) and survey (spatially diverse, temporally sparse, N = 1,912 over 17 months) data. Data from a global R S database were also used to examine interannual variability in comparable forests. The coefficient of variability of successive measurements generally varied from the minute (median coefficient of variation 16%) to hourly and daily (11–12%) time scales. Successive R S values measured at a given collar exhibited a strong hour‐to‐hour correlation ( r = 0.931) and a moderate correlation at a 2‐hr lag (0.289); day‐to‐day (i.e., 24 hr lag) hourly observations were uncorrelated. Daily R S means were well correlated at a 1‐day lag ( r = 0.856) but not at any further time lag. In a linear mixed‐effects model predicting R S , soil temperature and moisture exerted consistently strong effects regardless of time scale, and model coefficient of variability was generally high (>80%). These results provide new opportunities to explore the drivers and variability of R S fluxes, quantify sampling requirements, and improve error propagation.

Bond‐Lamberty, Ben↗

Analytical and numerical modeling of optical second harmonic generation in anisotropic crystals using #SHAARP package

Electric-dipole optical second harmonic generation (SHG) is a second-order nonlinear process that is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the SHG polarimetry of a nonlinear optical material is essential to extract its point group symmetry and the absolute nonlinear susceptibilities. Current literature on SHG analysis involves numerous approximations and a wide range of (in)accuracies. We have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (#SHAARP.si) which derives analytical and numerical solutions of reflection SHG polarimetry from a single interface (.si) for bulk homogeneous crystals with arbitrary symmetry group, arbitrary crystal orientation, complex and anisotropic linear dielectric tensor with frequency dispersion, a general SHG tensor and arbitrary light polarization. #SHAARP.si enables accurate modeling of polarimetry measurements in reflection geometry from highly absorbing crystals or wedge-shaped transparent crystals. The package is extendable to multiple interfaces.

36 MATERIALS SCIENCE↗

Observational signatures of anisotropic inflationary models

We study observational signatures of two classes of anisotropic inflationary models in which an inflaton field couples to (i) a vector kinetic term F{sub μν}F{sup μν} and (ii) a two-form kinetic term H{sub μνλ}H{sup μνλ}. We compute the corrections from the anisotropic sources to the power spectrum of gravitational waves as well as the two-point cross correlation between scalar and tensor perturbations. The signs of the anisotropic parameter g{sub *} are different depending on the vector and the two-form models, but the statistical anisotropies generally lead to a suppressed tensor-to-scalar ratio r and a smaller scalar spectral index n{sub s} in both models. In the light of the recent Planck bounds of n{sub s} and r, we place observational constraints on several different inflaton potentials such as those in chaotic and natural inflation in the presence of anisotropic interactions. In the two-form model we also find that there is no cross correlation between scalar and tensor perturbations, while in the vector model the cross correlation does not vanish. The non-linear estimator f{sub NL} of scalar non-Gaussianities in the two-form model is generally smaller than that in the vector model for the same orders of |g{sub *}|, so that the former is easier to be compatible with observational bounds of non-Gaussianities than the latter.

79 ASTRONOMY AND ASTROPHYSICS↗

Applicability/evaluation of flux based representations for linear/higher order elements for heat transfer in structures - Generalized gamma(T)-family

The applicability and evaluation of a generalized gamma(T) family of flux-based representations are examined for two different thermal analysis formulations for structures and materials which exhibit no phase change effects. The so-called H-theta and theta forms are demonstrated for numerous test models and linear and higher-order elements. The results show that the theta form with flux-based representations is generally superior to traditional approaches.

Namburu, R. R.↗

Machine-learning based prediction of injection rate and solenoid voltage characteristics in GDI injectors

We report that current state-of-the-art gasoline direct-injection (GDI) engines use multiple injections as one of the key technologies to improve exhaust emissions and fuel efficiency. For this technology to be successful, secured adequate control of fuel quantity for each injection is mandatory. However, nonlinearity and variations in the injection quantity can deteriorate the accuracy of fuel control, especially with small fuel injections. Therefore, it is necessary to understand the complex injection behavior and to develop a predictive model to be utilized in the development process. This study presents a methodology for rate of injection (ROI) and solenoid voltage modeling using artificial neural networks (ANNs) constructed from a set of Zeuch-style hydraulic experimental measurements conducted over a wide range of conditions. A quantitative comparison between the ANN model and the experimental data shows that the model is capable of predicting not only general features of the ROI trend, but also transient and non-linear behaviors at particular conditions. In addition, the end of injection (EOI) could be detected precisely with a virtually generated solenoid voltage signal and the signal processing method, which applies to an actual engine control unit. A correlation between the detected EOI timings calculated from the modeled signal and the measurement results showed a high coefficient of determination.

42 ENGINEERING↗

Towards enhancing and delaying disturbances in free shear flows

The family of shear flows comprising the jet, wake, and the mixing layer are subjected to perturbations in an inviscid incompressible fluid. By modeling the basic mean flows as parallel with piecewise linear variations for the velocities, complete and general solutions to the linearized equations of motion can be obtained in closed form as functions of all space variables and time when posed as an initial value problem. The results show that there is a continuous as well as the discrete spectrum that is more familiar in stability theory and therefore there can be both algebraic and exponential growth of disturbances in time. These bases make it feasible to consider control of such flows. To this end, the possibility of enhancing the disturbances in the mixing layer and delaying the onset in the jet and wake is investigated. It is found that growth of perturbations can be delayed to a considerable degree for the jet and the wake but, by comparison, cannot be enhanced in the mixing layer. By using moving coordinates, a method for demonstrating the predominant early and long time behavior of disturbances in these flows is given for continuous velocity profiles. It is shown that the early time transients are always algebraic whereas the asymptotic limit is that of an exponential normal mode. Numerical treatment of the new governing equations confirm the conclusions reached by use of the piecewise linear basic models. Although not pursued here, feedback mechanisms designed for control of the flow could be devised using the results of this work.

Criminale, W. O.↗

Elevation distributed micro-climatology data in a coastal glaciated watershed

Air temperature, ground temperature and relative humidity data were collected in a longitudinal transect of the Nooksack watershed at varying elevations from 500 to 1800 m above sea level. Data were collected by anchoring sensors from trees above winter snow levels and shaded from direct solar radiation. Paired sensors were also buried 3 cm under ground near each air temperature sensor to determine snow absence or presence. Select sites included relative humidity sensors to indicate whether precipitation was occurring. Data were collected every 3-4 h from December 2015 to Sept 2018 (with ongoing collection). Code for analysis of daily mean, minimum, maximum, and temperature change with elevation (lapse rates) are available on Github (https://doi.org/10.5281/zenodo.3239539). The sensor download and intermediate data products are available on HydroShare at (http://www.hydroshare.org/resource/222e832d3df24dea9bae9bbeb6f4219d) with publicly accessible visualization available from the Nooksack Observatory at data.cuahsi.org. Hydrologic models are generally structured with a single annual average lapse rate parameter which assumes a linear temperature gradient with elevation. The daily data (2016-2018) is used as part of ongoing studies on the non-linear dynamics and temporal variability of temperature with elevation to improve assessments of watershed function and salmon habitat.

54 ENVIRONMENTAL SCIENCES↗

Dynamic mesoscale model of dipolar fluids via fluctuating hydrodynamics

Fluctuating hydrodynamics (FHD) is a general framework of mesoscopic modeling and simulation based on conservational laws and constitutive equations of linear and nonlinear responses. However, explicit representation of electrical forces in FHD has yet to appear. In this work, we devised an Ansatz for the dynamics of dipole moment densities that is linked with the Poisson equation of the electrical potential ϕ in coupling to the other equations of FHD. The resulting ϕ-FHD equations then serve as a platform for integrating the essential forces, including electrostatics in addition to hydrodynamics, pressure-volume equation of state, surface tension, and solvent-particle interactions that govern the emergent behaviors of molecular systems at an intermediate scale. This unique merit of ϕ-FHD is illustrated by showing that the water dielectric function and ion hydration free energies in homogeneous and heterogenous systems can be captured accurately via the mesoscopic simulation. Furthermore, we show that the field variables of ϕ-FHD can be mapped from the trajectory of an all-atom molecular dynamics simulation such that model development and parametrization can be based on the information obtained at a finer-grained scale. With the aforementioned multiscale capabilities and a spatial resolution as high as 5 Å, the ϕ-FHD equations represent a useful semi-explicit solvent model for the modeling and simulation of complex systems, such as biomolecular machines and nanofluidics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing variation in faecal glucocorticoid concentrations in gray whales exposed to anthropogenic stressors

Understanding how individual animals respond to stressors behaviourally and physiologically is a critical step towards quantifying long-term population consequences and informing management efforts. Glucocorticoid (GC) metabolite accumulation in various matrices provides an integrated measure of adrenal activation in baleen whales and could thus be used to investigate physiological changes following exposure to stressors. In this study, we measured GC concentrations in faecal samples of Pacific Coast Feeding Group (PCFG) gray whales (Eschrichtius robustus) collected over seven consecutive years to assess the association between GC content and metrics of exposure to sound levels and vessel traffic at different temporal scales, while controlling for contextual variables such as sex, reproductive status, age, body condition, year, time of year and location. We develop a Bayesian Generalized Additive Modelling approach that accommodates the many complexities of these data, including non-linear variation in hormone concentrations, missing covariate values, repeated samples, sampling variability and some hormone concentrations below the limit of detection. Estimated relationships showed large variability, but emerging patterns indicate a strong context-dependency of physiological variation, depending on sex, body condition and proximity to a port. Our results highlight the need to control for baseline hormone variation related to context, which otherwise can obscure the functional relationship between faecal GCs and stressor exposure. Therefore, extensive data collection to determine sources of baseline variation in well-studied populations, such as PCFG gray whales, could shed light on cetacean stress physiology and be used to extend applicability to less-well-studied taxa. GC analyses may offer greatest utility when employed as part of a suite of markers that, in aggregate, provide a multivariate measure of physiological status, better informing estimates of individuals’ health and ultimately the consequences of anthropogenic stressors on populations.

59 BASIC BIOLOGICAL SCIENCES↗

KAYENTA: Theory and User's Guide

The physical foundations and domain of applicability of the Kayenta constitutive model are presented along with descriptions of the source code and user instructions. Kayenta, which is an outgrowth of the Sandia GeoModel, includes features and fitting functions appropriate to a broad class of materials including rocks, rock-like engineered materials (such as concretes and ceramics), and metals. Fundamentally, Kayenta is a computational framework for generalized plasticity models. As such, it includes a yield surface, but the term (3z(Byield(3y (Bis generalized to include any form of inelastic material response (including microcrack growth and pore collapse) that can result in non-recovered strain upon removal of loads on a material element. Kayenta supports optional anisotropic elasticity associated with joint sets, as well as optional deformation-induced anisotropy through kinematic hardening (in which the initially isotropic yield surface is permitted to translate in deviatoric stress space to model Bauschinger effects). The governing equations are otherwise isotropic. Because Kayenta is a unification and generalization of simpler models, it can be run using as few as 2 parameters (for linear elasticity) to as many as 40 material and control parameters in the exceptionally rare case when all features are used. For high-strain-rate applications, Kayenta supports rate dependence through an overstress model. Isotropic damage is modeled through loss of stiffness and strength.

58 GEOSCIENCES↗

Numerical analysis of flow non-uniformity in the hot gas manifold of the Space Shuttle main engine

Three-dimensional viscous flow in a conceptual hot gas manifold (HGM) for the Space Shuttle Main Engine High Pressure Fuel Turbopump (SSME HPFTP) was numerically analyzed. A finite difference scheme was used to solve the Navier-Stokes equations. The exact geometry of the SSME HGM was modeled using boundary fitted curvilinear coordinates and the General Interpolants Method (GIM) code. Slight compressibility of the subsonic flow was modeled using a linearized equation of state with artificial compressibility. A time relaxation method was used to obtain a steady state solution. The feasibility and potential usefulness of computational methods in assisting the design of SSME components which involves the flow of fluids within complex geometrical shapes is demonstrated.

Thoenes, J.↗