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 847 records · Page 47

Determining the general circulation of the ocean - A nonlinear inverse problem

The usual inversion procedures for determining the oceanic circulation take the density field as given. This does not explicitly account for the mismatch between the data (which are typically a single realization of the density field) and the circulation being estimated (which may be a time mean or some other smoothed quantity), but explicit consideration of the uncertainty of the density data renders the problem nonlinear, and new estimation methods have to be used. As an example, a nonlinear inversion is successfully applied here to the simple analytical model of Stommel and Veronis (1981). It is also shown that their linear procedure, if properly used, yields a much better result than they obtained.

Mercier, H.↗

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Numerical Results - 2

We present numerical results of the diffusion coefficients (DCs) in the coupled diffusion model derived in the preceding paper for a semiconductor quantum well. These include self and mutual DCs in the general two-component case, as well as density- and temperature-related DCs under the single-component approximation. The results are analyzed from the viewpoint of free Fermi gas theory with many-body effects incorporated. We discuss in detail the dependence of these DCs on densities and temperatures in order to identify different roles played by the free carrier contributions including carrier statistics and carrier-LO phonon scattering, and many-body corrections including bandgap renormalization and electron-hole (e-h) scattering. In the general two-component case, it is found that the self- and mutual- diffusion coefficients are determined mainly by the free carrier contributions, but with significant many-body corrections near the critical density. Carrier-LO phonon scattering is dominant at low density, but e-h scattering becomes important in determining their density dependence above the critical electron density. In the single-component case, it is found that many-body effects suppress the density coefficients but enhance the temperature coefficients. The modification is of the order of 10% and reaches a maximum of over 20% for the density coefficients. Overall, temperature elevation enhances the diffusive capability or DCs of carriers linearly, and such an enhancement grows with density. Finally, the complete dataset of various DCs as functions of carrier densities and temperatures provides necessary ingredients for future applications of the model to various spatially inhomogeneous optoelectronic devices.

Li, Jianzhong↗

A model of the human observer in failure detection tasks

A model for the human observer in failure detection tasks is proposed which consists of two stages: a linear estimator and a decision mechanism. The estimator is a Kalman filter, and the decision mechanism, which is based on Wald's sequential analysis, leads to a decision function which is the integration of the filter residuals. The final result is a simple detection system which depends on only three parameters, and the sensitivity of the model to these parameters is analyzed. The results of an experiment designed to test the validity of the model are reported. The question of open and closed decision intervals as well as the generalization of the model to more complicated cases is discussed.

Gai, E. G.↗

The attenuation of sound in a randomly lined duct

The kinetic theory of wave propagation in random media is applied to assess the relative merits of the peripherally segmented, axially segmented and 'checkerboard' configured acoustic duct treatments proposed by Mani (1976). It is argued that, in a first approximation, a multisegmented liner involving many discrete and distinct sections may be modelled as a liner whose impedance varies randomly with position. This results in a relatively simple analytical description of the acoustic properties of the liner in terms of a system of linear, first order ordinary differential equations, whose solution gives the distribution with axial distance of the acoustic modal intensities. The theory is expected to be valid provided that the variation in the impedance from segment to segment is small, so that the conclusions obtained here are of a preliminary nature, and must be regarded as giving a general indication of features which should be incorporated into the design of efficient acoustic liners. The principal such conclusion is that a checkerboard liner offers the best means of increasing the attenuation over and above that of a uniformly lined duct. On the other hand, only minimal gains are to be derived from the use of strictly peripherally segmented liners.

Howe, M. S.↗

Static shape determination and control for large space structures. I - The flexible beam. II - A large space antenna

A method for determining and controlling the shape of large, continuous space structures by means of discrete or pointwise observations and control devices is presented. The general linear boundary value problem satisfied by a one-dimensional shape function is defined, and the existence of solutions is studied. The static shape control problems for one-dimensional systems with and without rigid body modes and the static shape estimation problem are presented and solved. Eigenfunction expansions are presented which provide approximations to the algorithm terms when the associated Green's function is not known. An integral operator approach is applied to the multidimensional static problem, and the results are illustrated with a finite element model of the disk of a large space antenna which assumes no rigid body modes. It is shown that the shape control algorithm must be modified for systems with rigid body modes.

Weeks, C. J.↗

Crack Growth Prediction Methodology for Multi-Site Damage: Layered Analysis and Growth During Plasticity

A finite element program has been developed to perform quasi-static, elastic-plastic crack growth simulations. The model provides a general framework for mixed-mode I/II elastic-plastic fracture analysis using small strain assumptions and plane stress, plane strain, and axisymmetric finite elements. Cracks are modeled explicitly in the mesh. As the cracks propagate, automatic remeshing algorithms delete the mesh local to the crack tip, extend the crack, and build a new mesh around the new tip. State variable mapping algorithms transfer stresses and displacements from the old mesh to the new mesh. The von Mises material model is implemented in the context of a non-linear Newton solution scheme. The fracture criterion is the critical crack tip opening displacement, and crack direction is predicted by the maximum tensile stress criterion at the crack tip. The implementation can accommodate multiple curving and interacting cracks. An additional fracture algorithm based on nodal release can be used to simulate fracture along a horizontal plane of symmetry. A core of plane strain elements can be used with the nodal release algorithm to simulate the triaxial state of stress near the crack tip. Verification and validation studies compare analysis results with experimental data and published three-dimensional analysis results. Fracture predictions using nodal release for compact tension, middle-crack tension, and multi-site damage test specimens produced accurate results for residual strength and link-up loads. Curving crack predictions using remeshing/mapping were compared with experimental data for an Arcan mixed-mode specimen. Loading angles from 0 degrees to 90 degrees were analyzed. The maximum tensile stress criterion was able to predict the crack direction and path for all loading angles in which the material failed in tension. Residual strength was also accurately predicted for these cases.

James, Mark Anthony↗

Modeling the spatiotemporal organization of velocity storage in the vestibuloocular reflex by optokinetic studies

1. A generalized three-dimensional state space model of visual vestibular interaction was developed. Matrix and dynamical system operators associated with inputs from the semicircular canals, otolith velocity estimator, and the visual system have been incorporated into the model, which focus on their relationship to the velocity storage integrator. 2. A relationship was postulated between the eigenvalues and the direction of the eigenvectors of the system matrix and the orientation of the spatial vertical. It was assumed that the system matrix for a tilted position was a composition of two linear transformations of the system matrix for the upright position. One transformation modifies the eigenvalues of the system matrix, whereas another rotates the eigenvectors of the system matrix. The pitch and roll eigenvectors rotate with the head, whereas the yaw axis eigenvector remains approximately spatially invariant. 3. Based on the three-dimensional model, a computational procedure was formulated to identify the eigenvalues and eigenvectors of the system matrix with the use of a modification of the marquardt algorithm. With the use of data obtained from a monkey, it was shown that the three-dimensional behavior of velocity storage cannot be predicted solely in terms of its time constants, i.e., the inverse of its eigenvalues. With the use of the same eigenvalues the data could either be fit or not fit, depending on the eigenvector directions. Therefore, it is necessary to specify eigenvector directions when characterizing velocity storage in three dimensions. 4. Parameters found with the use of the Marquardt algorithm were incorporated into the model. Diagonal matrices in a head coordinate frame were introduced for coupling the visual system to the integrator and to the direct optokinetic pathway. Simulations of optokinetic nystagmus (OKN) and optokinetic after-nystagmus (OKAN) were run. The model predicted the behavior of yaw and pitch OKN and OKAN when the animal is upright. It also predicted the cross-coupling in the side down position. The trajectories in velocity space were also accurately simulated. 5. One of the predictions of the model is that when the stimulus direction is along an eigenvector, the trajectory in velocity space is a straight line. Using the "spectral width" of the residuals from a straight line sequence during OKAN, we developed a methodology to estimate how close the OKAN decay was to an eigenvector trajectory. 6. Thus we have developed a model-based approach for studying and interpreting the response characteristics of velocity storage in three dimensions.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Field Observations of Soil Moisture Variability across Scales

In this study, over 36,000 ground-based soil moisture measurements collected during the SGP97, SGP99, SMEX02, and SMEX03 field campaigns were analyzed to characterize the behavior of soil moisture variability across scales. The field campaigns were conducted in Oklahoma and Iowa in the central USA. The Oklahoma study region is sub-humid with moderately rolling topography, while the Iowa study region is humid with low-relief topography. The relationship of soil moisture standard deviation, skewness and the coefficient of variation versus mean moisture content was explored at six distinct extent scales, ranging from 2.5 m to 50 km. Results showed that variability generally increases with extent scale. The standard deviation increased from 0.036 cm3/cm3 at the 2.5-m scale to 0.071 cm3/cm3 at the 50-km scale. The log standard deviation of soil moisture increased linearly with the log extent scale, from 16 m to 1.6 km, indicative of fractal scaling. The soil moisture standard deviation versus mean moisture content exhibited a convex upward relationship at the 800-m and 50-km scales, with maximum values at mean moisture contents of roughly 0.17 cm3/cm3 and 0.19 cm3/cm3, respectively. An empirical model derived from the observed behavior of soil moisture variability was used to estimate uncertainty in the mean moisture content for a fixed number of samples at the 800-m and 50-km scales, as well as the number of ground-truth samples needed to achieve 0.05 cm3/cm3 and 0.03 cm3/cm3 accuracies. The empirical relationships can also be used to parameterize surface soil moisture variations in land surface and hydrological models across a range of scales. To our knowledge, this is the first study to document the behavior of soil moisture variability over this range of extent scales using ground-based measurements. Our results will contribute not only to efficient and reliable satellite validation, but also to better utilization of remotely sensed soil moisture products for enhanced modeling and prediction.

Famiglietti, James S.↗

Jet Surface Interaction-Scrubbing Noise

Generation of sound due to scrubbing of a jet flow past a nearby solid surface is investigated within the framework of the generalized acoustic analogy theory. The analysis applies to the boundary layer noise generated at and near a wall, and excludes the scattered noise component that is produced at the leading or the trailing edge. While compressibility effects are relatively unimportant at very low Mach numbers, frictional heat generation and thermal gradient normal to the surface could play important roles in generation and propagation of sound in high speed jets of practical interest. A general expression is given for the spectral density of the far-field sound as governed by the variable density Pridmore- Brown equation. The propagation Green's function should be solved numerically starting with the boundary conditions on the surface and subject to specified mean velocity and temperature profiles between the surface and the observer. The equivalent sources of aerodynamic sound are associated with non-linear momentum flux and enthalpy flux terms that appear in the linearized Navier-Stokes equations. These multi-pole sources should be modeled and evaluated with input from a Reynolds-Averaged Navier-Stokes (RANS) solver with an appropriate turbulence model.

Khavaran, Abbas↗

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗

General Rotorcraft Aeromechanical Stability Program (GRASP): Theory manual

The general rotorcraft aeromechanical stability program (GRASP) was developed to calculate aeroelastic stability for rotorcraft in hovering flight, vertical flight, and ground contact conditions. GRASP is described in terms of its capabilities and its philosophy of modeling. The equations of motion that govern the physical system are described, as well as the analytical approximations used to derive them. The equations include the kinematical equation, the element equations, and the constraint equations. In addition, the solution procedures used by GRASP are described. GRASP is capable of treating the nonlinear static and linearized dynamic behavior of structures represented by arbitrary collections of rigid-body and beam elements. These elements may be connected in an arbitrary fashion, and are permitted to have large relative motions. The main limitation of this analysis is that periodic coefficient effects are not treated, restricting rotorcraft flight conditions to hover, axial flight, and ground contact. Instead of following the methods employed in other rotorcraft programs. GRASP is designed to be a hybrid of the finite-element method and the multibody methods used in spacecraft analysis. GRASP differs from traditional finite-element programs by allowing multiple levels of substructure in which the substructures can move and/or rotate relative to others with no small-angle approximations. This capability facilitates the modeling of rotorcraft structures, including the rotating/nonrotating interface and the details of the blade/root kinematics for various types. GRASP differs from traditional multibody programs by considering aeroelastic effects, including inflow dynamics (simple unsteady aerodynamics) and nonlinear aerodynamic coefficients.

Hodges, Dewey H.↗

Dynamic Reconstruction and Multivariable Control for Force-Actuated, Thin Facesheet Adaptive Optics

The Multiple Mirror Telescope (MMT) under development at the University of Arizona takes a new approach in adaptive optics placing a large (0.65 m) force-actuated, thin facesheet deformable mirror at the secondary of an astronomical telescope, thus reducing the effects of emissivity which are important in IR astronomy. However, The large size of the mirror and low stiffness actuators used drive the natural frequencies of the mirror down into the bandwidth of the atmospheric distortion. Conventional adaptive optics takes a quasi-static approach to controlling the, deformable mirror. However, flexibility within the control bandwidth calls for a new approach to adaptive optics. Dynamic influence functions are used to characterize the influence of each actuator on the surface of the deformable mirror. A linearized model of atmospheric distortion is combined with dynamic influence functions to produce a dynamic reconstructor. This dynamic reconstructor is recognized as an optimal control problem. Solving the optimal control problem for a system with hundreds of actuators and sensors is formidable. Exploiting the circularly symmetric geometry of the mirror, and a suitable model of atmospheric distortion, the control problem is divided into a number of smaller decoupled control problems using circulant matrix theory. A hierarchic control scheme which seeks to emulate the quasi-static control approach that is generally used in adaptive optics is compared to the proposed dynamic reconstruction technique. Although dynamic reconstruction requires somewhat more computational power to implement, it achieves better performance with less power usage, and is less sensitive than the hierarchic technique.

Grocott, Simon C. O.↗

Development of the Transport Class Model (TCM) Aircraft Simulation From a Sub-Scale Generic Transport Model (GTM) Simulation

A six degree-of-freedom, flat-earth dynamics, non-linear, and non-proprietary aircraft simulation was developed that is representative of a generic mid-sized twin-jet transport aircraft. The simulation was developed from a non-proprietary, publicly available, subscale twin-jet transport aircraft simulation using scaling relationships and a modified aerodynamic database. The simulation has an extended aerodynamics database with aero data outside the normal transport-operating envelope (large angle-of-attack and sideslip values). The simulation has representative transport aircraft surface actuator models with variable rate-limits and generally fixed position limits. The simulation contains a generic 40,000 lb sea level thrust engine model. The engine model is a first order dynamic model with a variable time constant that changes according to simulation conditions. The simulation provides a means for interfacing a flight control system to use the simulation sensor variables and to command the surface actuators and throttle position of the engine model.

Hueschen, Richard M.↗

Bending of Light in Modified Gravity at Large Distances

We discuss the bending of light in a recent model for gravity at large distances containing a Rindler type acceleration proposed by Grumiller. We consider the static, spherically symmetric metric with cosmological constant and Rindler-like term 2ar presented in this model, and we use the procedure by Rindler and Ishak. to obtain the bending angle of light in this metric. Earlier work on light bending in this model by Carloni, Grumiller, and Preis, using the method normally employed for asymptotically flat space-times, led to a conflicting result (caused by the Rindler-like term in the metric) of a bending angle that increases with the distance of closest approach r(sub 0) of the light ray from the centrally concentrated spherically symmetric matter distribution. However, when using the alternative approach for light bending in nonasymptotically flat space-times, we show that the linear Rindler-like term produces a small correction to the general relativistic result that is inversely proportional to r(sub 0). This will in turn affect the bounds on Rindler acceleration obtained earlier from light bending and casts doubts on the nature of the linear term 2ar in the metric

Grumiller↗

An experimental and theoretical study of the ice accretion process during artificial and natural icing conditions

Real-time measurements of ice growth during artificial and natural icing conditions were conducted using an ultrasonic pulse-echo technique. This technique allows ice thickness to be measured with an accuracy of + or - 0.5 mm; in addition, the ultrasonic signal characteristics may be used to detect the presence of liquid on the ice surface and hence discern wet and dry ice growth behavior. Ice growth was measured on the stagnation line of a cylinder exposed to artificial icing conditions in the NASA Lewis Icing Research Tunnel (IRT), and similarly for a cylinder exposed in flight to natural icing conditions. Ice thickness was observed to increase approximately linearly with exposure time during the initial icing period. The ice accretion rate was found to vary with cloud temperature during wet ice growth, and liquid runback from the stagnation region was inferred. A steady-state energy balance model for the icing surface was used to compare heat transfer characteristics for IRT and natural icing conditions. Ultrasonic measurements of wet and dry ice growth observed in the IRT and in flight were compared with icing regimes predicted by a series of heat transfer coefficients. The heat transfer magnitude was generally inferred to be higher for the IRT than for the natural icing conditions encountered in flight. An apparent variation in the heat transfer magnitude was also observed for flights conducted through different natural icing-cloud formations.

Kirby, Mark S.↗

Passive dynamic controllers for non-linear mechanical systems

A methodology for model-independent controller design for controlling large angular motion of multi-body dynamic systems is outlined.The controlled system may consist of rigid and flexible components that undergo large rigid body motion and small elastic deformations. Control forces/torques are applied to drive the system, and at the same time suppress the vibrations due to flexibility of the components. The proposed controller consists of passive second-order systems which may be designed with little knowledge of the system parameters, even if the controlled system is non-linear. Under rather general assumptions, the passive design assures that the closed loop system has guaranteed stability properties. Unlike positive real controller design, stabilization can be accomplished without direct velocity feedback. In addition, the second-order passive design allows dynamic feedback controllers with considerable freedom to tune for desired system response, and to avoid actuator saturation. After developing the basic mathematical formulation of the design methodology, simulation results are presented to illustrate the proposed approach applied to a flexible six-degree-of-freedom manipulator.

Juang, Jer-Nan↗

Linear regression in astronomy. II

A wide variety of least-squares linear regression procedures used in observational astronomy, particularly investigations of the cosmic distance scale, are presented and discussed. The classes of linear models considered are (1) unweighted regression lines, with bootstrap and jackknife resampling; (2) regression solutions when measurement error, in one or both variables, dominates the scatter; (3) methods to apply a calibration line to new data; (4) truncated regression models, which apply to flux-limited data sets; and (5) censored regression models, which apply when nondetections are present. For the calibration problem we develop two new procedures: a formula for the intercept offset between two parallel data sets, which propagates slope errors from one regression to the other; and a generalization of the Working-Hotelling confidence bands to nonstandard least-squares lines. They can provide improved error analysis for Faber-Jackson, Tully-Fisher, and similar cosmic distance scale relations.

Feigelson, Eric D.↗

High Reynolds number analysis of flat plate and separated afterbody flow using non-linear turbulence models

The ability of the three-dimensional Navier-Stokes method, PAB3D, to simulate the effect of Reynolds number variation using non-linear explicit algebraic Reynolds stress turbulence modeling was assessed. Subsonic flat plate boundary-layer flow parameters such as normalized velocity distributions, local and average skin friction, and shape factor were compared with DNS calculations and classical theory at various local Reynolds numbers up to 180 million. Additionally, surface pressure coefficient distributions and integrated drag predictions on an axisymmetric nozzle afterbody were compared with experimental data from 10 to 130 million Reynolds number. The high Reynolds data was obtained from the NASA Langley 0.3m Transonic Cryogenic Tunnel. There was generally good agreement of surface static pressure coefficients between the CFD and measurement. The change in pressure coefficient distributions with varying Reynolds number was similar to the experimental data trends, though slightly over-predicting the effect. The computational sensitivity of viscous modeling and turbulence modeling are shown. Integrated afterbody pressure drag was typically slightly lower than the experimental data. The change in afterbody pressure drag with Reynolds number was small both experimentally and computationally, even though the shape of the distribution was somewhat modified with Reynolds number.

Carlson, John R.↗