Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials image simulations”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2000

TE/TM Simulations of Interferometric Measurements

Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.

Houshmand, Bijan↗
DOE OSTI2024

Implosion simulations and optimum times of radiographic Images in OMEGA high-density-carbon cylinder experiments

The Los Alamos Ignition Threshold team is planning direct-drive cylinder implosion experiments (CylDRT24B) at the OMEGA laser, scheduled for Thursday 23 May 2024, to investigate the effect of granular microstructure in high density carbon (HDC) shells. Granular microstructure is believed to play a crucial role in seeding perturbations, reducing compression, and enhancing mix in the implosion of HDC shells at the National Ignition Facility (NIF). But the ICF community lacks a detailed understanding of how grains induce these effects, owing partly to the practical difficulty of numerically simulating the behavior of tiny nanometer-scale granular structures in laboratory-scale experiments. Our planned experiments are aimed at acquiring data to help constrain simulations and calibrate reduced models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multi frame radiography of supersonic water jets interacting with a foil target

Pulsed-power-driven underwater electrical explosion of cylindrical or conical wire arrays produces supersonic water jets that emerge from a bath, propagating through the air above it. Interaction of these jets with solid targets may represent a new platform for attaining materials at high pressure (>1010 Pa) conditions in a university-scale laboratory. However, measurements of the internal structure of such jets and how they interact with targets are difficult optically due to large densities and density contrasts involved. We utilized multi-frame x-ray radiographic imaging capabilities of the ID19 beamline at the European Synchrotron Radiation Facility to explore the water jet and its interaction with a 50 μm thick copper foil placed a few mm from the surface of water. The jet was generated with a ∼130 kA-amplitude current pulse of ∼450 ns rise time applied to a conical wire array. X-ray imaging revealed a droplet-type structure of the jet with an average density of <400 kg/m3 propagating with a velocity of ∼1400 m/s. Measurements of deformation and subsequent perforation of the target by the jet suggested pressures at the jet–target interface of ∼5 × 109 Pa. The results were compared to hydrodynamic simulations for better understanding of the jet parameters and their interaction with the foil target. These results can be used in future research to optimize the platform, and extend it to larger jet velocities in the case of higher driving currents supplied to the wire array.

Physics↗
NASA NTRS2008

NASA Tech Briefs, June 2008

Topics covered include: Charge-Control Unit for Testing Lithium-Ion Cells; Measuring Positions of Objects Using Two or More Cameras; Lidar System for Airborne Measurement of Clouds and Aerosols; Radiation-Insensitive Inverse Majority Gates; Reduced-Order Kalman Filtering for Processing Relative Measurements; Spaceborne Processor Array; Instrumentation System Diagnoses a Thermocouple; Chromatic Modulator for a High-Resolution CCD or APS; Commercial Product Activation Using RFID; Cup Cylindrical Waveguide Antenna; Aerobraking Maneuver (ABM) Report Generator; ABM Drag_Pass Report Generator; Transformation of OODT CAS to Perform Larger Tasks; Visualization Component of Vehicle Health Decision Support System; Mars Reconnaissance Orbiter Uplink Analysis Tool; Problem Reporting System; G-Guidance Interface Design for Small Body Mission Simulation; DSN Scheduling Engine; Replacement Sequence of Events Generator; Force-Control Algorithm for Surface Sampling; Tool for Merging Proposals Into DSN Schedules; Micromachined Slits for Imaging Spectrometers; Fabricating Nanodots Using Lift-Off of a Nanopore Template; Making Complex Electrically Conductive Patterns on Cloth; Special Polymer/Carbon Composite Films for Detecting SO2; Nickel-Based Superalloy Resists Embrittlement by Hydrogen; Chemical Passivation of Li+-Conducting Solid Electrolytes; Organic/Inorganic Polymeric Composites for Heat-Transfer Reduction; Composite Cathodes for Dual-Rate Li-Ion Batteries; Improved Descent-Rate Limiting Mechanism; Alignment-Insensitive Lower-Cost Telescope Architecture; Micro-Resistojet for Small Satellites; Using Piezoelectric Devices to Transmit Power through Walls; Miniature Latching Valve; Apparatus for Sampling Surface Contamination; Novel Species of Non-Spore-Forming Bacteria; Chamber for Aerosol Deposition of Bioparticles; Hyperspectral Sun Photometer for Atmospheric Characterization and Vicarious Calibrations; Dynamic Stability and Gravitational Balancing of Multiple Extended Bodies; Simulation of Stochastic Processes by Coupled ODE-PDE; Cluster Inter-Spacecraft Communications; Genetic Algorithm Optimizes Q-LAW Control Parameters; Low-Impact Mating System for Docking Spacecraft; Non-Destructive Evaluation of Materials via Ultraviolet Spectroscopy; Gold-on-Polymer-Based Sensing Films for Detection of Organic and Inorganic Analytes in the Air; and Quantum-Inspired Maximizer.

Source record↗
DOE OSTI2020

Smart Scattering Scanning Near-Field Optical Microscopy

Scattering scanning near-field optical microscopy (s-SNOM) provides spectroscopic imaging from molecular to quantum materials with few nanometer deep subdiffraction limited spatial resolution. However, in its conventional implementation s-SNOM is slow to effectively acquire a series of spatio-spectral images, especially with large fields of view. This problem is further exacerbated for weak resonance contrast or when using light sources with limited spectral irradiance. Indeed, the generally limited signal-to-noise ratio prevents sampling a weak signal at the Nyquist sampling rate. Here, we demonstrate how acquisition time and sampling rate can be significantly reduced by using compressed sampling, matrix completion, and adaptive random sampling, while maintaining or even enhancing the physical or chemical image content. We use fully sampled real data sets of molecular, biological, and quantum materials as ground-truth physical data and show how deep under-sampling with a corresponding reduction of acquisition time by 1 order of magnitude or more retains the core s-SNOM image information. We demonstrate that a sampling rate of up to 6× smaller than the Nyquist criterion can be applied, which would provide a 30-fold reduction in the data required under typical experimental conditions. Furthermore, our smart s-SNOM approach is generally applicable and provides systematic full spatio-spectral s-SNOM imaging with a large field of view at high spectral resolution and reduced acquisition time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2007

NASA Tech Briefs, September 2007

Topics covered include; Rapid Fabrication of Carbide Matrix/Carbon Fiber Composites; Coating Thermoelectric Devices To Suppress Sublimation; Ultrahigh-Temperature Ceramics; Improved C/SiC Ceramic Composites Made Using PIP; Coating Carbon Fibers With Platinum; Two-Band, Low-Loss Microwave Window; MCM Polarimetric Radiometers for Planar Arrays; Aperture-Coupled Thin-Membrane L-Band Antenna; WGM-Based Photonic Local Oscillators and Modulators; Focal-Plane Arrays of Quantum-Dot Infrared Photodetectors; Laser Range and Bearing Finder With No Moving Parts; Microrectenna: A Terahertz Antenna and Rectifier on a Chip; Miniature L-Band Radar Transceiver; Robotic Vision-Based Localization in an Urban Environment; Programs for Testing an SSME-Monitoring System; Cathodoluminescent Source of Intense White Light; Displaying and Analyzing Antenna Radiation Patterns; Payload Operations Support Team Tools; Space-Shuttle Emulator Software; Soft Real-Time PID Control on a VME Computer; Analyzing Radio-Frequency Coverage for the ISS; Nanorod-Based Fast-Response Pressure-Sensitive Paints; Capacitors Would Help Protect Against Hypervelocity Impacts; Diaphragm Pump With Resonant Piezoelectric Drive; Improved Quick-Release Pin Mechanism; Designing Rolling-Element Bearings; Reverse-Tangent Injection in a Centrifugal Compressor; Inertial Measurements for Aero-assisted Navigation (IMAN); Analysis of Complex Valve and Feed Systems; Improved Path Planning Onboard the Mars Exploration Rovers; Robust, Flexible Motion Control for the Mars Explorer Rovers; Solar Sail Spaceflight Simulation; Fluorine-Based DRIE of Fused Silica; Mechanical Alloying for Making Thermoelectric Compounds; Process for High-Rate Fabrication of Alumina Nanotemplates; Electroform/Plasma-Spray Laminates for X-Ray Optics; An Automated Flying-Insect Detection System; Calligraphic Poling of Ferroelectric Material; Blackbody Cavity for Calibrations at 200 to 273 K; KML Super Overlay to WMS Translator; High-Performance Tiled WMS and KML Web Server; Modeling of Radiative Transfer in Protostellar Disks; Composite Pulse Tube; Photometric Calibration of Consumer Video Cameras; Criterion for Identifying Vortices in High- Pressure Flows; Amplified Thermionic Cooling Using Arrays of Nanowires; Delamination-Indicating Thermal Barrier Coatings; Preventing Raman Lasing in High-Q WGM Resonators; Procedures for Tuning a Multiresonator Photonic Filter; Robust Mapping of Incoherent Fiber-Optic Bundles; Extended-Range Ultrarefractive 1D Photonic Crystal Prisms; Rapid Analysis of Mass Distribution of Radiation Shielding; Modeling Magnetic Properties in EZTB; Deep Space Network Antenna Logic Controller; Modeling Carbon and Hydrocarbon Molecular Structures in EZTB; BigView Image Viewing on Tiled Displays; and Imaging Sensor Flight and Test Equipment Software.

Source record↗
DOE OSTI2024

Identification of Defects and the Origins of Surface Noise on Hydrogen–Terminated (100) Diamond

Near-surface nitrogen vacancy centres are critical to many diamond-based quantum technologies such as information processors and nanosensors. Surface defects play an important role in the design and performance of these devices. The targeted creation of defects is central to proposed bottom-up approaches to nanofabrication of quantum diamond processors, and uncontrolled surface defects may generate noise and charge trapping which degrade shallow NV device performance. Surface preparation protocols may be able to control the production of desired defects and eliminate unwanted defects, but only if their atomic structure can first be conclusively identified. This work uses a combination of scanning tunnelling microscopy (STM) imaging and first-principles simulations to identify several surface defects on H:C(100)—2 × 1 surfaces prepared using chemical vapour deposition (CVD). The atomic structure of these defects is elucidated, from which the microscopic origins of magnetic noise and charge trapping are determined based on the modeling of their paramagnetic properties and acceptor states. Rudimentary control of these deleterious properties is demonstrated through STM tip-induced manipulation of the defect structure. Furthermore, the results validate accepted models for CVD diamond growth by identifying key adsorbates responsible for the nucleation of new layers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Controlling solute channel formation using magnetic fields

Solute channel formation introduces compositional and microstructural variations in a range of processes, from metallic alloy solidification, to salt fingers in ocean and water reservoir flows. Applying an external magnetic field interacts with thermoelectric currents at solid/liquid interfaces generating additional flow fields. This thermoelectric (TE) magnetohydrodynamic (TEMHD) effect can impact on solute channel formation, via a mechanism recently drawing increasing attention. To investigate this phenomenon, we combined in situ synchrotron X-ray imaging and Parallel-Cellular-Automata-Lattice-Boltzmann based numerical simulations to study the characteristics of flow and solute transport under TEMHD. Observations suggest the macroscopic TEMHD flow appearing ahead of the solidification front, coupled with the microscopic TEMHD flow arising within the mushy zone are the primary mechanisms controlling plume migration and channel bias. Two TE regimes were revealed, each with distinctive mechanisms that dominate the flow. Further, we show that grain orientation modifies solute flow through anisotropic permeability. These insights led to a proposed strategy for producing solute channel-free solidification using a time-modulated magnetic field.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topology-dependent stability of vortex-antivortex structures

The non-trivial topology of magnetic structures such as vortices and skyrmions is considered as a key concept to explain the stability of those structures. The stability, dictated by non-trivial topology, provides great potential for device applications. Although it is a very critical scientific and technological issue, it is elusive to experimentally study the topology-dependent stability owing to the difficulties in establishing stably formed magnetic structures with different topologies. Here, we establish a platform for vortex-antivortex structures with different topological charges within Ni 80 Fe 20 rectangular elements thick enough to stabilize a unique three-dimensional magnetic structure with non-uniform magnetization along the thickness of the elements. The detailed magnetization configurations of the three-dimensional vortex-antivortex structures and their annihilations during their field-driven motions are investigated by utilizing magnetic transmission soft x-ray microscopy and micromagnetic simulation. We demonstrate that the stability of vortex-antivortex structures significantly depends on their topologies and the topology-dependent stability is associated with their different annihilation mechanisms. Furthermore, we believe that this work provides in-depth insight into the stability of magnetic structures and its topology dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Impact of ferromagnetic electrode length and thickness on Magnetic Tunnel Junction-Based Molecular Spintronic Devices (MTJMSD)

A knowledge gap exists about the impact of variation in length and thickness of ferromagnetic(FM) electrodes on molecular spintronics devices’ ability to manifest molecular coupling impact. Magnetic Tunnel Junction-Based Molecular Spintronic Devices (MTJMSDs) are potential candidates for inventing highly correlated materials and devices and investigating the fundamental science of organic spintronics. This paper reports our experi- mental observations providing the dramatic impact of variation in ferromagnetic electrode length and thickness on paramagnetic molecule-based MTJMSD. Room temperature transport studies were performed to investigate the effect of FM electrode thickness. On the other hand, magnetic force microscopy measurements were con- ducted to understand the effect of FM electrode length extending beyond the molecular junction area, i.e., the site where paramagnetic molecules bridged between two FM. In the strong molecular coupling regime, transport study suggested thickness variation caused ~1000 to million-fold differences in junction conductivity. MFM study revealed near-zero magnetic contrast for pillar-shaped MTJMSD without any extended FM electrode. However, MFM images showed a multitude of microscopic magnetic phases on cross junction shaped MTJMSD where FM electrodes extended beyond the junction area. To understand the intriguing experimental results, we conducted an in-depth theoretical study using Monte Carlo Simulation (MCS) approach. MCS study utilized a Heisenberg atomic model of cross junction shaped MTJMSD to gain insights about room temperature transport and MFM experimental observations of microscopic MTJMSD. To make this study applicable for a wide variety of MTJMSDs, we systematically studied the effect of variation in molecular coupling strength between magnetic molecules and ferromagnetic (FM) electrodes of various dimensions. Here, our theoretical results suggest 25 atoms thick FM electrode show very weak molecular coupling impact and of the same order of experimentally studied FM electrode dimensions showing weak effect. Also, MCS results agree with experimentally observed multiple magnetic phases on lengthy FM electrodes of an MTJMSD.

42 ENGINEERING↗
DOE OSTI2024

Mitigation of distortion of Al/steel part under simulated paint baking condition: Experiment and numerical model studies

Multi-material joining of lightweight structures is essential to reduce vehicle weight for more energy savings and less greenhouse gas emission. However, mismatch of thermal expansion coefficient for dissimilar materials during the paint baking process can induce part distortion and joint failure for adhesive bonding. Here, in the present work, a thermomechanical model based on contact mechanics and large deformation theory was developed for dissimilar high-strength Al alloy and steel components to study the distortion mechanism and influential factors of the residual gap. The established model was used to optimize joint conditions, such as pitch distance and part geometry. When a weld pitch is shorter than 100 mm, the maximum gap between Al and steel part can be greatly reduced to 0.1 mm, and the local stress and plastic strain around the joint during the oven heating and cooling cycle are also substantially reduced compared with the long pitch case (900 mm). The numerical modeling results revealed that a comparable bending stiffness ratio between the steel and Al cross sections is critical to the minimization of gap and distortion under paint baking condition. Digital image correlation technique was used to measure the overall part distortion and local strain distribution that were used to validate the model prediction. Weld bonding (adhesive bonding with friction bit joining) process was successfully employed to join Al to steel component without gap opening in adhesive after the paint baking and cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Z-Target Radiography Postprocessing With A Deep Convolution Neural Network

Analyzing X-ray radiographs is crucial for understanding target behavior in Inertial Confinement Fusion (ICF) and High Energy Density (HED) platforms. However, the density of Magneto Raleigh Taylor (MRT) bands and limitations of target materials often obscure relevant spike growth and density information. To address this issue, machine learning postprocessing techniques can be applied to remove darkened regions in radiography images. In this study, a novel method is presented for removing MRT darkened regions from z-target radiographs using a convolutional neural network (CNN). The CNN, consisting of six layers, treats the darkened regions as noise and employs a mixed loss function and end-to-end frameworks to suppress them while preserving sharpness. The six-layer architecture is designed to effectively learn features when provided with a larger volume of learning space. Each layer is optimized using a mixed loss function that combines a standard loss pixel approach with a multi-scaled structural similarity index loss, which considers luminance, contrast, and structure in local neighborhoods. This approach is particularly beneficial for capturing the stochastic structure of MRT limbs. Due to the limited availability of experimental data, training is conducted using synthetic target radiography from 3D Alegra simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dynamic Focusing with High-Quality-Factor Metalenses

Metasurface lenses provide an ultrathin platform in which to focus light, but weak light–matter interactions limit their dynamic tunability. Here we design submicron-thick, ultrahigh quality factor (high-Q) metalenses that enable dynamic modulation of the focal length and intensity. Using full-field simulations, we show that quality factors exceeding 5000 can be generated by including subtle, periodic perturbations within the constituent Si nanoantennas. Such high-Q resonances enable lens modulation based on the nonlinear Kerr effect, with focal lengths varying from 4 to 6.5 μm and focal intensities decreasing by half as input intensity increases from 0.1 to 1 mW/μm 2 . We also show how multiple high-Q resonances can be embedded in the lens response through judicious placement of the perturbations. Our high-Q lens design, with quality factors 2 orders of magnitude higher than existing lens designs, provides a foundation for reconfigurable, multiplexed, and hyperspectral metasurface imaging platforms.

36 MATERIALS SCIENCE↗
NASA NTRS2008

New Insights into the Composition and Texture of Lunar Regolith Using Ultrafast Automated Electron-Beam Analysis

Sieved grain mounts of Apollo 16 drive tube samples have been examined using QEMSCAN - an innovative electron beam technology. By combining multiple energy-dispersive X-ray detectors, fully automated control, and off-line image processing, to produce digital mineral maps of particles exposed on polished surfaces, the result is an unprecedented quantity of mineralogical and petrographic data, on a particle-by-particle basis. Experimental analysis of four size fractions (500-250 microns, 150-90 microns, 75-45 microns and < 20 microns), prepared from two samples (64002,374 and 64002,262), has produced a robust and uniform dataset which allows for the quantification of mineralogy; texture; particle shape, size and density; and the digital classification of distinct particle types in each measured sample. These preliminary data show that there is a decrease in plagioclase modal content and an opposing increase in glass modal content, with decreasing particle size. These findings, together with data on trace phases (metals, sulphides, phosphates, and oxides), provide not only new insights into the make-up of lunar regolith at the Apollo 16 landing site, but also key physical parameters which can be used to design lunar simulants, and compute Figures of Merit for each material produced.

Rickman, Doug↗
DOE OSTI2021

Improved Gamma-Ray Point Source Quantification in Three Dimensions by Modeling Attenuation in the Scene

We discuss using a series of detector measurements taken at different locations to localize a source of radiation as a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2D and even 3D, but few have studied the effect of intervening material on the problem. In this work we examine gamma-ray data taken from a freely moving system and develop voxelized 3-D models of the scene using data from the onboard LiDAR. Ray casting is used to compute the distance each gamma ray travels through the scene material, which is then used to calculate attenuation assuming a single attenuation coefficient for solids within the geometry. Parameter estimation using maximum likelihood is performed to simultaneously find the attenuation coefficient, source activity, and source position that best match the data. Using a simulation, we validate the ability of this method to reconstruct the true location and activity of a source, along with the true attenuation coefficient of the structure it is inside, and then we apply the method to measured data with sources and find good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2017

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

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

Inverse theory↗