Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density matrix method”

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 181 records · Page 10

Bayesian Approach to the Joint Inversion of Gravity and Magnetic Data, with Application to the Ismenius Area of Mars

This viewgraph presentation reviews a Bayesian approach to the inversion of gravity and magnetic data with specific application to the Ismenius Area of Mars. Many inverse problems encountered in geophysics and planetary science are well known to be non-unique (i.e. inversion of gravity the density structure of a body). In hopes of reducing the non-uniqueness of solutions, there has been interest in the joint analysis of data. An example is the joint inversion of gravity and magnetic data, with the assumption that the same physical anomalies generate both the observed magnetic and gravitational anomalies. In this talk, we formulate the joint analysis of different types of data in a Bayesian framework and apply the formalism to the inference of the density and remanent magnetization structure for a local region in the Ismenius area of Mars. The Bayesian approach allows prior information or constraints in the solutions to be incorporated in the inversion, with the "best" solutions those whose forward predictions most closely match the data while remaining consistent with assumed constraints. The application of this framework to the inversion of gravity and magnetic data on Mars reveals two typical challenges - the forward predictions of the data have a linear dependence on some of the quantities of interest, and non-linear dependence on others (termed the "linear" and "non-linear" variables, respectively). For observations with Gaussian noise, a Bayesian approach to inversion for "linear" variables reduces to a linear filtering problem, with an explicitly computable "error" matrix. However, for models whose forward predictions have non-linear dependencies, inference is no longer given by such a simple linear problem, and moreover, the uncertainty in the solution is no longer completely specified by a computable "error matrix". It is therefore important to develop methods for sampling from the full Bayesian posterior to provide a complete and statistically consistent picture of model uncertainty, and what has been learned from observations. We will discuss advanced numerical techniques, including Monte Carlo Markov

data analysis↗

Advancing Development of Environmental Barrier Coatings Resistant to Attack by Molten Calcium-Magnesium-Aluminosilicate (CMAS)

Ceramic matrix composites (CMCs) are a leading material system to replace metal-based parts in the hot-section of air-breathing turbine engines to improve fuel efficiency in aircraft engines. CMCs have higher temperature capabilities and lower density compared with traditional metallic structural materials. However, silicon-based CMCs are susceptible to oxidation in the harsh combustion environment encountered in turbine engines. Consequently, environmental barrier coatings (EBCs) are being developed to protect CMC components to improve durability and extend service life of CMCs. Sand, volcanic ash and other particulate debris, which are generally comprised of calcium-magnesium-aluminosilicate (CMAS) and other trace oxides, are routinely ingested by aircraft engines. At temperatures above 1200°C, CMAS particulates melt. Near target operating temperatures (~1500°C) of future CMC-based aircraft engines, molten CMAS behaves like a viscous melt that can infiltrate and chemically interact with protective coatings. These interactions can cause premature failure of the EBC system and ultimately the overall CMC engine component. Degradation of candidate EBC materials by molten CMAS will be presented with a focus on recent work, as well as methods of evaluating the complex high-temperature materials interactions, underway at NASA Glenn Research Center.

Wiesner, Valerie L.↗

Properties of silicon carbide fiber-reinforced silicon nitride matrix composites

The mechanical properties of NASA Lewis developed SiC/RBSN composites and their thermal and environmental stability havd been studied. The composites consist of nearly 30 vol pct of aligned 142 micron diameter chemically vapor-deposited SiC fibers in a relatively porous silicon nitride matrix. In the as-fabricated condition, the unidirectional and 2-D composites exhibited metal-like stress-strain behavior, graceful failure, and showed improved properties when compared with unreinforced matrix of comparable density. Furthermore, the measured room temperature tensile properties were relativley independent of tested volume and were unaffected by artifical notches normal to the loading direction or by thermal shocking from temperatures up to 800 C. The four-point bend strength data measured as a function of temperature to 1400 C in air showed that as-fabricated strength was maintained to 1200 C. At 1400 C, however, nearly 15 pct loss in strength was observed. Measurement of room temperature tensile strength after 100 hr exposure at temperatures to 1400 C in a nitrogen environment indicated no loss from the as-fabricated composite strength. On the other hand, after 100 hr exposure in flowing oxygen at 1200 and 1400 C, the composites showed approximately 40 pct loss from their as-fabricated ultimate tensile strength. Those exposed between 400 to 1200 C showed nearly 60 pct strength loss. Oxidation of the fiber/matrix interface as well as internal oxidation of the porous Si3N4 matrix are likely mechanisms for strength degradation. The excellent strength reproducibility, notch insensitivity, and high temperature strength of the composite makes it an ideal candidate for advanced heat engine applications provided coating or densification methods are developed to avoid internal oxidation attack.

Bhatt, Ramakrishna T.↗

An Investigation of Reliability Models for Ceramic Matrix Composites and their Implementation into Finite Element Codes

The development of modeling approaches for the failure analysis of ceramic-based material systems used in high temperature environments was the primary objective of this research effort. These materials have the potential to support many key engineering technologies related to the design of aeropropulsion systems. Monolithic ceramics exhibit a number of useful properties such as retention of strength at high temperatures, chemical inertness, and low density. However, the use of monolithic ceramics has been limited by their inherent brittleness and a large variation in strength. This behavior has motivated material scientists to reinforce the monolithic material with a ceramic fiber. The addition of a second ceramic phase with an optimized interface increases toughness and marginally increases strength. The primary purpose of the fiber is to arrest crack growth, not to increase strength. The material systems of interest in this research effort were laminated ceramic matrix composites, as well as two- and three- dimensional fabric reinforced ceramic composites. These emerging composite systems can compete with metals in many demanding applications. However, the ongoing metamorphosis of ceramic composite material systems, and the lack of standardized design data has in the past tended to minimize research efforts related to structural analysis. Many structural components fabricated from ceramic matrix composites (CMC) have been designed by "trial and error." The justification for this approach lies in the fact that during the initial developmental phases for a material system fabrication issues are paramount. Emphasis is placed on demonstrating feasibility rather than fully understanding the processes controlling mechanical behavior. This is understandable during periods of rapid improvements in material properties for any composite system. But to avoid the ad hoc approach, the analytical methods developed under this effort can be used to develop rational structural design protocols.

Duffy, Stephen F.↗

Efficient Implementation of an Optimal Interpolator for Large Spatial Data Sets

Scattered data interpolation is a problem of interest in numerous areas such as electronic imaging, smooth surface modeling, and computational geometry. Our motivation arises from applications in geology and mining, which often involve large scattered data sets and a demand for high accuracy. The method of choice is ordinary kriging. This is because it is a best unbiased estimator. Unfortunately, this interpolant is computationally very expensive to compute exactly. For n scattered data points, computing the value of a single interpolant involves solving a dense linear system of size roughly n x n. This is infeasible for large n. In practice, kriging is solved approximately by local approaches that are based on considering only a relatively small'number of points that lie close to the query point. There are many problems with this local approach, however. The first is that determining the proper neighborhood size is tricky, and is usually solved by ad hoc methods such as selecting a fixed number of nearest neighbors or all the points lying within a fixed radius. Such fixed neighborhood sizes may not work well for all query points, depending on local density of the point distribution. Local methods also suffer from the problem that the resulting interpolant is not continuous. Meyer showed that while kriging produces smooth continues surfaces, it has zero order continuity along its borders. Thus, at interface boundaries where the neighborhood changes, the interpolant behaves discontinuously. Therefore, it is important to consider and solve the global system for each interpolant. However, solving such large dense systems for each query point is impractical. Recently a more principled approach to approximating kriging has been proposed based on a technique called covariance tapering. The problems arise from the fact that the covariance functions that are used in kriging have global support. Our implementations combine, utilize, and enhance a number of different approaches that have been introduced in literature for solving large linear systems for interpolation of scattered data points. For very large systems, exact methods such as Gaussian elimination are impractical since they require 0(n(exp 3)) time and 0(n(exp 2)) storage. As Billings et al. suggested, we use an iterative approach. In particular, we use the SYMMLQ method, for solving the large but sparse ordinary kriging systems that result from tapering. The main technical issue that need to be overcome in our algorithmic solution is that the points' covariance matrix for kriging should be symmetric positive definite. The goal of tapering is to obtain a sparse approximate representation of the covariance matrix while maintaining its positive definiteness. Furrer et al. used tapering to obtain a sparse linear system of the form Ax = b, where A is the tapered symmetric positive definite covariance matrix. Thus, Cholesky factorization could be used to solve their linear systems. They implemented an efficient sparse Cholesky decomposition method. They also showed if these tapers are used for a limited class of covariance models, the solution of the system converges to the solution of the original system. Matrix A in the ordinary kriging system, while symmetric, is not positive definite. Thus, their approach is not applicable to the ordinary kriging system. Therefore, we use tapering only to obtain a sparse linear system. Then, we use SYMMLQ to solve the ordinary kriging system. We show that solving large kriging systems becomes practical via tapering and iterative methods, and results in lower estimation errors compared to traditional local approaches, and significant memory savings compared to the original global system. We also developed a more efficient variant of the sparse SYMMLQ method for large ordinary kriging systems. This approach adaptively finds the correct local neighborhood for each query point in the interpolation process.

Memarsadeghi, Nargess↗

Calculation of molecular excitation rates

State-to-state collisional excitation rates for interstellar molecules observed by radio astronomers continue to be required to interpret observed line intensities in terms of local temperatures and densities. A problem of particular interest is collisional excitation of water which is important for modeling the observed interstellar masers. In earlier work supported by a different NASA Grant, excitation of water in collisions with He atoms was studied; after many years of successively more refined calculations that problem now seems to be well understood, and discrepancies with earlier experimental data for related (pressure broadening) phenomena are believed to reflect experimental errors. Because of interstellar abundances, excitation by H2, the dominant interstellar species, is much more important than excitation by He, although it has been argued that rates for excitation by these are similar. Under the current grant theoretical study of this problem has begun which is greatly complicated by the additional degrees of freedom which must be included both in determining the interaction potential and also in the molecular scattering calculation. We have now computed the interaction forces for nearly a thousand molecular geometries and are close to having an acceptable global fit to these points which is necessary for the molecular dynamics calculations. Also, extensive modifications have been made to the molecular scattering code, MOLSCAT. These included coding the rotational basis sets and coupling matrix elements required for collisions of an asymmetric top with a linear rotor. A new method for numerical solution of the coupled equations has been incorporated. Because of the long-ranged nature of the water-hydrogen interaction it is necessary to integrate the equations to rather large intermolecular separations, and the integration methods previously available in MOLSCAT are not ideal for such cases. However, the method used by Alexander in his HIBRIDON code is particularly suited for such cases. We have obtained this code and incorporated that part which solves the coupled differential equations as an option in the MOLSCAT program.

Flynn, George↗

Acousto-ultrasonics to Assess Material and Structural Properties

This report was created to serve as a manual for applying the Acousto-Ultrasonic NDE method, as practiced at NASA Glenn, to the study of materials and structures for a wide range of applications. Three state of the art acousto-ultrasonic (A-U) analysis parameters, ultrasonic decay (UD) rate, mean time (or skewing factor, "s"), and the centroid of the power spectrum, "f(sub c)," have been studied and applied at GRC for NDE interrogation of various materials and structures of aerospace interest. In addition to this, a unique application of Lamb wave analysis is shown. An appendix gives a brief overview of Lamb Wave analysis. This paper presents the analysis employed to calculate these parameters and the development and reasoning behind their use. It also discusses the planning of A-U measurements for materials and structures to be studied. Types of transducer coupling are discussed including contact and non-contact via laser and air. Experimental planning includes matching transducer frequency range to material and geometry of the specimen to be studied. The effect on results of initially zeroing the DC component of the ultrasonic waveform is compared with not doing so. A wide range of interrogation problems are addressed via the application of these analysis parameters to real specimens is shown for five cases: Case 1: Differences in density in [0] SiC/RBSN ceramic matrix composite. Case 2: Effect of tensile fatigue cycling in [+/-45] SiC/SiC ceramic matrix composite. Case 3: Detecting creep life, and failure, in Udimet 520 Nickel-Based Super Alloy. Case 4: Detecting Surface Layer Formation in T-650-35/PMR-15 Polymer Matrix Composites Panels due to Thermal Aging. Case 5: Detecting Spin Test Degradation in PMC Flywheels. Among these cases a wide range of materials and geometries are studied.

Kautz, Harold E.↗

Effects of Fiber Reinforcement Architecture on the Hygrothermal-Mechanical Performance of Polyimide Matrix Composites for Aeropropulsion Applications

A lightweight sandwich support structure, for the combustor chamber of a new generation liquid propellant rocket engine, was designed and fabricated using a polymer matrix composite (PMC) facesheet on a Ti honeycomb core. The PMC facesheet consisted of high stiffness carbon fiber, M40JB, and high temperature Polyimides, such as PMR-II-50 and HFPE-II-52. Six different fiber architectures; four harness satin (4HS) woven fabric, uni-tape, woven-uni hybrid, stitched woven fabric, stitched uni-tape and triaxial braided structures have been investigated for optimum stiffness-thickness-weight-hygrothermal performance design criteria for the hygrothermal-mechanical propulsion service exposure conditions including rapid heating up to 200 F/sec, maximum operating temperature of 600 F, internal pressure up to 100 psi. One of the specific objectives in this study is to improve composite blistering resistance in z-direction at minimum expense of in-plane mechanical properties. An extensive property-performance database including dry-wet mechanical properties at various temperatures, thermal-physical properties, such as blistering onset condition was generated for fiber architecture down-selection and design guidelines. Various optimized process methods such as vacuum bag compression molding, solvent assistant resin transfer molding (SaRTM), resin film infusion (RFI) and autoclaving were utilized for PMC panel fabrication depending on the architecture type. In the case of stitched woven fabric architecture, the stitch pattern in terms of stitch density and yarn size was optimized based on both in-plane mechanical properties and blistering performance. Potential reduction of the in-plane properties transverse to the line of stitching was also evaluated. Efforts have been made to correlate the experimental results with theoretical micro-mechanics predictions. Changes in deformation mechanism and failure sequences in terms of fiber architecture will be discussed.

Shin, E. Eugene↗

Method and apparatus for the collection, storage, and real time analysis of blood and other bodily fluids

The present invention provides a simple, portable, relatively inexpensive apparatus for accurately and efficiently collecting, separating, testing, and even storing between about 1-20 ml, preferably about 1-10 ml, of blood or other bodily fluid in situ. The apparatus includes a collection chamber bounded on its sides by an opening in a sheet of material, preferably clear plastic, abutting a filter card. The filter card is made of fibrous material, preferably less than about a millimeter thick, having an average pore size of less than about 3 microns. Preferably, the fibers are glass and the fibrous material has an average pore size of about 1 micron. The fibrous material is treated with a carbohydrate/protein mixture which contains between about 1-40 percent wt/vol carbohydrate and about 0.1-15 percent wt/vol nonspecific protein, preferably between about 10-20 percent carbohydrate and about 5-8 percent protein. A preferred carbohydrate/protein mixture comprises about 10 percent mannitol and about 6 percent albumin. The blood or other fluid moves through the filter card by capillary action aided by an absorbent matrix with a high Klemm factor which abuts the filter card. The absorbent matrix and/or filter card can be treated with a wide spectrum of test reagents. The speed, cleanliness, and efficiency of the separation process can be altered by: (a) changing the absolute concentration of the carbohydrate/protein mixture; (b) applying positive or negative pressure to one side of the filter; and/or (c) varying the relative density and pore size of the filter card and absorbent matrix.

Whitson, Peggy A.↗

Sintering Behaviors of ZrC, NbC, and TaC Mono-and Binary Carbides

Nuclear Thermal Propulsion (NTP) has undergone development as an alternate in-space propulsion system to traditional chemical propulsion methods since the 1950s. In an NTP system, the energy released from fission in the core is utilized as the heat source to directly heat a propellant for propulsion, rather than chemical combustion in a traditional rocket. NTP has many desirable capabilities including flexible mission launch dates and reduced transit times due to it’s capability for high specific impulse. One of the main challenges with NTP systems is the structural integrity of the fuel. The fuel form required in the reactor core must withstand temperatures above 2700 Kelvin. Ceramic-metallic matrix fuel, ceramic-ceramic matrix fuel, and solid solution carbide fuels are the three strongest candidates for the extreme environments. Solid solution carbides have the potential to exhibit promising behavior as a fuel form in an NTP system. Of the multiple refractory metal carbides of interest, here we focus on zirconium carbide (ZrC), niobium carbide (NbC), and tantalum carbide (TaC). ZrC, NbC, and TaC powders were consolidated in monocarbide (ZrC, NbC, and TaC) and bi-carbide (ZrC-NbC, ZrC-TaC, NbC-TaC) forms using spark plasma sintering (SPS). In addition to the pure endpoint carbides, the examined compositions of the various bi-carbides ranged from 25-75 mol%. The sintering temperatures, pressures, and hold times were varied to determine the ideal sintering conditions. Grain size analysis, Archimedes’ density, scanning electron microscopy (SEM), X-ray diffraction (XRD), and energy dispersion spectroscopy (EDS) were used to determine and characterize the grain size, ideal density, porosity, phase stability, and chemical composition of each sample. The data from each sample was then used to generate a Master Sintering Curve (MSC) unique to each monocarbide or bi-carbide.

nuclear thermal propulsion↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Use of Design of Experiments and Rule-Based Inference in Determining Neural Network Architectures for Loss of Control Detection

In this work, we describe methods for selecting the neural network architectures and input spaces to implement belief state inference on generic commercial transport aircraft. First, we highlight a case study on the planning, execution, and analysis of a set of experiments to determine the configurations of a conditional variational autoencoder (CVAE). We present a structured method that can be used in a number of aerospace applications, to optimize the structure and training parameters of the CVAE for belief state inference, using Design of Experiments (DOE) statistical methodologies. The motivation for this specific DOE was to identify the appropriate hyperparameters for measuring the CVAE reconstruction probability and latent space, such that the measurements can be used to infer qualitative state changes for the aircraft. We demonstrate that this process yields information about a trained neural network’s utility for this specific application, along with a quantifiable range of certainty. We execute 84 experiments using loss-of-control flight maneuver data from a NASA T-2 aircraft, demonstrating that this empirical process allows us to construct cheap and simple models with specific attributes amenable to belief state inference in aerospace applications. While theoretically, we could create a single CVAE with an input space the size of all measurable flight variables and environmental dynamics, it becomes intractable to use such a neural network in an in-situ intelligent multi-agent system. Using the recommendations from our case study, we introduce a technical approach for feasibly describing the belief space by (1) identifying significant statistical relationships among flight variables using rule induction, (2) using a set of rules that cover all features to define the input space of multiple CVAEs, and (3) forming a belief space based on the joint probability density of their collective latent spaces. This results in a series of relatively small matrix multiplications that can be performed in real time, as opposed to large matrix computations in a single CVAE. We demonstrate the application of this approach on the T-2 flight loss-of control experiments, using the architecture and hyperparameter recommendations from the case study. We compare the utilities of an individual CVAE trained on all flight variables and multiple CVAEs defined on subsets of flight variables for detecting qualitative changes in flight. We demonstrate that the use of multiple CVAEs with smaller input spaces permits the CVAE to capture more granular relationships in the latent space, permitting better state space characterization and loss-of-control detection.

Design of experiments↗

Hot-wire calibration in subsonic/transonic flow regimes

A different approach for calibrating hot-wires, which simplifies the calibration procedure and reduces the tunnel run-time by an order of magnitude was sought. In general, it is accepted that the directly measurable quantities in any flow are velocity, density, and total temperature. Very few facilities have the capability of varying the total temperature over an adequate range. However, if the overheat temperature parameter, a(sub w), is used to calibrate the hot-wire then the directly measurable quantity, voltage, will be a function of the flow variables and the overheat parameter i.e., E = f(u,p,a(sub w), T(sub w)) where a(sub w) will contain the needed total temperature information. In this report, various methods of evaluating sensitivities with different dependent and independent variables to calibrate a 3-Wire hot-wire probe using a constant temperature anemometer (CTA) in subsonic/transonic flow regimes is presented. The advantage of using a(sub w) as the independent variable instead of total temperature, t(sub o), or overheat temperature parameter, tau, is that while running a calibration test it is not necessary to know the recovery factor, the coefficients in a wire resistance to temperature relationship for a given probe. It was deduced that the method employing the relationship E = f (u,p,a(sub w)) should result in the most accurate calibration of hot wire probes. Any other method would require additional measurements. Also this method will allow calibration and determination of accurate temperature fluctuation information even in atmospheric wind tunnels where there is no ability to obtain any temperature sensitivity information at present. This technique greatly simplifies the calibration process for hot-wires, provides the required calibration information needed in obtaining temperature fluctuations, and reduces both the tunnel run-time and the test matrix required to calibrate hotwires. Some of the results using the above techniques are presented in an appendix.

Nagabushana, K. A.↗

Demonstration of Minimally Machined Honeycomb Silicon Carbide Mirrors

Honeycomb silicon carbide composite mirrors are made from a carbon fiber preform that is molded into a honeycomb shape using a rigid mold. The carbon fiber honeycomb is densified by using polymer infiltration pyrolysis, or through a reaction with liquid silicon. A chemical vapor deposit, or chemical vapor composite (CVC), process is used to deposit a polishable silicon or silicon carbide cladding on the honeycomb structure. Alternatively, the cladding may be replaced by a freestanding, replicated CVC SiC facesheet that is bonded to the honeycomb. The resulting carbon fiber-reinforced silicon carbide honeycomb structure is a ceramic matrix composite material with high stiffness and mechanical strength, high thermal conductivity, and low CTE (coefficient of thermal expansion). This innovation enables rapid, inexpensive manufacturing. The web thickness of the new material is less than 1 millimeter, and core geometries tailored. These parameters are based on precursor carbon-carbon honeycomb material made and patented by Ultracor. It is estimated at the time of this reporting that the HoneySiC(Trademark) will have a net production cost on the order of $38,000 per square meter. This includes an Ultracor raw material cost of about $97,000 per square meter, and a Trex silicon carbide deposition cost of $27,000 per square meter. Even at double this price, HoneySiC would beat NASA's goal of $100,000 per square meter. Cost savings are estimated to be 40 to 100 times that of current mirror technologies. The organic, rich prepreg material has a density of 56 kilograms per cubic meter. A charred carbon-carbon panel (volatile organics burnt off) has a density of 270 kilograms per cubic meter. Therefore, it is estimated that a HoneySiC panel would have a density of no more than 900 kilograms per cubic meter, which is about half that of beryllium and about onethird the density of bulk silicon carbide. It is also estimated that larger mirrors could be produced in a matter of weeks. Each cell is completely uniform, maintaining the shape of the inserted mandrel. Furthermore, the layup creates pressure that insures node bond strength. Each node is a composite laminate using only the inherent resin system to form the bond. This contrasts starkly with the other known method of producing composite honeycomb, in which individual corrugations are formed, cured, and then bonded together in a secondary process. By varying the size of the mandrels within the layup, varying degrees of density can be achieved. Typical sizes are 3/8 and 3/16 in. (approximately 10 and 5 millimeters). Cell sizes up to 1 in. (approximately 25 millimeters) have been manufactured. Similarly, the shape of the core can be altered for a flexible honeycomb structure.

Goodman, William↗

Assessment of Internal Fabric of Particulate Materials

Particle arrangement and distribution within a soil matrix has long been recognized as having significant influence on the mechanical behavior of cohesionless soils. It is well known that two soil specimens having the same grain type (e.g., quartz, feldspar, etc.), same grain size distribution and relative density (or void ratio) can display completely different mechanical behavior. Because of the different fabric configurations in the otherwise similar specimens, they are likely to have different mechanical properties such as stress-strain response, anisotropy, dilatancy, etc. Soil Fabric is defined as the arrangement of particles, particle groups and associated pore space. In the literature, fabric analysis techniques are mainly classified as destructive (e.g., specimen stabilization, thin-sectioning, and microscopy), and nondestructive techniques (e.g., magnetic resonance imaging, ultrasonic testing, x-ray radiography, and computed tomography). Quantifying the void ratio and its distribution is the main parameter used to describe the fabric of particulate materials. This paper presents a comprehensive literature review of fabric analysis techniques applied to particulate materials. In addition, the results of a comprehensive investigation to quantify void ratio of sand specimens will be presented and discussed. The sand used in the experiments in a natural, uniform rounded to sub-rounded silica sand known as F-75 banding sand with mean particle size of 0.22 mm. Uniform specimens and specimens subjected to different axial-strain levels tested under triaxial and biaxial conditions are examined to evaluate void ratio evolution and its distribution using destructive thin-sectioning and nondestructive Computed Tomography (CT) techniques. Details of a new innovative polygon generation technique called Voronoi tessellation used to quantify void ratio of microscopic images of sand grains will be presented and compared to classical Oda's method. Finally, frequency distribution of void ratio calculated using CT analysis will be presented and discussed.

Alshibi, Khalid A.↗

The Digital Astronaut Project Computational Bone Remodeling Model (Beta Version) Bone Summit Summary Report

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur [1-3]. The most commonly used countermeasure against bone loss in microgravity has been prescribed exercise [4]. However, data has shown that existing exercise countermeasures are not as effective as desired for preventing bone loss in long duration, 4 to 6 months, spaceflight [1,3,5,6]. This spaceflight related bone loss may cause early onset of osteoporosis to place the astronauts at greater risk of fracture later in their lives. Consequently, NASA seeks to have improved understanding of the mechanisms of bone demineralization in microgravity in order to appropriately quantify this risk, and to establish appropriate countermeasures [7]. In this light, NASA's Digital Astronaut Project (DAP) is working with the NASA Bone Discipline Lead to implement well-validated computational models to help predict and assess bone loss during spaceflight, and enhance exercise countermeasure development. More specifically, computational modeling is proposed as a way to augment bone research and exercise countermeasure development to target weight-bearing skeletal sites that are most susceptible to bone loss in microgravity, and thus at higher risk for fracture. Given that hip fractures can be debilitating, the initial model development focused on the femoral neck. Future efforts will focus on including other key load bearing bone sites such as the greater trochanter, lower lumbar, proximal femur and calcaneus. The DAP has currently established an initial model (Beta Version) of bone loss due to skeletal unloading in femoral neck region. The model calculates changes in mineralized volume fraction of bone in this segment and relates it to changes in bone mineral density (vBMD) measured by Quantitative Computed Tomography (QCT). The model is governed by equations describing changes in bone volume fraction (BVF), and rates of changes in bone cell populations that remove and replace bone in packets within the bone region. The DAP bone model is unique in several respects. In particular in takes former models of volume fraction changes one step higher in fidelity and separates BVF into separate equations for mineralized and osteoid volume fractions governed by a mineralization rate. This more closely follows the physiology of the remodeling unit cycles where bone is first resorbed and then followed by the action of osteoblasts to lay down collagen matrix which eventually becomes mineralized. In another respect, the modules allow the functional description of the time rate of change of other parameters and variables in the model during a computational simulation. More detailed description of the model, preliminary validation results, current limitation and caveats, and planned advancements are provided in sections 2 through 5. The DAP bone model is being developed primarily as a research tool, and not as a clinical tool like QCT. Even if it transitions to a clinical tool, it is not intended to replace QCT or any other clinical tool. Moreover, the DAP bone model does not predict bone fracture. Its purpose is to provide valuable additional data via "forward prediction" simulations for during and after spaceflight missions to gain insight on, (1) mechanisms of bone demineralization in microgravity, and (2) the volumetric changes at the various bone sites in response to in-flight and post-flight exercise countermeasures. This data can then be used as input to the Keyak [8] (or equivalent) FE analysis method to gain insight on how bone strength may change during and after flight. This information can also be useful to help optimize exercise countermeasure protocols to minimize changes in bone strength during flight, and improve regain of bone strength post-flight. To achieve this goal, the bone model will be integrated with DAP's exercise countermeasure models to simulate the effect of exercise prescriptions on preserving bone. More specifically, the model will accept loading history due to muscle and joint force on bone and produce quantified remodeling within the bone region under influence of the applied stress. Furthermore, because they tend to respond differently, the bone remodeling model includes both trabecular bone and cortical bone.

Pennline, James↗

Radiation dosimetry using three-dimensional optical random access memories

Three-dimensional optical random access memories (3D ORAMs) are a new generation of high-density data storage devices. Binary information is stored and retrieved via a light induced reversible transformation of an ensemble of bistable photochromic molecules embedded in a polymer matrix. This paper describes the application of 3D ORAM materials to radiation dosimetry. It is shown both theoretically and experimentally, that ionizing radiation in the form of heavy charged particles is capable of changing the information originally stored on the ORAM material. The magnitude and spatial distribution of these changes are used as a measure of the absorbed dose, particle type and energy. The effects of exposure on 3D ORAM materials have been investigated for a variety of particle types and energies, including protons, alpha particles and 12C ions. The exposed materials are observed to fluoresce when exposed to laser light. The intensity and the depth of the fluorescence is dependent on the type and energy of the particle to which the materials were exposed. It is shown that these effects can be modeled using Monte Carlo calculations. The model provides a better understanding of the properties of these materials. which should prove useful for developing systems for charged particle and neutron dosimetry/detector applications. c2001 Published by Elsevier Science B.V.

Non-NASA Center↗

Vibration Characteristics Determined for Stainless Steel Sandwich Panels With a Metal Foam Core for Lightweight Fan Blade Design

The goal of this project at the NASA Glenn Research Center is to provide fan materials that are safer, weigh less, and cost less than the currently used titanium alloy or polymer matrix composite fans. The proposed material system is a sandwich fan construction made up of thin solid face sheets and a lightweight metal foam core. The stiffness of the sandwich structure is increased by separating the two face sheets by the foam layer. The resulting structure has a high stiffness and lighter weight in comparison to the solid facesheet material alone. The face sheets carry the applied in-plane and bending loads (ref. 1). The metal foam core must resist the transverse shear and transverse normal loads, as well as keep the facings supported and working as a single unit. Metal foams have ranges of mechanical properties, such as light weight, impact resistance, and vibration suppression (ref. 2), which makes them more suitable for use in lightweight fan structures. Metal foams have been available for decades (refs. 3 and 4), but the difficulties in the original processes and high costs have prevented their widespread use. However, advances in production techniques and cost reduction have created a new interest in this class of materials (ref. 5). The material chosen for the face sheet and the metal foam for this study was the aerospace-grade stainless steel 17-4PH. This steel was chosen because of its attractive mechanical properties and the ease with which it can be made through the powder metallurgy process (ref. 6). The advantages of a metal foam core, in comparison to a typical honeycomb core, are material isotropy and the ease of forming complex geometries, such as fan blades. A section of a 17-4PH sandwich structure is shown in the following photograph. Part of process of designing any blade is to determine the natural frequencies of the particular blade shape. A designer needs to predict the resonance frequencies of a new blade design to properly identify a useful operating range. Operating a blade at or near the resonance frequencies leads to high-cycle fatigue, which ultimately limits the blade's durability and life. So the aim of this study is to determine the variation of the resonance frequencies for an idealized sandwich blade as a function of its face-sheet thickness, core thickness, and foam density. The finite element method is used to determine the natural frequencies for an idealized rectangular sandwich blade. The proven Lanczos method (ref. 7) is used in the study to extract the natural frequency.

Ghosn, Louis J.↗