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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.

Lee-Rausch, Elizabeth M.↗
NASA NTRS1991

Improved large perturbation propulsion models for control system design (1988-1989) and large perturbation models of high velocity propulsion systems (1989-1990) and reduced order propulsion models for control system design (1990-1991)

Methods for modeling high speed propulsion systems will be discussed. Included in this category are internal flow propulsion systems without rotating machinery, such as inlets, ramjets, and scramjets. Among the modeling topics discussed are modeling of linear isentropic flow, heat exchange, gasdynamics, lumped parameter systems, and infinite dimensional systems. Furthermore, a generalized overview of modeling high speed propulsion systems is presented in this collection of papers.

Hartley, Tom T.↗
NASA NTRS1983

PAN AIR modeling studies

PAN AIR is a computer program that predicts subsonic or supersonic linear potential flow about arbitrary configurations. The code's versatility and generality afford numerous possibilities for modeling flow problems. Although this generality provides great flexibility, it also means that studies are required to establish the dos and don'ts of modeling. The purpose of this paper is to describe and evaluate a variety of methods for modeling flows with PAN AIR. The areas discussed are effects of panel density, internal flow modeling, forebody modeling in subsonic flow, propeller slipstream modeling, effect of wake length, wing-tail-wake interaction, effect of trailing-edge paneling on the Kutta condition, well- and ill-posed boundary-value problems, and induced-drag calculations. These nine topics address problems that are of practical interest to the users of PAN AIR.

Towne, M. C.↗
NASA NTRS1959

Static Stability Characteristics of Three Thick Wing Models with Parabolic Plan Forms at a Mach Number of 3.11

An experimental investigation has been made to determine the static stability characteristics of three thick wing models with parabolic plan forms at a Mach number of 3.11 for angles of attack from about -6 to 16 deg. The primary variable was aspect ratio, with the plan-form area and the ratio of base height to span kept the same for all three models. All models had stable, linear pitching-moment curves about the quarter chord of the wing mean aerodynamic chord. The model with the lowest aspect ratio attained a maximum untrimmed lift-drag ratio of about 5.0 at an angle of attack of about 8 deg. Increasing the aspect ratio (which was accompanied by an increase in base area because the ratio of the base height to span was kept constant) caused a decrease in maximum lift-drag ratio. All models were directionally stable for the range of angle of attack of the tests. Addition of a vertical tail to the models caused an increase in the directional stability over the angle-of-attack range. In general, the lateral aerodynamic characteristics of the models were not linear functions of angle of attack over any appreciable angle-of-attack range.

Queijo, M. J.↗
NASA NTRS1993

A Regional Analysis of Cloudy Mean Spherical Albedo over the Marine Stratocumulus Region and the Tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. This study examines the relationships between cloud properties and cloud fraction in order to supplement grid scale parameterizations. The satellite data used in this study is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.50 x 2.50 latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakleys (1991) theory of reflection for uniform and broken layered cloud and to Kedem, et al.(1990) findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 in) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al. theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗
NASA NTRS1993

A regional analysis of cloudy mean spherical albedo over the marine stratocumulus region and the tropical Atlantic Ocean

Since clouds are the largest variable in Earth's radiation budget, it is critical to determine both the spatial and temporal characteristics of their radiative properties. The relationships between cloud properties and cloud fraction are studied in order to supplement grid scale parameterizations. The satellite data used is from three hourly ISCCP (International Satellite Cloud Climatology Project) and monthly ERBE (Earth Radiation Budget Experiment) data on a 2.5 deg x 2.5 deg latitude-longitude grid. Mean cloud spherical albedo, the mean optical depth distribution, and cloud fraction are examined and compared off the coast of California and the mid-tropical Atlantic for July 1987 and 1988. Individual grid boxes and spatial averages over several grid boxes are correlated to Coakley's theory of reflection for uniform and broken layered cloud and to Kedem, et al.'s findings that rainfall volume and fractional area of rain in convective systems is linear. Kedem's hypothesis can be expressed in terms of cloud properties. That is, the total volume of liquid in a box is a linear function of cloud fraction. Results for the marine stratocumulus regime indicate that albedo is often invariant for cloud fractions of 20% to 80%. Coakley's satellite model of small and large clouds with cores (1 km) and edges (100 m) is consistent with this observation. The cores maintain high liquid water concentrations and large droplets while the edges contain low liquid water concentrations and small droplets. Large clouds are just a collection of cores. The mean optical depth (TAU) distributions support the above observation with TAU values of 3.55 to 9.38 favored across all cloud fractions. From these results, a method based upon Kedem, et al's theory is proposed to separate the cloud fraction and liquid water path (LWP) calculations in a general circulation model (GCM). In terms of spatial averaging, a linear relationship between albedo and cloud fraction is observed. For tropical locations outside the Intertropical Convergence Zone (ITCZ), results of cloud fraction and albedo spatial averaging followed that of the stratus boxes containing few overcast scenes. Both the ideas of Coakley and Kedem, et al. apply. Within the ITCZ, the grid boxes tended to have the same statistical properties as stratus boxes containing many overcast scenes. Because different dynamical forcing mechanisms are present, it is difficult to devise a method for determining subgrid scale variations. Neither of the theories proposed by Kedem, et al. or Coakley works well for the boxes with numerous overcast scenes.

Ginger, Kathryn M.↗
NASA NTRS1993

Computation of selected eigenvalues of generalized eigenvalue problems

The Arnoldi (1951) and Lanczos (1950) iterative algorithms are used to obtain a few selected eigenvalues of the generalized eigenvalue problem A(x) = lambdaB(x), where A and B are n x n nonsymmetric banded complex matrices. A shift-and-invert strategy is used to increase the rate of convergence toward the desired eigenvalues, and the two approaches are evaluated in view of their performance in a model problem concerning the linear stability of compressible boundary layers. A general scheme for computing the eigenvalues that lie within a 'box' in the complex plane is presented.

Nayar, Narinder↗
NASA NTRS1992

Pinatubo eruption winter climate effects: Model versus observations

Large volcanic eruptions, in addition to the well-known effect of producing global cooling for a year or two, have been observed to produce shorter-term responses in the climate system involving non-linear dynamical processes. In this paper, we use the ECHAM2 general circulation model forced with stratospheric aerosols to test some of these ideas. Run in a perpetual-January mode, with tropical stratospheric heating from the volcanic aerosols typical of the 1982 El Chichon eruption or the 1991 Pinatubo eruption, we find a dynamical response with an increased polar night jet in the Northern Hemisphere (NH) and stronger zonal winds which extended down into the troposphere. The Azores High shifts northward with increased tropospheric westerlies at 60N and increased easterlies at 30N. Surface temperatures are higher both in northern Eurasia and North America, in agreement with observations for the NH winters or 1982-83 and 1991-92 as well as the winters following the other 10 largest volcanic eruptions since 1883.

Graf, HANS-F.↗
NASA NTRS1972

Hamiltonian cosmology.

The study of cosmological models by means of equations of motion in Hamiltonian form is considered. Hamiltonian methods applied to gravity seem to go back to Rosenfeld (1930), who constructed a quantum-mechanical Hamiltonian for linearized general relativity theory. The first to notice that cosmologies provided a simple model in which to demonstrate features of Hamiltonian formulation was DeWitt (1967). Applications of the ADM formalism to homogeneous cosmologies are discussed together with applications of the Hamiltonian formulation, giving attention also to Bianchi-type universes. Problems involving the concept of superspace and techniques of quantization are investigated.

Ryan, M.↗
NASA NTRS1987

Fuel-air mixing and combustion in a two-dimensional Wankel engine

A two-equation turbulence model, an algebraic grid generalization method, and an approximate factorization time-linearized numerical technique are used to study the effects of mixture stratification at the intake port and gaseous fuel injection on the flow field and fuel-air mixing in a two-dimensional rotary engine model. The fuel distribution in the combustion chamber is found to be a function of the air-fuel mixture fluctuations at the intake port. It is shown that the fuel is advected by the flow field induced by the rotor and is concentrated near the leading apex during the intake stroke, while during compression, the fuel concentration is highest near the trailing apex and is lowest near the rotor. It is also found that the fuel concentration near the trailing apex and rotor is small except at high injection velocities.

Shih, T. I.-P.↗
DOE OSTI2020

Computational modelling and experimental tank testing of the multi float WaveSub under regular wave forcing

A submerged wave device generates energy from the relative motion of floating bodies. In WaveSub, three floats are joined to a reactor; each connected to a spring and generator. Electricity generated damps the orbital movements of the floats. The forces are non-linear and each float interacts with the others. Tuning to the wave climate is achieved by changing the line lengths, so there is a need to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design. to understand the performance trade-offs for a large number of configurations. This requires an efficient, large displacement, multidirectional, multi-body numerical scheme. Results from a 1/25 scale wave basin experiment are described. Here, we show that a time domain linear potential flow formulation (Nemoh, WEC-Sim) can match the tank testing provided that suitably tuned drag coefficients are employed. Inviscid linear potential models can match some wave device experiments; however, additional viscous terms generally provide better accuracy. Scale experiments are also prone to mechanical friction, and we estimate friction terms to improve the correlation further. The resulting error in mean power between numerical and physical models is approximately 10%. Predicted device movement shows a good match. Overall, drag terms in time domain wave energy modelling will improve simulation accuracy in wave renewable energy device design.

16 TIDAL AND WAVE POWER↗
NASA NTRS1989

Vibration analysis of rotor systems using reduced subsystem models

A general impedance method using reduced submodels has been developed for the linear dynamic analysis of rotor systems. Formulated in terms of either modal or physical coordinates of the subsystems, the method enables imbalance responses at specific locations of the rotor systems to be efficiently determined from a small number of 'master' degrees of freedom. To demonstrate the capability of this impedance approach, the Space Shuttle Main Engine high-pressure oxygen turbopump has been investigated to determine the bearing loads due to imbalance. Based on the same formulation, an eigenvalue analysis has been performed to study the system stability. A small 5-DOF model has been utilized to illustrate the application of the method to eigenvalue analysis. Because of its inherent characteristics of allowing formulation of reduced submodels, the impedance method can significantly increase the computational speed.

Fan, Uei-Jiun↗
NASA NTRS1990

Stability analysis of spacecraft power systems

The problems in applying standard electric utility models, analyses, and algorithms to the study of the stability of spacecraft power conditioning and distribution systems are discussed. Both single-phase and three-phase systems are considered. Of particular concern are the load and generator models that are used in terrestrial power system studies, as well as the standard assumptions of load and topological balance that lead to the use of the positive sequence network. The standard assumptions regarding relative speeds of subsystem dynamic responses that are made in the classical transient stability algorithm, which forms the backbone of utility-based studies, are examined. The applicability of these assumptions to a spacecraft power system stability study is discussed in detail. In addition to the classical indirect method, the applicability of Liapunov's direct methods to the stability determination of spacecraft power systems is discussed. It is pointed out that while the proposed method uses a solution process similar to the classical algorithm, the models used for the sources, loads, and networks are, in general, more accurate. Some preliminary results are given for a linear-graph, state-variable-based modeling approach to the study of the stability of space-based power distribution networks.

Halpin, S. M.↗
DOE OSTI2022

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Generalized multifractality at spin quantum Hall transition

Generalized multifractality characterizes scaling of eigenstate observables at Anderson-localization critical points. We explore generalized multifractality in 2D systems, with the main focus on the spin quantum Hall (SQH) transition in superconductors of symmetry class C. Relations and differences with the conventional integer quantum Hall (IQH) transition are also studied. Using the field-theoretical formalism of non-linear sigma-model, we derive the pure-scaling operators representing generalizing multifractality and then “translate” them to the language of eigenstate observables. Performing numerical simulations on network models for SQH and IQH transitions, we confirm the analytical predictions for scaling observables and determine the corresponding exponents. Remarkably, the generalized-multifractality exponents at the SQH critical point strongly violate the generalized parabolicity of the spectrum, which implies violation of the local conformal invariance at this critical point.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Analysis of a Computational Framework for Bayesian Inverse Problems: Ensemble Kalman Updates and MAP Estimators under Mesh Refinement

This paper analyzes a popular computational framework to solve infinite-dimensional Bayesian inverse problems, discretizing the prior and the forward model in a finite-dimensional weighted inner product space. We demonstrate the benefit of working on a weighted space by establishing operator-norm bounds for finite element and graph-based discretizations of Matérn-type priors and deconvolution forward models. For linear-Gaussian inverse problems, we develop a general theory to characterize the error in the approximation to the posterior. We also embed the computational framework into ensemble Kalman methods and MAP estimators for nonlinear inverse problems. Furthermore, our operator-norm bounds for prior discretizations guarantee the scalability and accuracy of these algorithms under mesh refinement.

Bayesian inverse problem↗
NASA NTRS1991

On Markov parameters in system identification

A detailed discussion of Markov parameters in system identification is given. Different forms of input-output representation of linear discrete-time systems are reviewed and discussed. Interpretation of sampled response data as Markov parameters is presented. Relations between the state-space model and particular linear difference models via the Markov parameters are formulated. A generalization of Markov parameters to observer and Kalman filter Markov parameters for system identification is explained. These extended Markov parameters play an important role in providing not only a state-space realization, but also an observer/Kalman filter for the system of interest.

Phan, Minh↗
NASA NTRS1994

Transfer function modeling of damping mechanisms in distributed parameter models

This work formulates a method for the modeling of material damping characteristics in distributed parameter models which may be easily applied to models such as rod, plate, and beam equations. The general linear boundary value vibration equation is modified to incorporate hysteresis effects represented by complex stiffness using the transfer function approach proposed by Golla and Hughes. The governing characteristic equations are decoupled through separation of variables yielding solutions similar to those of undamped classical theory, allowing solution of the steady state as well as transient response. Example problems and solutions are provided demonstrating the similarity of the solutions to those of the classical theories and transient responses of nonviscous systems.

Slater, J. C.↗