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 685 records · Page 38

A discourse on sensitivity analysis for discretely-modeled structures

A descriptive review is presented of the most recent methods for performing sensitivity analysis of the structural behavior of discretely-modeled systems. The methods are generally but not exclusively aimed at finite element modeled structures. Topics included are: selections of finite difference step sizes; special consideration for finite difference sensitivity of iteratively-solved response problems; first and second derivatives of static structural response; sensitivity of stresses; nonlinear static response sensitivity; eigenvalue and eigenvector sensitivities for both distinct and repeated eigenvalues; and sensitivity of transient response for both linear and nonlinear structural response.

Adelman, Howard M.↗

Momentum-Based Dynamics for Spacecraft with Chained Revolute Appendages

An efficient formulation is presented for a sub-class of multi-body dynamics problems that involve a six degree-of-freedom base body and a chain of N rigid linkages connected in series by single degree-of-freedom revolute joints. This general method is particularly well suited for simulations of spacecraft dynamics and control that include the modeling of an orbiting platform with or without internal degrees of freedom such as reaction wheels, dampers, and/or booms. In the present work, particular emphasis is placed on dynamic simulation of multi-linkage robotic manipulators. The differential equations of motion are explicitly given in terms of linear and angular momentum states, which can be evaluated recursively along a serial chain of linkages for an efficient real-time solution on par with the best of the O(N3) methods.

Queen, Steven↗

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.

Huang, T. C.↗

Recursive partitioned inversion of large (1500 x 1500) symmetric matrices

A recursive algorithm was designed to invert large, dense, symmetric, positive definite matrices using small amounts of computer core, i.e., a small fraction of the core needed to store the complete matrix. The described algorithm is a generalized Gaussian elimination technique. Other algorithms are also discussed for the Cholesky decomposition and step inversion techniques. The purpose of the inversion algorithm is to solve large linear systems of normal equations generated by working geodetic problems. The algorithm was incorporated into a computer program called SOLVE. In the past the SOLVE program has been used in obtaining solutions published as the Goddard earth models.

Putney, B. H.↗

An analysis of the two-dimensional acoustic field in a nonuniform duct carrying compressible flow

An analytical/numerical linear acoustic solution in a nonuniform two-dimensional duct carrying a compressible mean flow is developed. A quasi-one dimensional mean flow model is employed together with a consistent expression for the cross-flow velocity. The acoustic solution is obtained using the Wave Envelope Method. Numerical results are compared with those of an existing Wave Envelope solution which includes a more general mean flow model. In addition, the singular behavior of the solutions derived in the current work for nearly sonic flow conditions is examined explicitly.

Uenishi, K.↗

Computer studies of baroclinic flow

Programs necessary for computing the transition curve on the regime diagram for the atmospheric general circulation experiment (AGOE) were completed and used to determine the regime diagram for the rotating annulus and some axisymmetric flows for one possible AGOE configuration. The effect of geometrical constraints on the size of eddies developing from a basic state is being examined. In AGOE, the geometric constraint should be the width of the shear zone or the baroclinic zone. Linear and nonlinear models are to be used to examine both barotropic and baroclinic flows. The results should help explain the scale selection mechanism of baroclinic eddies in the atmosphere experimental models such as AGOE, and the multiple vortex phenomenon in tornadoes.

Gall, R.↗

Modelling of photodegradation in solar cell modules of substrate and superstrate design made with ethylene-vinyl acetate as pottant material

A computer model was developed which can generate realistic concentration versus time profiles of the chemical species formed during photooxidation of hydrocarbon polymers using as input data a set of elementary reactions with corresponding rate constants and initial conditions. The results of computer simulation have been shown to be consistent with the general experimental observations of the photooxidation of polyethylene exposed to sunlight at ambient temperatures. The useful lifetime (5% oxidation) of the unstabilized polyethylene is predicted to vary from a few months in hot weather (100 F) to almost two years in cool weather (45 F) with an apparent net activation energy of 10 kcal/mol. Modelling studies of alternate mechanisms for stabilization of clear, amorphous, linear polyethylene suggest that the optimum stabilizer would be a molecularly dispensed additive in very low concentration which can trap peroxy radicals and also decompose hydroperoxides.

Somersall, A. C.↗

A parametric heat transfer study for cryogenic ball bearings in SSME HPOTP

A numerical modeling is to examine the effects of coolant convective heat transfer coefficient and frictional heating on the local temperature characteristics of a ball element in Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump (HPOTP) bearing. The present modeling uses a control-volume based, finite-difference method to solve the non-dimensionalized heat conduction equation in spherical coordinate system. The dimensionless temperature is found as a function of Biot number, heat flux ratio between the two race contacts, and location in the ball. The current results show that, for a given cooling capability, the ball temperature generally increases almost linearly with the heat input from the race-contacts. This increase is always very high at one of the two contacts. An increase in heat transfer coefficient generally reduces the ball temperature and alleviates the temperature gradient, except for the regions very close to the race contacts. For a 10-fold increase of heat transfer coefficient, temperature decrease is 35 percent for the average over entire ball, and 10 percent at the inner-race contact. The corresponding change of temperature gradient displays opposing trends between the regions immediately adjacent to the contacts and the remaining portion of the ball. The average temperature gradient in the vicinity of both contacts increases approximately 70 to 100 percent. A higher temperature gradient produces excessive thermal stress locally which may be detrimental to the material integrity. This, however, is the only unfavorable issue for an increase of heat transfer coefficient.

Chyu, Mingking K.↗

Laboratory Measurements of Solar-Wind/Comet X-Ray Emission and Charge Exchange Cross Sections

The detection of X-rays from comets such as Hyakutake, Hale-Bopp, d Arrest, and Linear as they approach the Sun has been unexpected and exciting. This phenomenon, moreover, should be quite general, occurring wherever a fast solar or stellar wind interacts with neutrals in a comet, a planetary atmosphere, or a circumstellar cloud. The process is, O(+8) + H2O --> O(+7*) + H2O(+), where the excited O(+7*) ions are the source of the X-ray emissions. Detailed modeling has been carried out of X-ray emissions in charge-transfer collisions of heavy solar-wind Highly Charged Ions (HCIs) and interstellar/interplanetary neutral clouds. In the interplanetary medium the solar wind ions, including protons, can charge exchange with interstellar H and He. This can give rise to a soft X-ray background that could be correlated with the long-term enhancements seen in the low-energy X-ray spectrum of ROSAT. Approximately 40% of the soft X-ray background detected by Exosat, ROSAT, Chandra, etc. is due to Charge Exchange (CXE): our whole heliosphere is glowing in the soft X-ray due to CXE.

Chutjian, A.↗

Damage to cellular DNA from particulate radiations, the efficacy of its processing and the radiosensitivity of mammalian cells. Emphasis on DNA double strand breaks and chromatin breaks

For several years, it has been evident that cellular radiation biology is in a necessary period of consolidation and transition (Lett 1987, 1990; Lett et al. 1986, 1987). Both changes are moving apace, and have been stimulated by studies with heavy charged particles. From the standpoint of radiation chemistry, there is now a consensus of opinion that the DNA hydration shell must be distinguished from bulk water in the cell nucleus and treated as an integral part of DNA (chromatin) (Lett 1987). Concomitantly, sentiment is strengthening for the abandonment of the classical notions of "direct" and "indirect" action (Fielden and O'Neill 1991; O'Neill 1991; O'Neill et al. 1991; Schulte-Frohlinde and Bothe 1991 and references therein). A layer of water molecules outside, or in the outer edge of, the DNA (chromatin) hydration shell influences cellular radiosensitivity in ways not fully understood. Charge and energy transfer processes facilitated by, or involving, DNA hydration must be considered in rigorous theories of radiation action on cells. The induction and processing of double stand breaks (DSBs) in DNA (chromatin) seem to be the predominant determinants of the radiotoxicity of normally radioresistant mammalian cells, the survival curves of which reflect the patterns of damage induced and the damage present after processing ceases, and can be modelled in formal terms by the use of reaction (enzyme) kinetics. Incongruities such as sublethal damage are neither scientifically sound nor relevant to cellular radiation biology (Calkins 1991; Lett 1990; Lett et al. 1987a). Increases in linear energy transfer (LET infinity) up to 100-200 keV micron-1 cause increases in the extents of neighboring chemical and physical damage in DNA denoted by the general term DSB. Those changes are accompanied by decreasing abilities of cells normally radioresistant to sparsely ionizing radiations to process DSBs in DNA and chromatin and to recover from radiation exposure, so they make significant contributions to the relative biological effectiveness (RBE) of a given radiation.(ABSTRACT TRUNCATED AT 400 WORDS).

Non-NASA Center↗

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. 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). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, 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 were 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 accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying 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.

Nemec, Marian↗

The Analysis of the Patterns of Radiation-Induced DNA Damage Foci by a Stochastic Monte Carlo Model of DNA Double Strand Breaks Induction by Heavy Ions and Image Segmentation Software

To create a generalized mechanistic model of DNA damage in human cells that will generate analytical and image data corresponding to experimentally observed DNA damage foci and will help to improve the experimental foci yields by simulating spatial foci patterns and resolving problems with quantitative image analysis. Material and Methods: The analysis of patterns of RIFs (radiation-induced foci) produced by low- and high-LET (linear energy transfer) radiation was conducted by using a Monte Carlo model that combines the heavy ion track structure with characteristics of the human genome on the level of chromosomes. The foci patterns were also simulated in the maximum projection plane for flat nuclei. Some data analysis was done with the help of image segmentation software that identifies individual classes of RIFs and colocolized RIFs, which is of importance to some experimental assays that assign DNA damage a dual phosphorescent signal. Results: The model predicts the spatial and genomic distributions of DNA DSBs (double strand breaks) and associated RIFs in a human cell nucleus for a particular dose of either low- or high-LET radiation. We used the model to do analyses for different irradiation scenarios. In the beam-parallel-to-the-disk-of-a-flattened-nucleus scenario we found that the foci appeared to be merged due to their high density, while, in the perpendicular-beam scenario, the foci appeared as one bright spot per hit. The statistics and spatial distribution of regions of densely arranged foci, termed DNA foci chains, were predicted numerically using this model. Another analysis was done to evaluate the number of ion hits per nucleus, which were visible from streaks of closely located foci. In another analysis, our image segmentaiton software determined foci yields directly from images with single-class or colocolized foci. Conclusions: We showed that DSB clustering needs to be taken into account to determine the true DNA damage foci yield, which helps to determine the DSB yield. Using the model analysis, a researcher can refine the DSB yield per nucleus per particle. We showed that purely geometric artifacts, present in the experimental images, can be analytically resolved with the model, and that the quantization of track hits and DSB yields can be provided to the experimentalists who use enumeration of radiation-induced foci in immunofluorescence experiments using proteins that detect DNA damage. An automated image segmentaiton software can prove useful in a faster and more precise object counting for colocolized foci images.

Ponomarev, Artem↗

Extension of a nonlinear systems theory to general-frequency unsteady transonic aerodynamic responses

A methodology for modeling nonlinear unsteady aerodynamic responses, for subsequent use in aeroservoelastic analysis and design, using the Volterra-Wiener theory of nonlinear systems is presented. The methodology is extended to predict nonlinear unsteady aerodynamic responses of arbitrary frequency. The Volterra-Wiener theory uses multidimensional convolution integrals to predict the response of nonlinear systems to arbitrary inputs. The CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code is used to generate linear and nonlinear unit impulse responses that correspond to each of the integrals for a rectangular wing with a NACA 0012 section with pitch and plunge degrees of freedom. The computed kernels then are used to predict linear and nonlinear unsteady aerodynamic responses via convolution and compared to responses obtained using the CAP-TSD code directly. The results indicate that the approach can be used to predict linear unsteady aerodynamic responses exactly for any input amplitude or frequency at a significant cost savings. Convolution of the nonlinear terms results in nonlinear unsteady aerodynamic responses that compare reasonably well with those computed using the CAP-TSD code directly but at significant computational cost savings.

Silva, Walter A.↗

Dynamic analysis of flexible multibody structures

A system composed of several interconnected elastic components that may experience large angular motion relative to each other during operation is referred to as a flexible multibody structure. Several formulations were proposed for the determination of the dynamic response of controlled flexible multibody structures. In general, these formulations consist of superposing elastic deformations of the component body onto the large rigid body motion of the component. It was shown that this particular methodology for combining linear structural deformations with nonlinear kinematics can lead to erroneous response predictions when either the beam member is very flexible or the rotational speed is high. In addition, previous formulations introduce constraint equations to define the interrelations among system components. This approach increases the number of equations that must be solved, and may result in contraint violation when numerical error accumulates during the integration process. In order to overcome the difficulties, a new approach was suggested. The approach is essentially a finite element formulation which takes advantage of the fact that many multibody structures are joint dominated. The Large Angle Transient Dynamic Analysis (LATDYN) program for clarity of documentation, ease of use, user friendliness, modeling generality, and accuracy of results was evaluated. This required gaining a working familiarity with the code and performing several case studies.

Hernried, Alan G.↗

Technical Evaluation of the NASA Model for Cancer Risk to Astronauts Due to Space Radiation

At the request of NASA, the National Research Council's (NRC's) Committee for Evaluation of Space Radiation Cancer Risk Model reviewed a number of changes that NASA proposes to make to its model for estimating the risk of radiation-induced cancer in astronauts. The NASA model in current use was last updated in 2005, and the proposed model would incorporate recent research directed at improving the quantification and understanding of the health risks posed by the space radiation environment. NASA's proposed model is defined by the 2011 NASA report Space Radiation Cancer Risk Projections and Uncertainties 2010 (Cucinotta et al., 2011). The committee's evaluation is based primarily on this source, which is referred to hereafter as the 2011 NASA report, with mention of specific sections or tables cited more formally as Cucinotta et al. (2011). The overall process for estimating cancer risks due to low linear energy transfer (LET) radiation exposure has been fully described in reports by a number of organizations. They include, more recently: (1) The "BEIR VII Phase 2" report from the NRC's Committee on Biological Effects of Ionizing Radiation (BEIR) (NRC, 2006); (2) Studies of Radiation and Cancer from the United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR, 2006), (3) The 2007 Recommendations of the International Commission on Radiological Protection (ICRP), ICRP Publication 103 (ICRP, 2007); and (4) The Environmental Protection Agency s (EPA s) report EPA Radiogenic Cancer Risk Models and Projections for the U.S. Population (EPA, 2011). The approaches described in the reports from all of these expert groups are quite similar. NASA's proposed space radiation cancer risk assessment model calculates, as its main output, age- and gender-specific risk of exposure-induced death (REID) for use in the estimation of mission and astronaut-specific cancer risk. The model also calculates the associated uncertainties in REID. The general approach for estimating risk and uncertainty in the proposed model is broadly similar to that used for the current (2005) NASA model and is based on recommendations by the National Council on Radiation Protection and Measurements (NCRP, 2000, 2006). However, NASA's proposed model has significant changes with respect to the following: the integration of new findings and methods into its components by taking into account newer epidemiological data and analyses, new radiobiological data indicating that quality factors differ for leukemia and solid cancers, an improved method for specifying quality factors in terms of radiation track structure concepts as opposed to the previous approach based on linear energy transfer, the development of a new solar particle event (SPE) model, and the updates to galactic cosmic ray (GCR) and shielding transport models. The newer epidemiological information includes updates to the cancer incidence rates from the life span study (LSS) of the Japanese atomic bomb survivors (Preston et al., 2007), transferred to the U.S. population and converted to cancer mortality rates from U.S. population statistics. In addition, the proposed model provides an alternative analysis applicable to lifetime never-smokers (NSs). Details of the uncertainty analysis in the model have also been updated and revised. NASA's proposed model and associated uncertainties are complex in their formulation and as such require a very clear and precise set of descriptions. The committee found the 2011 NASA report challenging to review largely because of the lack of clarity in the model descriptions and derivation of the various parameters used. The committee requested some clarifications from NASA throughout its review and was able to resolve many, but not all, of the ambiguities in the written description.

Source record↗

The effects of curvature and viscosity on baroclinic instability: A two-layer model

A linear stability analysis of a baroclinic zonal current contained between two parallel rigid boundaries is presented. Curvature is included by performing the analysis on a beta b-plane and viscosity by allowing for the effects of Ekman layers on the rigid boundaries. A two-layer model is used. This calculation was carried out to assist in the design of a spherical model of the general circulation of the earth's atmosphere for Spacelab. In the low-gravity environment on an orbiting vehicle, a dominant radial dielectric body force, analogous to planetary gravity, can be achieved over a volume of liquid held between two concentric spheres. The results show the Eady short wavelength cutoff, and long wavelength cutoffs due to Ekman damping and curvature.

Fowlis, W. W.↗

New nonlinear control algorithms for multiple robot arms

Multiple coordinated robot arms are modeled by considering the arms as closed kinematic chains and as a force-constrained mechanical system working on the same object simultaneously. In both formulations, a novel dynamic control method is discussed. It is based on feedback linearization and simultaneous output decoupling technique. By applying a nonlinear feedback and a nonlinear coordinate transformation, the complicated model of the multiple robot arms in either formulation is converted into a linear and output decoupled system. The linear system control theory and optimal control theory are used to design robust controllers in the task space. The first formulation has the advantage of automatically handling the coordination and load distribution among the robot arms. In the second formulation, it was found that by choosing a general output equation it became possible simultaneously to superimpose the position and velocity error feedback with the force-torque error feedback in the task space.

Tarn, T. J.↗

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗