SEARCH · Engineering Papers

Results for “matrix approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

969 records · Page 54

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Collagenase-3 binds to a specific receptor and requires the low density lipoprotein receptor-related protein for internalization

We have previously identified a specific receptor for collagenase-3 that mediates the binding, internalization, and degradation of this ligand in UMR 106-01 rat osteoblastic osteosarcoma cells. In the present study, we show that collagenase-3 binding is calcium-dependent and occurs in a variety of cell types, including osteoblastic and fibroblastic cells. We also present evidence supporting a two-step mechanism of collagenase-3 binding and internalization involving both a specific collagenase-3 receptor and the low density lipoprotein receptor-related protein. Ligand blot analysis shows that (125)I-collagenase-3 binds specifically to two proteins ( approximately 170 kDa and approximately 600 kDa) present in UMR 106-01 cells. Western blotting identified the 600-kDa protein as the low density lipoprotein receptor-related protein. Our data suggest that the 170-kDa protein is a specific collagenase-3 receptor. Low density lipoprotein receptor-related protein-null mouse embryo fibroblasts bind but fail to internalize collagenase-3, whereas UMR 106-01 and wild-type mouse embryo fibroblasts bind and internalize collagenase-3. Internalization, but not binding, is inhibited by the 39-kDa receptor-associated protein. We conclude that the internalization of collagenase-3 requires the participation of the low density lipoprotein receptor-related protein and propose a model in which the cell surface interaction of this ligand requires a sequential contribution from two receptors, with the collagenase-3 receptor acting as a high affinity primary binding site and the low density lipoprotein receptor-related protein mediating internalization.

NASA Discipline Musculoskeletal↗

NASA NTRS2001

Environmental Barrier Coatings for Silicon-Based Ceramics

Silicon-based ceramics, such as SiC fiber-reinforced SiC (SiC/SiC ceramic matrix composites (CMC) and monolithic silicon nitride (Si3N4), are prime candidates for hot section structural components of next generation gas turbine engines. Silicon-based ceramics, however, suffer from rapid surface recession in combustion environments due to volatilization of the silica scale via reaction with water vapor, a major product of combustion. Therefore, application of silicon-based ceramic components in the hot section of advanced gas turbine engines requires development of a reliable method to protect the ceramic from environmental attack. An external environmental barrier coating (EBC) is considered a logical approach to achieve protection and CP long-term stability. The first generation EBC consisted of two layers, mullite (3Al2O3-2SiO2) bond coat and yttria-stabilized zirconia (YSZ, ZrO2-8 Wt.% Y2O3) top coat. Second generation EBCs, with substantially improved performance compared with the first generation EBC, were developed in the NASA High Speed Research-Enabling Propulsion Materials (HSR-EPM) Program. The first generation EBC consisted of two layers, mullite (3Al2O3-2SiO2) bond coat and yttria-stabilized zirconia (YSZ, ZrO2-8 wt.% Y2O3) top coat. Second generation EBCs, with substantially improved performance compared with the first generation EBC, were developed in the NASA High Speed Research-Enabling Propulsion Materials (HSR-EPM) Program (5). They consist of three layers, a silicon first bond coat, a mullite or a mullite + BSAS (BaO(1-x)-SrO(x)-Al2O3-2SiO2) second bond coat, and a BSAS top coat. The EPM EBCs were applied on SiC/SiC CMC combustor liners in three Solar Turbines (San Diego, CA) Centaur 50s gas turbine engines. The combined operation of the three engines has accumulated over 24,000 hours without failure (approximately 1,250 C maximum combustor liner temperature), with the engine in Texaco, Bakersfield, CA, accumulating about 14,000 hours. As the commercialization of Si-based ceramic components in gas turbines is on the horizon, a major emphasis is placed on EBCs for two reasons. First, they are absolute necessity for the protection of Si-based ceramics from water vapor. Second, they can enable a major enhancement in the performance of gas turbines by creating temperature gradients with the incorporation of a low thermal conductivity layer. Thorough understanding of current state-of-the-art EBCs will provide the foundation upon which development of future EBCs will be based. Phase stability and thermal conductivity of EPM EBCs are published elsewhere. This paper will discuss the chemical/environmental durability and silica volatility of EPM EBCs and their impact on the coating's upper temperature limit.

Lee, Kang N.↗

NASA NTRS2016

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗

NASA NTRS1995

Cratering and penetration experiments in Teflon targets at velocities from 1 to 7 km/s

Approximately 20 sq m of protective thermal blankets, largely composed of Teflon, were retrieved from the Long Duration Exposure Facility (LDEF) after the spacecraft had spent approximately 5.7 years in space. Examination of these blankets revealed that they contained thousands of hypervelocity impact features ranging from micron-sized craters to penetration holes several millimeters in diameter. We conducted impact experiments in an effort to reproduce such features and to -- hopefully -- understand the relationships between projectile size and the resulting crater or penetration-hole diameter over a wide range of impact velocity. Such relationships are needed to derive the size- and mass-frequency distribution and flux of natural and man-made particles in low-Earth orbit. Powder propellant and light-gas guns were used to launch soda-lime glass spheres of 3.175 mm (1/8 inch) nominal diameter (Dp) into pure Teflon FEP targets at velocities ranging from 1 to 7 km/s. Target thickness (T) was varied over more than three orders of magnitude from infinite halfspace targets (Dp/T less than 0.1) to very thin films (Dp/T greater than 100). Cratering and penetration of massive Teflon targets is dominated by brittle failure and the development of extensive spall zones at the target's front and, if penetrated, the target's rear side. Mass removal by spallation at the back side of Teflon targets may be so severe that the absolute penetration-hole diameter (Dh) can become larger than that of a standard crater (Dc) at relative target thicknesses of Dp/T = 0.6-0.9. The crater diameter is infinite halfspace Teflon targets increases -- at otherwise constant impact conditions -- with encounter velocity by a factor of V0.44. In contrast, the penetration-hole size is very thin foils (Dp/T greater than 50) is essentially unaffected by impact velocity. Penetrations at target thicknesses intermediate to these extremes will scale with variable exponents of V. Our experimental matrix is sufficiently systematic and complete, up to 7 km/s, to make reasonable recommendations for the velocity-scaling of Teflon craters and penetrations. We specifically suggest that cratering behavior and associated equations dominate all impacts in which the shock-pulse duration of the projectile (tp) is shorter than that of the target (tt). We also demonstrate that each penetration hole from space-retrieved surfaces may be assigned a unique projectile size, provided an impact velocity is known or assumed. This calibration seems superior to the traditional ballistic-limit approach.

Hoerz, Friedrich↗

NASA NTRS2024

Hybrid Thermally Efficient Core (HyTEC) HyTEC Phase 1 – Advanced Aerodynamics Final Report

The objective of the Hybrid Thermally Efficient Core (HyTEC) – Advanced High Pressure Turbine (HPT) Aerodynamics project is to develop technology for a compact core that contributes to significant fuel burn reductions of 5-10% over current generation technologies. To accomplish this, the HPT is incorporating a range of aerodynamic features and technologies to improve component efficiency and provide favorable systems level trades. In particular, this project explored low solidity airfoils, advanced tip treatments, platform contouring, and advanced ceramic matrix composite (CMC) Stage 2 Nozzle (S2N) airfoils to eliminate the need for post throat cooling in a compact core environment. The maturation of these technology areas is expected to provide a significant improvement in component efficiencies, and consequently reductions in fuel burn, over the current state of the art (SoA). To mature these technologies to Technology Readiness Level (TRL) 4, a test campaign was performed that consisted of four tests at three facilities. TRL 3 testing was performed in the CW22 linear cascade at NASA Glenn for blade and nozzle technologies, TRL 4 nozzle testing was performed at GE Aerospace (GEA) Test Cell A8, and TRL 4 blade testing was performed at the Notre Dame Turbomachinery Laboratory (NDTL) using the Transonic Research Turbine (TRT) rig. Low solidity was successfully demonstrated to a TRL 4 level. Low solidity nozzles showed benefits in line with pre-project expectations, while low solidity blades were shown to have an aerodynamic penalty. Crucially, this program only considered the aerodynamic losses, and systems trades such as reductions in cooling flows are expected to continue to make low solidity blades a net positive. By successfully quantifying the aerodynamic performance in this project, these trades can be conducted to determine where in the engine architecture low solidity blades will contribute positively to system operation. Platform contouring was demonstrated to a TRL 4 level, with performance in line with the lower end of the pre-project expected range. Advanced tip treatments performance levels were indeterminate, showing the expected improvements to flow physics but with a performance level confounded by several rig issues including whirl mode induced variation in tip clearance. The elimination of post-throat cooling on the S2N was successfully demonstrated to provide a performance benefit, however that benefit was approximately half the level that was expected in pre-project predictions. Overall, the technology maturation plan for HyTEC Phase 1 was successful, bringing the suite of technologies to TRL 4.

High pressure turbine↗

NASA NTRSDate not supplied

Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions

Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.

machine learning↗

NASA NTRS2012

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↗

NASA NTRS2012

Radiation Shielding Materials Containing Hydrogen, Boron, and Nitrogen: Systematic Computational and Experimental Study

The key objectives of this study are to investigate, both computationally and experimentally, which forms, compositions, and layerings of hydrogen, boron, and nitrogen containing materials will offer the greatest shielding in the most structurally robust combination against galactic cosmic radiation (GCR), secondary neutrons, and solar energetic particles (SEP). The objectives and expected significance of this research are to develop a space radiation shielding materials system that has high efficacy for shielding radiation and that also has high strength for load bearing primary structures. Such a materials system does not yet exist. The boron nitride nanotube (BNNT) can theoretically be processed into structural BNNT and used for load bearing structures. Furthermore, the BNNT can be incorporated into high hydrogen polymers and the combination used as matrix reinforcement for structural composites. BNNT's molecular structure is attractive for hydrogen storage and hydrogenation. There are two methods or techniques for introducing hydrogen into BNNT: (1) hydrogen storage in BNNT, and (2) hydrogenation of BNNT (hydrogenated BNNT). In the hydrogen storage method, nanotubes are favored to store hydrogen over particles and sheets because they have much larger surface areas and higher hydrogen binding energy. The carbon nanotube (CNT) and BNNT have been studied as potentially outstanding hydrogen storage materials since 1997. Our study of hydrogen storage in BNNT - as a function of temperature, pressure, and hydrogen gas concentration - will be performed with a hydrogen storage chamber equipped with a hydrogen generator. The second method of introducing hydrogen into BNNT is hydrogenation of BNNT, where hydrogen is covalently bonded onto boron, nitrogen, or both. Hydrogenation of BN and BNNT has been studied theoretically. Hyper-hydrogenated BNNT has been theoretically predicted with hydrogen coverage up to 100% of the individual atoms. This is a higher hydrogen content than possible with hydrogen storage; however, a systematic experimental hydrogenation study has not been reported. A combination of the two approaches may be explored to provide yet higher hydrogen content. The hydrogen containing BNNT produced in our study will be characterized for hydrogen content and thermal stability in simulated space service environments. These new materials systems will be tested for their radiation shielding effectiveness against high energy protons and high energy heavy ions at the HIMAC facility in Japan, or a comparable facility. These high energy particles simulate exposure to SEP and GCR environments. They will also be tested in the LaRC Neutron Exposure Laboratory for their neutron shielding effectiveness, an attribute that determines their capability to shield against the secondary neutrons found inside structures and on lunar and planetary surfaces. The potential significance is to produce a radiation protection enabling technology for future exploration missions. Crew on deep space human exploration missions greater than approximately 90 days cannot remain below current crew Permissible Exposure Limits without shielding and/or biological countermeasures. The intent of this research is to bring the Agency closer to extending space missions beyond the 90-day limit, with 1 year as a long-term goal. We are advocating a systems solution with a structural materials component. Our intent is to develop the best materials system for that materials component. In this Phase I study, we have shown, computationally, that hydrogen containing BNNT is effective for shielding against GCR, SEP, and neutrons over a wide range of energies. This is why we are focusing on hydrogen containing BNNT as an innovative advanced concept. In our future work, we plan to demonstrate, experimentally, that hydrogen, boron, and nitrogen based materials can provide mechanically strong, thermally stable, structural materials with effective radiation shielding against GCR, SEP, and neutrons.

Radiation↗

NASA NTRS2007

Improved Underwater Excitation-Emission Matrix Fluorometer

A compact, high-resolution, two-dimensional excitation-emission matrix fluorometer (EEMF) has been designed and built specifically for use in identifying and measuring the concentrations of organic compounds, including polluting hydrocarbons, in natural underwater settings. Heretofore, most EEMFs have been designed and built for installation in laboratories, where they are used to analyze the contents of samples collected in the field and brought to the laboratories. Because the present EEMF can be operated in the field, it is better suited to measurement of spatially and temporally varying concentrations of substances of interest. In excitation-emission matrix (EEM) fluorometry, fluorescence is excited by irradiating a sample at one or more wavelengths, and the fluorescent emission from the sample is measured at multiple wavelengths. When excitation is provided at only one wavelength, the technique is termed one-dimensional (1D) EEM fluorometry because the resulting matrix of fluorescence emission data (the EEM) contains only one row or column. When excitation is provided at multiple wavelengths, the technique is termed two-dimensional (2D) EEM fluorometry because the resulting EEM contains multiple rows and columns. EEM fluorometry - especially the 2D variety - is well established as a means of simultaneously detecting numerous dissolved and particulate compounds in water. Each compound or pool of compounds has a unique spectral fluorescence signature, and each EEM is rich in information content, in that it can contain multiple fluorescence signatures. By use of deconvolution and/or other mixture-analyses techniques, it is often possible to isolate the spectral signature of compounds of interest, even when their fluorescence spectra overlap. What distinguishes the present 2D EEMF over prior laboratory-type 2D EEMFs are several improvements in packaging (including a sealed housing) and other aspects of design that render it suitable for use in natural underwater settings. In addition, the design of the present 2D EEMF incorporates improvements over the one prior commercial underwater 2D EEMF, developed in 1994 by the same company that developed the present one. Notable advanced features of the present EEMF include the following: 1) High sensitivity and spectral resolution are achieved by use of an off-the-shelf grating spectrometer equipped with a sensor in the form of a commercial astronomical- grade 256 532-pixel charge-coupled-device (CCD) array. 2) All of the power supply, timing, control, and readout circuits for the illumination source and the CCD, ancillary environmental monitoring sensors, and circuitry for controlling a shutter or filter motor are custom-designed and mounted compactly on three circuit boards below a fourth circuit board that holds the CCD (see figure). 3) The compactness of the grating spectrometer, CCD, and circuit assembly makes it possible to fit the entire instrument into a compact package that is intended to be maneuverable underwater by one person. 4) In mass production, the cost of the complete instrument would be relatively low - estimated at approximately $30,000 at 2005 prices.

Moore, Casey↗

NASA NTRS2003

NASA Glenn/AADC-Rolls Royce Collaborated to Measure Erosion Resistance on Coated Polymer Matrix Composites

Polymer matrix composites (PMCs) are increasingly used in aerospace and automotive applications because of their light weight and high strength-to-weight ratio relative to metals. However, a major drawback of PMCs is poor abrasion resistance, which restricts their use, especially at high temperatures. Simply applying a hard coating on PMCs to improve abrasion and erosion resistance is not effective since coating durability is short lived (ref. 1). Generally, PMCs have higher coefficients of thermal expansion than metallic or ceramic coatings have, and coating adhesion suffers because of poor interfacial adhesion strength. One technique commonly used to improve coating adhesion or durability is the use of bond coats that are interleaved between a coating and a substrate with vastly different coefficients of thermal expansion. An example of this remedy is the use of bondcoats for ceramic thermal barrier coatings on metallic turbine components (ref. 2). Prior collaborative research between the NASA Glenn Research Center and the Allison Advanced Development Company (AADC) demonstrated that bond coats sandwiched between PMCs and high-quality plasma-sprayed, erosion-resistant coatings substantially improved the erosion resistance of PMCs (ref. 3). One unresolved problem in this earlier collaboration was that there was no easy, accurate way to measure the coating erosion wear scar. Coating wear was determined by both profilometry and optical microscopy. Both techniques are time consuming. Wear measurement by optical microscopy requires sample destruction and does not provide a comprehensive measure of the entire wear volume. An even more subtle, yet critical, problem is that these erosion coatings contain two or more materials with different densities. Therefore, simply measuring specimen mass loss before and after erosion will not provide an accurate gauge for coating and/or substrate volume loss. By using a noncontact technique called scanning optical interferometry, which was recently developed at Glenn, researchers can accurately determine the wear performance of erosion-coated PMCs while preserving the sample. An example of this interferometry technique is shown in the preceding figure for an erosion-coated inlet guide vane from a Rolls Royce AE3007 regional gas turbine jet engine. Erosion was conducted with coated and uncoated PMC vanes, with the abrasive material moving at a velocity of 229 m/s at impingement angles of 20 and 90 degrees. The coatings for PMCs remarkably reduced the erosion volume loss by a factor of approximately 10. Currently, several erosion coatings for PMCs are being compared and downselected for engine testing at Rolls Royce.

Miyoshi, Kazuhisa↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

NASA NTRS1992

Control-structure-thermal interactions in analysis of lunar telescopes

The lunar telescope project was an excellent model for the CSTI study because a telescope is a very sensitive instrument, and thermal expansion or mechanical vibration of the mirror assemblies will rapidly degrade the resolution of the device. Consequently, the interactions are strongly coupled. The lunar surface experiences very large temperature variations that range from approximately -180 C to over 100 C. Although the optical assemblies of the telescopes will be well insulated, the temperature of the mirrors will inevitably fluctuate in a similar cycle, but of much smaller magnitude. In order to obtain images of high quality and clarity, allowable thermal deformations of any point on a mirror must be less than 1 micron. Initial estimates indicate that this corresponds to a temperature variation of much less than 1 deg through the thickness of the mirror. Therefore, a lunar telescope design will most probably include active thermal control, a means of controlling the shape of the mirrors, or a combination of both systems. Historically, the design of a complex vehicle was primarily a sequential process in which the basic structure was defined without concurrent detailed analyses or other subsystems. The basic configuration was then passed to the different teams responsible for each subsystem, and their task was to produce a workable solution without requiring major alterations to any principal components or subsystems. Consequently, the final design of the vehicle was not always the most efficient, owing to the fact that each subsystem design was partially constrained by the previous work. This procedure was necessary at the time because the analysis process was extremely time-consuming and had to be started over with each significant alteration of the vehicle. With recent advances in the power and capacity of small computers, and the parallel development of powerful software in structural, thermal, and control system analysis, it is now possible to produce very detailed analyses of intermediate designs in a much shorter period of time. The subsystems can thus be designed concurrently, and alterations in the overall design can be quickly adopted into each analysis; the design becomes an iterative process in which it is much easier to experiment with new ideas, configurations, and components. Concurrent engineering has the potential to produce efficient, highly capable designs because the effect of one subystem on another can be assessed in much more detail at a very early point in the program. The research program consisted of several tasks: scale a prototype telescope assembly to a 1 m aperture, develop a model of the telescope assembly by using finite element (FEM) codes that are available on site, determine structural deflections of the mirror surfaces due to the temperature variations, develop a prototype control system to maintain the proper shape of the optical elements, and most important of all, demonstrate the concurrent engineering approach with this example. In addition, the software used for the finite element models and thermal analysis was relatively new within the Program Development Office and had yet to be applied to systems this large or complex; understanding the software and modifying it for use with this project was also required. The I-DEAS software by Structural Dynamics Research Corporation (SDRC) was used to build the finite element models, and TMG developed by Maya Heat Transfer Technologies, Ltd. (which runs as an I-DEAS module) was used for the thermal model calculations. All control system development was accomplished with MATRIX(sub X) by Integrated Systems, Inc.

Thompson, Roger C.↗

NASA NTRS2009

Modeling and Simulation of Variable Mass, Flexible Structures

The advent of the new Ares I launch vehicle has highlighted the need for advanced dynamic analysis tools for variable mass, flexible structures. This system is composed of interconnected flexible stages or components undergoing rapid mass depletion through the consumption of solid or liquid propellant. In addition to large rigid body configuration changes, the system simultaneously experiences elastic deformations. In most applications, the elastic deformations are compatible with linear strain-displacement relationships and are typically modeled using the assumed modes technique. The deformation of the system is approximated through the linear combination of the products of spatial shape functions and generalized time coordinates. Spatial shape functions are traditionally composed of normal mode shapes of the system or even constraint modes and static deformations derived from finite element models of the system. Equations of motion for systems undergoing coupled large rigid body motion and elastic deformation have previously been derived through a number of techniques [1]. However, in these derivations, the mode shapes or spatial shape functions of the system components were considered constant. But with the Ares I vehicle, the structural characteristics of the system are changing with the mass of the system. Previous approaches to solving this problem involve periodic updates to the spatial shape functions or interpolation between shape functions based on system mass or elapsed mission time. These solutions often introduce misleading or even unstable numerical transients into the system. Plus, interpolation on a shape function is not intuitive. This paper presents an approach in which the shape functions are held constant and operate on the changing mass and stiffness matrices of the vehicle components. Each vehicle stage or component finite element model is broken into dry structure and propellant models. A library of propellant models is used to describe the distribution of mass in the fuel tank or Solid Rocket Booster (SRB) case for various propellant levels. Based on the mass consumed by the liquid engine or SRB, the appropriate propellant model is coupled with the dry structure model for the stage. Then using vehicle configuration data, the integrated vehicle model is assembled and operated on by the constant system shape functions. The system mode shapes and frequencies can then be computed from the resulting generalized mass and stiffness matrices for that mass configuration. The rigid body mass properties of the vehicle are derived from the integrated vehicle model. The coupling terms between the vehicle rigid body motion and elastic deformation are also updated from the constant system shape functions and the integrated vehicle model. This approach was first used to analyze variable mass spinning beams and then prototyped into a generic dynamics simulation engine. The resulting code was tested against Crew Launch Vehicle (CLV-)class problems worked in the TREETOPS simulation package and by Wilson [2]. The Ares I System Integration Laboratory (SIL) is currently being developed at the Marshall Space Flight Center (MSFC) to test vehicle avionics hardware and software in a hardware-in-the-loop (HWIL) environment and certify that the integrated system is prepared for flight. The Ares I SIL utilizes the Ares Real-Time Environment for Modeling, Integration, and Simulation (ARTEMIS) tool to simulate the launch vehicle and stimulate avionics hardware. Due to the presence of vehicle control system filters and the thrust oscillation suppression system, which are tuned to the structural characteristics of the vehicle, ARTEMIS must incorporate accurate structural models of the Ares I launch vehicle. The ARTEMIS core dynamics simulation models the highly coupled nature of the vehicle flexible body dynamics, propellant slosh, and vehicle nozzle inertia effects combined with mass and flexible body properties that vary significant with time during the flight. All forces that act on the vehicle during flight must be simulated, including deflected engine thrust force, spatially distributed aerodynamic forces, gravity, and reaction control jet thrust forces. These forces are used to excite an integrated flexible vehicle, slosh, and nozzle dynamics model for the vehicle stack that simulates large rigid body translations and rotations along with small elastic deformations. Highly effective matrix math operations on a distributed, threaded high-performance simulation node allow ARTEMIS to retain up to 30 modes of flex for real-time simulation. Stage elements that separate from the stack during flight are propagated as independent rigid six degrees of freedom (6DOF) bodies. This paper will present the formulation of the resulting equations of motion, solutions to example problems, and describe the resulting dynamics simulation engine within ARTEMIS.

Tobbe, Patrick A.↗

NASA NTRSDate not supplied

Using X-Ray Computed Tomography to Catalog Rock Fragments in Apollo Drive Tube 73002

Overview: The Apollo missions collected 382 kg of rock, regolith, and core samples from six locations on the nearside of the Moon. Approximately 84% by mass of the Apollo collection remains in pristine condition within the curation facility at Johnson Space Center (i.e., never allocated, continuously stored in dry-N2 purged cabinets, exposure history restricted to Teflon, stainless steel, and Al-metal). Although most Apollo samples have been well characterized, there are several types of samples that have remained wholly or largely unstudied since their return, and/or that have been cu-rated under special conditions, e.g., frozen samples, samples stored in He-purged environment, and previously unopened drive tubes. NASA solicited proposals for the Apollo Next Generation Sample Analysis Pro-gram (ANGSA), and 9 teams were selected to study a subset of the unopened and frozen samples [1]. The first sample opened as part of the ANGSA pro-gram was drive tube 73002. This was originally a ~30 cm long, 4 cm diameter drive tube collected on a land-slide deposit near Lara Crater at the Apollo 17 landing site. It was part of a ~60 cm long double drive tube collected, and the bottom half of the tube (73001) was sealed under vacuum on the Moon [2]. Prior to opening sample 73002, the sample was imaged with a high resolution X-ray Computed Tomography (XCT) scan of the entire tube [3], which provided invaluable information during the dissection process [4]. In addition to the pre-dissection XCT scans, individual >4 mm particles were separated from the 73002 regolith during processing and scanned at high resolution by XCT. Here we pre-sent the initial lithologic classification of 134 individual >4 mm rock fragments separated from the 73002 core during the dissection process. Methodology: Drive tube 73002 was manually dissected in 0.5 cm depth intervals in three passes (Fig. 1) [4,5]. Each interval from pass 1 and 2 was sieved to <1 mm and >1 mm size fractions, and >1 mm particles were further manually subdivided into 1-2, 2-4, 4-10, and >10 mm size fractions. Pass 3 was not sieved, but >10 mm clasts were separated manually. Each 4-10 mm and >10 mm fragment was individually weighed, triply bagged in Teflon, and scanned by XCT. There are 60 rock fragments in Pass 1, 64 rock fragments in Pass 2 (from the 4-10 mm and >10 mm size fractions), and 8 rock fragments from Pass 3 (>10 mm size fraction). Each individually bagged rock fragment was scanned using the 180 kV nano-focus transmission source on the Nikon XTH 320 XCT system at NASA Johnson Space Center [6]. Scanning conditions varied considerably for individual particles, in large part be-cause of the large variation in size (0.008-19.623 g). All scans fell within the following range of scan conditions: 2.8-20.6 um voxel size; 90-155 kV voltage; 18-39 uA current; 1891-3141 projections; and 902-2000 slices. Results and Discussion: The 132 rock fragments from sample 73002 fall into the following general cat-gories: agglutinates (n = 6); basalts (13); impact melts (5); impact melt breccias (IMB; 42); regolith breccias (62); and soil breccias (4); see Figure 2 for representative examples of each lithology. Within most of these broad lithologic groups are recognizable subgroups. For example, a significant portion of regolith breccia fragments contain some agglutinate-like glass (n = 9) or are dilithologic (7) because they contain a single large clast (~50% or more by volume). Subgroups can be based on similarities to previously identified lunar lithologies, such as high-Ti basalts (9) and VLT basalts (4; Fig. 3), or based on commonly seen features within the fragments, such as poikilitic ilmenite IMB (10), ilmenite-lath IMB (12), or vesicular IMB (11) “groups”. Particles in the same “group” are not necessarily intended to be genetically related, but rather identify particles that are similar and that follow-up studies can classify in more detail [7]. Conclusion: Identification of lithologies based on XCT is a powerful tool, but a more absolute classification will sometimes require additional textural information from thin sections (e.g., glassy-matrix regolith vs. impact-melt vs. granulitic breccia) or quantitative mineral compositions (e.g., basalt vs. monomict breccia).

R A Zeigler↗