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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Model Resolution Sensitivity of the Simulation of North Atlantic Oscillation Teleconnections to Precipitation Extremes

We evaluate a high-resolution (0.25°), four-member ensemble simulation of the global climate (1979–2005) with the U.S. Department of Energy's Energy Exascale Earth System Model v0.3—forced with observed ocean surface temperatures and sea ice extent—for its ability to represent the North Atlantic Oscillation (NAO) teleconnections to winter precipitation extremes over western Europe. As compared to the low-resolution model (1°), it simulates a stronger impact of NAO on daily precipitation extremes over the western slopes of mountain ranges over southwestern Norway, northwestern United Kingdom, and the Western Balkan states. Precipitation extremes and their linear relationship with NAO are quantified using the generalized extreme value distribution. NAO-dependent large-scale (stratiform) precipitation intensity strengthens in the high-resolution model on seasonal time scales and plays a dominant role during simulated daily precipitation extremes. Improvements in the high-resolution model over these varied-topography regions largely appear to be due to finer resolved scales of motion that amplify NAO-dependent seasonal vertical moisture fluxes and enhance stable condensation. However, the high-resolution model simulates a weaker than observed impact of NAO on extratropical cyclone activity and total precipitable water, generally underperforming the low-resolution model These effects possibly offset the impact of enhanced vertical moisture fluxes on NAO-dependent precipitation extremes in the high-resolution model in these regions. Over the southwestern Iberian peninsula, the high-resolution model underperforms the low-resolution model simulating weaker than observed impact of NAO on precipitation extremes. Finally, this appears to be due to the reduction in total precipitable water despite an increase in NAO-dependent extratropical activity there.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1981

Analysis of a high-performance tubular solar collector

This article analyzes the thermal performance of a new vacuum tube solar collector. The assumptions and mathematical modeling are presented. The problem is reduced to the formulation of two simultaneous linear differential equations characterizing the collector thermal behavior. After applying the boundary conditions, a general solution is obtained which is found similar to the general Hottel, Whillier and Bliss form, but with a complex flow factor. The details of the two-dimensional thermal model of the solar collector at steady state is also presented to include the computer simulation and the performance parameterization. Comparison of the simulated performance with the manufacturer's test data showed good agreement at wide ranges of operating conditions.

Lansing, F. L.↗
NASA NTRS2006

Room-temperature Broadening and Pressure-shift Coefficients in the nu(exp 2) Band of CH3D-O2: Measurements and Semi-classical Calculations

We report measured Lorentz O2-broadening and O2-induced pressure-shift coefficients of CH3D in the nu(exp 2) fundamental band. Using a multispectrum fitting technique we have analyzed 11 laboratory absorption spectra recorded at 0.011 cm(exp 1) resolution using the McMath-Pierce Fourier transform spectrometer, Kitt Peak, Arizona. Two absorption cells with path lengths of 10.2 and 25 cm were used to record the spectra. The total sample pressures ranged from 0.98 to 339.85 Torr with CH3D volume mixing ratios of 0.012 in oxygen. We report measurements for O2 pressure-broadening coefficients of 320 nu(exp 2) transitions with quantum numbers as high as J0(sup w) = 17 and K = 14, where K(sup w) = K' is equivalent to K (for a parallel band). The measured O2-broadening coefficients range from 0.0153 to 0.0645 cm(exp -1) atm(exp -1) at 296 K. All the measured pressure-shifts are negative. The reported O2-induced pressure-shift coefficients vary from about -0.0017 to -0.0068 cm(exp -1) atm(exp -1). We have examined the dependence of the measured broadening and shift parameters on the J(sup W), and K quantum numbers and also developed empirical expressions to describe the broadening coefficients in terms of m (m = -J(sup W), J(sup W), and J(sup w) + 1 in the QP-, QQ-, and QR-branch, respectively) and K. On average, the empirical expressions reproduce the measured broadening coefficients to within 4.4%. The O2-broadening and pressure shift coefficients were calculated on the basis of a semiclassical model of interacting linear molecules performed by considering in addition to the electrostatic contributions the atom-atom Lennard-Jones potential. The theoretical results of the broadening coefficients are generally larger than the experimental data. Using for the trajectory model an isotropic Lennard-Jones potential derived from molecular parameters instead of the spherical average of the atom-atom model, a better agreement is obtained with these data, especially for |m| <= 12 values (11.3% for the first calculation and 8.1% for the second calculation). The O2-pressure shifts whose vibrational contribution are either derived from parameters fitted in the QQ-branch of selfinduced shifts of CH3D or those obtained from pressure shifts induced by Xe in the nu(sup 3) band of CH3D are in reasonable agreement with the scattered experimental data (17.0% for the first calculation and 18.7% for the second calculation).

Predoi-Cross, Adriana↗
DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes

XGC (X-point Gyrokinetic Code) is a whole-volume, total-f gyrokinetic particle-in-cell code developed for modelling tokamaks.In recent work, XGC has been extended to model more general 3D toroidal magnetic configurations, such as stellarators.These improvements have resulted in the XGC-S version.In this paper, XGC-S is benchmarked in the reduced delta-f limit for linear electrostatic ion temperature gradient-driven microinstabilities, which can underlie turbulent transport in stellarators.An initial benchmark of XGC-S in tokamak geometry shows good agreement with the XGC1, ORB5, and global GENE codes.A benchmark between XGC-S and the EUTERPE global gyrokinetic code for stellarators has also been performed, this time in geometry of the optimised stellarator Wendelstein 7-X.Good agreement has been found for the mode number spectrum, mode structure, and growth rate.

M, Cole↗
DOE OSTI2011

Dark degeneracy and interacting cosmic components

We study some properties of the dark degeneracy, which is the fact that what we measure in gravitational experiments is the energy-momentum tensor of the total dark sector, and any split into components (as in dark matter and dark energy) is arbitrary. In fact, just one dark fluid is necessary to obtain exactly the same cosmological and astrophysical phenomenology as the {Lambda}CDM model. We work explicitly the first-order perturbation theory and show that beyond the linear order the dark degeneracy is preserved under some general assumptions. Then we construct the dark fluid from a collection of interacting fluids. Finally, we try to break the degeneracy with a general class of couplings to baryonic matter. Nonetheless, we show that these interactions can also be understood in the context of the {Lambda}CDM model as between dark matter and baryons. For this last investigation we choose two independent parametrizations for the interactions, one inspired by electromagnetism and the other by chameleon theories. Then, we constrain them with a joint analysis of CMB and supernovae observational data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes

XGC (X-point Gyrokinetic Code) is a whole-volume, total-f gyrokinetic particle-in-cell code developed for modeling tokamaks. Recently, XGC has been extended to model more general 3D toroidal magnetic configurations, such as stellarators. These improvements have concluded in the XGC-S version. In this paper, XGC-S is benchmarked in the reduced delta-f limit for linear electrostatic ion temperature gradient-driven microinstabilities, which can underlie turbulent transport in stellarators. An initial benchmark of XGC-S in tokamak geometry shows good agreement with the XGC1, ORB5, and global GENE codes. A benchmark between XGC-S and the EUTERPE global gyrokinetic code for stellarators has also been performed, here in the geometry of the optimized stellarator Wendelstein 7-X. Good agreement has been discovered for the mode number spectrum, mode structure, and growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

The dynamics and control of large flexible space structures X, part 1

The effect of delay in the control system input on the stability of a continuously acting controller which is designed without considering the delay is studied. The stability analysis of a second order plant is studied analytically and verified numerically. For this example it is found that the system becomes unstable for a delay which is equivalent to only 16 percent of its natural period of motion. It is also observed that even a small amount of natural damping in the system can increase the amount of delay that can be tolerated before the onset of instability. The delay problem is formulated in the discrete time domain and an analysis procedure suggested. The maximum principle from optimal control theory is applied to minimize the time required for the slewing of a general rigid spacecraft. The slewing motion need not be restricted to a single axis maneuver. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. Numerical examples based on the rigidized in-orbit model of the SCOLE also include the more general reflector line-of-sight slewing maneuvers.

Bainum, Peter M.↗
NASA NTRS2002

A New Forced Oscillation Capability for the Transonic Dynamics Tunnel

A new forced oscillation system has been installed and tested at NASA Langley Research Center's Transonic Dynamics Tunnel (TDT). The system is known as the Oscillating Turntable (OTT) and has been designed for the purpose of oscillating, large semispan models in pitch at frequencies up to 40 Hz to acquire high-quality unsteady pressure and loads data. Precisely controlled motions of a wind-tunnel model on the OTT can yield unsteady aerodynamic phenomena associated with flutter, limit cycle oscillations, shock dynamics, and non-linear aerodynamic effects on many vehicle configurations. This paper will discuss general design and components of the OTT and will present test data from performance testing and from research tests on two rigid semispan wind-tunnel models. The research tests were designed to challenge the OTT over a wide range of operating conditions while acquiring unsteady pressure data on a small rectangular supercritical wing and a large supersonic transport wing. These results will be presented to illustrate the performance capabilities, consistency of oscillations, and usefulness of the OTT as a research tool.

Piatak, David J.↗
NASA NTRS1998

Current Limitations on VLBI Accuracy

The contribution of VLBI to geophysics and geodesy arises from its ability to measure distances between stations in a network and to determine the orientation of stations in the network as well as the orientation of the network with respect to the external reference frame of extragalactic radio objects. Integrating nearly two decades of observations provides useful information about station positions and velocities and the orientation of the Earth, but the complications of the real world and the limitations of observing, modeling and analysis prevent recovery of all effects. Of the factors that limit the accuracy of seemingly straightforward geodetic parameters, the neutral propagation medium has been subject to the greatest scrutiny, but the treatment of the mapping function, the wet component and spatial/temporal inhomogeneities is still improving. These affect both the terrestrial scale and consistency over time. The modeling of non-secular site motions (tides and loading) has increased in sophistication, but there are some differences between the models and the observations. VLBI antennas are massive objects, so their behavior is quite unlike GPS monuments, but antenna deformations add some (generally) unmodeled signal. Radio sources used in geodetic VLBI observations are selected for strength and (relative) absence of structure, but apparent changes in position can leak into geodetic parameters. A linear rate of change of baseline or site parameters is the simplest model and its error improves with time span. However, in most cases the VLBI data distribution is insufficient to look for real non-linear behavior that might affect the average rate. A few sites have multiple VLBI antennas, and some show small differences in rate. VLBI intrinsically measures relative positions and velocities, but individual site positions and velocities are generally more useful. The creation of the VLBI terrestrial reference frame, which transforms relative information into individual results, is an empirical process that has intrinsic errors. While UT1 is uniquely measured by VLBI, the geographical distribution and availability of VLBI stations, especially in the southern hemisphere, and the consistency of the VLBI terrestrial reference frame may limit the accuracy of Earth orientation measurements. The effects of particular error sources on geodetic and geophysical parameters derived from VLBI data will be illustrated.

Ma, Chopo↗
NASA NTRS1990

Climate Fluctuations and Climate Sensitivity

Some evidence is presented that the main part of the atmospheric climate system is such that small forcings in the heat balance lead to linear responses in the surface temperature field. By examining first a noise forced energy-balance climate model and then comparing it with a long run of a highly symmetrical general circulation model, one finds a remarkable connection between spatial autocorrelation statistics and the thermal influence function for a point heat source. These findings are brought together to indicate that this particular climatological field may be largely governed by linear processes.

North, Gerald R.↗
NASA NTRS2004

Cyclic Oxidation Modeling and Life Prediction

The cyclic oxidation process can be described as an iterative scale growth and spallation sequence by a number of similar models. Model input variable include oxide scale type and growth parameters, spalling geometry, spall constant, and cycle duration. Outputs include net weight change, the amounts of retained and spalled oxide, the total oxygen and metal consumed, and the terminal rates of weight loss and metal consumption. All models and their variations produce a number of similar characteristic features. In general, spalling and material consumption increase to a steady state rate, at which point the retained scale approaches a constant and the rate of weight loss becomes linear. For one model, this regularity was demonstrated as dimensionless, universal expressions, obtained by normalizing the variables by critical performance factors. These insights were enabled through the use of the COSP for Windows cyclic oxidation spalling program.

Smialek, James L.↗
NASA NTRS1985

Rift systems on Venus: An assessment of mechanical and thermal models

The formation and distribution of major tectonic features on Venus are closely linked to the dominant mechanism of lithospheric heat loss. Among the most spectacular and extensive of the major tectonic features on Venus are the Chasmata, deep linear valleys generally interpreted to be the products of lithospheric extension and rifting. Systems of chasmata and related features can be traced along several tectonic zones up to 20,000 km in linear extent. Mechanical and thermal models for terrestrial continental-rifting are applied to the rift systems of Venus. The models are tested against known topographic and tectonic characteristics of Venus chasmata as well as independent information on the physical properties of the Venus crust and lithosphere.

Solomon, S. C.↗
DOE OSTI2022

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Effects of Hall current and electron temperature anisotropy on proton fire-hose instabilities

The standard magnetohydrodynamic (MHD) theory predicts that the Alfvén wave may become fire-hose unstable for β{sub ∥}−β{sub ⊥}>2. In this study, we examine the proton fire-hose instability (FHI) based on the gyrotropic two-fluid model, which incorporates the ion inertial effects arising from the Hall current and electron temperature anisotropy but neglects the electron inertia in the generalized Ohm's law. The linear dispersion relation is derived and analyzed which in the long wavelength approximation, λ{sub i}k→0 or α{sub e}=μ{sub 0}(p{sub ∥,e}−p{sub ⊥,e})/B{sup 2}=1, recovers the ideal MHD model with separate temperature for ions and electrons. Here, λ{sub i} is the ion inertial length and k is the wave number. For parallel propagation, both ion cyclotron and whistler waves become propagating and growing for β{sub ∥}−β{sub ⊥}>2+λ{sub i}{sup 2}k{sup 2}(α{sub e}−1){sup 2}/2. For oblique propagation, the necessary condition for FHI remains to be β{sub ∥}−β{sub ⊥}>2 and there exist one or two unstable fire-hose modes, which can be propagating and growing or purely growing. For large λ{sub i}k values, there exists no nearly parallel FHI leaving only oblique FHI and the effect of α{sub e}>1 may greatly enhance the growth rate of parallel and oblique FHI. The magnetic field polarization of FHI may be reversed due to the sign change associated with (α{sub e}−1) and the purely growing FHI may possess linear polarization while the propagating and growing FHI may possess right-handed or left-handed polarization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Forward Analytic Model of Neutron Time-of-Flight Signals for Inferring Ion Temperatures from MagLIF Experiments

A forward analytic model is required to rapidly simulate the neutron time-of-flight (nToF) signals that result from magnetized liner inertial fusion (MagLIF) experiments at Sandia’s Z Pulsed Power Facility. Various experimental parameters, such as the burn-weighted fuel-ion temperature and liner areal density, determine the shape of the nToF signal and are important for characterizing any given MagLIF experiment. Extracting these parameters from measured nToF signals requires an appropriate analytic model that includes the primary deuterium-deuterium neutron peak, once-scattered neutrons in the beryllium liner of the MagLIF target, and direct beamline attenuation. Here, mathematical expressions for this model were derived from the general-geometry time- and energy-dependent neutron transport equation with anisotropic scattering. Assumptions consistent with the time-of-flight technique were used to simplify this linear Boltzmann transport equation into a more tractable form. Models of the uncollided and once-collided neutron scalar fluxes were developed for one of the five nToF detector locations at the Z-Machine. Numerical results from these models were produced for a representative MagLIF problem and found to be in good agreement with similar neutron transport simulations. Twenty experimental MagLIF data sets were analyzed using the forward models, which were determined to only be significantly sensitive to the ion temperature. The results of this work were also found to agree with values obtained separately using a zero scatter analytic model and a high-fidelity Monte Carlo simulation. Finally, inherent difficulties in this and similar techniques are identified, and a new approach forward is suggested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measurement of Fission Product Concentration Profiles in AGR-3/4 TRISO Fuel Graphitic Matrix and Nuclear Graphites

The third Advanced Gas Reactor (AGR) irradiation experiment, AGR-3/4, was designed to investigate the migration of fission products in fuel compact graphitic matrix and reactor graphite components. Using destructive methods, radial fission product concentration profiles were measured for gamma-emitting fission products (e.g., Ag-110m, Cs-134, and Eu-154) and beta-emitting Sr-90 in irradiated graphitic and graphite components from six different AGR-3/4 irradiation capsules. These new measured concentration profiles can now be compared to non-destructive measurements and fission product transport simulations and will be used to derive new diffusivities and sorptivities to support refinement of fission product transport models and high-temperature gas-cooled reactor (HTGR) source-term analyses. Each capsule in the AGR-3/4 experiment had four fuel compacts in the middle of two concentric rings of graphitic matrix material, PCEA graphite, or IG-110 graphite. In addition to the approximately 1898 tristructural isotropic (TRISO) coated particles in each compact, there were 20 designed to fail (DTF) particles coated only in pyrocarbon so that they released fission products into the surrounding cylindrical rings of carbonaceous materials. Destructive sampling of the rings involved machining/milling material from around the circumference of the rings, collecting that material, and performing radiochemical analyses on it. Milling operations were performed in multiple steps or segments, and each segment was generally 0.508 mm (0.020 in) thick. Knowing the radial position at which each segment was milled, the volume of the milled material at each segment, and the fission product content in each segment, the radial fission product concentrations were constructed for select isotopes in each ring. Ag-110m profiles had the most variation. Some profiles were peaked at an inner or outer surface. Some were peaked at the middle of the ring wall thickness. Some increased radially outward, and some decreased radially outward. These types of variations and the fact that the measured profiles do not generally compare favorably with the transport model employed for AGR-3/4 may adversely impact the ability to extract reasonable transport parameters for this isotope. In many cases, the Cs-134 profiles decreased somewhat linearly in the outward radial direction, and in cursory comparisons, the shapes of these profiles appeared similar to those from model predictions. The step changes in concentration across the inner-outer ring gap were generally consistent with the model predictions as well. In some cases, there were local maxima in concentration at the outer surface of the rings. This suggests that fission products could have transported in the small gaps between the inner ring and the outer ring and between the outer ring and the sink ring such that some portion of a given fission product can bypass diffusion through the ring itself. The analysis of the small nubs on the outer surfaces of some of the outer rings revealed fission product concentrations in the nubs that were often higher than in the outermost segments of the rings. This further supports the hypothesis that short-circuit, gap transport occurred, causing relatively high surface concentrations on the outer surfaces of the rings. Eu-154 and Sr-90 profiles tend to have very similar shapes, suggesting that they transport via the same mechanisms. The observed profiles were indicative of a transport process where the isotopes are sorbed on the inner surface of the ring, but diffusion into the ring from that surface is quite slow. Some elevated concentrations of Sr-90 (relative to Eu-154) on the outer surface of a ring suggested that rapid, gas-gap transport of gaseous precursor Kr-90 and volatile Rb-90 could have occurred prior to their decaying to Sr-90. Overall, the Eu 154 and Sr-90 profiles were still very similar, which indicates that the transport of short-lived Sr-90 precursors is not a major effect. In some capsules, the qualitative Sr-90 behavior across the ring gaps was consistent with the model (using the available legacy Sr-90 transport parameters), but in other capsules the model was inconsistent with the measurements and seems to underestimate the amount of Sr-90 in the outer rings. The total ring Sr-90 inventories were estimated for all the rings that were subject to physical sampling. These results will be used to adjust the predicted particle and/or compact releases used in the AGR-3/4 fission product transport model. Given the different irradiation temperatures among the capsules and the rings, it was not possible to discern fundamental differences in the transport of isotopes within the different carbon materials, i.e., graphitic matrix, IG-110, or PCEA. It may be possible to do this in the course of determining transport from the concentration profiles in future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1989

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↗