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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1999

SeaWiFS Postlaunch Technical Report Series: The 1997 Prelaunch Radiometric Calibration of SeaWiFS - Volume 4

The Sea-viewing Wide Field-of-view Sensor (SeaWiFS) was originally calibrated by the instrument's manufacturer, Santa Barbara Research Center (SBRC), in November 1993. In preparation for an August 1997 launch, the SeaWiFS Project and the National Institute of Standards and Technology (NIST) undertook a second calibration of SeaWiFS in January and April 1997 at the facility of the spacecraft integrator, Orbital Sciences Corporation (OSC). This calibration occurred in two phases, the first after the final thermal vacuum test, and the second after the final vibration test of the spacecraft. For the calibration, SeaWiFS observed an integrating sphere from the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) at four radiance levels. The spectral radiance of the sphere at these radiance levels was also measured by the SeaWiFS Transfer Radiometer (SXR). In addition, during the calibration, SeaWiFS and the SXR observed the sphere at 16 radiance levels to determine the linearity of the SeaWiFS response. As part of the calibration analysis, the GSFC sphere was also characterized using a GSFC spectroradiometer. The 1997 calibration agrees with the initial 1993 calibration to within +/- 4%. The new calibration coefficients, computed before and after the vibration test, agree to within 0.5%. The response of the SeaWiFS channels in each band is linear to better than 1%. In order to compare to previous and current methods, the SeaWiFS radiometric responses are presented in two ways: using the nominal center wave-lengths for the eight bands; and using band-averaged spectral radiances. The band-averaged values are used in the flight calibration table. An uncertainty analysis for the calibration coefficients is also presented.

Hooker, Stanford B.↗

NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗

NASA NTRSDate not supplied

Spaceflight Medical Evacuation Risk Assessment Principles - A Qualitative Investigation

BACKGROUND Future human space exploration beyond Low Earth Orbit (LEO) will require innovative solutions in many areas, primarily those that provide direct medical support to crew on long-duration missions and optimize their health and performance. The associated challenges therein will be numerous, including but not necessarily limited to extended one-way or "asynchronous" communication delays, minimal to non-existent resupply, and a prolonged transit time to "definitive care" ranging from 3-days to 9-months. At the same time, long-duration exploration spacecraft and crew will face restrictions on mass, power, volume, and data far more significant than that seen in current LEO settings. Given these limitations, medical risk assessment is of primary importance, especially evaluating the implications of a medical evacuation of an ill or injured crewmember. Such evacuations are complicated, potentially dangerous, and well may be impossible in certain phases of the mission. Regardless, such issues must be weighed against the risks of the injured crew remaining aboard a spacecraft with limited medical resources. OBJECTIVE This qualitative study drew from the experiences of subject matter experts (SMEs) in spaceflight and appropriate analog environments (i.e., military, disaster, and extreme environment fields) to identify unique principles common amongst medical evacuation considerations helpful in informing future risk assessment tools. Appropriate analog environments included austere operational settings where multiple factors (weather, logistics/limited resupply, denied/extreme environments, and patient condition) resulted in a limited ability to provide definitive local medical care. The fundamental principles in question revolved around scenarios where evacuation became a complicating yet necessary consideration and where life-threatening medical concerns had to be weighed against critical mission objective(s). The primary authors collected semi-structured data gathered through in-depth interviews with 16 subject matter experts (SMEs). Interview questions investigated how these SMEs consider and weigh the attendant risks present in medical evacuation scenarios. Among the critical questions posed were those that sought to understand how the SME balanced the challenges and requirements of medical evacuation (or keeping an injured patient "on-site" aka: "prolonged field care") against the evacuation operation's risks on impacting overarching mission success. The team analyzed interview transcripts for common themes and principles using the qualitative methods of thematic analysis based on consensus, co-occurrence, and comparison. As a result, nine primary risk consideration themes and nine contributing factor themes emerged, all of which will ideally inform future medical evacuation decision-making tools and operational decision-making for exploration class missions. Specific aims for this study included: 1. Identification of common principles used to assess risks and benefits of medical evacuations in extreme environments 2. Identification of common points of friction or complication and challenges in extreme environment evacuations

A T Almand↗

NASA NTRSDate not supplied

A New Era of H-O-C-S Magma Solubility Modeling: Better, Faster, Stronger

H 2 O, CO 2 , and S are the most abundant volatiles in magmatic systems and are critical to understanding magma storage, phase equilibria, and volcanic eruptions. Models that consider all three of these components, however, may not allow for critical examination and adjustment of assumptions underlying the model, or provide benchmark testing or extensible interfaces. Thus, understanding why models produce different results can be challenging. We have gathered authors of established (D-Compress) and recent (VolFe, EVo, Sulfur_X, MAGEC) H-O-C-S volatile solubility models to work together to understand how and why our models diverge. We present a series of benchmark basalt degassing scenarios revealing that often understated model assumptions such as fO 2 buffer equations, fO 2 -Fe 3+ /ΣFe relationships, and even major element normalization routines have outsized effects on model results. All models consider S 2- and S 6+ melt species but with different approaches to sulfate/sulfide capacities, partition coefficients, and species fugacities, leading to divergence in the evolution of modeled gas compositions, melt S and Fe speciation, and fO 2 , with the extent of divergence depending on melt composition. Such scenarios enable meaningful intercomparison of existing models and lay the groundwork for a user-friendly yet powerful solubility modeling framework. Given our wealth of existing solubility literature, we suggest that the field of magmatic volatiles should focus now on the creation of modern tools and the modular implementation of existing model equations or methods, and that the evaluation of code usability, transparency, and benchmarking should be codified pillars of the peer-review process. As an example of such an endeavor, we present early work coupling these sulfur solubility models with VESIcal, an extensible and rigorously tested python library containing seven existing H 2 O-CO 2 solubility models. VESIcal includes the ability to extract, edit, and even interchange assumptions underlying any model. For example, users may combine or swap separately published H 2 O, CO 2 , and S models, as well as underlying model choices, such as Equations of State and redox models.

volatiles in magmas↗

NASA NTRS2002

Microphysics, Radiation and Surface Processes in the Goddard Cumulus Ensemble (GCE) Model

One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and in-coming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a CRM, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. The GCE model has been used to understand the following: 1) water and energy cycles and their roles in the tropical climate system; 2) the vertical redistribution of ozone and trace constituents by individual clouds and well organized convective systems over various spatial scales; 3) the relationship between the vertical distribution of latent heating (phase change of water) and the large-scale (pre-storm) environment; 4) the validity of assumptions used in the representation of cloud processes in climate and global circulation models; and 5) the representation of cloud microphysical processes and their interaction with radiative forcing over tropical and midlatitude regions. Four-dimensional cloud and latent heating fields simulated from the GCE model have been provided to the TRMM Science Data and Information System (TSDIS) to develop and improve algorithms for retrieving rainfall and latent heating rates for TRMM and the NASA Earth Observing System (EOS). More than 90 referred papers using the GCE model have been published in the last two decades. Also, more than 10 national and international universities are currently using the GCE model for research and teaching. In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.

Tao, Wei-Kuo↗

NASA NTRS2008

Micro-Vibration Measurements on Thermally Loaded Multi-Layer Insulation Samples in Vacuum

Some scientific missions require to an extreme extent the absence of any on-board microvibration. Recent projects dedicated to measuring the Earth's gravity field and modeling the geoid with extremely high accuracy are examples. Their missions demand for extremely low micro-vibration environment on orbit for: (1) Not disturbing the measurement of earth gravity effects with the installed gradiometer or (2) Even not damaging the very high sensitive instruments. Based on evidence from ongoing missions multi-layer insulation (MLI) type thermal control blankets have been identified as a structural element of spacecrafts which might deform under temperature variations being caused by varying solar irradiation in orbit. Any such deformation exerts tiny forces which may cause small reactions resulting in micro-vibrations, in particular by exciting the spacecraft eigenmodes. The principle of the test set-up for the micro-vibration test was as follows. A real side wall panel of the spacecraft (size about 0.25 m2) was low-frequency suspended in a thermal vacuum chamber. On the one side of this panel, the MLI samples were fixed by using the standard methods. In front of the MLI, an IR-rig was installed which provided actively controlled IR-radiation power of about 6 kW/m2 in order to heat the MLI surface. The cooling was passive using the shroud temperature at a chamber pressure <1E-5mbar. The resulting micro-vibrations due to MLI motion in the heating and the cooling phase were measured via seismic accelerometers which were rigidly mounted to the panel. Video recording was used to correlate micro-vibration events to any visual MLI motion. Different MLI sample types were subjected to various thermal cycles in a temperature range between -60 C to +80 C. In this paper, the experience on these micro-vibration measurements will be presented and the conclusions for future applications will be discussed

Deutsch, Georg↗

NASA NTRS2005

NASA Tech Briefs, September 2005

Topivs include: Diamond-Coated Carbon Nanotubes for Efficient Field Emission; Improved Anode Coatings for Direct Methanol Fuel Cells; Advanced Ablative Insulators and Methods of Making Them; PETIs as High-Temperature Resin-Transfer-Molding Materials; Stable Polyimides for Terrestrial and Space Uses; Low-Density, Aerogel-Filled Thermal-Insulation Tiles; High-Performance Polymers Having Low Melt Viscosities; Nonflammable, Hydrophobic Aerogel Composites for Insulation; Front-Side Microstrip Line Feeding a Raised Antenna Patch; Medium-Frequency Pseudonoise Georadar; Facilitating Navigation Through Large Archives; Program for Weibull Analysis of Fatigue Data; Comprehensive Micromechanics-Analysis Code - Version 4.0; Component-Based Visualization System; Software for Engineering Simulations of a Spacecraft; LabVIEW Interface for PCI-SpaceWire Interface Card; Path Following with Slip Compensation for a Mars Rover; International Space Station Electric Power System Performance Code-SPACE; Software for Automation of Real-Time Agents, Version 2; Software for Optimizing Plans Involving Interdependent Goals; Computing Gravitational Fields of Finite-Sized Bodies; Custom Sky-Image Mosaics from NASA's Information Power Grid; ANTLR Tree Grammar Generator and Extensions; Generic Kalman Filter Software; Alignment Stage for a Cryogenic Dilatometer; Rugged Iris Mechanism; Treatments To Produce Stabilized Aluminum Mirrors for Cryogenic Uses; Making AlNx Tunnel Barriers Using a Low-Energy Nitrogen-Ion Beam; Making Wide-IF SIS Mixers with Suspended Metal-Beam Leads; Sol-Gel Glass Holographic Light-Shaping Diffusers; Automated Counting of Particles To Quantify Cleanliness; Phase Correction for GPS Antenna with Nonunique Phase Center; Compact Infrasonic Windscreen; Broadband External-Cavity Diode Laser; High-Efficiency Solar Cells Using Photonic-Bandgap Materials; Generating Solid Models from Topographical Data; Computationally Lightweight Air-Traffic-Control Simulation; Spool Valve for Switching Air Flows Between Two Beds; Partial Model of Insulator/ Insulator Contact Charging; Asymmetric Electrostatic Radiation Shielding for Spacecraft; and Reusable Hybrid Propellant Modules for Outer-Space Transport.

Source record↗

NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗

NASA NTRS2019

Fusion-Enabled Pluto Orbiter and Lander

Direct Fusion Drive (DFD) is a unique fusion engine concept based on the Princeton Field-Reversed Configuration (PFRC) fusion reactor conceived by Dr. Sam Cohen of the Princeton Plasma Physics Laboratory. DFD would enable the Pluto orbiter and lander context mission and more broadly enable true "rapid transit" to outer-planet and near interstellar space. The truly game-changing levels of thrust and power in a modestly sized package could integrate with our current launch infrastructure while radically expanding the science capability of these missions. Our Phase I was our first funded work on the DFD, with previous work at PSS occurring only under internal R&D. We established the feasibility of our Pluto mission trajectories using straight-line and planar models, including a departure spiral from Earth and insertion at Pluto. We developed our first thrust and specific impulse model using the results of the UEDGE multi-fluid code. Our specific power model was improved. During this Phase II effort, we continued our efforts to increase the fidelity of the designs for the RF, magnet, and shielding subsystems. Dedicated thrust augmentation experiments were run on the PFRC experiment, using a supersonic gas puffing valve. For the first time, we analyzed the design of a closed-loop operation mode and estimated the hardware that would be required for a dual-mode engine. In an exciting new development, we have invented a new thermophotovoltaic thermal conversion method that has the potential to have efficiencies of a Brayton or Stirling system. Our report presents details of these analyses. Our roadmap to bringing DFD to flight predicts that with sufficient support, a first flight unit could be built by 2040. We anticipate that three machine generations are required before this point: a ~1 T PFRC-3 machine hitting new plasma temperature and density levels, a ~5 T PFRC-4 machine with first demonstration of D-3He fusion, and a flight. (The current experiment, PFRC-2, is limited to about 0.1 T). In order to achieve a flight in 2035-2040, the TRL of the supporting systems must be increased in parallel, including low mass radiators, cryogenic propellant storage, and large (>100 kW) thermal conversion systems. Fortunately, many of these systems are dual-use and are required for other technologies including fission systems, so DFD would contribute to and benefit from those programs. This NIAC support and the results of our work have led to multiple follow-on contracts. We won two NASA STTRs focused on DFD subsystems, one on the RF system and one on the superconducting magnets, and our magnet STTR is now midway through a Phase II. We will be receiving a superconducting test magnet for experiments at PPPL this summer. In addition, we won an ARPA-E OPEN grant; this is the first time that OPEN has included fusion technologies, and we are part of a cohort of three alternative fusion companies now supported directly by DOE. We are very optimistic that if we are able to meet our experimental milestones in the next 12-18 months, we will be competitive for additional DOE grants to build PFRC-3.

Propulsion↗

NASA NTRSDate not supplied

Propulsive trajectory optimization to minimize surface contamination

MOTIVATION: We present an optimization technique for propulsive vehicles that autonomously minimizes contamination during surface approach and landing. In addition to short-range hoppers, the optimization technique is also fully applicable to traditional orbit-to-surface landers. This study addresses scenarios where surface alterations from propulsion events are counterproductive or hazardous to the mission objectives. This is of immediate interest for landers (whether human or robotic), that may rely on pristine soils collected in the immediate vicinity of landing sites to accomplish science investigations, mining, or ISRU surface operations. Such missions are averse to various surface-plume interactions such as thermal scoring, physical agitation, and contamination. The capability can be applied with minimal impact to the baseline mission concept. METHODS: Optimization algorithms have been developed to calculate descent trajectories and maneuvers, thrust magnitude, and attitude for various mission cases. These parameters are determined as an optimal solution when minimizing either fuel consumption, contamination deposited at the landing site, or some weighted combination of both. Among constraints imposed on the solution, we examined pitch rate, vertical takeoff and vertical landing (VTVL) requirements, size of the contamination zone, and minimum ground clearance during flight. This tool provides unique, non-intuitive solutions and can be a valuable resource for mission planners. RESULTS: A variety of agile trajectory solutions were obtained, each yielding different reductions in landing site contamination and corresponding to only modest increases in fuel consumption. Several optimal trajectories were obtained by varying the contamination weight in the fitness function. As expected, when the contamination weight is zero, the trajectory appears close to parabolic since the optimization scheme only attempts to minimize for fuel utilization, yielding essentially, the expected ballistic trajectory. Notably for contamination weights greater than zero, trajectory inflections are observed in the descent phase, which manifests as hovering or additional, mini “pseudo hops” before the final touchdown. A trajectory inflection is characterized by arresting the majority of the spacecraft vertical velocity component at a coordinate outside of the landing target, and without violating ground clearance constraints. FUTURE WORK: Our optimization technique is ready for laboratory or field demonstrations to validate the sophisticated maneuvering solutions obtained for fuel optimization and surface preservation. An appropriate testbed would validate the optimal guidance algorithms, the navigation system, and sensor suite by emulating vehicle flight in closed loop robotic tests. Critically, these algorithms could then be ported to flight software for implementation.

surface contamination↗

NASA NTRS2016

Development and Calibration of a Field-Deployable Microphone Phased Array for Propulsion and Airframe Noise Flyover Measurements

A new aeroacoustic measurement capability has been developed consisting of a large channelcount, field-deployable microphone phased array suitable for airframe noise flyover measurements for a range of aircraft types and scales. The array incorporates up to 185 hardened, weather-resistant sensors suitable for outdoor use. A custom 4-mA current loop receiver circuit with temperature compensation was developed to power the sensors over extended cable lengths with minimal degradation of the signal to noise ratio and frequency response. Extensive laboratory calibrations and environmental testing of the sensors were conducted to verify the design's performance specifications. A compact data system combining sensor power, signal conditioning, and digitization was assembled for use with the array. Complementing the data system is a robust analysis system capable of near real-time presentation of beamformed and deconvolved contour plots and integrated spectra obtained from array data acquired during flyover passes. Additional instrumentation systems needed to process the array data were also assembled. These include a commercial weather station and a video monitoring / recording system. A detailed mock-up of the instrumentation suite (phased array, weather station, and data processor) was performed in the NASA Langley Acoustic Development Laboratory to vet the system performance. The first deployment of the system occurred at Finnegan Airfield at Fort A.P. Hill where the array was utilized to measure the vehicle noise from a number of sUAS (small Unmanned Aerial System) aircraft. A unique in-situ calibration method for the array microphones using a hovering aerial sound source was attempted for the first time during the deployment.

Humphreys, William M., Jr.↗

NASA NTRS2014

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

The Solar Dynamics Observatory Helioseismic and Magnetic Imager (SDOHMI) filtergrams, taken at six wavelengths around the Fe I 6173.3 Å line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multiheight velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection-simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams.We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30 40 km above (line center) and 30 40 km below (average wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic-cutoff frequency.

velocity fields↗

NASA NTRS2000

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗

NASA NTRS2000

Recent Trends in the Analysis of Quasioptical Systems

The recent trend in microwave instruments is the use of multiple millimeter and submillimeter wavelength bands. These systems are typically analyzed by using physical optics, Gaussian beams or ray tracing techniques. Physical optics offers high accuracy at the expense of computation time. This trade-off becomes particularly apparent in the analysis of multiple reflector antennas, such as beam waveguide antennas, where physical optics is used to compute the current on each reflector from the current on the previous reflector. At the other end of the spectrum is ray tracing approaches that ignore diffraction effects entirely. These methods are fast but sacrifice the ability to predict some effects accurately. An intermediate approach is to use an appropriate set of expansion functions to model the field between the reflectors. If the set is chosen wisely only a few coefficients need to be determined from each reflector current. The field is then computed at the next reflector through the use of the expansion functions and their coefficients rather than by using the previous reflector current. For a beam waveguide system with no enclosing tubes an excellent set of expansion functions is the Gaussian beam mode set. In many cases a preliminary design which includes the effects on diffraction may be obtained by considering only the fundamental mode and a thin lens model for the reflectors. Higher-order modes are included to model the effects of the curved reflector, which include asymmetric distortion of the beam, cross polarization, and beam truncation. This paper describes a computer code implementing higher-order Gaussian beam scattering by multiple reflector systems. There are four essential steps in the algorithm. (1) Compute the current on the first reflector using physical optics using either a feed model or by an incident set of Gaussian beam modes. (2) Find the direction of propagation for the reflected Gaussian beam-set using ray tracing. (3) Determine the waist size and location for the output beam set by examining the amplitude and phase distribution of the current on the reflector. (4) Compute the amplitudes of the individual modes in the output mode set. These steps are then repeated for each addition reflector in the chain. In each of these cases the previous Gaussian beam set provides the input field for the current calculation. Details of the four steps discussed above will be discussed. Examples will compare results from the Gaussian beam approach to pure physical optics, illustrating both its merits and limitations. Hybrid approaches capable of eliminating some of the limitations will also be discussed.

Hoppe, D. J.↗

NASA NTRS2014

Open Circuit Resonant (SansEC) Sensor Technology for Lightning Mitigation and Damage Detection and Diagnosis for Composite Aircraft Applications

Traditional methods to protect composite aircraft from lightning strike damage rely on a conductive layer embedded on or within the surface of the aircraft composite skin. This method is effective at preventing major direct effect damage and minimizes indirect effects to aircraft systems from lightning strike attachment, but provides no additional benefit for the added parasitic weight from the conductive layer. When a known lightning strike occurs, the points of attachment and detachment on the aircraft surface are visually inspected and checked for damage by maintenance personnel to ensure continued safe flight operations. A new multi-functional lightning strike protection (LSP) method has been developed to provide aircraft lightning strike protection, damage detection and diagnosis for composite aircraft surfaces. The method incorporates a SansEC sensor array on the aircraft exterior surfaces forming a "Smart skin" surface for aircraft lightning zones certified to withstand strikes up to 100 kiloamperes peak current. SansEC sensors are open-circuit devices comprised of conductive trace spiral patterns sans (without) electrical connections. The SansEC sensor is an electromagnetic resonator having specific resonant parameters (frequency, amplitude, bandwidth & phase) which when electromagnetically coupled with a composite substrate will indicate the electrical impedance of the composite through a change in its resonant response. Any measureable shift in the resonant characteristics can be an indication of damage to the composite caused by a lightning strike or from other means. The SansEC sensor method is intended to diagnose damage for both in-situ health monitoring or ground inspections. In this paper, the theoretical mathematical framework is established for the use of open circuit sensors to perform damage detection and diagnosis on carbon fiber composites. Both computational and experimental analyses were conducted to validate this new method and system for aircraft composite damage detection and diagnosis. Experimental test results on seeded fault damage coupons and computational modeling simulation results are presented. This paper also presents the shielding effectiveness along with the lightning direct effect test results from several different SansEC LSP and baseline protected and unprotected carbon fiber reinforced polymer (CFRP) test panels struck at 40 and 100 kiloamperes following a universal common practice test procedure to enable damage comparisons between SansEC LSP configurations and common practice copper mesh LSP approaches. The SansEC test panels were mounted in a LSP test bed during the lightning test. Electrical, mechanical and thermal parameters were measured during lightning attachment and are presented with post test nondestructive inspection comparisons. The paper provides correlational results between the SansEC sensors computed electric field distribution and the location of the lightning attachment on the sensor trace and visual observations showing the SansEC sensor's affinity for dispersing the lightning attachment.

Szatkowski, George N.↗

NASA NTRS2021

A Study into Validating a Coupled Method of Characteristics and Direct Simulation Monte Carlo Method Against Empirical Data

The following will outline the methodology and results of validating a coupled Method of Characteristics (MOC) and Direct Simulation Monte Carlo (DSMC) method. This research focused specifically on modeling plume impingement, induced by Reaction Control System (RCS) thrusters that flew on the National Aeronautics and Space Administration’s (NASA’s) space shuttle Discovery. For each simulation, the continuum portion of the RCS thruster was simulated using MOC for solving hyperbolic Partial Differential Equations (PDEs) and computed with the NASA code, Reacting and Multi-phase Program (RAMP). The solution was then implemented as a starting condition into the NASA DSMC code, Direct Simulation and Monte Carlo Analysis Code (DAC). Typically, DSMC models rely on code-to-code validation for fidelity. The significance of this research is in its ability to validate its models against empirical data. Prior to computing solutions for these simulations, the mesh size and structure were optimized and many variants of DSMC input parameters were iterated on in order to acquire a reliable, mesh-independent, fully optimized numerical solution. This research will discuss the mathematical formulation of MOC for nozzle flow and DSMC for rarefied gases. Additionally, it will provide an explanation of how to implement these mathematical concepts into the two solvers: RAMP and DAC. Ultimately, this research will demonstrate that the overall process illustrated produces results in good agreement with empirical data. As a consequence, the methodology presented is granted an increased level of confidence and will greatly contribute to the aerospace industry and its effort in understanding and predicting rarefied flow fields.

Direct Simulation Monte Carlo Analysis Code↗

NASA NTRS2021

A Study into Validating A Coupled Method of Characteristics And Direct Simulation Monte Carlo Method Against Empirical Data

DAC↗

NASA NTRS2000

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 optomechanical 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.↗