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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Plasma Modeling Hierarchy and Verification Approach

This report reviews a hierarchy of formal mathematical models for describing plasma phenomena. Starting with the Boltzmann equation, a sequence of approximations and modeling assumptions can be made that progressively reduce to the equations for magnetohydrodynamics. Understanding the assumptions behind each of these models and their mathematical form is essential to appropriate use of each level of the hierarchy. A sequence of moment models of the Boltzmann equation are presented, then focused into a generalized three-fluid model for neutral species, electrons, and ions. This model is then further reduced to a two-fluid model, for which Braginskii described a useful closure. Further reduction of the two-fluid model yields a Generalized Ohm's Law model, which provides a connection to magnetohydrodynamic approaches. A verification approach based on linear plasma waves is presented alongside the model hierarchy, which is intended as an initial and necessary but not sufficient step for verification of plasma models within this hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1971

Space Trajectories Error Analysis (STEAP) Programs. Volume 1: Analytic manual, update

Manual revisions are presented for the modified and expanded STEAP series. The STEAP 2 is composed of three independent but related programs: NOMAL for the generation of n-body nominal trajectories performing a number of deterministic guidance events; ERRAN for the linear error analysis and generalized covariance analysis along specific targeted trajectories; and SIMUL for testing the mathematical models used in the navigation and guidance process. The analytic manual provides general problem description, formulation, and solution and the detailed analysis of subroutines. The programmers' manual gives descriptions of the overall structure of the programs as well as the computational flow and analysis of the individual subroutines. The user's manual provides information on the input and output quantities of the programs. These are updates to N69-36472 and N69-36473.

Source record↗
NASA NTRS1972

U-Th-Pb systematics in three Apollo 14 basalts and the problem of initial Pb in lunar rocks.

The isotopic composition of Pb and the elemental concentration of U, Th and Pb were measured on 'total' rock samples 14053, 14073 and 14310 and on mineral separates of 14310 and 14053. These are the first Pb-U isochrons obtained for lunar basalts and indicate a reasonable solution to the previous discrepancy between the different methods of 'absolute' age determination. The resulting U-Pb isochron ages are compatible with the Rb-Sr and K-Ar ages on the same rocks. However, it is not possible to establish a precise time of 'crystallization' from the Pb-U data because of the small angle of intersection between the linear arrays and the concordia curve. These data show that total rock model ages do not in general yield crystallization ages. The data on 14310 and 14053 show that these rocks were formed containing a highly radiogenic initial lead which accounts for the excessively high total rock model ages by the U-Th-Pb method. The data prove that at the time of extrusion of some basalts, unsupported lead with extremely high Pb-207/Pb-206 ratios was added to the lunar surface.

Tera, F.↗
NASA NTRS1992

The gravitational wave contribution to cosmic microwave background anisotropies and the amplitude of mass fluctuations from COBE results

A stochastic background of primordial gravitational waves may substantially contribute, via the Sachs-Wolfe effect, to the large-scale cosmic microwave background (CMB) anisotropies recently detected by COBE. This implies a bias in any resulting determination of the primordial amplitude of density fluctuations. We consider the constraints imposed on n is less than 1 ('tilted') power-law fluctuation spectra, taking into account the contribution from both scalar and tensor waves, as predicted by power-law inflation. The gravitational wave contribution to CMB anisotropies generally reduces the required rms level of mass fluctuation, thereby increasing the linear bias parameter, even in models where the spectral index is close to the Harrison-Zel'dovich value n = 1. This 'gravitational wave bias' helps to reconcile the predictions of CDM models with observations on pairwise galaxy velocity dispersion on small scales.

Lucchin, Francesco↗
NASA NTRS1996

Application of Hybrid Optimization-Expert System for Optimal Power Management on Board Space Power Station

The space power system has two sources of energy: photo-voltaic blankets and batteries. The optimal power management problem on-board has two broad operations: off-line power scheduling to determine the load allocation schedule of the next several hours based on the forecast of load and solar power availability. The nature of this study puts less emphasis on speed requirement for computation and more importance on the optimality of the solution. The second category problem, on-line power rescheduling, is needed in the event of occurrence of a contingency to optimally reschedule the loads to minimize the 'unused' or 'wasted' energy while keeping the priority on certain type of load and minimum disturbance of the original optimal schedule determined in the first-stage off-line study. The computational performance of the on-line 'rescheduler' is an important criterion and plays a critical role in the selection of the appropriate tool. The Howard University Center for Energy Systems and Control has developed a hybrid optimization-expert systems based power management program. The pre-scheduler has been developed using a non-linear multi-objective optimization technique called the Outer Approximation method and implemented using the General Algebraic Modeling System (GAMS). The optimization model has the capability of dealing with multiple conflicting objectives viz. maximizing energy utilization, minimizing the variation of load over a day, etc. and incorporates several complex interaction between the loads in a space system. The rescheduling is performed using an expert system developed in PROLOG which utilizes a rule-base for reallocation of the loads in an emergency condition viz. shortage of power due to solar array failure, increase of base load, addition of new activity, repetition of old activity etc. Both the modules handle decision making on battery charging and discharging and allocation of loads over a time-horizon of a day divided into intervals of 10 minutes. The models have been extensively tested using a case study for the Space Station Freedom and the results for the case study will be presented. Several future enhancements of the pre-scheduler and the 'rescheduler' have been outlined which include graphic analyzer for the on-line module, incorporating probabilistic considerations, including spatial location of the loads and the connectivity using a direct current (DC) load flow model.

Momoh, James↗
NASA NTRS2007

The Local Stellar Velocity Field via Vector Spherical Harmonics

We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism.We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) = (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) = (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star...

kinematics↗
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↗
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 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↗