Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field 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 649 records · Page 36

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

NASA Tech Briefs, July 2009

Topics covered include: Dual Cryogenic Capacitive Density Sensor; Hail Monitor Sensor; Miniature Six-Axis Load Sensor for Robotic Fingertip; Improved Blackbody Temperature Sensors for a Vacuum Furnace; Wrap-Around Out-the-Window Sensor Fusion System; Wide-Range Temperature Sensors with High-Level Pulse Train Output; Terminal Descent Sensor Simulation; A Robust Mechanical Sensing System for Unmanned Sea Surface Vehicles; Additive for Low-Temperature Operation of Li-(CF)n Cells; Li/CFx Cells Optimized for Low-Temperature Operation; Number Codes Readable by Magnetic-Field-Response Recorders; Determining Locations by Use of Networks of Passive Beacons; Superconducting Hot-Electron Submillimeter-Wave Detector; Large-Aperture Membrane Active Phased-Array Antennas; Optical Injection Locking of a VCSEL in an OEO; Measuring Multiple Resistances Using Single-Point Excitation; Improved-Bandwidth Transimpedance Amplifier; Inter-Symbol Guard Time for Synchronizing Optical PPM; Novel Materials Containing Single-Wall Carbon Nanotubes Wrapped in Polymer Molecules; Light-Curing Adhesive Repair Tapes; Thin-Film Solid Oxide Fuel Cells; Zinc Alloys for the Fabrication of Semiconductor Devices; Small, Lightweight, Collapsible Glove Box; Radial Halbach Magnetic Bearings; Aerial Deployment and Inflation System for Mars Helium Balloons; Steel Primer Chamber Assemblies for Dual Initiated Pyrovalves; Voice Coil Percussive Mechanism Concept for Hammer Drill; Inherently Ducted Propfans and Bi-Props; Silicon Nanowire Growth at Chosen Positions and Orientations; Detecting Airborne Mercury by Use of Gold Nanowires; Detecting Airborne Mercury by Use of Palladium Chloride; Micro Electron MicroProbe and Sample Analyzer; Nanowire Electron Scattering Spectroscopy; Electron-Spin Filters Would Offer Spin Polarization Greater than 1; Subcritical-Water Extraction of Organics from Solid Matrices; A Model for Predicting Thermoelectric Properties of Bi2Te3; Integrated Miniature Arrays of Optical Biomolecule Detectors; A Software Rejuvenation Framework for Distributed Computing; Kurtosis Approach to Solution of a Nonlinear ICA Problem; Robust Software Architecture for Robots; R4SA for Controlling Robots; Bio-Inspired Neural Model for Learning Dynamic Models; Evolutionary Computing Methods for Spectral Retrieval; Monitoring Disasters by Use of Instrumented Robotic Aircraft; Complexity for Survival of Living Systems; Using Drained Spacecraft Propellant Tanks for Habitation; Connecting Node; and Electrolytes for Low-Temperature Operation of Li-CFx Cells.

Source record↗

Alignment-Insensitive Lower-Cost Telescope Architecture

This architecture features an active wavefront sensing and control scheme along with methods for measuring the relative positions of the primary to aft optics, such as the secondary mirror, and should enable larger and cheaper telescope architectures needed for future applications. A wavefront source/sensor is placed at the center of curvature of the primary mirror. The system provides continuous light onto a primary mirror that is retro-reflected onto itself. This allows the wavefront controller to constantly update the positions of the primary mirror segments (or deformable mirror actuators). Another function of this innovation involves using a concave mirror on the back of the secondary mirror (or other aft optic) that has the same center-of-curvature location (in defocus) as the primary mirror. The two return beams can be aligned next to each other on a detector, or radially on top of each other. This provides a means with which to measure the relative position of the primary to the secondary (or other aft optics), thus allowing for the removal of misalignment of the center-of-curvature source/sensor (meaning it doesn't need precision placement) and also provides a means with which to monitor the relative alignment over time. This innovation does not require extremely good thermal stability on the primary mirror and can thus be used in any thermal environment and with cheaper materials. In addition to this, the architecture lets one phase (or align) the primary mirror independent of whether a star or scene is in the field. The segmented, spherical primary allows for cost-effective three-meter class (e.g. Midex and Discovery) missions as well as enabling 30-meter telescope solutions that can be manufactured in a reasonable amount of time. The continuous wavefront sensing and control architecture enables missions for low-Earth-orbit.

Feinberg, Lee↗

The AgMIP Coordinated Climate-Crop Modeling Project (C3MP): Methods and Protocols

Climate change is expected to alter a multitude of factors important to agricultural systems, including pests, diseases, weeds, extreme climate events, water resources, soil degradation, and socio-economic pressures. Changes to carbon dioxide concentration ([CO2]), temperature, and water (CTW) will be the primary drivers of change in crop growth and agricultural systems. Therefore, establishing the CTW-change sensitivity of crop yields is an urgent research need and warrants diverse methods of investigation. Crop models provide a biophysical, process-based tool to investigate crop responses across varying environmental conditions and farm management techniques, and have been applied in climate impact assessment by using a variety of methods (White et al., 2011, and references therein). However, there is a significant amount of divergence between various crop models' responses to CTW changes (Rotter et al., 2011). While the application of a site-based crop model is relatively simple, the coordination of such agricultural impact assessments on larger scales requires consistent and timely contributions from a large number of crop modelers, each time a new global climate model (GCM) scenario or downscaling technique is created. A coordinated, global effort to rapidly examine CTW sensitivity across multiple crops, crop models, and sites is needed to aid model development and enhance the assessment of climate impacts (Deser et al., 2012). To fulfill this need, the Coordinated Climate-Crop Modeling Project (C3MP) (Ruane et al., 2014) was initiated within the Agricultural Model Intercomparison and Improvement Project (AgMIP; Rosenzweig et al., 2013). The submitted results from C3MP Phase 1 (February 15, 2013-December 31, 2013) are currently being analyzed. This chapter serves to present and update the C3MP protocols, discuss the initial participation and general findings, comment on needed adjustments, and describe continued and future development. AgMIP aims to improve substantially the climate, crop, and economic simulation tools that are used to characterize the agricultural sector, to assess future world food security under changing climate conditions, and to enhance adaptation capacity both globally and regionally. To understand better and improve the modeled crop responses, AgMIP has conducted detailed crop model intercomparisons at closely observed field sites for wheat (Asseng et al., 2013), rice (Li et al., in review), maize (Bassu et al., 2014), and sugarcane (Singels et al., 2013). A coordinated modeling exercise was one of the original motivations for AgMIP, and C3MP provides rapid estimation of crop responses to CO2, water, and temperature (CTW) changes, adding dimension and insight into the crop model intercomparisons, while facilitating interactions within the global community of modelers. C3MP also contributes a fast-track, multi-model climate sensitivity assessment for the AgMIP climate and crop modeling teams on Research Track 2 (Fig. 1), which seeks to understand the impact of projected climatic changes on crop production and food security (Rosenzweig et al., 2013; Ruane et al., 2014).

agriculture↗

Construction of Polarimetric Radar-Based Reference Rain Maps for the Iowa Flood Studies Campaign

The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.

Radar↗

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗

Analysis of One-Way Laser Ranging Data to LRO, Time Transfer and Clock Characterization

We processed and analyzed one-way laser ranging data from International Laser Ranging Service ground stations to NASA's Lunar Reconnaissance Orbiter (LRO), obtained from June 13, 2009 until September 30, 2014. We pair and analyze the one-way range observables from station laser fire and spacecraft laser arrival times by using nominal LRO orbit models based on the GRAIL gravity field. We apply corrections for instrument range walk, as well as for atmospheric and relativistic effects. In total we derived a tracking data volume of approximately 3000 hours featuring 64 million Full Rate and 1.5 million Normal Point observations. From a statistical analysis of the dataset we evaluate the experiment and the ground station performance. We observe a laser ranging measurement precision of 12.3 centimeters in case of the Full Rate data which surpasses the LOLA (Lunar Orbiting Laser Altimeter) timestamp precision of 15 centimeters. The averaging to Normal Point data further reduces the measurement precision to 5.6 centimeters. We characterized the LRO clock with fits throughout the mission time and estimated the rate to 6.9 times10 (sup -8), the aging to 1.6 times 10 (sup -12) per day and the change of aging to 2.3 times 10 (sup -14) per day squared over all mission phases. The fits also provide referencing of onboard time to the TDB (Barycentric Dynamical Time) time scale at a precision of 166 nanoseconds over two and 256 nanoseconds over all mission phases, representing ground to space time transfer. Furthermore we measure ground station clock differences from the fits as well as from simultaneous passes which we use for ground to ground time transfer from common view observations. We observed relative offsets ranging from 33 to 560 nanoseconds and relative rates ranging from 2 times 10 (sup -13) to 6 times 10 (sup -12) between the ground station clocks during selected mission phases. We study the results from the different methods and discuss their applicability for time transfer.

Bauer, S.↗

A New and General Formulation of the Parametric HFGMC Micromechanical Method for Three-Dimensional Multi-Phase Composites

The recent two-dimensional (2-D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3-D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2-D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3-D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2-D parametric and orthogonal HFGMC are special cases of the present 3-D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3-D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3-D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, discontinuous aligned short fiber. A 3-D repeating unit-cell for foam material composite is simulated.

Haj-Ali, Rami↗

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗

Heating of the nighttime D region by very low frequency transmitters

VLF signals propagating in the Earth-ionosphere waveguide are used to probe the heated nighttime D region over three U.S. Navy very low frequency (VLF,3-30 kHz) transmitters. Ionospheric cooling and heating are observed when a transmitter turns off and on in the course of normal operations. Heating by the 24.0-kHz NAA transmitter in Cutler, Maine, (1000 kW radiated power) was observed by this method in 41 of 52 off/on episodes during December 1992, increasing the amplitude and retarding the phase of the 21.4-kHz NSS probe wave propagating from Annapolis, Maryland, to Gander, Newfoundland, by as much as 0.84 dB and 5.3 deg, respectively. In 6 of these 41 episodes, the amplitude of the 28.5-kHz NAU probe wave propagating from Puerto Rico to Gander was also perturbed by as much as 0.29 dB. The latter observations were unexpected due to the greater than 770 km distance between NAA and the NAU-Gander great circle path. Heating by the NSS (21.4 kHz, 265 kW) and NLK (24.8 kHz, 850 kW) transmitters was observed serendipitously in data from earlier measurements of the amplitudes of VLF signals propagating in the Earth-ionosphere waveguide. A three-dimensional model of wave absorption and electron heating in a magnetized, weakly ionized plasma is used to calculate the extent nad shape of the collision frequency (i.e., electron temperature) enhancement above a VLF transmitter. The enhancements are annular, with a geomagnetic north-south asymmetry and a radius at the outer half-maximum of the collision frequency enhancement of about 150 km. Heating by the NAA transmitter is predicted to increase the nighttime D region electron temperature by as much as a factor of 3. The calculated changes in the D region conductivity are used in a three-dimensional model of propagation in the Earth-ionosphere wavelength to predict the effect of the heated patch on a subionospheric VLF probe wave. The range of predicted scattered field amplitudes is in general consistent with the observed signal perturbations. Discrepanices in the predictions are attributed to lack of knowledge of the D region electron density profile along the probe wave great circle paths.

Rodriguez, Juan V.↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Solar System Escape Architecture for Revolutionary Science

Our most distance spacecraft, Voyager 1 and Voyager 2, recently reached the boundary of the solar system known as the Heliopause. Beyond this boundary lies interstellar space, and, at tremendous distance, the stars. The Voyager spacecraft have shown us that the Heliopause is a dynamic, rapidly changing environment, varying greatly as the fields and effects emanating from the sun are met by a host of forces pushing in from our nearby galactic neighborhood. This is the region we aim to explore. Specifically we aim to develop a mission architecture and spacecraft concept capable of reaching the Heliopause (approximately100 AU (Astronomical Units)) region in multiple directions of interest, with a variety of scientific instrument suites, within a reasonable timeframe (about 15 years.) This report details our Phase 1 work including science goal definition, trajectory planning, technology research, mission planning, instrument selection, spacecraft design, and more. We present an architecture employing a 250-meter by 250-meter solar sail with a 175 kilogram spacecraft. This spacecraft could carry a variety of instrument suites depending on the destination, time, and other factors. The spin-stabilized spacecraft would powered by a small radioisotope power system (about 20 watts), makes use of an optical communication system, and carries 3 extensible booms for instrument accommodation as visible in the sketch and rendering below. We conclude that it would be reasonable to implement the architecture described herein in the 2020-2030 timeframe based on likely near-term technological and material advances. We are confident that cruise times of around 15 years could be achieved to science-rich destinations. This report details the approach we used to come to these conclusions and identifies some of the steps along the way. We utilized a variety of methodologies involving our full capability at JPL (Jet Propulsion Laboratory) - from small point studies to large collaborative engineering processes. This paper also describes potential industry partnerships and costing information, and includes the methods used to account for our project’s significant deviation from any prior mission in the NASA costing tool database. Outreach and public engagement is an important part of NASA activities. We engaged in several different forms of this outreach and describe them in this report. We conclude with some open questions and a summary of activities we hope to engage in with Phase 2 funding. The entire team would like to thank the NIAC (NASA Innovative Advanced Concepts) program for the opportunity to investigate this fascinating and challenging concept. The Voyager program gives us tantalizing clues as to the nature and behavior of the Heliopause – the boundary between the sun’s influence and the interstellar medium. This information comes from forty-year old instruments designed to study the outer planets. A targeted Heliopause investigation would give insight into the formation of the solar system, the role of the sun in the local interstellar neighborhood, and contribute to human exploration planning by helping to predict periods of low galactic cosmic ray (GCR) penetration into the inner solar system.

Spacecraft↗

Numerical Simulation of One- and Two-Phase Flows in Propulsion Systems

In this report, we present some results of problems investigated during joint research between the Hampton University Fluid Mechanics and Acoustics Laboratory (HU/FM&AL), NASA GRC, and the LaRC Hyper-X Program. This work is supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The main areas of current scientific interest of the HU/FM&AL include an investigation of the proposed and patented advanced methods for aircraft engine thrust and noise benefits. These methods are based on nontraditional 3D corrugated and composite nozzle, inlet, propeller and screw designs such as a Bluebell and Telescope nozzle, Mobius-shaped screw, etc. This is the main subject of our other projects, of which one is presented at the current conference. Here we analyze additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. Such additional approaches are: (1) to add some solid, fluid, or gas mass at discrete locations to the main supersonic gas stream to minimize the negative influence of strong shock waves formed in propulsion systems. This mass addition may be accompanied by heat addition to the main stream as a result of the fuel combustion or by cooling of this stream as a result of the liquid mass evaporation and boiling; (2) Use of porous or permeable nozzles and additional shells at the nozzle exit for preliminary cooling of the hot jet exhaust and pressure compensation for off-design conditions (so-called continuous ejector with small mass flow rate); and (3) to propose and analyze new effective methods of fuel injection into the flow stream in air-breathing engines. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations, and experimental tests at the Hampton University campus, NASA, and IM/MSU. The main outcomes during this reporting period are: (1) Publications: The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01-3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference being held on 08-11 July, in Salt Lake City, UT; (2) Grants and proposals: The HU/FM&AL was awarded the NASA grant NAG3-2495 in October 2000 and the CRDF award was granted to the NASA GRC-HU FM&AL and IM/MSU (Russia) in July 2000. A solicited proposal was submitted for the NASA NRA-01GRC-02 competition and two unsolicited proposals to NASA are in preparation; (3) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc. Such investigations have been conducted for three main problems: (a) Combustion efficiency optimization in the half-duct combustor system; (b) Drag reduction effects for blunt bodies with solid needles; and (c) Solid particle injection from the butt-end against a supersonic flow. The NASA CFL3D, HU/FM&AL, and IM/MSU GODUNOV codes were used and modified for solution of these problems. The codes are based on full Euler and Navier-Stokes solvers with and without nonequilibrium oxygen-nitrogen and air-hydrogen chemical reactions in laminar and turbulent gas flow regimes; (4) Experimental tests: In the IM/MSU supersonic wind tunnel, experimental tests were conducted with different number, location, and geometric parameters of solid needles mounted at the front of the butt-end in supersonic flow. Optimal parameters were determined that provide minimal butt-end drag in the stationary flow regime. Experimental and numerical simulation results are in good agreement; (5) Student Research Activity: Involvement of one graduate and two undergraduate students as research assistants in the current project.

Gilinsky, Mikhail↗

Laboratory Investigation of the Effect of Venusian Weathering on Mineral Spectra

Introduction: The recent selection of two missions to Venus has renewed the importance of deter-mining weathering reactions between minerals and the Venusian atmosphere, and the spectral signatures of minerals before, during, and after these reactions. The rate at which weathering reactions progress also constrain show long unstable minerals will be present on the surface (e.g., [1]), and enables the use of mineralogy as a constraint on surface age(e.g., [2-3]).In order to gain an understanding of how mineral compositions and spectra change with weathering, we have begun conducting experiments in a 1 atm experimental setup at Wesleyan University. This setup exposes minerals to the temperature and most abundant gases of the Venus atmosphere (CO2, SO2, N2).We conducted initial experiments using biotite, calcite, and montmorillonite in order to test our methodology with minerals that may be relevant to recording the history of water on Venus. Methods: Experiments were conducted in Thermo Fisher Scientific Lindberg/Blue M Mini-Mite horizontal tube furnaces at Wesleyan University. Experiments used natural mineral chips and powders, and were conducted at 1 atmosphere and 460 °C under pre-mixed gases provided by Air Gas(Table 1). The furnaces are set up in a flow through configuration so that solid samples are exposed to a fixed gas composition, and quartz glass process tubes were used in all experiments (1/4” diameter for experiment V4, ½” diameter for all others). These conditions were maintained for the du-rations listed in Table 1, at which point the furnace was turned off with gas flowing until the sample was cool enough to extract under N2and be placed in a desiccator for storage. Run products were carbon coated and examined using a Hitachi SU5000 Field Emission Gun Scanning Electron Microscope (SEM) equipped with an EDAX Octane Pro EDS detector located at Wesley-an University. Visible-Near Infrared Spectroscopy (VNIR) analysis was performed on powdered samples under a nitrogen atmosphere using an ASD Fieldspec Proover the 350-2500 nm range. Powdered mineral samples were milled to a particle size of < 45μm and were spiked with an internal standard (Al2O3, corundum) to obtain quantitative mineralogy. Samples were analyzed using a Panalytical X’Pert pro X-ray Diffractometer (XRD), with an X’Celerator high speed detector and Co Kα radiation, with data collected at a step size of 0.02 ̊/minute step counting rate from 2 to 80 degrees 2θ at 45 mA/40kV in the X-ray Diffraction Laboratory located at NASA Johnson Space Center. Materials Data Inc (MDI) software suite, JadeTMv9 was used for Rietveld refinement to determine phase abundances and mineral identification by comparing XRD patterns to International Center for Diffraction Data (ICDD) database patterns. Table 1: Experimental matrix. Experiment Name Duration (Days)Gas Com-position (trace gas)MineralsV487SO2/N2(1.4%)Montmorillonite, biotiteV519SO2/N2(1.4%)CalciteV619CO2/SO2/N2(1.4% SO2, 2.1% N2)Montmorillonite, biotiteV828SO2/N2(1.4%)Biotite, calcite Results: Calcite. In both experiments, calcite was exposed to the SO2/N2gas mixture, and in both experimental run products, XRD analysis detected anhydrite, which is consistent with EDS measurements conducted in the SEM. The XRD analyses show greater amounts of anhydrite present after 28 days than 19, suggesting the calcite reaction progressed further given longer duration. Grain surface morphology as seen in the SEM also shows secondary mineral growth (Fig. 1). VNIR spectra show no change, as expected since anhydrite lacks spectral features in this wavelength range. Montmorillonite. Montmorillonite was exposed to two different gas mixtures, the SO2/N2mix and CO2/SO2/N2mix over different durations (87 and 19 days, respectively). VNIR analyses of both run products show a reduction in the 1441and1910 nm water features as well as a shift of the 1411 and 2011 nm features to shorter wavelengths that may indicate re-structuring in the crystal lattice and the production of amorphous phases [4]. XRD results are consistent with this, showing a shift of the 001 peak from 15Åto 10Åin both run products. The amount of X-ray amorphous material in the montmorillonite run products was greater than that present in the unreacted clay, being the greatest in the 87 day V4 experiment. No other secondary phases were detected in the run products, however sulfur was present in EDS analyses of powder samples.

A R Santos↗

A Method for Calculating the Probability of Successfully Completing a Rocket Propulsion Ground Test

Propulsion ground test facilities face the daily challenges of scheduling multiple customers into limited facility space and successfully completing their propulsion test projects. Due to budgetary and schedule constraints, NASA and industry customers are pushing to test more components, for less money, in a shorter period of time. As these new rocket engine component test programs are undertaken, the lack of technology maturity in the test articles, combined with pushing the test facilities capabilities to their limits, tends to lead to an increase in facility breakdowns and unsuccessful tests. Over the last five years Stennis Space Center's propulsion test facilities have performed hundreds of tests, collected thousands of seconds of test data, and broken numerous test facility and test article parts. While various initiatives have been implemented to provide better propulsion test techniques and improve the quality, reliability, and maintainability of goods and parts used in the propulsion test facilities, unexpected failures during testing still occur quite regularly due to the harsh environment in which the propulsion test facilities operate. Previous attempts at modeling the lifecycle of a propulsion component test project have met with little success. Each of the attempts suffered form incomplete or inconsistent data on which to base the models. By focusing on the actual test phase of the tests project rather than the formulation, design or construction phases of the test project, the quality and quantity of available data increases dramatically. A logistic regression model has been developed form the data collected over the last five years, allowing the probability of successfully completing a rocket propulsion component test to be calculated. A logistic regression model is a mathematical modeling approach that can be used to describe the relationship of several independent predictor variables X(sub 1), X(sub 2),..,X(sub k) to a binary or dichotomous dependent variable Y, where Y can only be one of two possible outcomes, in this case Success or Failure. Logistic regression has primarily been used in the fields of epidemiology and biomedical research, but lends itself to many other applications. As indicated the use of logistic regression is not new, however, modeling propulsion ground test facilities using logistic regression is both a new and unique application of the statistical technique. Results from the models provide project managers with insight and confidence into the affectivity of rocket engine component ground test projects. The initial success in modeling rocket propulsion ground test projects clears the way for more complex models to be developed in this area.

Messer, Bradley P.↗

A Human Factors Approach to Bridging Systems and Introducing New Technologies

The application of human factors in aviation has grown to cover a wide range of disciplines and methods capable of assessing human-systems integration at many levels. For example, at the individual level, pilot workload may be studied while at the team level, coordinated workload distribution may be the focal point. At the organizational level, the way in which individuals and teams are supported by training and standards, policies and procedures may introduce additional, relevant topics. A consideration of human factors at each level contributes to our understanding of successes and failures in pilot performance, but this system focused on the flight deck alone -- is only one part of the airspace system. In the FAA's NextGen plan to overhaul the National Airspace System (NAS), new capabilities will enhance flightdeck systems (pilots), flight operations centers (dispatchers) and air traffic control systems (controllers and air traffic managers). At a minimum, the current roles and responsibilities of these three systems are likely to change. Since increased automation will be central to many of the enhancements, the role of automation is also likely to change. Using NextGen examples, a human factors approach for bridging complex airspace systems will be the main focus of this presentation. It is still crucial to consider the human factors within each system, but the successful implementation of new technologies in the NAS requires an understanding of the collaborations that occur when these systems intersect. This human factors approach to studying collaborative systems begins with detailed task descriptions within each system to establish a baseline of the current operations. The collaborative content and context are delineated through the review of regulatory and advisory materials, letters of agreement, policies, procedures and documented practices. Field observations and interviews also help to fill out the picture. Key collaborative functions across systems are identified and placed on a phase-of-flight timeline including information requirements, decision authority and use of automation, as well as level of frequency and criticality.

Kanki, Barbara G.↗

NASA Tech Briefs, September 2006

Topics covered include: Improving Thermomechanical Properties of SiC/SiC Composites; Aerogel/Particle Composites for Thermoelectric Devices; Patches for Repairing Ceramics and Ceramic- Matrix Composites; Lower-Conductivity Ceramic Materials for Thermal-Barrier Coatings; An Alternative for Emergency Preemption of Traffic Lights; Vehicle Transponder for Preemption of Traffic Lights; Automated Announcements of Approaching Emergency Vehicles; Intersection Monitor for Traffic-Light-Preemption System; Full-Duplex Digital Communication on a Single Laser Beam; Stabilizing Microwave Frequency of a Photonic Oscillator; Microwave Oscillators Based on Nonlinear WGM Resonators; Pointing Reference Scheme for Free-Space Optical Communications Systems; High-Level Performance Modeling of SAR Systems; Spectral Analysis Tool 6.2 for Windows; Multi-Platform Avionics Simulator; Silicon-Based Optical Modulator with Ferroelectric Layer; Multiplexing Transducers Based on Tunnel-Diode Oscillators; Scheduling with Automated Resolution of Conflicts; Symbolic Constraint Maintenance Grid; Discerning Trends in Performance Across Multiple Events; Magnetic Field Solver; Computing for Aiming a Spaceborne Bistatic- Radar Transmitter; 4-Vinyl-1,3-Dioxolane-2-One as an Additive for Li-Ion Cells; Probabilistic Prediction of Lifetimes of Ceramic Parts; STRANAL-PMC Version 2.0; Micromechanics and Piezo Enhancements of HyperSizer; Single-Phase Rare-Earth Oxide/Aluminum Oxide Glasses; Tilt/Tip/Piston Manipulator with Base-Mounted Actuators; Measurement of Model Noise in a Hard-Wall Wind Tunnel; Loci-STREAM Version 0.9; The Synergistic Engineering Environment; Reconfigurable Software for Controlling Formation Flying; More About the Tetrahedral Unstructured Software System; Computing Flows Using Chimera and Unstructured Grids; Avoiding Obstructions in Aiming a High-Gain Antenna; Analyzing Aeroelastic Stability of a Tilt-Rotor Aircraft; Tracking Positions and Attitudes of Mars Rovers; Stochastic Evolutionary Algorithms for Planning Robot Paths; Compressible Flow Toolbox; Rapid Aeroelastic Analysis of Blade Flutter in Turbomachines; General Flow-Solver Code for Turbomachinery Applications; Code for Multiblock CFD and Heat-Transfer Computations; Rotating-Pump Design Code; Covering a Crucible with Metal Containing Channels; Repairing Fractured Bones by Use of Bioabsorbable Composites; Kalman Filter for Calibrating a Telescope Focal Plane; Electronic Absolute Cartesian Autocollimator; Fiber-Optic Gratings for Lidar Measurements of Water Vapor; Simulating Responses of Gravitational-Wave Instrumentation; SOFTC: A Software Correlator for VLBI; Progress in Computational Simulation of Earthquakes; Database of Properties of Meteors; Computing Spacecraft Solar-Cell Damage by Charged Particles; Thermal Model of a Current-Carrying Wire in a Vacuum; Program for Analyzing Flows in a Complex Network; Program Predicts Performance of Optical Parametric Oscillators; Processing TES Level-1B Data; Automated Camera Calibration; Tracking the Martian CO2 Polar Ice Caps in Infrared Images; Processing TES Level-2 Data; SmaggIce Version 1.8; Solving the Swath Segment Selection Problem; The Spatial Standard Observer; Less-Complex Method of Classifying MPSK; Improvement in Recursive Hierarchical Segmentation of Data; Using Heaps in Recursive Hierarchical Segmentation of Data; Tool for Statistical Analysis and Display of Landing Sites; Automated Assignment of Proposals to Reviewers; Array-Pattern-Match Compiler for Opportunistic Data Analysis; Pre-Processor for Compression of Multispectral Image Data; Compressing Image Data While Limiting the Effects of Data Losses; Flight Operations Analysis Tool; Improvement in Visual Target Tracking for a Mobile Robot; Software for Simulating Air Traffic; Automated Vectorization of Decision-Based Algorithms; Grayscale Optical Correlator Workbench; "One-Stop Shopping" for Ocean Remote-Sensing and Model Data; State Analysis Database Tool; Generating CAHV and CAHVOmages with Shadows in ROAMS; Improving UDP/IP Transmission Without Increasing Congestion; FORTRAN Versions of Reformulated HFGMC Codes; Program for Editing Spacecraft Command Sequences; Flight-Tested Prototype of BEAM Software; Mission Scenario Development Workbench; Marsviewer; Tool for Analysis and Reduction of Scientific Data; ASPEN Version 3.0; Secure Display of Space-Exploration Images; Digital Front End for Wide-Band VLBI Science Receiver; Multifunctional Tanks for Spacecraft; Lightweight, Segmented, Mostly Silicon Telescope Mirror; Assistant for Analyzing Tropical-Rain-Mapping Radar Data; and Anion-Intercalating Cathodes for High-Energy- Density Cells.

Source record↗