Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fast Region Evaluation”

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 73 records · Page 4

Chance-Constrained Guidance With Non-Convex Constraints

Missions to small bodies, such as comets or asteroids, require autonomous guidance for descent to these small bodies. Such guidance is made challenging by uncertainty in the position and velocity of the spacecraft, as well as the uncertainty in the gravitational field around the small body. In addition, the requirement to avoid collision with the asteroid represents a non-convex constraint that means finding the optimal guidance trajectory, in general, is intractable. In this innovation, a new approach is proposed for chance-constrained optimal guidance with non-convex constraints. Chance-constrained guidance takes into account uncertainty so that the probability of collision is below a specified threshold. In this approach, a new bounding method has been developed to obtain a set of decomposed chance constraints that is a sufficient condition of the original chance constraint. The decomposition of the chance constraint enables its efficient evaluation, as well as the application of the branch and bound method. Branch and bound enables non-convex problems to be solved efficiently to global optimality. Considering the problem of finite-horizon robust optimal control of dynamic systems under Gaussian-distributed stochastic uncertainty, with state and control constraints, a discrete-time, continuous-state linear dynamics model is assumed. Gaussian-distributed stochastic uncertainty is a more natural model for exogenous disturbances such as wind gusts and turbulence than the previously studied set-bounded models. However, with stochastic uncertainty, it is often impossible to guarantee that state constraints are satisfied, because there is typically a non-zero probability of having a disturbance that is large enough to push the state out of the feasible region. An effective framework to address robustness with stochastic uncertainty is optimization with chance constraints. These require that the probability of violating the state constraints (i.e., the probability of failure) is below a user-specified bound known as the risk bound. An example problem is to drive a car to a destination as fast as possible while limiting the probability of an accident to 10(exp -7). This framework allows users to trade conservatism against performance by choosing the risk bound. The more risk the user accepts, the better performance they can expect.

FROM↗

Implementation of a Discrete Dipole Approximation Scattering Database into Community Radiative Transfer Model

The Community Radiative Transfer Model (CRTM) is a fast model that requires bulk optical properties of hydrometeors in the form of lookup tables to simulate all-sky satellite radiances. Current cloud scattering lookup tables of CRTM were generated using the Mie-Lorenz theory thus assuming spherical shapes for all frozen habits, while actual clouds contain frozen hydrometeors with different shapes. The Discrete Dipole Approximation (DDA) technique is an effective technique for simulating the optical properties of non-spherical hydrometeors in the microwave region. This paper discusses the implementation and validation of a comprehensive DDA cloud scattering database into CRTM for the microwave frequencies. The original DDA database assumes total random orientation in the calculation of single scattering properties. The mass scattering parameters required by CRTM were then computed from single scattering properties and water content dependent particle size distributions. The new lookup tables eliminate the requirement for providing the effective radius as input to CRTM by using the cloud water content for the mass dimension. A collocated dataset of short-term forecasts from Integrated Forecast System of the European Centre for Medium-Range Weather Forecasts and satellite microwave data was used for the evaluation of results. The results overall showed that the DDA lookup tables, in comparison with the Mie tables, greatly reduce the differences among simulated and observed values. The Mie lookup tables especially introduce excessive scattering for the channels operating below 90 GHz and low scattering for the channels above 90 GHz.

Discrete dipole↗

Solutions Network Formulation Report. Integration of OMI and TES Aerosol Products into the EPA Regional Planning Organizations' FASTNET Aerosol Tracking and Analysis Tool

Every year, more than 280 million visitors tour our Nation s most treasured parks and wilderness areas. Unfortunately, many visitors are unable to see the spectacular vistas they expect because of white or brown haze in the air. Most of this haze is not natural; it is air pollution, carried by the wind often hundreds of miles from its origin. Some of the pollutants have been linked to serious health problems, such as asthma and other lung disorders, and even premature death. In addition, nitrates and sulfates contribute to acid rain formation, which contaminates rivers and lakes and erodes buildings and historical monuments. The U.S. Environmental Protection Agency RPOs (Regional Planning Organizations) have been tasked with monitoring and determining the nature and origin of haze in Class I scenic areas, and finding ways to reduce haze in order to improve visibility in these areas. The RPOs have developed an Internet-based air quality DST (Decision Support Tool) called FASTNET (Fast Aerosol Sensing Tools for Natural Event Tracking). While FASTNET incorporates a few satellite datasets, most of the data utilized by this DST comes from ground-based instrument networks. The problem is that in many areas the sensors are sparsely located, with long distances between them, causing difficulties in tracking haze over the United States, determining its source, and analyzing its content. Satellite data could help to fill in the data gaps and to supplement and verify ground-recorded air quality data. Although satellite data are now being used for air quality research applications, such data are not routinely used for environmental decision support, in part because of limited resources, difficulties with interdisciplinary data interpretation, and the need for advanced inter-agency partnerships. As a result, the validation and verification of satellite data for air quality operational system applications has been limited This candidate solution evaluates the usefulness of OMI (Ozone Monitoring Instrument) and TES (Tropospheric Emission Spectrometer) air quality data for the RPOs by comparing OMI and TES data with ground-based data that are acquired during identified episodes of air pollution. The air quality data from OMI and TES are of different spectral ranges than data from satellites currently included in FASTNET, giving them potential advantages over the existing satellites. If the OMI and TES data are shown to be useful to the RPOs, they would then be integrated into the FASTNET DST for use on an operational basis.

Knowlton, Kelly↗

A description and evaluation of TIMER: A time-based terminal flow-control concept

A description of a time-based ATC concept called TIMER (Traffic Intelligence for the Management of Efficient Runway-scheduling) and the results of a fast time and real time computer evaluation are presented. The concept was designed to improve the efficiency of extended terminal area operations (en route approach, transition, and terminal flight to the runway). TIMER integrates en route metering, fuel efficient cruise and profile descents, terminal sequencing and spacing together with computer-generated controller aids, in order to fully use runway capacity and improve efficiency of delay absorption. The concept, by using simplified aircraft models, accommodates both 4-D and non 4-D equipped aircraft and is designed for integration into the manual, voice linked ATC system in an evolutionary manner and still be able to accommodate proposed system upgrade features such as data link and further ground automation. Fast time and real time computer simulation results identify and show the effects and interactions of such key variables as horizon of control, metering fix and final approach delivery time errors, aircraft separation requirements, delay discounting, wind, flight technical error, and knowledge of aircraft final approach speed. The current ATC system has a runway interarrival-error standard deviation of approx. 26 seconds. Simulation results indicate that, with computer aiding, the runway interarrival-error standard deviation for non 4-D equipped traffic can be reduced to the region of 8 to 12 seconds if expected final approach speed is known; however, the reduction is only in the region of 16 to 20 seconds if expected final approach speed is unknown. Another major finding is that en route metering fix delivery-error standard deviation should be kept to less than a number somewhere between 35 to 45 seconds to achieve full runway capacity. This requirement implies the need for either airborne automation or assistance to the controller since the current manual performance in today's en route metering environment is in the order of 1.5 minutes.

Credeur, Leonard↗

Planetary Nebulae with Supporting IR Data

We present new HST/WFPC2 imagery for the planetary nebula (PN) NGC 6818. Observations were made in line filters F437N, F487N, F502N, and F656N plus continuum filter F547M. The primary goal was to develop a high spatial resolution (approx. 0.1 in.) map of the intrinsic line ratio [O III] 4363/5007 and thereby evaluate the electron temperature (T(sub e)) and the mean-square T(sub e) variation (t(sup 2) across the PN. In this process we developed an extinction map from the F487N (H(beta)) and F656N (H(alpha)) images by comparing the observed line ratios in each pixel to the theoretical ratio and computing a c(H(beta)) map which was used to correct the observed 4363/5007 ratios for reddening. We also adjust for the continuum contribution to the line filter data. We present color-coded pictures of the reddening (c(H(beta))) map, the [O III] T(sub e) map, as well as our determinations of t(sup 2). The T(sub e) map shows a decline from approx. 14000 K in the inner regions to approx. 11000 K at the outer edge. Such a radial T(sub e) gradient is expected for a high-excitation nebula with a prominent He(++) zone such as NGC 6818. A composite of images taken in 3 filters (F656N, red; F487N, blue; and F502N, [O III] 5007, green) shows a roughly spherical outer envelope as well as a brighter vase-shaped interior "bubble". There is a prominent orifice to the North and a smaller one to the South, along the major axis, possibly caused by a blow-out from a fast wind. This nebula has an appearance remarkably similar to that of the PN NGC 3918 previously imaged with HST by H. Bond. We note from the continuum images (F547M) two stars in the nebular field that are fainter than the prominent central star; these are roughly 2-4 sec. N and NE of the central star. Further study is needed to establish whether or not there may be a physical association of either star with the central star.

Harrington, J. Patrick↗

The Thermal Infrared Optical Depth of Mineral Dust Retrieved from Integrated CALIOP and IIR Observations

Recent studies reveal that a higher fraction of coarse mineral dust particles than that estimated by climate model simulations has been observed in the atmosphere, leading to a more significant positive (i.e., warming) longwave (LW) thermal infrared (TIR) direct radiative effect (DRE). However, the magnitude of this DRE remains highly uncertain because our understanding of the radiative properties, quantitatively represented by the optical depth of dust, especially information on the TIR, remains limited. This study presents a simple approach to retrieve the thermal infrared dust aerosol optical depth (DAODTIR) over oceans during nighttime using the observations from the Infrared Imaging Radiometer (IIR) and the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard CALIPSO. For each cloud-free dust-laden profile identified by the IIR-CALIOP observation, a Lookup-Table (LUT) of the 10.6 μm IIR band brightness temperatures difference (dBT) under different DAODTIR with respect to their dust-free BTs is constructed based on the CALIOP retrieved dust vertical profile and pre-assumed dust scattering properties using a fast radiative transfer model. Then the DAODTIR is retrieved by projecting the IIR-observed dBT on the LUT. Sensitivity studies show that the DAODTIR retrieval at 10.6 μm is more susceptible to the dust particle size distribution (PSD) assumption than dust refractive indices. To estimate the uncertainty caused by PSD assumption, two DAODTIR retrieval products, one based on the dust PSD from the AERONET at Cape Verde and the other on an in situ measured PSD from the recent Fennec campaign, are provided. The retrieval uncertainty is mainly contributed by the BT difference between the observation and simulation using auxiliary atmospheric data. The climatology of the retrieval from 2013 to 2019 shows confident spatiotemporal variations of DAODTIR with the global-averaged value of 0.006 and 0.008 based on different pre-assumed dust PSDs. Climatological results agree reasonably well with two independent DAODTIR retrieval products based on the Infrared Atmospheric Sounding Interferometer (IASI) over the active dust transport regions, such as North and Tropical Atlantic (r = 0.904 and 0.819) and Indian Ocean (r = 0.832). The seasonal and interannual variation is also well-compared (r = 0.758) with AERONET coarse-mode AOD at 97 selected sites. The synergic CALIOP observation allows the retrieved DAODTIR to directly compare with the extrapolated DAODTIR from DAOD in visible (i.e., 532 nm), which helps evaluate the observational constraints on DAODTIR. This study offers a unique prospect of collocating active lidar and passive IR observations for retrieving dust DAODTIR.

aerosol↗

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗

Lunar and Planetary Science XXXVI, Part 13

Contents include the following: A Fast, Non-Destructive Method for Classifying Ordinary Chondrite Falls Using Density and Magnetic Susceptibility. An Update on Results from the Magnetic Properties Experiments on the Mars Exploration Rovers, Spirit and Opportunity. Measurement Protocols for In Situ Analysis of Organic Compounds at Mars and Comets. Piping Structures on Earth and Possibly Mars: Astrobiological Implications. Uranium and Lead in the Early Planetary Core Formation: New Insights Given by High Pressure and Temperature Experiments. The Mast Cameras and Mars Descent Imager (MARDI) for the 2009 Mars Science Laboratory. MGS MOC: First Views of Mars at Sub-Meter Resolution from Orbit. Analysis of Candor Chasma Interior Layered Deposits from OMEGA/MEX Spectra. Analysis of Valley Networks on Valles Marineris Plateau Using HRSC/MEX Data. Solar Abundance of Elements from Neutron-Capture Cross Sections. Preliminary Evaluation of the Secondary Ion/Accelerator Mass Spectrometer, MegaSIMS. Equilibrium Landforms in the Dry Valleys of Antarctica: Implications for Landscape Evolution and Climate Change on Mars. Continued Study of Ba Isotopic Compositions of Presolar Silicon Carbide Grains from Supernovae. Paleoenviromental Evolution of the Holden-Uzboi Area. Stability of Magnesium Sulfate Minerals in Martian Environments. Tungsten Isotopic Constraints on the Formation and Evolution of Iron Meteorite Parent Bodies. Migration of Dust Particles and Volatiles Delivery to the Inner Planets. On the Sitting of Trapped Noble Gases in Insoluble Organic Matter of Primitive Meteorites. Trapping of Xenon Upon Evaporation-Condensation of Organic Matter Under UV Irradiation: Isotopic Fractionation and Electron Paramagnetic Resonance Analysis. Stability of Water on Mars. A Didactic Activity. Analysis of Coronae in the Parga Chasma Region, Venus. Photometric and Compositional Surface Properties of the Gusev Crater Region, Mars, as Derived from Multi-Angle, Multi-Spectral Investigation of Mars Express HRSC Data. Mapping Compositional Diversity on Mars: Spatial Distribution and Geological Implications. A New Simulation Chamber for Studying Planetary Environments. Folded Structure in Terra Sirenum. Mars. Nitrogen-Noble Gas Static Mass Spectrometry of Genesis Collector Materials. Neon Isotope Heterogeneity in the Terrestrial Mantle: Implication for the Acquisition of Volatile Elements in Terrestrial Planets. The Cosmic Clock, the Cycle of Terrestrial Mass Extinctions.

Source record↗

Evaluation of the Interference of Tenax®TA Adsorbent With Dimethylformamide Dimethyl Acetal Reagent for Gas Chromatography-Dragonfly Mass Spectrometry and Future Gas Chromatography-Mass Spectrometry in Situ Analysis

Among future space missions, national aeronautics and space administration (NASA) selected two of them to analyze the diversity in organic content within Martian and Titan soil samples using a gas chromatograph – mass spectrometer (GC–MS) instrument. The Dragonfly space mission is planned to be launched in 2027 to Titan's surface and explore the Shangri-La surface region for years. One of the main goals of this mission is to understand the past and actual abundant prebiotic chemistry on Titan, which is not well characterized yet. The ExoMars space mission is planned to be launched in 2028 to Mars’ surface and explore the Oxia Planum and Mawrth Vallis region for years. The main objectives focus on the exploration of the subsurface soil samples, potentially richer in organics, that might be relevant for the search of past life traces on Mars where irradiation does not impact the matrices and organics. One recently used sample pre-treatment for gas chromatography – mass spectrometry analysis is planned on both space missions to detect refractory organic molecules of interest for astrobiology. This pre-treatment is called derivatization and uses a chemical reagent – called dimethylformamide dimethyl acetal (DMF-DMA) – to sublimate organic compounds keeping them safe from thermal degradation and conserving the chirality of the molecules extracted from Titan or Mars’ matrices. Indeed, the detection of building blocks of life or enantiomeric excess of some organics (e.g. amino acids) after DMF-DMA pre-treatment and GC–MS analyses would be both bioindicators. The main results highlighted by our work on DMF-DMA and Tenax®TA interaction and efficiency to detect organic compounds at ppb levels in a fast and single preparation are first that Tenax®TA did not show the onset of degradation until after 150 experiments – a 120 h at 300 °C experiment – which greatly exceeds the experimental lifetimes for the DraMS and GC-space in situ investigations. Tenax®TA polymer and DMF-DMA produce many by-products (about 70 and 46, respectively, depending on the activation temperature). Further, the interaction between the two leads to the production of 22 additional by-products from DMF-DMA degradation, but these listed by-products do not prevent the detection of trace-level organic molecules after their efficient derivatization and volatilization by DMF-DMA in the oven ahead the GC–MS trap and column.

DraMS-Dragonfly mission↗

Formation of a Tropopause Cirrus Layer Observed over Florida during CRYSTAL-FACE

On July 13, 2002 a widespread, subvisible tropopause cirrus layer occurred over the Florida region. This cloud was observed in great detail with the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) instrumentation, including in situ measurements with the WB-57 aircraft. In this paper, we use the 13 July cloud as a case study to evaluate the physical processes controlling the formation and evolution of tropopause cirrus layers. Microphysics measurements indicate that ice crystal diameters in the cloud layer ranged from about 7 to 50 microns, and the peak number mode was about 10-25 microns. In situ water vapor and temperature measurements in the cloud indicated supersaturation with respect to ice throughout, with ice saturation ratios as large as 1.8. Even when the ice surface area density was as high as about 500 sq microns/cu cm, ice supersaturations of 20-30% were observed. Trajectory analysis shows that the air sampled near the tropopause on this day generally came from the north and cooled considerably during the previous few days. Examination of infrared satellite imagery along air parcel back trajectories from the WB-57 flight track indicates that the tropopause cloud layer formation was, in general, not simply left over ice from recently generated anvil cirrus. Simulations of cloud formation using time-height curtains of temperature along the trajectory paths show that the cloud could have formed in situ near the tropopause as the air was advected into the south Florida region and cooled to unusually low temperatures. If we assume a high threshold for ice nucleation via homogeneous freezing of aqueous sulfate aerosols, the model reproduces the observed cloud structure, ice crystal size distributions, and ice supersaturation statistics. Inclusion of observed gravity wave temperature perturbations in the simulations is essential to reproduce the observed cloud properties. Without waves, crystal number densities are too low, crystal sizes are too large, and the crystals fall out too fast, leaving very little cloud persisting at the end of the simulations. In the cloud simulations, coincidence of high supersaturations and high surface areas can be produced by either recent nucleation or sedimentation of crystals into supersaturated layers. The agreement between model results and observed supersaturations is improved somewhat if we assume that the steady state relative humidity within cirrus at T<200 K is enhanced by about 30%. The WB-57 measurements and the model results suggest that the cloud layer irreversibly dehydrated air near the tropopause.

Jensen, Eric↗

Constraints on Nonlinear and Stochastic Growth Theories for Type 3 Solar Radio Bursts from the Corona to 1 AU

Existing, competing theories for coronal and interplanetary type III solar radio bursts appeal to one or more of modulational instability, electrostatic (ES) decay processes, or stochastic growth physics to preserve the electron beam, limit the levels of Langmuir-like waves driven by the beam, and produce wave spectra capable of coupling nonlinearly to generate the observed radio emission. Theoretical constraints exist on the wavenumbers and relative sizes of the wave bandwidth and nonlinear growth rate for which Langmuir waves are subject to modulational instability and the parametric and random phase versions of ES decay. A constraint also exists on whether stochastic growth theory (SGT) is appropriate. These constraints are evaluated here using the beam, plasma, and wave properties (1) observed in specific interplanetary type III sources, (2) predicted nominally for the corona, and (3) predicted at heliocentric distances greater than a few solar radii by power-law models based on interplanetary observations. It is found that the Langmuir waves driven directly by the beam have wavenumbers that are almost always too large for modulational instability but are appropriate to ES decay. Even for waves scattered to lower wavenumbers (by ES decay, for instance), the wave bandwidths are predicted to be too large and the nonlinear growth rates too small for modulational instability to occur for the specific interplanetary events studied or the great majority of Langmuir wave packets in type III sources at arbitrary heliocentric distances. Possible exceptions are for very rare, unusually intense, narrowband wave packets, predominantly close to the Sun, and for the front portion of very fast beams traveling through unusually dilute, cold solar wind plasmas. Similar arguments demonstrate that the ES decay should proceed almost always as a random phase process rather than a parametric process, with similar exceptions. These results imply that it is extremely rare for modulational instability or parametric decay to proceed in type III sources at any heliocentric distance: theories for type III bursts based on modulational instability or parametric decay are therefore not viable in general. In contrast, the constraint on SGT can be satisfied and random phase ES decay can proceed at all heliocentric distances under almost all circumstances. (The contrary circumstances involve unusually slow, broad beams moving through unusually hot regions of the Corona.) The analyses presented here strongly justify extending the existing SGT-based model for interplanetary type III bursts (which includes SGT physics, random phase ES decay, and specific electromagnetic emission mechanisms) into a general theory for type III bursts from the corona to beyond 1 AU. This extended theory enjoys strong theoretical support, explains the characteristics of specific interplanetary type III bursts very well, and can account for the detailed dynamic spectra of type III bursts from the lower corona and solar wind.

Cairns, Iver H.↗

Engineering Evaluation of International Low Impact Docking System Latch Hooks

The international Low Impact Docking System (iLIDS) provides a structural arrangement that allows for visiting vehicles to dock with the International Space Station (ISS) (Fig 1). The iLIDS docking units are mechanically joined together by a series of active and passive latch hooks. In order to preserve docking capability at the existing Russian docking interfaces, the iLIDS latch hooks are required to conform to the existing Russian design. The latch hooks are classified as being fail-safe. Since the latch hooks are fail-safe, the hooks are not fracture critical and a fatigue based service life assessment will satisfy the structural integrity requirements. Constant amplitude fatigue testing to failure on four sets of active/passive iLIDS latch hooks was performed at load magnitudes of 10, 11, and 12 kips. Failure analysis of the hook fatigue failures identified multi-site fatigue initiation that was effectively centered about the hook mid-plane (consistent with the 3D model results). The fatigue crack initiation distribution implies that the fatigue damage accumulation effectively results in a very low aspect ratio surface crack (which can be simulated as thru-thickness crack). Fatigue damage progression resulted in numerous close proximity fatigue crack initiation sites. It was not possible to determine if fatigue crack coalescence occurs during cyclic loading or as result of the fast fracture response. The presence of multiple fatigue crack initiation sites on different planes will result in the formation of ratchet marks as the cracks coalesce. Once the stable fatigue crack becomes unstable and the fast fracture advances across the remaining ligament and the plane stress condition at a free-surface will result in failure along a 45 deg. shear plane (slant fracture) and the resulting inclined edge is called a shear lip. The hook thickness on the plane of fatigue crack initiation is 0.787". The distance between the shear lips on this plane was on the order of 0.48" and it was effectively centered about the mid-plane of the section. The numerous ratchet marks between the shear lips on the fracture initiation plane are indicative of multiple fatigue initiation sites within this region. The distribution of the fatigue damage about the centerline of the hook is consistent with the analytical results that demonstrate peak stress/strain response at the mid-plane that decreases in the direction of the hook outer surfaces. Scanning electron microscope images of the failed sections detected fatigue crack striations in close proximity to the free surface of the hook radius. These findings were documented at three locations on the fracture surface : 1) adjacent to the left shear lip, 2) adjacent to the right shear lip, and 3) near the centerline of the section. The features of the titanium fracture surface did not allow for a determination of a critical crack size via identification of the region where the fatigue crack propagation became unstable. The fracture based service life projections where benchmarked with strain-life analyses. The strainrange response in the hook radius was defined via the correlated finite element models and the modified method of universal slopes was incorporated to define the strain-life equation for the titanium alloy. The strain-life assessment confirmed that the fracture based projections were reasonable for the loading range of interest. Based upon the analysis and component level fatigue test data a preliminary service life capability for the iLIDS active and passive hooks of 2 lifetimes is projected (includes a scatter factor of 4).

Martinez, J.↗

Fission Surface Power Technology Development Testing at NASA's Early Flight Fission Test Facility

Fission surface power (FSP) systems could be used to provide power anytime, anywhere on the surface of the Moon or Mars. FSP systems could be used at polar locations, at locations away from the poles, or in permanently shaded regions, with excellent performance at all sites. A potential reference 40 kWe option has been devised that is cost-competitive with alternatives while providing more power for less mass anywhere on the lunar surface. The reference FSP system (FSPS) is also readily extensible for use on Mars. At Mars the system would be capable of operating through global dust storms and providing year-round power at any Martian latitude. Under the NASA Exploration Technology Development Program (ETDP), NASA and the Department of Energy (DOE) have begun technology development on Fission Surface Power (FSP). The primary customer for this technology is the NASA Constellation Program which is responsible for the development of surface systems to support human exploration on the moon and Mars. The objectives of the FSP technology project are: 1) Develop FSP concepts that meet expected surface power requirements at reasonable cost with added benefits over other options. 2) Establish a hardware-based technical foundation for FSP design concepts and reduce overall development risk. 3) Reduce the cost uncertainties for FSP and establish greater credibility for flight system cost estimates. 4) Generate the key products to allow Agency decision-makers to consider FSP as a viable option for flight development. To be mass efficient, FSP systems must operate at higher coolant temperatures and use different types of power conversion than typical terrestrial systems. The primary reason is the difficulty in rejecting excess heat to space. Although many options exist, NASA s current reference FSP system uses a fast spectrum, pumped-NaK cooled reactor coupled to a Stirling power conversion subsystem. The reference system uses technology with significant terrestrial heritage while still providing excellent performance on the surface of the moon or Mars. Recent testing at NASA s Early Flight Fission Test Facility (EFF-TF) has helped assess the viability of the reference FSP system, and has helped evaluate methods for system integration. In June, 2009, a representative pumped NaK loop (provided by Marshall Space Flight Center) was coupled to a Stirling power converter (provided by Glenn Research Center) and tested at various conditions representative of those that would be seen during actual FSP system operation. In all areas, performance of the integrated system exceeded project goals. High-temperature NaK pump testing has also been performed at the EFF-TF, as has testing of methods for providing long-duration NaK purity.

Houts. Michael G.↗

Simulation of a Periodic Jet in a Crossflow with a RANS Solver Using an Unstructured Grid

A second-order unstructured-grid code, developed and used primarily for steady aerodynamic simulations, is applied to the synthetic jet in a cross flow. The code, FUN3D, is a vertex-centered finite-volume method originally developed by Anderson[1, 2], and is currently supported by members of the Fast Adaptive Aerospace Tools team at NASA Langley. Used primarily for design[3] and analysis[4] of steady aerodynamic configurations, FUN3D incorporates a discrete adjoint capability, and supports parallel computations using MPI. A detailed description of the FUN3D code can be found in the references given above. The code is under continuous development and contains a variety of flux splitting algorithms for the inviscid terms, two methods for computing gradients, several turbulence models, and several solution methodologies; all in varying states of development. Only the most robust and reliable components, based on experiences with steady aerodynamic simulations, were employed in this work. As applied in this work, FUN3D solves the Reynolds averaged Navier-Stokes equations using the one equation turbulence model of Spalart and Allmaras[5]. The spatial discretization is formed on unstructured meshes using a vertex-centered approach. The inviscid terms are evaluated by a flux-difference splitting formulation using least-squares reconstruction and Roe-type approximate Riemann fluxes. Green-Gauss gradient evaluations are used for viscous and turbulence modeling terms. The discrete spatial operator is combined with a backward time operator which is then solved iteratively using point or line Gauss-Seidel and local time stepping in a pseudo time. For steady flows, the physical time step is set to infinity and the pseudo time step is ramped up with the iteration count. A second-order backward in time operator is used for time accurate flows with 20 to 50 steps in the pseudo time applied at each physical time step. For this effort, FUN3D was modified to support spatially varying boundary and initial conditions, and unsteady boundary conditions. Also, a specialized in/out flow boundary condition was implemented to model the action of the diaphragm. This boundary condition is described below in more detail. The grids were generated using the internally developed codes GridEX[6] for meshing the surfaces and inviscid regions of the domain, and for CAD access; and MesherX[7] for meshing the viscous regions. Grid spacing in on the surfaces and in the inviscid regions are indirectly controlled by specifying sources. The viscous layers are generated using an advancing layer technique. MeshersX allows the user to control the spatial variation of the first step off the surface, growth rates, and the termination criterion by providing small problem dependent subroutines.

Atkins, H. L.↗

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew↗

Platform for Postprocessing Waveform-Based NDE

Taking advantage of the similarities that exist among all waveform-based non-destructive evaluation (NDE) methods, a common software platform has been developed containing multiple- signal and image-processing techniques for waveforms and images. The NASA NDE Signal and Image Processing software has been developed using the latest versions of LabVIEW, and its associated Advanced Signal Processing and Vision Toolkits. The software is useable on a PC with Windows XP and Windows Vista. The software has been designed with a commercial grade interface in which two main windows, Waveform Window and Image Window, are displayed if the user chooses a waveform file to display. Within these two main windows, most actions are chosen through logically conceived run-time menus. The Waveform Window has plots for both the raw time-domain waves and their frequency- domain transformations (fast Fourier transform and power spectral density). The Image Window shows the C-scan image formed from information of the time-domain waveform (such as peak amplitude) or its frequency-domain transformation at each scan location. The user also has the ability to open an image, or series of images, or a simple set of X-Y paired data set in text format. Each of the Waveform and Image Windows contains menus from which to perform many user actions. An option exists to use raw waves obtained directly from scan, or waves after deconvolution if system wave response is provided. Two types of deconvolution, time-based subtraction or inverse-filter, can be performed to arrive at a deconvolved wave set. Additionally, the menu on the Waveform Window allows preprocessing of waveforms prior to image formation, scaling and display of waveforms, formation of different types of images (including non-standard types such as velocity), gating of portions of waves prior to image formation, and several other miscellaneous and specialized operations. The menu available on the Image Window allows many further image processing and analysis operations, some of which are found in commercially-available image-processing software programs (such as Adobe Photoshop), and some that are not (removing outliers, Bscan information, region-of-interest analysis, line profiles, and precision feature measurements).

Roth, Don↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗