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 577 records · Page 32

Parallel processing for scientific computations

The scope of this project dealt with the investigation of the requirements to support distributed computing of scientific computations over a cluster of cooperative workstations. Various experiments on computations for the solution of simultaneous linear equations were performed in the early phase of the project to gain experience in the general nature and requirements of scientific applications. A specification of a distributed integrated computing environment, DICE, based on a distributed shared memory communication paradigm has been developed and evaluated. The distributed shared memory model facilitates porting existing parallel algorithms that have been designed for shared memory multiprocessor systems to the new environment. The potential of this new environment is to provide supercomputing capability through the utilization of the aggregate power of workstations cooperating in a cluster interconnected via a local area network. Workstations, generally, do not have the computing power to tackle complex scientific applications, making them primarily useful for visualization, data reduction, and filtering as far as complex scientific applications are concerned. There is a tremendous amount of computing power that is left unused in a network of workstations. Very often a workstation is simply sitting idle on a desk. A set of tools can be developed to take advantage of this potential computing power to create a platform suitable for large scientific computations. The integration of several workstations into a logical cluster of distributed, cooperative, computing stations presents an alternative to shared memory multiprocessor systems. In this project we designed and evaluated such a system.

Alkhatib, Hasan S.↗

Computational Aeroelasticity: Success, Progress, Challenge

The formal term Computational Aeroelasticity (CAE) has only been recently adopted to describe aeroelastic analysis methods coupling high-level computational fluid dynamics codes with structural dynamics techniques. However, the general field of aeroelastic computations has enjoyed a rich history of development and application since the first hand-calculations performed in the mid 1930 s. This paper portrays a much broader definition of Computational Aeroelasticity; one that encompasses all levels of aeroelastic computation from the simplest linear aerodynamic modeling to the highest levels of viscous unsteady aerodynamics, from the most basic linear beam structural models to state-of-the-art Finite Element Model (FEM) structural analysis. This paper is not written as a comprehensive history of CAE, but rather serves to review the development and application of aeroelastic analysis methods. It describes techniques and example applications that are viewed as relatively mature and accepted, the "successes" of CAE. Cases where CAE has been successfully applied to unique or emerging problems, but the resulting techniques have proven to be one-of-a-kind analyses or areas where the techniques have yet to evolve into a routinely applied methodology are covered as "progress" in CAE. Finally the true value of this paper is rooted in the description of problems where CAE falls short in its ability to provide relevant tools for industry, the so-called "challenges" to CAE.

Schuster, David M.↗

Real-time Adaptive Control Using Neural Generalized Predictive Control

The objective of this paper is to demonstrate the feasibility of a Nonlinear Generalized Predictive Control algorithm by showing real-time adaptive control on a plant with relatively fast time-constants. Generalized Predictive Control has classically been used in process control where linear control laws were formulated for plants with relatively slow time-constants. The plant of interest for this paper is a magnetic levitation device that is nonlinear and open-loop unstable. In this application, the reference model of the plant is a neural network that has an embedded nominal linear model in the network weights. The control based on the linear model provides initial stability at the beginning of network training. In using a neural network the control laws are nonlinear and online adaptation of the model is possible to capture unmodeled or time-varying dynamics. Newton-Raphson is the minimization algorithm. Newton-Raphson requires the calculation of the Hessian, but even with this computational expense the low iteration rate make this a viable algorithm for real-time control.

Haley, Pam↗

Shape and Stress Sensing of Multilayered Composite and Sandwich Structures Using an Inverse Finite Element Method

The marked increase in the use of composite and sandwich material systems in aerospace, civil, and marine structures leads to the need for integrated Structural Health Management systems. A key capability to enable such systems is the real-time reconstruction of structural deformations, stresses, and failure criteria that are inferred from in-situ, discrete-location strain measurements. This technology is commonly referred to as shape- and stress-sensing. Presented herein is a computationally efficient shape- and stress-sensing methodology that is ideally suited for applications to laminated composite and sandwich structures. The new approach employs the inverse Finite Element Method (iFEM) as a general framework and the Refined Zigzag Theory (RZT) as the underlying plate theory. A three-node inverse plate finite element is formulated. The element formulation enables robust and efficient modeling of plate structures instrumented with strain sensors that have arbitrary positions. The methodology leads to a set of linear algebraic equations that are solved efficiently for the unknown nodal displacements. These displacements are then used at the finite element level to compute full-field strains, stresses, and failure criteria that are in turn used to assess structural integrity. Numerical results for multilayered, highly heterogeneous laminates demonstrate the unique capability of this new formulation for shape- and stress-sensing.

Cerracchio, Priscilla↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic (LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖, 𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, there are some instances where the fitting needs to be done programmatically. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic(LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖,𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, it is often convenient to do the fitting programmatically, especially when a large number of datasets should be analyzed. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

A Quasi-Linear Behavioral Model and an Application to Self-Directed Learning

A model is presented that describes the relationship between one's knowledge of the world and the concomitant personal behaviors that serve as a mechanism to obtain desired outcomes. Integrated within this model are the differing roles that outcomes serve as motivators and as modifiers to one's worldview. The model is dichotomized between general and contextual applications. Because learner self-directedness (a personal characteristic) involves cognition and affection while self-directed learning (a pedagogic process) encompasses conation, behavior and introspection, the model can be dichotomized again in another direction. Presented also are the roles that cognitive motivation theories play in moving an individual through this behavioral model and the roles of wishes, self-efficacy, opportunity and self-influence.

Ponton, Michael K.↗

Hydrostatic calculations of axisymmetric flow and its stability for the AGCE model

Baroclinic waves in the atmospherics general circulation experiment (AGCE) apparatus by the use of numerical hydrostatic primitive equation models were determined. The calculation is accomplished by using an axisymmetric primitive equation model to compute, for a given set of experimental parameters, a steady state axisymmetric flow and then testing this axisymmetric flow for stability using a linear primitive equation model. Some axisymmetric flows are presented together with preliminary stability calculations.

Miller, T. L.↗

Aircraft control surface failure detection and isolation using the OSGLR test

The performance of the orthogonal series generalized likelihood ratio (OSGLR) test in detecting and isolating commercial aircraft control surface and actuator failures is evaluated. A modification to incorporate age-weighting which significantly reduces the sensitivity of the algorithm to modeling errors is presented. The steady-state implementation of the algorithm based on a single linear model valid for a cruise flight condition is tested using a nonlinear aircraft simulation. A number of off-nominal no-failure flight conditions including maneuvers, nonzero flap deflections, different turbulence levels and steady winds were tested. Based on the no-failure decision functions produced by off-nominal flight conditions, the failure detection and isolation performance at the nominal flight condition was determined. The extension of the algorithm to a wider flight envelope by scheduling on dynamic pressure and flap deflection is examined. Based on this testing, the OSGLR algorithm should be capable of detecting control surface failures that would affect the safe operation of a commercial aircraft. Isolation may be difficult if there are several surfaces which produce similar effects on the aircraft. Extending the algorithm over the entire operating envelope of a commercial aircraft appears feasible.

Bonnice, W. F.↗

Unsteady aerodynamics of blade rows

The requirements placed on an unsteady aerodynamic theory intended for turbomachinery aeroelastic applications are discussed along with a brief description of the various theoretical models that are available to address these requirements. The main emphasis is placed on the description of a linearized inviscid theory which fully accounts for the effects of a nonuniform mean or steady flow on unsteady aerodynamic response. Although this theory has been developed primarily for blade flutter prediction, more general equations are presented which account for unsteady excitations due to incident external aerodynamic disturbances as well as those due to blade motions. The resulting equations consist of a system of three field equations along with conditions imposed at blade, wake and shock surfaces and in the far field. These equations can be solved to determine the fluctuations in all fluid dynamic properties throughout the required solution domain. Example solutions are presented to demonstrate several effects associated with nonuniform steady flows on the linearized unsteady aerodynamic response to prescribed blade motions.

Verdon, J. M.↗

Steady induction effects in geomagnetism. Part 1B: Geomagnetic estimation of steady surficial core motions: A non-linear inverse problem

The problem of estimating a steady fluid velocity field near the top of Earth's core which induces the secular variation (SV) indicated by models of the observed geomagnetic field is examined in the source-free mantle/frozen-flux core (SFI/VFFC) approximation. This inverse problem is non-linear because solutions of the forward problem are deterministically chaotic. The SFM/FFC approximation is inexact, and neither the models nor the observations they represent are either complete or perfect. A method is developed for solving the non-linear inverse motional induction problem posed by the hypothesis of (piecewise, statistically) steady core surface flow and the supposition of a complete initial geomagnetic condition. The method features iterative solution of the weighted, linearized least-squares problem and admits optional biases favoring surficially geostrophic flow and/or spatially simple flow. Two types of weights are advanced radial field weights for fitting the evolution of the broad-scale portion of the radial field component near Earth's surface implied by the models, and generalized weights for fitting the evolution of the broad-scale portion of the scalar potential specified by the models.

Voorhies, Coerte V.↗

Symmetry properties in polarimetric remote sensing

This paper presents the relations among polarimetric backscattering coefficients from the viewpoint of symmetry groups. Symmetry of geophysical media encountered in remote sensing due to reflection, rotation, azimuthal, and centrical symmetry groups is considered for both reciprocal and nonreciprocal cases. On the basis of the invariance under symmetry transformations in the linear polarization basis, the scattering coefficients are related by a set of equations which restrict the number of independent parameters in the polarimetric covariance matrix. The properties derived under these transformations are general and apply to all scattering mechanisms in a given symmetrical configuration. The scattering coefficients calculated from theoretical models for layer random media and rough surfaces are shown to obey the derived symmetry relations. Use of symmetry properties in remote sensing of structural and environmental responses of scattering media is discussed. As a practical application, the results from this paper provide new methods for the external calibration of polarimetric radars without the deployment of man-made calibration targets.

Nghiem, S. V.↗

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.

Samtaney, Ravi↗

Modifications and Modeling of the Fission Surface Power Primary Test Circuit (FSP-PTC)

An actively pumped alkali metal flow circuit, designed and fabricated at the NASA Marshall Space Flight Center, underwent a range of tests at MSFC in early 2007. During this period, system transient responses and the performance of the liquid metal pump were evaluated. In May of2007, the circuit was drained and cleaned to prepare for multiple modifications: the addition of larger upper and lower reservoirs, the installation of an annular linear induction pump (ALIP), and the inclusion of a closeable orifice in the test section. Performance of the ALIP, provided by Idaho National Laboratory (INL), will be evaluated when testing resumes. Data from the first round of testing has been used to refine the working system model, developed using the Generalized Fluid System Simulation Program (GFSSP). This paper covers the modifications of the FSP-PTC and the updated GFSSP system model.

Garber, Anne E.↗

High-Fidelity Generalization Method of Cells for Inelastic Periodic Multiphase Materials

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress said inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model, and with finite-element results for transverse loading.

Aboudi, Jacob↗

On the applicability of integrated circuit technology to general aviation orientation estimation

The criteria of the significant value of the panel instruments used in general aviation were examined and kinematic equations were added for comparison. An instrument survey was performed to establish the present state of the art in linear and angular accelerometers, pressure transducers, and magnetometers. A very preliminary evaluation was done of the computers available for data evaluation and estimator mechanization. The mathematical model of a light twin aircraft employed in the evaluation was documented, the results of the sensor survey and the results of the design studies were presented.

Debra, D. B.↗

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.↗

Zircon Lu-Hf systematics: Evidence for the episodic development of Archaean greenstone belts

A combined U-Th-Pb and Lu-Hf isotopic study of zircons was undertaken in order to determine the provenance and age of an Archean granite-greenstone terrain and to test the detailed application of the Lu-Hf system in various Archean zircons. The eastern Wawa subprovince of the Superior province consists of the low grade Michipicoten and Gamitagama greenstone belts and the granitic terrain. The Hf isotopic data indicate that the typical lithological features of a greenstone belt cycle could be accommodated in a crustal growth model that involved decreasing depth of melting in three isotopically distinct reservoirs: mantle, lower crust and upper crust. The model age of the sources given by the intersection of the lower crustal curve with the bulk earth evolution curve is about 2900 My, in good agreement with the zircon U-Pb basement age. This linear array also has a similar intersection age to that of Proterozoic carbonatite complexes. The general convergence of the other reservoir vectors around this age suggests that mantle depletion, crustal extraction and intracrustal differentiation were all part of the same episodic event. It is also apparent that recycling of older basement was important in the formation of many of the later greenstone belt rocks.

Smith, P. E.↗