Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Matrix”

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 307 records · Page 17

Finite element applications to explore the effects of partial bonding on metal matrix composite properties

The mechanics of materials approach (definition of E, G, Nu, and Alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E sub (111) and alpha sub (111) significantly.

Caruso, J. J.↗

Finite element applications to explore the effects of partial bonding on metal matrix composite properties

The mechanics of materials approach (definition of E, G, nu, and alpha) and the finite element method are used to explore the effects of partial bonding and fiber fracture on the behavior of high temperature metal matrix composites. Composite ply properties are calculated for various degrees of disbonding to evaluate the sensitivity of these properties to the presence of fiber/matrix disbonding and fiber fracture. The mechanics of materials approach allows for the determination of the basic ply material properties needed for design/analysis of composites. The finite element method provides the necessary structural response (forces and displacements) for the mechanics of materials equations. Results show that disbonding of fractured fibers affect only E-l(11) and alpha-l(11) significantly.

Caruso, J. J.↗

Multiresponse Imaging For Improved Resolution

Images restored with resolution finer than sampling lattice. Multiresponse imaging overcomes sampling-passband constraint and critically constrained only by optical response and sensitivity of image-gathering device. Process consists of multiresponse image gathering and Wiener-matrix restoration. Image-gathering process acquires sequence of images, each with different optical response so within-pass-band and aliased signal components weighted differently. Wiener-matrix filter, in turn, unscrambles within-passband and aliased frequency components of undersampled signal and restores them up to cutoff frequency of optical response.

Fales, Carl L., Jr.↗

Electron density diagnostics for gaseous nebulae involving the O 4 intercombination lines near 1400 A

Theoretical O IV electron density sensitive emission line ratios, determined using electron impact excitation rates calculated with the R-matrix code, are presented for R(sub 1) = I(1407.4 A)/I(1401.2 A), R(sub 2) = I(1404.8 A)/I(1401.2A), R(sub 3) = I(1399.8 A)/(1401.2 A), and R(sub 4) = I(1397.2 A)/I(1401.2 A). The observed values of R(sub 1)-R(sub 4), measured from high resolution spectra obtained with the International Ultraviolet Explorer (IUE) satellite, lead to electron densities that are compatible, and which are also in good agreement with those deduced from line ratios in other species. This provides observational support for the accuracy of the atomic data adopted in the present calculations.

Keenan, F. P.↗

Uncertainty Modeling for Structural Control Analysis and Synthesis

The development of an accurate model of uncertainties for the control of structures that undergo a change in operational environment, based solely on modeling and experimentation in the original environment is studied. The application used throughout this work is the development of an on-orbit uncertainty model based on ground modeling and experimentation. A ground based uncertainty model consisting of mean errors and bounds on critical structural parameters is developed. The uncertainty model is created using multiple data sets to observe all relevant uncertainties in the system. The Discrete Extended Kalman Filter is used as an identification/parameter estimation method for each data set, in addition to providing a covariance matrix which aids in the development of the uncertainty model. Once ground based modal uncertainties have been developed, they are localized to specific degrees of freedom in the form of mass and stiffness uncertainties. Two techniques are presented: a matrix method which develops the mass and stiffness uncertainties in a mathematical manner; and a sensitivity method which assumes a form for the mass and stiffness uncertainties in macroelements and scaling factors. This form allows the derivation of mass and stiffness uncertainties in a more physical manner. The mass and stiffness uncertainties of the ground based system are then mapped onto the on-orbit system, and projected to create an analogous on-orbit uncertainty model in the form of mean errors and bounds on critical parameters. The Middeck Active Control Experiment is introduced as experimental verification for the localization and projection methods developed. In addition, closed loop results from on-orbit operations of the experiment verify the use of the uncertainty model for control analysis and synthesis in space.

Campbell, Mark E.↗

Energetic Consistency and Coupling of the Mean and Covariance Dynamics

The dynamical state of the ocean and atmosphere is taken to be a large dimensional random vector in a range of large-scale computational applications, including data assimilation, ensemble prediction, sensitivity analysis, and predictability studies. In each of these applications, numerical evolution of the covariance matrix of the random state plays a central role, because this matrix is used to quantify uncertainty in the state of the dynamical system. Since atmospheric and ocean dynamics are nonlinear, there is no closed evolution equation for the covariance matrix, nor for the mean state. Therefore approximate evolution equations must be used. This article studies theoretical properties of the evolution equations for the mean state and covariance matrix that arise in the second-moment closure approximation (third- and higher-order moment discard). This approximation was introduced by EPSTEIN [1969] in an early effort to introduce a stochastic element into deterministic weather forecasting, and was studied further by FLEMING [1971a,b], EPSTEIN and PITCHER [1972], and PITCHER [1977], also in the context of atmospheric predictability. It has since fallen into disuse, with a simpler one being used in current large-scale applications. The theoretical results of this article make a case that this approximation should be reconsidered for use in large-scale applications, however, because the second moment closure equations possess a property of energetic consistency that the approximate equations now in common use do not possess. A number of properties of solutions of the second-moment closure equations that result from this energetic consistency will be established.

Cohn, Stephen E.↗

Optimization of Aerospace Structure Subject to Damage Tolerance Criteria

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Impacts of Crystalline Host Rock on Repository Barrier Materials at 250 °C: Hydrothermal Co-Alteration of Wyoming Bentonite and Steel in the Presence of Grimsel Granodiorite

Direct disposal of dual-purpose canisters (DPC) has been proposed to streamline the disposal of spent nuclear fuel. However, there are scenarios where direct disposal of DPCs may result in temperatures in excess of the specified upper temperature limits for some engineered barrier system (EBS) materials, which may cause alteration within EBS materials dependent on local conditions such as host rock composition, chemistry of the saturating groundwaters, and interactions between barrier materials themselves. Here we report the results of hydrothermal experiments reacting EBS materials—bentonite buffer and steel—with an analogue crystalline host rock and groundwater at 250 °C. Experiment series explored the effect of reaction time on the final products and the effects of the mineral and fluid reactants on different steel types. Post-mortem X-ray diffraction, electron microprobe, and scanning electron microscopy analyses showed characteristic alteration of both bentonite and steel, including the formation of secondary zeolite and calcium silicate hydrate minerals within the bentonite matrix and the formation of iron-bearing clays and metal oxides at the steel surfaces. Swelling clays in the bentonite matrix were not quantitatively altered to non-swelling clay species by the hydrothermal conditions. The combined results of the solution chemistry over time and post-mortem mineralogy suggest that EBS alteration is more sensitive to initial groundwater chemistry than the presence of host rock, where limited potassium concentration in the solution prohibits conversion of the smectite minerals in the bentonite matrix to non-swelling clay species.

54 ENVIRONMENTAL SCIENCES↗

Experiences on p-Version Time-Discontinuous Galerkin's Method for Nonlinear Heat Transfer Analysis and Sensitivity Analysis

The focus of this research is on the development of analysis and sensitivity analysis equations for nonlinear, transient heat transfer problems modeled by p-version, time discontinuous finite element approximation. The resulting matrix equation of the state equation is simply in the form ofA(x)x = c, representing a single step, time marching scheme. The Newton-Raphson's method is used to solve the nonlinear equation. Examples are first provided to demonstrate the accuracy characteristics of the resultant finite element approximation. A direct differentiation approach is then used to compute the thermal sensitivities of a nonlinear heat transfer problem. The report shows that only minimal coding effort is required to enhance the analysis code with the sensitivity analysis capability.

Hou, Gene↗

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

Comparison of the Baseline USL Calculation Methods for Loosely-Coupled and Novel Neutronic Systems [Slides]

Current work includes reconstructing 187-group ENDF B/VII.1 covariance matrix for comparison study to 44-group ENDF B/VII.1 matrix to determine how covariance matrix structures affect the USL calculations and investigating the exact domination mechanism of region-wise sensitivities in a loosely-coupled system. Future work will involve seeking to understand how bias distributions change with reactivity and how to calculate propagate distribution error into USL calculations, as well as how cross section perturbation studies can be extended to other types of calculations, such as shielding calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Failure of morphology of (0 deg)8 graphite/epoxy as influenced by environments and processing

Optical and scanning electron microscopy were used to investigate the failure morphology of graphite/epoxy specimens which had been tested until tensile failure. Failure morphology was studied as a function of the quality control variables of specimen preparation technique, prepreg batch, and cure condition, and also as a function of the environmental parameters of temperature and moisture content. Defective specimens were found to exhibit a low energy failure morphology. Poor specimen edge preparation and one batch of prepreg when tested at elevated temperature or moisture content also exhibited energy failure morphology. Postcuring had no effect on strength but did slightly alter failure morphology. Temperature or moisture appeared to decrease flaw sensitivity and thus increase strength; however, moisture also appeared to increase interfacial debonding between filament and matrix. When combined moisture and temperature increased interfacial debonding and made the epoxy matrix more prone to fracture.

Clements, L. L.↗

Concepts for improving the damage tolerance of composite compression panels

The results of an experimental evaluation of graphite-epoxy composite compression panel impact damage tolerance and damage propagation arrest concepts are reported. The tests were conducted on flat plate specimens and blade-stiffened structural panels such as those used in commercial aircraft wings, and the residual strength of damaged specimens and their sensitivity to damage while subjected to in-plane compression loading were determined. Results suggest that matrix materials that fail by delamination have the lowest damage tolerance, and it is concluded that alternative matrix materials with transverse reinforcement to suppress the delamination failure mode and yield the higher-strain value transverse shear crippling mode should be developed.

Rhodes, M. D.↗

A solution of Mx(double dot) + Cx(dot) + Kx = 0 applicable to the design of active dampers

A solution is presented for the equations of motion for the damped linear oscillator, Mx(double dot) + Cx(dot) + Kx = 0. The algorithm solves a transformed set of equations in terms of the modal variables of the undamped system and, at the same time, solves the adjoint equation of the transformed problem. The adjoint solution is normalized to give the inverse of the solution matrix of the transformed problem. The normalized inverse is useful in design for direct computation of sensitivity derivatives of damping ratios with respect to damping rates. The algorithm is programmed to reduce storage requirements by a factor of three-fourths compared to standard complex eigenvalue subroutines. A numerical example is included.

Thurston, G. A.↗

N IV emission lines in the ultraviolet spectra of gaseous nebulae

Theoretical electron density sensitive emission-line ratios, determined using electron impact excitation rates calculated with the R-matrix code, are presented for R = I(2s(sup 2) (1)S-2s2p 3P(sub 2))/I(2s(sup 2) (1)S-2s2p 3P(sub 1) = 1(1483 A)/I(1486 A) in N IV. These are found to be up to an order of magnitude different from those deduced by previous authors, principally due to the inclusion of excitation rates for transitions among the 2s2p (3)p fine-structure levels. The observed values of R for several planetary nebulae, symbiotic stars and the Cygnus Loop supernova remnant, measured from spectra obtained with the International Ultraviolet Explorer (IUE) satellite and the Hopkins Ultraviolet Explorer (HUT), lead to electron densities which are in excellent agreement with those deduced from line ratios in other species. This provides observational support for the accuracy of the atomic data adopted in the present calculations.

Keenan, F. P.↗

Simultaneous Luminescence Pressure and Temperature Measurement System for Hypersonic Wind Tunnels

Surface pressures and temperatures are determined from visible emission brightness and green-to-red color ratioing of induced luminescence from a ceramic surface with an organic dye coating. A ceramic-dye matrix of porous silica ceramic with an adsorbed dye is developed for high-temperature pressure sensitivity and stability (up to 150 C). Induced luminescence may be excited using a broad range of incident radiation from visible blue light (488-nm wavelength) to the near ultraviolet (365 nm). Ceramic research models and test samples are fabricated using net-form slip-casting and sintering techniques. Methods of preparation and effects of adsorption film thickness on measurement sensitivity are discussed. With the present 8-bit imaging system a 10% pressure measurement uncertainty from 50 to 760 torr is estimated, with an improvement to 5% from 3 to 1500 torr with a 12-bit imaging system.

Buck, Gregory M.↗

SPARROW: Steam Propelled Autonomous Retrieval Robot for Ocean Worlds

The Steam Propelled Autonomous Retrieval Robot (SPARROW) for Ocean Worlds was a Phase I mission concept study funded under the NASA NIAC program. This report represents the findings of that study and recommendations for future work. SPARROW, envisioned as a soccer ball-sized payload to a primary lander mission, is a propulsively hopping robot for the exploration of Europa's rugged, icy surface. A multi-thruster, passively gimballed robot within a protective, spherical shell, SPARROW is able to freely rotate, self-right, and tumble over chaotic terrains. Europa's abundant surface ice would be harvested as an in situ propellant source. The principal objective of SPARROW is to increase the science return of a Europa landed asset by enabling access to distal, spatially distributed geologic units. The design of mobility systems for Europa is challenging, due in part to its almost entirely unconstrained surface topography and strength. Images returned by Voyager and Galileo yielded resolutions on the order of hundreds of meters per pixel, with localized regions reaching 6 meters per pixel—still far larger than a typical rover. A key benefit of SPARROW's hopping, impact-tolerant design, is that it eliminates the need for a priori information regarding terrain topography and surface strength; no surface reaction forces are required for motion. In this context, SPARROW is believed to be entirely terrain agnostic. In this report we detail the results of three study objectives: i) to quantify the energy required to collect surface ice, change its phase, and maintain propellant temperature, ii) to identify control and estimation strategies that enable SPARROW to successfully reach, and return from, regions of scientific interest, and iii) to characterize the impact of SPARROW's range on likely science return. Five water-based propellant architectures are presented alongside their mass, power, and volume requirements. Monte Carlo simulations of SPARROW hopping and tumbling over 1 km of glacial ice are summarized, characterizing SPARROW's sensitivity to uncertainty in: initial pose, thrust profile, and vehicle-terrain interaction. A science traceability matrix is presented, which details the effect of sortie range on three science goals: constraining Europa's evolutionary morphology, assessing sub-surface ocean habitability, and searching for life and/or biosignatures.

Autonomous↗