Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesh data”

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 1,531 records · Page 85

Coupled aerodynamic and acoustical predictions for turboprops

To predict the noise fields for proposed turboprop airplanes, an existing turboprop noise code by Farassat has been modified to accept blade pressure inputs from a three-dimensional aerodynamic code. A Euler-type code can handle the nonlinear transonic flow of these high-speed, highly swept blades. This turbofan code was modified to allow the calculation mesh to extend to about twice the blade radius and to apply circumferential periodicity rather than solid-wall boundary conditions on the blade in the region between the blade tip and the outer shroud. Outputs were added for input to the noise prediction program and for color contour plots of various flow variables. The Farassat input subroutines were modified to read files of blade coordinates and predicted surface pressures. Aerodynamic and acoustic results are shown for the SR-3 model blade. Comparison of the acoustic predicted results with measured data show good agreement.

Clark, Bruce J.↗

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗

Toward Time Synchronization in Delay Tolerant Network based Solar System Internetworking

The expanding presence in space will place an increased dependency on networked communications– a scalable communications infrastructure; that is, the Solar System Internet (SSI). Upcoming developments towards a SSI include NASA’s upcoming LunaNet, or lunar Internet, which provides multi-hop multi-path communications using Delay Tolerant Networking (DTN). DTN has been an active area of research and development, particularly in routing, security, and optimization. DTNs are marked by mobility, disconnection, and a wide variance of latencies (propagation and processing delays). In this paper, we outline progress towards a theory of time synchronization across such a network. An underlying assumption of DTN is that the network is time synchronized already, rather than synchronization being provided as a service. While this is necessary for schedule-based routing, which is necessarily prevalent in DTNs, it is so deeply ingrained as to be built into the primary unit of data in DTNs– the bundle. Indeed, a bundle’s creation timestamp and its time to live (called the lifetime) are based on time, and there are special recommendations for systems that lack accurate clocks. The assumption of time synchronization makes sense when limiting considerations to smaller-scale and more traditional space communication. However, just as end-to-end connectivity cannot be guaranteed in DTNs, neither can access to a reference or authoritative clock. In this more general case, it might be necessary to synchronize over time-varying meshes, and perhaps even to consider relativistic effects. Moreover, by imposing synchronization restrictions in order to sustain a network, the effectiveness of the network to achieve scalability will be necessarily muted. To work towards a time synchronization theory for DTNs, we build upon past successes in modeling DTNs using time-varying graphs and sheaves. This includes error and limitation estimation, which allows one to define domains over which schedule-based routing is possible, up to some threshold sensitivity. Despite the theoretical nature of these results, the approaches taken are also algorithmic, and hence lend themselves to practical implementations. The paper concludes with comparisons of the various methods along with suggestions for future work.

Delay Tolerant Networking↗

Development of a Multiphase Adjoint Capability in OpenFOAM

Thermal hydraulic predictions of flow regimes within nuclear reactor cores require significant computational resources and accuracy to ensure that the core design does not violate the thermal limits of materials. Multiphase computational fluid dynamics (CFD) has become a valuable computational capability for predicting flow characteristics as they relate to applications such as nuclear reactor safety and mixing phenomena around spacer grids within fuel assemblies. Multiphase CFD requires closure relationships for both the effects of turbulence as well as the transfer of mass, momentum, and energy between phases. Understanding the sensitivity of the parameters used in closure relationships is important for improving accuracy, but performing rigorous sensitivity analysis is computationally expensive. Adjoint mathematics have been utilized for reactor physics and fluid dynamics applications with respect to verification, sensitivity analysis, adaptive mesh refinement, and optimization. The usage of adjoint methods to predict the sensitivities of a variety of closure relationship parameters within a given problem has been shown to return sensitivity responses with significantly fewer computations than normal 'large-scale' methods. As the complexity of multiphase CFD methods increases, the demand for accurate verification tools and improvements in computational efficiency also increases, making the development of adjoint capabilities an attractive option. Closure relationships that describe interphase momentum physics in multiphase CFD modelling rely heavily on prescribed constants determined from experimental data. Due to the computational efficiency of an adjoint approach, the time to derive these closure relationships could be dramatically reduced, and data assimilation techniques could be used to return probability density functions of closure models derived using experimental data. This work develops a multiphase adjoint capability for the Eulerian-Eulerian two field, two phase equations. Specifically, the adjoint capability has been developed for the sub-cooled nucleate boiling solver boilEulerFoam developed by Dr. Nam Dinh of NC State University within the open-source CFD library OpenFOAM.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Estimating uncertainties in statistics computed from direct numerical simulation

Rigorous assessment of uncertainty is crucial to the utility of direct numerical simulation (DNS) results. Uncertainties in the computed statistics arise from two sources: finite statistical sampling and the discretization of the Navier–Stokes equations. Due to the presence of non-trivial sampling error, standard techniques for estimating discretization error (such as Richardson extrapolation) fail or are unreliable. This work provides a systematic and unified approach for estimating these errors. First, a sampling error estimator that accounts for correlation in the input data is developed. Then, this sampling error estimate is used as part of a Bayesian extension of Richardson extrapolation in order to characterize the discretization error. These methods are tested using the Lorenz equations and are shown to perform well. These techniques are then used to investigate the sampling and discretization errors in the DNS of a wall-bounded turbulent flow at Reτ ≈ 180. Both small (Lx/δ × Lz/δ = 4π × 2π) and large (Lx/δ × Lz/δ = 12π × 4π) domain sizes are investigated. For each case, a sequence of meshes was generated by first designing a “nominal” mesh using standard heuristics for wall-bounded simulations. These nominal meshes were then coarsened to generate a sequence of grid resolutions appropriate for the Bayesian Richardson extrapolation method. In addition, the small box case is computationally inexpensive enough to allow simulation on a finer mesh, enabling the results of the extrapolation to be validated in a weak sense. For both cases, it is found that while the sampling uncertainty is large enough to make the order of accuracy difficult to determine, the estimated discretization errors are quite small. This indicates that the commonly used heuristics provide adequate resolution for this class of problems. However, it is also found that, for some quantities, the discretization error is not small relative to sampling error, indicating that the conventional wisdom that sampling error dominates discretization error for this class of simulations needs to be reevaluated.

Oliver, Todd A.↗

SU-E-T-565: RAdiation Resistance of Cancer CElls Using GEANT4 DNA: RACE

Purpose: The objective of the RACE project is to develop a comparison between Monte Carlo simulation using the Geant4-DNA toolkit and measurements of radiation damage on 3D melanoma and chondrosarcoma culture cells coupled with gadolinium nanoparticles. We currently expose the status of the developments regarding simulations. Methods: Monte Carlo studies are driven using the Geant4 toolkit and the Geant4-DNA extension. In order to model the geometry of a cell population, the opensource CPOP++ program is being developed for the geometrical representation of 3D cell populations including a specific cell mesh coupled with a multi-agent system. Each cell includes cytoplasm and nucleus. The correct modeling of the cell population has been validated with confocal microscopy images of spheroids. The Geant4 Livermore physics models are used to simulate the interactions of a 250 keV X-ray beam and the production of secondaries from gadolinium nanoparticles supposed to be fixed on the cell membranes. Geant4-DNA processes are used to simulate the interactions of charged particles with the cells. An atomistic description of the DNA molecule, from PDB (Protein Data Bank) files, is provided by the so-called PDB4DNA Geant4 user application we developed to score energy depositions in DNA base pairs and sugar-phosphate groups. Results: At the microscopic level, our simulations enable assessing microscopic energy distribution in each cell compartment of a realistic 3D cell population. Dose enhancement factors due to the presence of gadolinium nanoparticles can be estimated. At the nanometer scale, direct damages on nuclear DNA are also estimated. Conclusion: We successfully simulated the impact of direct radiations on a realistic 3D cell population model compatible with microdosimetry calculations using the Geant4-DNA toolkit. Upcoming validation and the future integration of the radiochemistry module of Geant4-DNA will propose to correlate clusters of ionizations with in vitro experiments. All those developments will be released publicly. This work was supported by grants from Plan Cancer 2009-2013 French national initiative managed by INSERM (Institut National de la Sante et de la Recherche Medicale)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Parallel High-Order Anisotropic Meshing Using Discrete Metric Tensors

This paper presents a metric-aligned meshing algorithm that relies on the Lp-Centroidal Voronoi Tesselation approach. A prototype of this algorithm was first presented at the Scitech conference of 2018 and this work is an extension to that paper. At the end of the previously presented work, a set of problems were mentioned which we are trying to address in this paper. First, we show a significant improvement in code performance since we were limited to present relatively benign (analytical) test cases. Second, we demonstrate here that we are able to rely on discrete metric data that is delivered by a Computational Fluid Dynamics (CFD) solver. Third, we demonstrate how to generate high-order curved elements that are aligned with the underlying discrete metric field.

Ekelschot, Dirk↗

WE-EF-BRA-01: A Dual-Use Optical Tomography System for Small Animal Radiation Research Platform (SARRP)

Purpose: We develop a novel dual-use configuration for a tri-modality, CBCT/bioluminescence tomography(BLT)/fluorescence tomography(FT), imaging system with the SARRP that can function as a standalone system for longitudinal imaging research and on-board the SARRP to guide irradiation. BLT provides radiation guidance for soft tissue target, while FT offers functional information allowing mechanistic investigations. Methods: The optical assembly includes CCD camera, lens, filter wheel, 3-way mirrors, scanning fiber system and light-tight enclosure. The rotating mirror system directs the optical signal from the animal surface to the camera at multiple projection over 180 degree. The fiber-laser system serves as the external light source for the FT application. Multiple filters are used for multispectral imaging to enhance localization accuracy using BLT. SARRP CBCT provides anatomical information and geometric mesh for BLT/FT reconstruction. To facilitate dual use, the 3-way mirror system is cantilevered in front of the camera. The entire optical assembly is driven by a 1D linear stage to dock onto an independent mouse support bed for standalone application. After completion of on-board optical imaging, the system is retracted from the SARRP to allow irradiation of the mouse. Results: A tissue-simulating phantom and a mouse model with a luminescence light source are used to demonstrate the function of the dual-use optical system. Feasibility data have been obtained based on a manual-docking prototype. The center of mass of light source determined in living mouse with on-board BLT is within 1±0.2mm of that with CBCT. The performance of the motorized system is expected to be the same and will be presented. Conclusion: We anticipate the motorized dual use system provide significant efficiency gain over our manual docking and off-line system. By also supporting off-line longitudinal studies independent of the SARRP, the dual-use system is a highly efficient and cost-effective platform to facilitate optical imaging for pre-clinical radiation research. The work is supported by NIH R01CA158100 and Xstrahl Ltd. Drs. John Wong and Iulian Iordachita receive royalty payment from a licensing agreement between Xstrahl Ltd and Johns Hopkins University. John Wong also has a consultant agreement with Xstrahl Ltd.

60 APPLIED LIFE SCIENCES↗

Selected computations of transonic cavity flows

An efficient diagonal scheme implemented in an overset mesh framework has permitted the analysis of geometrically complex cavity flows via the Reynolds averaged Navier-Stokes equations. Use of rapid hyperbolic and algebraic grid methods has allowed simple specification of critical turbulent regions with an algebraic turbulence model. Comparisons between numerical and experimental results are made in two dimensions for the following problems: a backward-facing step; a resonating cavity; and two quieted cavity configurations. In three-dimensions the flow about three early concepts of the stratospheric Observatory For Infrared Astronomy (SOFIA) are compared to wind-tunnel data. Shedding frequencies of resolved shear layer structures are compared against experiment for the quieted cavities. The results demonstrate the progress of computational assessment of configuration safety and performance.

Atwood, Christopher A.↗

Computer Aided Modeling to Determine the Effectiveness of Resistive Exercises as Countermeasures for Bone Mineral Density Loss

Due to the loss of gravitational loading, astronauts have a tendency to lose bone mineral density in their lumbar spine and lower extremities on orbit. NASA requires astronauts to perform exercises during space flight to help reduce the amount of demineralization. To test these exercises on earth, 17 week bed rest studies are conducted that consist of specific diet and exercise regimes. Developing a finite element model of these exercises will help to quantify the stress distribution imposed by of each of these exercises. To help develop this model, MRI images are acquired from individuals participating in the bed rest studies. The MRIs can be used to create a subject specific model of each individual for testing. The MRIs are processed in the Magnetic Resonance Imaging Data Transfer System program to develop a three-dimensional finite element model of the femur for evaluation. Modifications were made to the MRIDTS that simplified the model creation process. These modifications made it possible to construct two separate models of different portions of a bone simultaneously and then later connect them manually. This helped alleviate the warping problem associated with the drastic changes in geometry found in some body parts, such as the joints. The code was also modified to incorporate material properties of various bone components into the model. Interior meshing was also incorporated into the program to allow for both the cortical shell and the entire bone to be modeled. A prototype model of the right femur of an adult female is being constructed and tested to determine the feasibility of finite element analysis as a tool for evaluating exercise effectiveness. The model is being run through the ANSYS finite element program on the Alabama Super Computer Network. After the model is validated, models of bedrest subjects can be generated to investigate exercise countermeasures.

Murphy, Benjamin M.↗

How much can we learn from a merging cold front cluster? Insights from X-ray temperature and radio maps of A3667

The galaxy cluster A3667 is an ideal laboratory to study the plasma processes in the intracluster medium. High-resolution Chandra X-ray observations show a cold front in A3667. At radio wavelengths, A3667 reveals a double radio-relic feature in the outskirts of the cluster. These suggest multiple merger events in this cluster. In this paper, we analyze the substantial archival X-ray observations of A3667 from the Chandra X-ray Observatory and compare these with existing radio observations as well as state-of-the-art adaptive mesh refinement MHD cosmological simulations using Enzo. We have used two temperature map making techniques, weighted Voronoi tessellation and adaptive circular binning, to produce the high-resolution and largest field-of-view temperature maps of A3667. These high-fidelity temperature maps allow us to study the X-ray shocks in the cluster using a new two-dimensional shock-finding algorithm. We have also estimated the Mach numbers from the shocks inferred from previous ATCA radio observations. The combined shock statistics from the X-ray and radio data are in agreement with the shock statistics in a simulated MHD cluster. We have also studied the profiles of the thermodynamic properties across the cold front using ∼447 ks from the combined Chandra observations on A3667. Our results show that the stability of the cold front in A3667 can be attributed to the suppression of the thermal conduction across the cold front by a factor of ∼100-700 compared to the classical Spitzer value.

79 ASTRONOMY AND ASTROPHYSICS↗

Shortening the Design and Certification Cycle for Additively Manufactured Materials by Improved Mesoscale Simulations and Validation Experiments: Fiscal Year 2019 Status Report.

This report outlines the fiscal year (FY) 2019 status of an ongoing multi-year effort to develop a general, microstructurally-aware, continuum-level model for representing the dynamic response of material with complex microstructures. This work has focused on accurately representing the response of both conventionally wrought processed and additively manufactured (AM) 304L stainless steel (SS) as a test case. Additive manufacturing, or 3D printing, is an emerging technology capable of enabling shortened design and certification cycles for stockpile components through rapid prototyping. However, there is not an understanding of how the complex and unique microstructures of AM materials affect their mechanical response at high strain rates. To achieve our project goal, an upscaling technique was developed to bridge the gap between the microstructural and continuum scales to represent AM microstructures on a Finite Element (FE) mesh. This process involves the simulations of the additive process using the Sandia developed kinetic Monte Carlo (KMC) code SPPARKS. These SPPARKS microstructures are characterized using clustering algorithms from machine learning and used to populate the quadrature points of a FE mesh. Additionally, a spall kinetic model (SKM) was developed to more accurately represent the dynamic failure of AM materials. Validation experiments were performed using both pulsed power machines and projectile launchers. These experiments have provided equation of state (EOS) and flow strength measurements of both wrought and AM 304L SS to above Mbar pressures. In some experiments, multi-point interferometry was used to quantify the variation is observed material response of the AM 304L SS. Analysis of these experiments is ongoing, but preliminary comparisons of our upscaling technique and SKM to experimental data were performed as a validation exercise. Moving forward, this project will advance and further validate our computational framework, using advanced theory and additional high-fidelity experiments.

36 MATERIALS SCIENCE↗

Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow

The proposed paper reports advances in developing a method for high Reynolds number compressible viscous flow simulations using a Cartesian cut-cell method with embedded boundaries. This preliminary work focuses on accuracy of the discretization near solid wall boundaries. A model problem is used to investigate the accuracy of various difference stencils for second derivatives and to guide development of the discretization of the viscous terms in the Navier-Stokes equations. Near walls, quadratic reconstruction in the wall-normal direction is used to mitigate mesh irregularity and yields smooth skin friction distributions along the body. Multigrid performance is demonstrated using second-order coarse grid operators combined with second-order restriction and prolongation operators. Preliminary verification and validation for the method is demonstrated using flat-plate and airfoil examples at compressible Mach numbers. Simulations of flow on laminar and turbulent flat plates show skin friction and velocity profiles compared with those from boundary-layer theory. Airfoil simulations are performed at laminar and turbulent Reynolds numbers with results compared to both other simulations and experimental data

Berger, Marsha↗

Serpent and MCNP Calculations of the Energy Deposition in the Transformational Challenge Reactor

This paper focuses on the calculation of the energy deposition in the Transformational Challenge Reactor by two major Monte Carlo codes: Serpent and MCNP. The first software computation relies on Kinetic Energy Released per unit Mass (KERMA) factors while the second one relies on Q-values. The results from these two independent computation methodologies are in very good agreement; however, Serpent runs much faster than MCNP (for the same computational model) and allows for a detailed energy deposition distribution from a 1-mm-side square mesh with a relative statistical error between 0.5% and 1%. This detailed energy deposition is suitable for multiphysics analyses aimed at design optimizations. In order to calculate the energy deposition, Serpent needs enhanced ACE files (distributed by the software developers). Unlike other Monte Carlo software that uses inputs based on Python or Java languages, the Serpent input syntax is very similar to that of MCNP; a Python script can convert a MCNP input to a Serpent input in seconds. For simulations not requiring the calculation of the energy deposition, Serpent can also read nuclear data from MCNP ACE files, which eventually improves the comparison of the results of the two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Localized Evaluation for Constructing Discrete Vector Fields

Topological abstractions offer a method to summarize the behavior of vector fields, but computing them robustly can be challenging due to numerical precision issues. One alternative is to represent the vector field using a discrete approach, which constructs a collection of pairs of simplices in the input mesh that satisfies criteria introduced by Forman's discrete Morse theory. While numerous approaches exist to compute pairs in the restricted case of the gradient of a scalar field, state-of-the-art algorithms for the general case of vector fields require expensive optimization procedures. This paper introduces a fast, novel approach for pairing simplices of two-dimensional, triangulated vector fields that do not vary in time. The key insight of our approach is that we can employ a local evaluation, inspired by the approach used to construct a discrete gradient field, where every simplex in a mesh is considered by no more than one of its vertices. Specifically, we observe that for any edge in the input mesh, we can uniquely assign an outward direction of flow. We can further expand this consistent notion of outward flow at each vertex, which corresponds to the concept of a downhill flow in the case of scalar fields. Working with outward flow enables a linear-time algorithm that processes the (outward) neighborhoods of each vertex one-by-one, similar to the approach used for scalar fields. Here, we couple our approach to constructing discrete vector fields with a method to extract, simplify, and visualize topological features. Empirical results on analytic and simulation data demonstrate drastic improvements in running time, produce features similar to the current state-of-the-art, and show the application of simplification to large, complex flows.

97 MATHEMATICS AND COMPUTING↗

Reduction of rotor-turbulence interaction noise in static fan noise testing

Three effective inlet cleanup methods are described. A flared reverse cone inlet is used to eliminate wakes from the fan casings and/or probe supports. Boundary layer suction is employed ahead of the fan rotor and on the outer flare of the cone to reduce the boundary layer turbulence and remove any residual wakes. To reduce the midstream turbulent intensity and length scales, a turbulence control structure is used which is constructed with both a layer of honeycomb and a fine mesh screen. The effects of these cleanup methods are quantified by measuring the far-field noise in an anechoic chamber, using a high-speed 20-in.-diam fan of the current high passband type. The changes in the turbulent field impinging on this rotor are quantified by mapping the streamwise and transverse turbulent properties (spectra, intensity and length scale) with crossed hot film probes. It is convincingly shown that it is possible to clean up the inlet flow of a static fan noise test facility to a point where the static acoustic data simulate flight data.

Kantola, R. A.↗

Medium-range, objective predictions of thunderstorm location and severity for aviation

This paper presents a computerized technique for medium-range (12-48h) prediction of both the location and severity of thunderstorms utilizing atmospheric predictions from the National Meteorological Center's limited-area fine-mesh model (LFM). A regional-scale analysis scheme is first used to examine the spatial and temporal distributions of forecasted variables associated with the structure and dynamics of mesoscale systems over an area of approximately 10 to the 6th sq km. The final prediction of thunderstorm location and severity is based upon an objective combination of these regionally analyzed variables. Medium-range thunderstorm predictions are presented for the late afternoon period of April 10, 1979, the day of the Wichita Falls, Texas tornado. Conventional medium-range thunderstorm forecasts, made from observed data, are presented with the case study to demonstrate the possible application of this objective technique in improving 12-48 h thunderstorm forecasts for aviation.

Wilson, G. S.↗

Remote Sensing of Volcanic Clouds: Sulfur Gases and Plume Top Topography

New absorption line parameters for H2S were published and submitted to the Gestion et Etude des Informations Spectroscopiques Atmospheriques (GEISA) and high resolution transmission molecular absorption (HITRAN) databases. These new absorption line parameters will make it possible to use observations from the future Tropospheric Emission Spectrometer (TES) instrument [Earth Observing System (EOS) Chemistry Mission (CHEM) platform] to make more accurate H2S measurements if it observes an H2S-rich volcanic cloud. H2S is the second most abundant volcanic sulfur gas, and like SO2, it also converts to H2SO4 aerosols and can have a climate impact. A paper on the Moderate-resolution Imaging-Spectroradiometer (MODIS) SO2 alert is being revised. New aspects in the revision include verification of the SO2 alert during the EOS mission; factors affecting SO2 detection at thermal infrared, ultraviolet, and microwave wavelengths; radiative transfer tests; more description of satellite instruments; and thermal surface alert installed for MODIS. Her research involves the use of remote sensing to generate maps of plume top altitude. This parameter is important for models of volcanic eruption, aircraft hazards, and climate impact. The topographic shape of the top surface of a volcanic plume can provide information necessary to understand the physics controlling the injection and dispersal of a volcanic plume in the atmosphere. Glaze et al. describe the application of a photoclinometric technique to volcanic plumes. The software algorithm has been improved to account for more general plume and illumination geometries and for easily extracting position information directly from Advanced Very High-Resolution Radiometer (AVHRR) level 1B data. Testing of the algorithm has focused on acquiring AVHRR data for a variety of volcanic plumes in an effort to identify problems with the software as well as model sensitivities. The plumes chosen were erupted from volcanoes at a variety of latitudes, at different times of day, and in different seasons. Another important issue that is being investigated is the importance of plume opacity. A diagram illustrates the method for deriving phi, the angle a wind-blown plume makes with the horizontal, in the direction perpendicular to the remote sensing scan. A figure shows a sample output of the algorithm, a wire mesh map of plume top topography. Additional information is contained in the original.

Crisp, Joy A.↗