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 901 records · Page 50

A comparative analysis of two land surface heterogeneity representations

Two contrasting representations of land surface variability used in general circulation models (GCMs) are compared through an analysis of their corresponding surface energy balance equations. In one representation (the 'mixture' approach), different vegetation types are assumed to be homogeneously mixed over a grid square, so that the GCM atmosphere sees near-surface conditions pertaining to the mixture only. In the second representation (the 'mosaic' approach), different vegetation types are viewed as separate 'tiles' of a grid-square 'mosaic', and each tile interacts with the atmosphere independently. The mosaic approach is computationally simpler and in many ways more flexible than the mixture approach. Analytical solutions to the linearized energy balance equations and numerical solutions to the nonlinear equations both demonstrate that the mixture strategy, when applied to two coexisting vegetation types that differ only in canopy transpiration resistance, promotes both total turbulent flux and latent heat flux relative to the mosaic strategy. The effective differences between the strategies, however, are small over a wide range of conditions. In particular, the strategies are effectively equivalent when the transpiration resistances of the different vegetation types are of the same order of magnitude.

Koster, Randal D.↗

3D Higher Order Modeling in the BEM/FEM Hybrid Formulation

Higher order divergence- and curl-conforming bases have been shown to provide significant benefits, in both convergence rate and accuracy, in the 2D hybrid finite element/boundary element formulation (P. Fink and D. Wilton, National Radio Science Meeting, Boulder, CO, Jan. 2000). A critical issue in achieving the potential for accuracy of the approach is the accurate evaluation of all matrix elements. These involve products of high order polynomials and, in some instances, singular Green's functions. In the 2D formulation, the use of a generalized Gaussian quadrature method was found to greatly facilitate the computation and to improve the accuracy of the boundary integral equation self-terms. In this paper, a 3D, hybrid electric field formulation employing higher order bases and higher order elements is presented. The improvements in convergence rate and accuracy, compared to those resulting from lower order modeling, are established. Techniques developed to facilitate the computation of the boundary integral self-terms are also shown to improve the accuracy of these terms. Finally, simple preconditioning techniques are used in conjunction with iterative solution procedures to solve the resulting linear system efficiently. In order to handle the boundary integral singularities in the 3D formulation, the parent element- either a triangle or rectangle-is subdivided into a set of sub-triangles with a common vertex at the singularity. The contribution to the integral from each of the sub-triangles is computed using the Duffy transformation to remove the singularity. This method is shown to greatly facilitate t'pe self-term computation when the bases are of higher order. In addition, the sub-triangles can be further divided to achieve near arbitrary accuracy in the self-term computation. An efficient method for subdividing the parent element is presented. The accuracy obtained using higher order bases is compared to that obtained using lower order bases when the number of unknowns is approximately equal. Also, convergence rates obtained using higher order bases are compared to those obtained with lower order bases for selected sample

Fink, P. W.↗

X-59 Sonic Boom Test Results from the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel

A wind tunnel test was conducted to investigate near-field sonic boom pressure signatures of the X-59 Low-Boom Flight Demonstrator aircraft. A 1.62%-scale model of the aircraft was fabricated for the wind tunnel test, which took place in the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel in September and October 2021. The model had provisions for being mounted by a swept blade strut that attached at top of model ahead of the inlet, or by a rear-entry sting that held the model at the location of the nacelle. The model had alternate parts for ±0.5° deflections of the flaps, ailerons, and stabilator, and ±1° deflections of the T-tail. Off-body static pressure measurements of the flow field below the model were made on a pressure rail which had 420 orifices along its tip. The model was positioned at various heights from the rail by vertical movement of the tunnel strut, and at various longitudinal stations relative to the rail by means of a linear actuator mounted between the wind tunnel strut and the balance. Spatial averaging of model pressure signatures acquired over a range of longitudinal positions reduced the effects of tunnel flow distortions and the interference of the rail flow field and shocks on the model pressure signatures. The test was run at approximate Mach numbers of 1.36, 1.4, and 1.47, and the model was set at various angles of attack and roll relative to the rail. Plots of the model signatures for representative variations of Mach number, model angles, control deflections, and height relative to the rail are provided throughout the report. Repeatability was generally very good and gave confidence in the quality of the measurements. The signatures measured at various heights from the rail provided insight into the aging of the model shocks as they propagated from 1.2 to 3 body lengths from the model. Off-track signatures up to 45° from centerline obtained by rolling the model gave indications of the shock flow fields across the width of the sonic boom carpet. The deflections of the various control surfaces allowed assessment of the boom sensitivity to the control surface movements.

Sonic boom↗

X-59 Sonic Boom Test Results from the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel

A wind tunnel test was conducted to investigate near-field sonic boom pressure signatures of the X-59 Low-Boom Flight Demonstrator aircraft. A 1.62%-scale model of the aircraft was fabricated for the wind tunnel test, which took place in the NASA Glenn 8- by 6-Foot Supersonic Wind Tunnel in September and October 2021. The model had provisions for being mounted by a swept blade strut that attached at top of model ahead of the inlet, or by a rear-entry sting that held the model at the location of the nacelle. The model had alternate parts for ±0.5° deflections of the flaps, ailerons, and stabilator, and ±1° deflections of the T-tail. Off-body static pressure measurements of the flow field below the model were made on a pressure rail which had 420 orifices along its tip. The model was positioned at various heights from the rail by vertical movement of the tunnel strut, and at various longitudinal stations relative to the rail by means of a linear actuator mounted between the wind tunnel strut and the balance. Spatial averaging of model pressure signatures acquired over a range of longitudinal positions reduced the effects of tunnel flow distortions and the interference of the rail flow field and shocks on the model pressure signatures. The test was run at approximate Mach numbers of 1.36, 1.4, and 1.47, and the model was set at various angles of attack and roll relative to the rail. Plots of the model signatures for representative variations of Mach number, model angles, control deflections, and height relative to the rail are provided throughout the report. Repeatability was generally very good and gave confidence in the quality of the measurements. The signatures measured at various heights from the rail provided insight into the aging of the model shocks as they propagated from 1.2 to 3 body lengths from the model. Off-track signatures up to 45° from centerline obtained by rolling the model gave indications of the shock flow fields across the width of the sonic boom carpet. The deflections of the various control surfaces allowed assessment of the boom sensitivity to the control surface movements.

Sonic boom↗

Simulation and control problems in elastic robots

Computational issues associated with modeling and control of robots with revolute joints and elastic arms are considered. A manipulator with one arm and pinned at one end is considered to investigate various aspects of the modeling procedure and the model, and the effect of coupling between the rigid-body and the elastic motions. The rigid-body motion of a manipulator arm is described by means of a reference frame attached to the shadow beam, and the linear elastic operator denoting flexibility is defined with respect to this reference frame. The small elastic motion assumption coupled with the method of assumed modes is used to model the elasticity in the arm. It is shown that only terms up to quadratic in these model amplitudes need to be retained. An important aspect of the coupling between the rigid-body and the elastic motion is the centrifugal stiffening effect. This effect stiffens the elastic structure, as to be expected on physical grounds, gives rise to a time-varying inertia term for the rigid-body motion, and, in general, results in an effective inertia term smaller than the rigid-body inertia term. Simulation results are presented for an elastic beam pinned at one end and free at the other, and rotating in a horizontal plane, and control issues such as the order of the model, number of sensors, and modal extraction are examined within this context.

Tadikonda, S. S. K.↗

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗

Hypersonic vehicle model and control law development using H(infinity) and micron synthesis

The control system design for a Single Stage To Orbit (SSTO) air breathing vehicle will be central to a successful mission because a precise ascent trajectory will preserve narrow payload margins. The air breathing propulsion system requires the vehicle to fly roughly halfway around the Earth through atmospheric turbulence. The turbulence, the high sensitivity of the propulsion system to inlet flow conditions, the relatively large uncertainty of the parameters characterizing the vehicle, and continuous acceleration make the problem especially challenging. Adequate stability margins must be provided without sacrificing payload mass since payload margins are critical. Therefore, a multivariable control theory capable of explicitly including both uncertainty and performance is needed. The H(infinity) controller in general provides good robustness but can result in conservative solutions for practical problems involving structured uncertainty. Structured singular value mu framework for analysis and synthesis is potentially much less conservative and hence more appropriate for problems with tight margins. An SSTO control system requires: highly accurate tracking of velocity and altitude commands while limiting angle-of-attack oscillations, minimized control power usage, and a stabilized vehicle when atmospheric turbulence and system uncertainty are present. The controller designs using H(infinity) and mu-synthesis procedures were compared. An integrated flight/propulsion dynamic mathematical model of a conical accelerator vehicle was linearized as the vehicle accelerated through Mach 8. Vehicle acceleration through the selected flight condition gives rise to parametric variation that was modeled as a structured uncertainty. The mu-analysis approach was used in the frequency domain to conduct controller analysis and was confirmed by time history plots. Results demonstrate the inherent advantages of the mu framework for this class of problems.

Gregory, Irene M.↗

Validity of the Temperature Reconstruction from Water Isotopes in Ice Cores

Well-documented present-day distributions of stable water isotopes (HDO and others) show the existence, in middle and high latitudes, of a linear relationship between the mean annual isotope content of precipitation (SD and 51"0) and the mean annual temperature at the precipitation site. Paleoclimatologists have used this relationship, which is particularly well obeyed over Greenland and Antarctica, to infer paleotemperatures from ice core data. There is, however, growing evidence that spatial and temporal isotope/ surface temperature slopes differ, thus complicating the use of stable water isotopes as paleothermometers. In this paper we review empirical estimates of temporal slopes in polar regions and relevant information that can be inferred from isotope models: simple, Rayleigh-type distillation models and (particularly over Greenland) general circulation models (GCMS) fitted with isotope tracer diagnostics. Empirical estimates of temporal slopes appear consistently lower than present-day spatial slopes and are dependent on the timescale considered. This difference is most probably due to changes in the evaporative origins of moisture, changes in the seasonality of the precipitation, changes in the strength of the inversion layer, or some combination of these changes. Isotope models have not yet been used to evaluate the relative influences of these different factors. The apparent disagreement in the temporal and spatial slopes clearly makes calibrating the isotope paleothermometer difficult. Nevertheless, the use of a (calibrated) isotope paleothermometer appears justified; empirical estimates and most (though not all) GCM results support the practice of interpreting ice core isotope records in terms of local temperature changes.

Jouzel, J.↗

Roll Damping Derivatives from Generalized Lifting-Surface Theory and Wind Tunnel Forced-Oscillation Tests

Improving aerodynamic models for adverse loss-of-control conditions in flight is an area being researched under the NASA Aviation Safety Program. Aerodynamic models appropriate for loss of control conditions require a more general mathematical representation to predict nonlinear unsteady behaviors. As more general aerodynamic models are studied that include nonlinear higher order effects, the possibility of measurements that confound aerodynamic and structural responses are probable. In this study an initial step is taken to look at including structural flexibility in analysis of rigid-body forced-oscillation testing that accounts for dynamic rig, sting and balance flexibility. Because of the significant testing required and associated costs in a general study, it makes sense to capitalize on low cost analytical methods where possible, especially where structural flexibility can be accounted for by a low cost method. This paper provides an initial look at using linear lifting surface theory applied to rigid-body aircraft roll forced-oscillation tests.

Pototzky, Anthony S↗

Modeling for CO poisoning of a fuel cell anode

Poisoning losses in a half-cell in the 110-190 C temperature range have been measured in 100 wt pct H3PO4 for various mixtures of H2, CO, and CO2 gases in order to investigate the polarization loss due to poisoning by CO of a porous fuel cell Pt anode. At a fixed current density, the poisoning loss was found to vary linearly with ln of the CO/H2 concentration ratio, although deviations from linearity were noted at lower temperatures and higher current densities for high CO/H2 concentration ratios. The surface coverages of CO were also found to vary linearly with ln of the CO/H2 concentration ratio. A general adsorption relationship is derived. Standard free energies for CO adsorption were found to vary from -14.5 to -12.1 kcal/mol in the 130-190 C temperature range. The standard entropy for CO adsorption was found to be -39 cal/mol per deg K.

Dhar, H. P.↗

Three-dimensional baroclinic instability of a Hadley cell for small Richardson number

For the case of a baroclinic flow whose Richardson number, Ri, is of order unity, a three-dimensional linear stability analysis is conducted on the basis of a model for a thin, horizontal, rotating fluid layer which is subjected to horizontal and vertical temperature gradients. The Hadley cell basic state and stability analysis are both based on the Navier-Stokes and energy equations, and perturbations possessing zonal, meridional, and vertical structures are considered. An attempt is made to extend the previous theoretical work on three-dimensional baroclinic instability for small Ri to a more realistic model involving the Prandtl and Ekman numbers, as well as to finite growth rates and a wider range of the zonal wavenumber. In general, it is found that the symmetric modes of maximum growth are not purely symmetric, but have a weak zonal structure.

Antar, B. N.↗

Shape design sensitivity analysis and optimization of three dimensional elastic solids using geometric modeling and automatic regridding

An automatic regridding method and a three dimensional shape design parameterization technique were constructed and integrated into a unified theory of shape design sensitivity analysis. An algorithm was developed for general shape design sensitivity analysis of three dimensional eleastic solids. Numerical implementation of this shape design sensitivity analysis method was carried out using the finite element code ANSYS. The unified theory of shape design sensitivity analysis uses the material derivative of continuum mechanics with a design velocity field that represents shape change effects over the structural design. Automatic regridding methods were developed by generating a domain velocity field with boundary displacement method. Shape design parameterization for three dimensional surface design problems was illustrated using a Bezier surface with boundary perturbations that depend linearly on the perturbation of design parameters. A linearization method of optimization, LINRM, was used to obtain optimum shapes. Three examples from different engineering disciplines were investigated to demonstrate the accuracy and versatility of this shape design sensitivity analysis method.

Yao, Tse-Min↗

Numerical Evaluation of Mode 1 Stress Intensity Factor as a Function of Material Orientation For BX-265 Foam Insulation Material

Foam; a cellular material, is found all around us. Bone and cork are examples of biological cell materials. Many forms of man-made foam have found practical applications as insulating materials. NASA uses the BX-265 foam insulation material on the external tank (ET) for the Space Shuttle. This is a type of Spray-on Foam Insulation (SOFI), similar to the material used to insulate attics in residential construction. This foam material is a good insulator and is very lightweight, making it suitable for space applications. Breakup of segments of this foam insulation on the shuttle ET impacting the shuttle thermal protection tiles during liftoff is believed to have caused the space shuttle Columbia failure during re-entry. NASA engineers are very interested in understanding the processes that govern the breakup/fracture of this complex material from the shuttle ET. The foam is anisotropic in nature and the required stress and fracture mechanics analysis must include the effects of the direction dependence on material properties. Material testing at NASA MSFC has indicated that the foam can be modeled as a transversely isotropic material. As a first step toward understanding the fracture mechanics of this material, we present a general theoretical and numerical framework for computing stress intensity factors (SIFs), under mixed-mode loading conditions, taking into account the material anisotropy. We present mode I SIFs for middle tension - M(T) - test specimens, using 3D finite element stress analysis (ANSYS) and FRANC3D fracture analysis software, developed by the Cornel1 Fracture Group. Mode I SIF values are presented for a range of foam material orientations. Also, NASA has recorded the failure load for various M(T) specimens. For a linear analysis, the mode I SIF will scale with the far-field load. This allows us to numerically estimate the mode I fracture toughness for this material. The results represent a quantitative basis for evaluating the strength and fracture properties of anisotropic foam insulation material.

Knudsen, Erik↗

Shape Memory Alloy (SMA)-Based Launch Lock

Most NASA missions require the use of a launch lock for securing moving components during the launch or securing the payload before release. A launch lock is a device used to prevent unwanted motion and secure the controlled components. The current launch locks are based on pyrotechnic, electro mechanically or NiTi driven pin pullers and they are mostly one time use mechanisms that are usually bulky and involve a relatively high mass. Generally, the use of piezoelectric actuation provides high precession nanometer accuracy but it relies on friction to generate displacement. During launch, the generated vibrations can release the normal force between the actuator components allowing shaft's free motion which could result in damage to the actuated structures or instruments. This problem is common to other linear actuators that consist of a ball screw mechanism. The authors are exploring the development of a novel launch lock mechanism that is activated by a shape memory alloy (SMA) material ring, a rigid element and an SMA ring holding flexure. The proposed design and analytical model will be described and discussed in this paper.

launch lock↗

Data Simulation Using VLIDORT and PCRTM and its Use for Radiometric Calibration of Hyperspectral UV-NIR-SWIR Sensors

Using the vector linearized discrete ordinate radiative transfer model (VLIDORT), a multiple-scattering multi-layer discrete ordinate scattering codes, coupled with a fast principal component (PC) radiative transfer model (PCRTM), we can simulate the radiation fields for hyperspectral UV-NIR-SWIR sensors from 250nm to 2500nm with a very high speed. Since the accurate gas absorptions, molecules and aerosol/cloud particle scatterings, and the surface reflectance have been considered in these simulations, these model data can be used for satellite or airborne sensor design, characterization, or end-to-end system level uncertainty assessment. Comparison of the modelled data with measured results has been a general practice for operational satellite sensor calibrations, such as the Ozone Mapping and Profiler Suite (OMPS) onboard on SNPP and NOAA-20. This presentation will give more details of this simulation system and models. Some examples of the generation of simulations data for CLARREO PathFinder (CPF) and the Surface Biology and Geology (SBG), as well as some examples to use these data for CPF and OMPS radiometric calibration will be presented.

Xiaozhen Xiong↗

The simulated Indian monsoon: A GCM sensitivity study

A series of sensitivity experiments are conducted in an attempt to understand and correct deficiencies in the simulation of the seasonal mean Indian monsoon with a global atmospheric general circulation model. The seasonal mean precipitation is less than half that observed. This poor simulation in seasonal integrations is independent of the choice of initial conditions and global sea surface temperature data used. Experiments are performed to test the sensitivity of the Indian monsoon simulation to changes in orography, vegetation, soil, wetness, and cloudiness. The authors find that the deficiency of the model precipitation simulation may be attributed to the use of an enhanced orography in the integrations. Replacement of this orography with a mean orography results in a much more realistic simulation of Indian monsoon circulation and rainfall. Experiments with a linear primitive equation model on the sphere suggest that this striking improvement is due to modulations of the orographically forced waves in the lower troposphere. This improvement in the monsoon simulation is due to the kinematic and dynamical effects of changing the topography, rather than the thermal effects, which were minimal. The magnitude of the impact on the Indian monsoon of the other sensitivity experiments varied considerably, but was consistently less than the impact of using the mean orography. However, results from the soil moisture sensitivity experiments suggest a possibly important role for soil moisture in simulating tropical precipitation, including that associated with the Indian monsoon.

Fennessy, M. J.↗

Equilibrium structure of solar magnetic flux tubes: Energy transport with multistream radiative transfer

We examine the equilibrium structure of vertical intense magnetic flux tubes on the Sun. Assuming cylindrical geometry, we solve the magnetohydrostatic equations in the thin flux-tube approximation, allowing for energy transport by radiation and convection. The radiative transfer equation is solved in the six-stream approximation, assuming gray opacity and local thermodynamic equilibrium. This constitutes a significant improvement over a previous study, in which the transfer was solved using the multidimensional generalization of the Eddington approximation. Convection in the flux tube is treated using mixing-length theory, with an additional parameter alpha, characterizing the suppression of convective energy transport in the tube by the strong magnetic field. The equations are solved using the method of partial linearization. We present results for tubes with different values of the magnetic field strength and radius at a fixed depth in the atmosphere. In general, we find that, at equal geometric heights, the temperature on the tube axis, compared to the ambient medium, is higher in the photosphere and lower in the convection zone, with the difference becoming larger for thicker tubes. At equal optical depths the tubes are generally hotter than their surroundings. The results are comparatively insensitive to alpha but depend upon whether radiative and convective energy transport operate simultaneously or in separate layers. A comparison of our results with semiempirical models shows that the temperature and intensity contrast are in broad agreement. However, the field strengths of the flux-tube models are somewhat lower than the values inferred from observations.

Hasan, S. S.↗

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗