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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

The Use of RGPS Kinematic Data to Estimate Nonlinear Sea Ice Motion

In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motion In current simulations of the interaction between sea ice and its environment, large significance is placed on the deformation of the sea ice. Sea ice deformation is an important process in determining the sea ice thickness distribution across a wide range of space and time scales. Changes in the sea ice thickness distribution affect energy and mass fluxes between the atmosphere and ocean and also the strength of the ice. While most current ice models assume linear variation in the ice motion field to calculate strain, deformation of sea ice occurs through the opening, closing and shearing of ice along discrete linear features. New numerical models are being developed which explicitly account for discontinuities in ice motion, and the need for requisite data sets for model validation has emerged. Multiple buoy data sets, as well as satellite data, have been used to examine the movement and deformation of sea ice. Generally it has been found that the ice motion field has been represented well by buoy data, as well as satellite data over a broad range of scales. However, the underlying deformation (spatial variation in displacement) as represented by different data sets may vary. For the work presented here, sea ice motio

Pruis, M.↗
NASA NTRS1995

An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximation of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudomeasurements are examined, including velocity, density, and volume transport at localized arrays and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform an empirical nudging scheme used in a previous study. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

Fukumori, Ichiro↗
NASA NTRS1973

General atmospheric circulation driven by polar and diurnal surface temperature variations.

Described is a global circulation model for the Venus atmosphere that includes the effects of both polar cooling and diurnal temperature variation. It is based on a linearized Boussinesq approximation and boundary conditions derived from theoretical and empirical considerations. The time-dependent, three-dimensional flow field is deduced without any a priori assumptions about its configuration. Results show that the mean atmospheric motions are essentially zonal in a narrow belt near the equator and change to become meridional over most of the globe. The circulation pattern is not symmetrical and rotates about the polar axis of the planet with the period of the solar day.

Bohachevsky, I. O.↗
NASA NTRS1994

An Approximate Kalman Filter for Ocean Data Assimilation; An Example with an Idealized Gulf Stream Model

A practical method of data assimilation for use with large, nonlinear, ocean general circulation models is explored. A Kalman filter based on approximations of the state error covariance matrix is presented, employing a reduction of the effective model dimension, the error's asymptotic steady-state limit, and a time-invariant linearization of the dynamic model for the error integration. The approximations lead to dramatic computational savings in applying estimation theory to large complex systems. We examine the utility of the approximate filter in assimilating different measurement types using a twin experiment of an idealized Gulf Stream. A nonlinear primitive equation model of an unstable east-west jet is studied with a state dimension exceeding 170,000 elements. Assimilation of various pseudo measurements is examined, including velocity, density, and volume transport at localized arrays, and realistic distributions of satellite altimetry and acoustic tomography observations. Results are compared in terms of their effects on the accuracies of the estimation. The approximate filter is shown to outperform a previous study that used an empirical nudging scheme. The examples demonstrate that useful approximate estimation errors can be computed in a practical manner for general circulation models.

['Kalman↗
NASA NTRS1976

High angle-of-attack stability-and-control analysis

Methods of linear systems analysis were applied to mathematical models of aircraft flying at high angle of attack and maneuver rate. First order longitudinal and lateral directional coupling is obtained by linearizing the complete nonlinear equations of motion about a generalized (quasi steady) trim point. Open loop stability boundaries are defined using the linear dynamic equations, and pilot in the loop effects are presented. Stability augmentation structures for maneuvering flight conditions are shown to be defined readily using optimal control theory.

Stengel, R. F.↗
NASA NTRS1990

Physical interpretation and application of principles of ultrasonic nondestructive evaluation of high-performance materials

An ultrasonic measurement system employed in the experimental interrogation of the anisotropic properties (through the measurement of the elastic stiffness constants) of the uniaxial graphite-epoxy composites is presented. The continuing effort for the development of improved visualization techniques for physical parameters is discussed. The background is set for the understanding and visualization of the relationship between the phase and energy/group velocity for propagation in high-performance anisotropic materials by investigating the general requirements imposed by the classical wave equation. The consequences are considered when the physical parameters of the anisotropic material are inserted into the classical wave equation by a linear elastic model. The relationship is described between the phase velocity and the energy/group velocity three dimensional surfaces through graphical techniques.

Miller, James G.↗
NASA NTRS1993

Multiple degree of freedom force-state component identification

Force-state mapping is extended to allow the quasistatic characterization of realistic multiple degree of freedom systems whose constitutive relations are nonlinear, dissipative, coupled and depend on memory of past states. A general constitutive force-state model is formulated, which includes dynamic hysteresis phenomena typical of that found in deployable structures. A two step parameter identification approach is developed, which employs a linear extended Kalman filter and a nonlinear least squares fit. A six degree of freedom force-state testing device was constructed, and used to obtain data on one bay of a simple truss and on one bay of the MODE deployable truss.

Masters, Brett P.↗
NASA NTRS1999

A Multipole Expansion Method for Analyzing Lightning Field Changes

Changes in the surface electric field are frequently used to infer the locations and magnitudes of lightning-caused changes in thundercloud charge distributions. The traditional procedure is to assume that the charges that are effectively deposited by the flash can be modeled either as a single point charge (the Q model) or a point dipole (the P model). The Q model has four unknown parameters and provides a good description of many cloud-to-ground (CG) flashes. The P model has six unknown parameters and describes many intracloud (IC) discharges. In this paper we introduce a new analysis method that assumes that the change in the cloud charge can be described by a truncated multipole expansion, i.e., there are both monopole and dipole terms in the unknown source distribution, and both terms are applied simultaneously. This method can be used to analyze CG flashes that are accompanied by large changes in the cloud dipole moment and complex IC discharges. If there is enough information content in the measurements, the model can also be generalized to include quadrupole and higher order terms. The parameters of the charge moments are determined using a dme-dimensional grid search in combination with a linear inversion, and because of this, local minima in the error function and the associated solution ambiguities are avoided. The multipole method has been tested on computer-simulated sources and on natural lightning at the NASA Kennedy Space Center and U.S. Air Force Eastern Range.

Koshak, William J.↗
NASA NTRS2002

Reply by the Authors to C. K. W. Tam

The prediction of noise generation and radiation by turbulence has been the subject of continuous research for over fifty years. The essential problem is how to model the noise sources when one s knowledge of the detailed space-time properties of the turbulence is limited. We attempted to provide a comparison of models based on acoustic analogies and recent alternative models. Our goal was to demonstrate that the predictive capabilities of any model are based on the choice of the turbulence property that is modeled as a source of noise. Our general definition of an acoustic analogy is a rearrangement of the equations of motion into the form L(u) = Q, where L is a linear operator that reduces to an acoustic propagation operator outside a region upsilon; u is a variable that reduces to acoustic pressure (or a related linear acoustic variable) outside upsilon; and Q is a source term that can be meaningfully estimated without knowing u and tends to zero outside upsilon.

Morris, Philip J.↗
NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗
NASA NTRS2014

Assimilation of MODIS Snow Cover Through the Data Assimilation Research Testbed and the Community Land Model Version 4

To improve snowpack estimates in Community Land Model version 4 (CLM4), the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) was assimilated into the Community Land Model version 4 (CLM4) via the Data Assimilation Research Testbed (DART). The interface between CLM4 and DART is a flexible, extensible approach to land surface data assimilation. This data assimilation system has a large ensemble (80-member) atmospheric forcing that facilitates ensemble-based land data assimilation. We use 40 randomly chosen forcing members to drive 40 CLM members as a compromise between computational cost and the data assimilation performance. The localization distance, a parameter in DART, was tuned to optimize the data assimilation performance at the global scale. Snow water equivalent (SWE) and snow depth are adjusted via the ensemble adjustment Kalman filter, particularly in regions with large SCF variability. The root-mean-square error of the forecast SCF against MODIS SCF is largely reduced. In DJF (December-January-February), the discrepancy between MODIS and CLM4 is broadly ameliorated in the lower-middle latitudes (2345N). Only minimal modifications are made in the higher-middle (4566N) and high latitudes, part of which is due to the agreement between model and observation when snow cover is nearly 100. In some regions it also reveals that CLM4-modeled snow cover lacks heterogeneous features compared to MODIS. In MAM (March-April-May), adjustments to snowmove poleward mainly due to the northward movement of the snowline (i.e., where largest SCF uncertainty is and SCF assimilation has the greatest impact). The effectiveness of data assimilation also varies with vegetation types, with mixed performance over forest regions and consistently good performance over grass, which can partly be explained by the linearity of the relationship between SCF and SWE in the model ensembles. The updated snow depth was compared to the Canadian Meteorological Center (CMC) data. Differences between CMC and CLM4 are generally reduced in densely monitored regions.

data assimilation↗
NASA NTRS1990

Reynolds stress closure in jet flows using wave models

A collection of papers is presented. The outline of this report is as follows. Chapter three contains a description of a weakly nonlinear turbulence model that was developed. An essential part of the application of such a closure scheme to general geometry jets is the solution of the local hydrodynamic stability equation for a given jet cross-section. Chapter four describes the conformal mapping schemes used to map such geometries onto a simple computational domain. Chapter five describes a solution of a stability problem for circular, elliptic, and rectangular geometries. In chapter six linear models for the shock shell structure in non-circular jets is given. The appendices contain reprints of papers also published during this study including the following topics: (1) instability of elliptic jets; (2) a technique for predicting the shock cell structure in non-circular jets using a vortex sheet model; and (3) the resonant interaction between twin supersonic jets.

Morris, Philip J.↗
NASA NTRS1986

Near millimeter wave imaging/multi-beam integrated antennas

Some preliminary results on a mixer design which is suitable for integration with tapered slot antennas have been obtained and published. This mixer design was tested both in a 4 to 10 GHz model, and (slightly modified) at 94 GHz. The latter utilized the same Hewlett-Packard beam-lead diodes which were used as detector diodes in the linearly tapered slot antennas (LTSA) arrays. These diodes are the most rugged to be found, and generally survive well on the flexible Kapton substrates. The 4 to 10 GHz version of this mixer has less than 6 dB conversion loss over an octave bandwidth. It uses a slot ring in a balanced configuration, and requires the LO to be fed through a separate port. A different design for a mixer which may be integrated with an LTSA antenna element is discussed. This mixer was tested at 38 GHz with the same HP beam-lead diodes, and has less than 10 dB conversion loss. Further work on mixers has emphasized theoretical modeling, using a computer program, which takes into account the effect of excess noise of Schottky-barrier diodes for the first time. Calculated results agree quantitatively with measured results on millimeter wave mixers.

Yngvesson, K. Sigfrid↗
NASA NTRS1991

On the formulation of a minimal uncertainty model for robust control with structured uncertainty

In the design and analysis of robust control systems for uncertain plants, representing the system transfer matrix in the form of what has come to be termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents a transfer function matrix M(s) of the nominal closed loop system, and the delta represents an uncertainty matrix acting on M(s). The nominal closed loop system M(s) results from closing the feedback control system, K(s), around a nominal plant interconnection structure P(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unsaturated uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, but for real parameter variations delta is a diagonal matrix of real elements. Conceptually, the M-delta structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the currently available literature addresses computational methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty, where the term minimal refers to the dimension of the delta matrix. Since having a minimally dimensioned delta matrix would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta would be useful. Hence, a method of obtaining the interconnection system P(s) is required. A generalized procedure for obtaining a minimal P-delta structure for systems with real parameter variations is presented. Using this model, the minimal M-delta model can then be easily obtained by closing the feedback loop. The procedure involves representing the system in a cascade-form state-space realization, determining the minimal uncertainty matrix, delta, and constructing the state-space representation of P(s). Three examples are presented to illustrate the procedure.

Belcastro, Christine M.↗
NASA NTRS1976

Thermal analyses of the International Ultraviolet Explorer (IUE) scientific instrument using the NASTRAN thermal analyzer (NTA): A general purpose summary

The NTA Level 15.5.2/3, was used to provide non-linear steady-state (NLSS) and non-linear transient (NLTR) thermal predictions for the International Ultraviolet Explorer (IUE) Scientific Instrument (SI). NASTRAN structural models were used as the basis for the thermal models, which were produced by a straight forward conversion procedure. The accuracy of this technique was sub-sequently demonstrated by a comparison of NTA predicts with the results of a thermal vacuum test of the IUE Engineering Test Unit (ETU). Completion of these tasks was aided by the use of NTA subroutines.

Jackson, C. E., Jr.↗
NASA NTRS2003

Transient Simulation of the Integrated Powerhead Demonstrator (IPD) Rocket Engine

The IPD rocket engine is the product of a joint Air Force/NASA program to demonstrate the concept of a full-flow staged combustion power cycle. In this type of rocket engine cycle, both the entire fuel flow and the entire oxidizer flow are combusted in one of two pre-bumers and used to drive two turbopumps, thus utilizing the flow work of the total propellant flow. The basis of the IPD engine program relies on integrating newly developed and pre-existing hardware to demonstrate the component and material technologies to make this concept feasible, while simultaneously saving development time and costs. To provide insight to the project team and contractors during engine development and test phases, the engine system was modeled at MSFC using the Rocket Engine Transient Simulation (ROCETS) software to analyze system performance and determine component integration issues. The ROCETS software is used extensively a MSFC to perform steady-state power-balances and transient simulations of thermodynamic power and general fluid systems. The software is favored for its capability to solve large systems of non-linear equations, its librarie of fluid properties and flow devices, its flexibility to modify existing code to improve the physics-derived approximations of real fluid thermodynamic behavior, and the ability to add unique system constraints. The purpose of this paper is to present the methodology used to model the IPD engine system, detail the pitfalls encountered with the software, and explain the approximations made to more accurately represent engine component and fluid behavior. Engine system performance output from the model will be presented and explained in comparison with real fluid and component behavior.

Guidos, Mike↗
NASA NTRS2020

The HadGEM3-GA7.1 radiative kernel: the importance of a well-resolved stratosphere

We present top-of-atmosphere and surface radiative kernels based on the atmospheric component (GA7.1) of the HadGEM3 general circulation model developed by the UK Met Office. We show that the utility of radiative kernels for forcing adjustments in idealised CO2 perturbation experiments is greatest where there is sufficiently high resolution in the stratosphere in both the target climate model and the radiative kernel. This is because stratospheric cooling to a CO2 perturbation continues to increase with height, and low-resolution or low-top kernels or climate model output are unable to fully resolve the full stratospheric temperature adjustment. In the sixth phase of the Coupled Model Intercomparison Project (CMIP6), standard atmospheric model data are available up to 1 hPa on 19 pressure levels, which is a substantial advantage compared to CMIP5. We show in the IPSL-CM6A-LR model where a full set of climate diagnostics are available that the HadGEM3-GA7.1 kernel exhibits linear behaviour and the residual error term is small, as well as from a survey of kernels available in the literature that in general low-top radiative kernels underestimate the stratospheric temperature response. The HadGEM3-GA7.1 radiative kernels are available at https://doi.org/10.5281/zenodo.3594673 (Smith, 2019).

radiative kernels↗
NASA NTRS1976

Aspects of job scheduling

A mathematical model for job scheduling in a specified context is presented. The model uses both linear programming and combinatorial methods. While designed with a view toward optimization of scheduling of facility and plant operations at the Deep Space Communications Complex, the context is sufficiently general to be widely applicable. The general scheduling problem including options for scheduling objectives is discussed and fundamental parameters identified. Mathematical algorithms for partitioning problems germane to scheduling are presented.

Phillips, K.↗