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

The Cloud Feedback Model Intercomparison Project (CFMIP) contribution to CMIP6.

The primary objective of CFMIP is to inform future assessments of cloud feedbacks through improved understanding of cloud-climate feedback mechanisms and better evaluation of cloud processes and cloud feedbacks in climate models. However, the CFMIP approach is also increasingly being used to understand other aspects of climate change, and so a second objective has now been introduced, to improve understanding of circulation, regional-scale precipitation, and non-linear changes. CFMIP is supporting ongoing model inter-comparison activities by coordinating a hierarchy of targeted experiments for CMIP6, along with a set of cloud-related output diagnostics. CFMIP contributes primarily to addressing the CMIP6 questions 'How does the Earth system respond to forcing?' and 'What are the origins and consequences of systematic model biases?' and supports the activities of the WCRP Grand Challenge on Clouds, Circulation and Climate Sensitivity. A compact set of Tier 1 experiments is proposed for CMIP6 to address this question: (1) what are the physical mechanisms underlying the range of cloud feedbacks and cloud adjustments predicted by climate models, and which models have the most credible cloud feedbacks? Additional Tier 2 experiments are proposed to address the following questions. (2) Are cloud feedbacks consistent for climate cooling and warming, and if not, why? (3) How do cloud-radiative effects impact the structure, the strength and the variability of the general atmospheric circulation in present and future climates? (4) How do responses in the climate system due to changes in solar forcing differ from changes due to CO2, and is the response sensitive to the sign of the forcing? (5) To what extent is regional climate change per CO2 doubling state-dependent (non-linear), and why? (6) Are climate feedbacks during the 20th century different to those acting on long-term climate change and climate sensitivity? (7) How do regional climate responses (e.g. in precipitation) and their uncertainties in coupled models arise from the combination of different aspects of CO2 forcing and sea surface warming? CFMIP also proposes a number of additional model outputs in the CMIP DECK, CMIP6 Historical and CMIP6 CFMIP experiments, including COSP simulator outputs and process diagnostics to address the following questions. 1. How well do clouds and other relevant variables simulated by models agree with observations? 2. What physical processes and mechanisms are important for a credible simulation of clouds, cloud feedbacks and cloud adjustments in climate models? 3. Which models have the most credible representations of processes relevant to the simulation of clouds? 4. How do clouds and their changes interact with other elements of the climate system?

carbon dioxide↗

Lunar electromagnetic scattering. I - Propagation parallel to the diamagnetic cavity axis.

A general analytic solution is obtained for the interaction of the moon and its downstream cavity with a linearly polarized plane electromagnetic wave propagating parallel to the cavity axis. The solution is formulated in terms of a spherical moon model with arbitrary radially dependent electromagnetic parameters and a nonconducting cylindrical downstream cavity. Use is made of a number of approximations that are consistent with the physical nature of the interaction between the moon and the solar wind.

Schwartz, K.↗

Hovering impulsive noise: Some measured and calculated results

In-plane impulsive noise radiating from a hovering model rotor was measured in an anechoic environment. The hover acoustic signature was compared with existing theoretical prediction models with previous forward flight experiments using the same model rotor. These hover tests showed good experimental consistency with forward flight measurements, both in pressure level, and waveform character, over the range of Mach numbers tested (0.8 to 1.0). Generally poor correlation, however, was confirmed with current linear theory prediction efforts. Failure to predict both the peak pressure levels and the shape was reported, especially with increasing tip Mach number.

Boxwell, D. A.↗

Viscous relaxation of impact crater relief on Venus - Constraints on crustal thickness and thermal gradient

Models for the viscous relaxation of impact crater topography are used to constrain the crustal thickness (H) and the mean lithospheric thermal gradient beneath the craters on Venus. A general formulation for gravity-driven flow in a linearly viscous fluid has been obtained which incorporates the densities and temperature-dependent effective viscosities of distinct crust and mantle layers. An upper limit to the crustal volume of Venus of 10 to the 10th cu km is obtained which implies either that the average rate of crustal generation has been much smaller on Venus than on earth or that some form of crustal recycling has occurred on Venus.

Grimm, Robert E.↗

Shade images of forested areas obtained from LANDSAT MSS data

The pixel size in the present day Remote Sensing systems is large enough to include different types of land cover. Depending upon the target area, several components may be present within the pixel. In forested areas, generally, three main components are present: tree canopy, soil (understory), and shadow. The objective is to generate a shade (shadow) image of forested areas from multispectral measurements of LANDSAT MSS (Multispectral Scanner) data by implementing a linear mixing model, where shadow is considered as one of the primary components in a pixel. The shade images are related to the observed variation in forest structure, i.e., the proportion of inferred shadow in a pixel is related to different forest ages, forest types, and tree crown cover. The Constrained Least Squares (CLS) method is used to generate shade images for forest of eucalyptus and vegetation of cerrado using LANDSAT MSS imagery over Itapeva study area in Brazil. The resulted shade images may explain the difference on ages for forest of eucalyptus and the difference on three crown cover for vegetation of cerrado.

Shimabukuro, Yosio Edemir↗

Improvements in Iterative Convergence of FUN3D Solutions

This paper presents a hierarchical adaptive nonlinear iteration method (HANIM) implemented in NASA computational fluid dynamics code, FUN3D, to improve robustness and computational efficiency of FUN3Dsolutions. In contrast to the baseline iterative solver that relies on an approximate Jacobian, a simple multicolor Gauss-Seidel point-implicit iteration scheme, and linear CFL ramping, HANIM is based upon a hierarchy of modules including pre conditioner, generalized conjugate residual, realizability check, nonlinear control,and CFL adaption modules. HANIM performance is systematically compared with the performance of the baseline solver. The iterative solutions are compared for three aerodynamic benchmark cases: a subsonic separated flow around a hemisphere cylinder, a supersonic flow through a long duct, and a subsonic flow over the NASA wing-body juncture model. Two Reynolds-averaged Navier-Stokes turbulence models are used in these computations, namely, the negative variant of the linear one-equation Spalart-Allmar as model and its nonlinear extension based on quadratic constitutive relations.

Li Wang↗

Equicontrollability and the model following problem

Equicontrollability and its application to the linear time-invariant model-following problem are discussed. The problem is presented in the form of two systems, the plant and the model. The requirement is to find a controller to apply to the plant so that the resultant compensated plant behaves, in an input-output sense, the same as the model. All systems are assumed to be linear and time-invariant. The basic approach is to find suitable equicontrollable realizations of the plant and model and to utilize feedback so as to produce a controller of minimal state dimension. The concept of equicontrollability is a generalization of control canonical (phase variable) form applied to multivariable systems. It allows one to visualize clearly the effects of feedback and to pinpoint the parameters of a multivariable system which are invariant under feedback. The basic contributions are the development of equicontrollable form; solution of the model-following problem in an entirely algorithmic way, suitable for computer programming; and resolution of questions on system decoupling.

Curran, R. T.↗

Critical issues in viscous flow computations

In developing computer programs to numerically solve the Navier-Stokes equations, the purpose of the computation must be clearly kept in mind. In the Air Force, the purpose is to provide design information on non-linear aerodynamic phenomenon for aircraft that perform throughout the flight corridor. This translates into the requirement for a computer program which can solve the time averaged compressible Navier-Stokes equations (with a turbulence model) in three dimensions for generalized geometries. The intended application of the results then controls the priorities in addressing critical issues. Recurrent problem areas encountered in the study of viscous flow include: (1) grid generation for arbitrary geometry; (2) numerical difficulties; (3) turbulence models; (4) accuracy and efficiency; and (5) smearing of discontinuities.

Hankey, W. L., Jr.↗

A frost formation model and its validation under various experimental conditions

A numerical model that was used to calculate the frost properties for all regimes of frost growth is described. In the first regime of frost growth, the initial frost density and thickness was modeled from the theories of crystal growth. The 'frost point' temperature was modeled as a linear interpolation between the dew point temperature and the fog point temperature, based upon the nucleating capability of the particular condensing surfaces. For a second regime of frost growth, the diffusion model was adopted with the following enhancements: the generalized correlation of the water frost thermal conductivity was applied to practically all water frost layers being careful to ensure that the calculated heat and mass transfer coefficients agreed with experimental measurements of the same coefficients.

Dietenberger, M. A.↗

Derivation and definition of a linear aircraft model

A linear aircraft model for a rigid aircraft of constant mass flying over a flat, nonrotating earth is derived and defined. The derivation makes no assumptions of reference trajectory or vehicle symmetry. The linear system equations are derived and evaluated along a general trajectory and include both aircraft dynamics and observation variables.

Duke, Eugene L.↗

Microgravity vibration isolation: Optimal preview and feedback control

In order to achieve adequate low-frequency vibration isolation for certain space experiments an active control is needed, due to inherent passive-isolator limitations. Proposed here are five possible state-space models for a one-dimensional vibration isolation system with a quadratic performance index. The five models are subsets of a general set of nonhomogeneous state space equations which includes disturbance terms. An optimal control is determined, using a differential equations approach, for this class of problems. This control is expressed in terms of constant, Linear Quadratic Regulator (LQR) feedback gains and constant feedforward (preview) gains. The gains can be easily determined numerically. They result in a robust controller and offers substantial improvements over a control that uses standard LQR feedback alone.

Hampton, R. D.↗

Probabilistic nonlinear finite element analysis of composite structures

A probabilistic finite element analysis procedure for laminated composite shells is developed. A total Lagrangian finite element formulation, employing a degenerated three-dimensional laminated composite shell element with the full Green-Lagrange strains and first-order shear deformable kinematics, is used. The first-order second-moment technique for probabilistic finite element analysis of random fields is employed, and results are presented in the form of mean and variance of the structural response. Reliability calculations are made by using the first-order reliability method combined with sensitivity derivatives from the finite element analysis. Both ply-level and micromechanics-level random variables are incorporated, the latter by means of the Aboudi micromechanics model. Two sample problems are solved to verify the accuracy of the procedures developed and to quantify the variability of certain material type/structure combinations. In general, the procedure is quite effective in determining the response statistics and reliability for linear and geometric nonlinear behavior of laminated composite shells.

Engelstad, S. P.↗

Analytic MHD Theory for Earth's Bow Shock at Low Mach Numbers

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven M(A) and sonic M(s) Mach numbers are both much greater than 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first-order perturbation solution is identical to the solution found in the earlier theory. The second-order perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle theta between the normal to the bow shock and magnetic field is not close to 0 deg or 180 deg (the solutions are symmetric about 90 deg). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 deg. In this exceptional case the new solution does not improve on the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta approx. = 90 deg another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) is less than or = 2. MHD and gasdynamic simulations have produced empirical models in which the shock's standoff distance a(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between a(s) and X obtained from MHD simulations, a(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots of X and a(s) predicted for the Earth's bow shock. The plots show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the corresponding phenomenological relations in comparisons with data. Furthermore, significant differences exist between the standoff distances predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.

Grabbe, Crockett L.↗

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗

Nonlinear/linear unified thermal stress formulations - Transfinite element approach

A new unified computational approach for applicability to nonlinear/linear thermal-structural problems is presented. Basic concepts of the approach including applicability to nonlinear and linear thermal structural mechanics are first described via general formulations. Therein, the approach is demonstrated for thermal stress and thermal-structural dynamic applications. The proposed transfinite element approach focuses on providing a viable hybrid computational methodology by combining the modeling versatility of contemporary finite element schemes in conjunction with transform techniques and the classical Bubnov-Galerkin schemes. Comparative samples of numerical test cases highlight the capabilities of the proposed concepts.

Tamma, Kumar K.↗

Proteus three-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A computer code called Proteus 3D was developed to solve the three-dimensional, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The objective in this effort was to develop a code for aerospace propulsion applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The governing equations are solved in generalized nonorthogonal body fitted coordinates, by marching in time using a fully-coupled ADI solution procedure. The boundary conditions are treated implicitly. All terms, including the diffusion terms, are linearized using second-order Taylor series expansions. Turbulence is modeled using either an algebraic or two-equation eddy viscosity model. The thin-layer or Euler equations may also be solved. The energy equation may be eliminated by the assumption of constant total enthalpy. Explicit and implicit artificial viscosity may be used. Several time step options are available for convergence acceleration. The documentation is divided into three volumes. The Programmer's Reference contains detailed information useful when modifying the program. The program structure, the Fortran variables stored in common blocks, and the details of each subprogram are described.

Towne, Charles E.↗

Possible isotopic fractionation effects in sputtered minerals

A model which makes definite predictions for the fractionation of isotopes in sputtered material is discussed. The fractionation patterns are nonlinear, and the pattern for a particular set of isotopes depends on the chemical matrix within which those isotopes are contained. Calculations are presented for all nonmonoisotopic elements contained in the minerals perovskite, anorthite, ackermanite, enstatite, and troilite. All isotopes are fractionated at the level of approximately 4-6 deg/o per atomic mass unit. Oxygen is always positively fractionated (heavier isotopes sputtered preferentially), and heavier elements are generally negatively fractioned (light isotopes sputtered preferentially). The value of Delta (O-18:O-16) is always less by about 1.8 deg/o than a linear extrapolation based upon the calculated delta (O-17:O-16) value would suggest. The phenomenon of both negative and positive fractionation patterns from a single target mineral are used to make an experimental test of the proposed model.

Haff, P. K.↗