Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “discrete element modeling”

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 523 records · Page 29

Thermal mechanical analysis of sprag clutches

Work done at Case Western Reserve University on the Thermal Mechanical analysis of sprag helicopter clutches is reported. The report is presented in two parts. The first part is a description of a test rig for the measurement of the heat generated by high speed sprag clutch assemblies during cyclic torsional loading. The second part describes a finite element modeling procedure for sliding contact. The test rig provides a cyclic torsional load of 756 inch-pounds at 5000 rpm using a four-square arrangement. The sprag clutch test unit was placed between the high speed pinions of the circulating power loop. The test unit was designed to have replaceable inner ad outer races, which contain the instrumentation to monitor the sprag clutch. The torque loading device was chosen to be a water cooled magnetic clutch, which is controlled either manually or through a computer. In the second part, a Generalized Eulerian-Lagrangian formulation for non-linear dynamic problems is developed for solid materials. This formulation is derived from the basic laws and axioms of continuum mechanics. The novel aspect of this method is that we are able to investigate the physics in the spatial region of interest as material flows through it without having to follow material points. A finite element approximation to the governing equations is developed. Iterative Methods for the solution of the discrete finite element equations are explored. A FORTRAN program to implement this formulation is developed and a number of solutions to problems of sliding contact are presented.

Mullen, Robert L.↗

Simulating Matrix Crack and Delamination Interaction in a Clamped Tapered Beam

Blind predictions were conducted to validate a discrete crack methodology based on the Floating Node Method to simulate matrix-crack/delamination interaction. The main novel aspects of the approach are: (1) the implementation of the floating node method via an 'extended interface element' to represent delaminations, matrix-cracks and their interaction, (2) application of directional cohesive elements to infer overall delamination direction, and (3) use of delamination direction and stress state at the delamination front to determine migration onset. Overall, good agreement was obtained between simulations and experiments. However, the validation exercise revealed the strong dependence of the simulation of matrix-crack/delamination interaction on the strength data (in this case transverse interlaminar strength, YT) used within the cohesive zone approach applied in this work. This strength value, YT, is itself dependent on the test geometry from which the strength measurement is taken. Thus, choosing an appropriate strength value becomes an ad-hoc step. As a consequence, further work is needed to adequately characterize and assess the accuracy and adequacy of cohesive zone approaches to model small crack growth and crack onset. Additionally, often when simulating damage progression with cohesive zone elements, the strength is lowered while keeping the fracture toughness constant to enable the use of coarser meshes. Results from the present study suggest that this approach is not recommended for any problem involving crack initiation, small crack growth or multiple crack interaction.

De Carvalho, N. V.↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

To understand the behavior of charged lunar regolith when perturbed by lunar landers, it is important to couple the grain dynamics with mechanical and electrical particle interactions. To accomplish this, improvements have been made to a discrete element method (DEM) software package to include both short- and long-range interactions between spherical particles. Short-range interactions rely on contact between the particles, such as electrical conduction and triboelectric charge transfer driven by work functions. Long-range interactions act at a distance between every pairing of particles, such as electrostatic forces and gravitational forces. Results from simulations between a few particles are compared with theory to verify these added behaviors prior to scaling up to more complex scenarios. The radii, initial charges, electrical conductivities, work functions, and separation of the particles are varied and the resultant charges as well as the time required to reach the final state are determined.

Electrostatics↗

Comparisons of Kinematics and Dynamics Simulation Software Tools

Kinematic and dynamic analyses for moving bodies are essential to system engineers and designers in the process of design and validations. 3D visualization and motion simulation plus finite element analysis (FEA) give engineers a better way to present ideas and results. Marshall Space Flight Center (MSFC) system engineering researchers are currently using IGRIP from DELMIA Inc. as a kinematic simulation tool for discrete bodies motion simulations. Although IGRIP is an excellent tool for kinematic simulation with some dynamic analysis capabilities in robotic control, explorations of other alternatives with more powerful dynamic analysis and FEA capabilities are necessary. Kinematics analysis will only examine the displacement, velocity, and acceleration of the mechanism without considering effects from masses of components. With dynamic analysis and FEA, effects such as the forces or torques at the joint due to mass and inertia of components can be identified. With keen market competition, ALGOR Mechanical Event Simulation (MES), MSC visualNastran 4D, Unigraphics Motion+, and Pro/MECHANICA were chosen for explorations. In this study, comparisons between software tools were presented in terms of following categories: graphical user interface (GUI), import capability, tutorial availability, ease of use, kinematic simulation capability, dynamic simulation capability, FEA capability, graphical output, technical support, and cost. Propulsion Test Article (PTA) with Fastrac engine model exported from IGRIP and an office chair mechanism were used as examples for simulations.

Yeu-Sheng Paul Shiue↗

Using Piezoelectric Devices to Transmit Power through Walls

A method denoted wireless acoustic-electric feed-through (WAEF) has been conceived for transmitting power and/or data signals through walls or other solid objects made of a variety of elastic materials that could be electrically conductive or nonconductive. WAEF would make it unnecessary to use wires, optical fibers, tubes, or other discrete wall-penetrating signal-transmitting components, thereby eliminating the potential for structural weakening or leakage at such penetrations. Avoidance of such penetrations could be essential in some applications in which maintenance of pressure, vacuum, or chemical or biological isolation is required. In a basic WAEF setup, a transmitting piezoelectric transducer on one side of a wall would be driven at resonance to excite ultrasonic vibrations in the wall. A receiving piezoelectric transducer on the opposite side of the wall would convert the vibrations back to an ultrasonic AC electric signal, which would then be detected and otherwise processed in a manner that would depend on the modulation (if any) applied to the signal and whether the signal was used to transmit power, data, or both. An electromechanical-network model has been derived as a computationally efficient means of analyzing and designing a WAEF system. This model is a variant of a prior model, known in the piezoelectric-transducer art as Mason's equivalent-circuit model, in which the electrical and mechanical dynamics, including electromechanical couplings, are expressed as electrical circuit elements that can include inductors, capacitors, and lumped-parameter complex impedances. The real parts of the complex impedances are used to account for dielectric, mechanical, and coupling losses in all components (including all piezoelectric-transducer, wall, and intermediate material layers). In an application to a three-layer piezoelectric structure, this model was shown to yield the same results as do solutions of the wave equations of piezoelectricity and acoustic propagation in their full complexity.

Sherrit, Stewart↗

Operator procedure verification with a rapidly reconfigurable simulator

Generating and testing procedures for controlling spacecraft subsystems composed of electro-mechanical and computationally realized elements has become a very difficult task. Before a spacecraft can be flown, mission controllers must envision a great variety of situations the flight crew may encounter during a mission and carefully construct procedures for operating the spacecraft in each possible situation. If, despite extensive pre-compilation of control procedures, an unforeseen situation arises during a mission, the mission controller must generate a new procedure for the flight crew in a limited amount of time. In such situations, the mission controller cannot systematically consider and test alternative procedures against models of the system being controlled, because the available simulator is too large and complex to reconfigure, run, and analyze quickly. A rapidly reconfigurable simulation environment that can execute a control procedure and show its effects on system behavior would greatly facilitate generation and testing of control procedures both before and during a mission. The How Things Work project at Stanford University has developed a system called DME (Device Modeling Environment) for modeling and simulating the behavior of electromechanical devices. DME was designed to facilitate model formulation and behavior simulation of device behavior including both continuous and discrete phenomena. We are currently extending DME for use in testing operator procedures, and we have built a knowledge base for modeling the Reaction Control System (RCS) of the space shuttle as a testbed. We believe that DME can facilitate design of operator procedures by providing mission controllers with a simulation environment that meets all these requirements.

Iwasaki, Yumi↗

Implicit Formulation of Muscle Dynamics in OpenSim

Astronauts lose bone and muscle mass during spaceflight. Exercise countermeasure is the primary method for counteracting bone and muscle mass loss in space. New spacecraft exercise device concepts are currently being developed for the NASAs new crew exploration vehicle. The NASA Digital Astronaut Project (DAP) uses computational modeling to help determine if the new exercise devices will be effective as countermeasures. The NASA Digital Astronaut Project is developing the ability to utilize predictive simulation to provide insight into the change in kinematics and kinetics with a change in device and gravitational environment (1-g versus 0-g). For example, in space exercise the subject's body weight is applied in addition to the loads prescribed for musculoskeletal maintenance. How and where these loads are applied obviously directly impacts bone and tissue loads. Additionally, due to space vehicle structural requirements, exercise devices are often placed on vibration isolation systems. This changes the apparent impedance or stiffness of the device as seen by the user. Data collection under these conditions is often impractical and limited. Predictive modeling provides a means to have a virtual subject to test hypotheses. Predictive simulation provides a virtual subject for which we are able to perform studies such as sensitivity to device loading and vibration isolation without the need for laboratory kinematic or kinetic test data.Direct Collocation optimization provides an efficient means to perform task based optimization and predictive modeling. It is relatively straight forward to structure a physical exercise task in a Direct Collocation mathematical formulation: perform a motion such that you start at an initial pose, achieve a given amount of deflection i.e a squat, return to the initial pose, and minimize muscle activation cost. Direct Collocation is advantageous in that it does not require numerical integration to evaluate the objective function. Instead, the system dynamics are transformed to discrete time and the optimizer is constrained such that the solution is not considered to be a valid unless the dynamic equations are satisfied at all time points. The simulation and optimization are effectively done simultaneously. Due to the implicit integration, time steps can be more coarse than in a differential equation solver. In a gait scenario this means that that the model constraints and cost function are evaluated at 100 nodes in the gait cycle versus 10,000 integration steps in a variable-step forward dynamic simulation. Furthermore, no time is wasted on accurate simulations of movements that are far from the optimum. Constrained optimization algorithms require a Jacobian matrix that contains the partial derivatives of each of the dynamic constraints with respect to of each of the state and control variables at all time points. This is a large but sparse matrix. An implicit dynamics formulation requires computation of the dynamic residuals f as a function of the states x and their derivatives, and controls u:f(x, dxdt, u) 0If the dynamics of musculoskeletal system are formulated implicitly, the Jacobian elements are often available analytically, eliminating the need for numerical differentiation; this is obviously computationally advantageous. Additionally, implicit formulation of musculoskeletal dynamics do not suffer from singularities from low mass bodies, zero muscle activation, or other stiff system or

physical exercise↗

Incorporation of Damage and Failure into an Orthotropic Elasto-Plastic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

A material model which incorporates several key capabilities which have been identified by the aerospace community as lacking in the composite impact models currently available in LS-DYNA(Registered Trademark) is under development. In particular, the material model, which is being implemented as MAT 213 into a tailored version of LS-DYNA being jointly developed by the FAA and NASA, incorporates both plasticity and damage within the material model, utilizes experimentally based tabulated input to define the evolution of plasticity and damage as opposed to specifying discrete input parameters (such as modulus and strength), and is able to analyze the response of composites composed with a variety of fiber architectures. The plasticity portion of the orthotropic, three-dimensional, macroscopic composite constitutive model is based on an extension of the Tsai-Wu composite failure model into a generalized yield function with a non-associative flow rule. The capability to account for the rate and temperature dependent deformation response of composites has also been incorporated into the material model. For the damage model, a strain equivalent formulation is utilized to allow for the uncoupling of the deformation and damage analyses. In the damage model, a diagonal damage tensor is defined to account for the directionally dependent variation of damage. However, in composites it has been found that loading in one direction can lead to damage in multiple coordinate directions. To account for this phenomena, the terms in the damage matrix are semi-coupled such that the damage in a particular coordinate direction is a function of the stresses and plastic strains in all of the coordinate directions. The onset of material failure, and thus element deletion, is being developed to be a function of the stresses and plastic strains in the various coordinate directions. Systematic procedures are being developed to generate the required input parameters based on the results of experimental tests.

Finite Element Method↗

CFD Analysis of Spray Combustion and Radiation in OMV Thrust Chamber

The Variable Thrust Engine (VTE), developed by TRW, for the Orbit Maneuvering Vehicle (OMV) uses a hypergolic propellant combination of Monomethyl Hydrazine (MMH) and Nitrogen Tetroxide (NTO) as fuel and oxidizer, respectively. The propellants are pressure fed into the combustion chamber through a single pintle injection element. The performance of this engine is dependent on the pintle geometry and a number of complex physical phenomena and their mutual interactions. The most important among these are (1) atomization of the liquid jets into fine droplets; (2) the motion of these droplets in the gas field; (3) vaporization of the droplets (4) turbulent mixing of the fuel and oxidizer; and (5) hypergolic reaction between MMH and NTO. Each of the above phenomena by itself poses a considerable challenge to the technical community. In a reactive flow field of the kind occurring inside the VTE, the mutual interactions between these physical processes tend to further complicate the analysis. The objective of this work is to develop a comprehensive mathematical modeling methodology to analyze the flow field within the VTE. Using this model, the effect of flow parameters on various physical processes such as atomization, spray dynamics, combustion, and radiation is studied. This information can then be used to optimize design parameters and thus improve the performance of the engine. The REFLEQS CFD Code is used for solving the fluid dynamic equations. The spray dynamics is modeled using the Eulerian-Lagrangian approach. The discrete ordinate method with 12 ordinate directions is used to predict the radiative heat transfer in the OMV combustion chamber, nozzle, and the heat shield. The hypergolic reaction between MMH and NTO is predicted using an equilibrium chemistry model with 13 species. The results indicate that mixing and combustion is very sensitive to the droplet size. Smaller droplets evaporate faster than bigger droplets, leading to a well mixed zone in the combustion chamber. The radiative heat flux at combustion chamber and nozzle walls are an order of negligible less than the conductive heat flux. Simulations performed with the heat shield show that a negligible amount of fluid is entrained into the heat shield region. However, the heat shield is shown to be effective in protecting the OMV structure surrounding the engine from the radiated heat.

Giridharan, M. G.↗

Integrated Aeromechanics with Three-Dimensional Solid-Multibody Structures

A full three-dimensional finite element-multibody structural dynamic solver is coupled to a three-dimensional Reynolds-averaged Navier-Stokes solver for the prediction of integrated aeromechanical stresses and strains on a rotor blade in forward flight. The objective is to lay the foundations of all major pieces of an integrated three-dimensional rotor dynamic analysis - from model construction to aeromechanical solution to stress/strain calculation. The primary focus is on the aeromechanical solution. Two types of three-dimensional CFD/CSD interfaces are constructed for this purpose with an emphasis on resolving errors from geometry mis-match so that initial-stage approximate structural geometries can also be effectively analyzed. A three-dimensional structural model is constructed as an approximation to a UH-60A-like fully articulated rotor. The aerodynamic model is identical to the UH-60A rotor. For preliminary validation measurements from a UH-60A high speed flight is used where CFD coupling is essential to capture the advancing side tip transonic effects. The key conclusion is that an integrated aeromechanical analysis is indeed possible with three-dimensional structural dynamics but requires a careful description of its geometry and discretization of its parts.

Three-Dimensional↗

NASA Tech Briefs, November 2009

Topics covered include: Cryogenic Chamber for Servo-Hydraulic Materials Testing; Apparatus Measures Thermal Conductance Through a Thin Sample from Cryogenic to Room Temperature; Rover Attitude and Pointing System Simulation Testbed; Desktop Application Program to Simulate Cargo-Air-Drop Tests; Multimodal Friction Ignition Tester; Small-Bolt Torque-Tension Tester; Integrated Spacesuit Audio System Enhances Speech Quality and Reduces Noise; Hardware Implementation of a Bilateral Subtraction Filter; Simple Optoelectronic Feedback in Microwave Oscillators; Small X-Band Oscillator Antennas; Free-Space Optical Interconnect Employing VCSEL Diodes; Discrete Fourier Transform Analysis in a Complex Vector Space; Miniature Scroll Pumps Fabricated by LIGA; Self-Assembling, Flexible, Pre-Ceramic Composite Preforms; Flight-speed Integral Image Analysis Toolkit; Work Coordination Engine; Multi-Mission Automated Task Invocation Subsystem; Autonomously Calibrating a Quadrupole Mass Spectrometer; Determining Spacecraft Reaction Wheel Friction Parameters; Composite Silica Aerogels Opacified with Titania; Multiplexed Colorimetric Solid-Phase Extraction; Detecting Airborne Mercury by Use of Polymer/Carbon Films; Lattice-Matched Semiconductor Layers on Single Crystalline Sapphire Substrate; Pressure-Energized Seal Rings to Better Withstand Flows; Rollerjaw Rock Crusher; Microwave Sterilization and Depyrogenation System; Quantifying Therapeutic and Diagnostic Efficacy in 2D Microvascular Images; NiF2/NaF:CaF2/Ca Solid-State High-Temperature Battery Cells; Critical Coupling Between Optical Fibers and WGM Resonators; Microwave Temperature Profiler Mounted in a Standard Airborne Research Canister; Alternative Determination of Density of the Titan Atmosphere; Solar Rejection Filter for Large Telescopes; Automated CFD for Generation of Airfoil Performance Tables; Progressive Classification Using Support Vector Machines; Active Learning with Irrelevant Examples; A Data Matrix Method for Improving the Quantification of Element Percentages of SEM/EDX Analysis; Deployable Shroud for the International X-Ray Observatory; Improved Model of a Mercury Ring Damper; Optoelectronic pH Meter: Further Details; X-38 Advanced Sublimator; and Solar Simulator Represents the Mars Surface Solar Environment.

Source record↗

Transmission X-ray Diffraction (XRD) Patterns Relevant to the MSL Chemin Amorphous Component: Sulfates And Silicates

The Mars Science Laboratory (MSL) CheMin instrument on the Curiosity rover is a transmission X-ray diffractometer (Co-Kalpha radiation source and a approx.5deg to approx.52deg 2theta range) where the analyzed powder samples are constrained to have discrete particle diameters <150 microns by a sieve. To date, diffraction patterns have been obtained for one basaltic soil (Rocknest (RN)) and four drill fines of coherent rock (John Klein (JK), Cumberland (CB), Windjana (WJ), and Confidence Hills (CH)). The CheMin instrument has detected and quantified the abundance of both primary igneous (e.g., feldspar, olivine, and pyroxene) and secondary (e.g., Ca-sulfates, hematite, akaganeite, and Fe-saponite) minerals. The diffraction patterns of all CheMin samples are also characterized by a broad diffraction band centered near 30deg 2theta and by increasing diffraction intensity (scattering continuum) from approx.15deg to approx.5deg, the 2theta minimum. Both the broad band and the scattering continuum are attributed to the presence of an XRD amorphous component. Estimates of amorphous component abundance, based on the XRD data itself and on mass-balance calculations using APXS data crystalline component chemistry derived from XRD data, martian meteorites, and/or stoichiometry [e.g., 6-9], range from approx.20 wt.% to approx.50 wt.% of bulk sample. The APXSbased calculations show that the amorphous component is rich in volatile elements (esp. SO3) and is not simply primary basaltic glass, which was used as a surrogate to model the broad band in the RN CheMin pattern. For RN, the entire volatile inventory (except minor anhydrite) is assigned to the amorphous component because no volatile-bearing crystalline phases were reported within detection limits [2]. For JK and CB, Fesaponite, basanite, and akaganeite are volatile-bearing crystalline components. Here we report transmission XRD patterns for sulfate and silicate phases relevant to interpretation of MSL-CheMin XRD amorphous components.

Morris, R. V.↗

Pterodactyl: Integrated Control Design for Precision Targeting of Deployable Entry Vehicles

Deployable Entry Vehicles (DEVs) enable in-situ scientific exploration at destinations with atmospheres across the solar system. Because they stow in a compact form and deploy only when ready to enter the atmosphere, DEVs relax the volume constraint imposed by rigid aeroshells. This work seeks to do for a DEV what the Wright Brothers did to propel modern day aviation: develop the guidance and control (G&C) methods that will make maneuvering and precision landing of DEV a reality. The Pterodactyl project objective is to deliver an integrated G&C methodology for a DEV, based on a detailed analysis that utilizes a Multi-disciplinary, Design, Analysis and Optimization (MDAO) framework. The current state-of-the-art for blunt body entry, G&C is rooted in the precision landing of vehicles such as Mars Science Laboratory (MSL) and Apollo, which used a propulsive reaction control system (RCS) to steer. Recent research has taken a particular interest in non-propulsive control for DEVs, including direct force control (angle of attack modulation via control surfaces or mass movement) and drag modulation (discrete change in ballistic coefficient). Using the MDAO framework that includes a guidance and control model to explore multiple control concepts for a DEV will shed light on the best design approach for these vehicles. In Pterodactyl, we will complete this study for a novel DEV concept, and then we will fabricate a functional prototype to help validate the design. The project is expected to down-select to a final control architecture by the end of 2018, and complete fabrication of the prototype by the end of 2019.The DEV chosen for detailed study in this project is the Adaptable Deployable Entry and Placement Technology (ADEPT). ADEPT uses a revolutionary 3D-woven carbon fabric that is foldable, can serve as primary structure, and can survive the extreme heating environment of atmospheric entry. The specific configuration of ADEPT under investigation is called Lifting Nano-ADEPT (LNA). LNA is designed for secondary payloads missions that require precision landing either for scientific objectives at a target destination or for payload recovery at Earth.The MDAO framework being created through this research, called COBRA-Pt (Composite Beam Roll-Up Solar Array-Prototype), will combine three critical elements of the system design: a guidance algorithm with Monte Carlo, a parametric control model, and vehicle geometry details. Novel control models being studied are deployable aerodynamic surfaces as well as shape morphing. These concepts will be compared at the system level with a more traditional propulsive RCS by comparing several key performance parameters. Upon completion of the design study, a functional prototype of LNA will be fabricated that will include the integration of guidance software and relevant control actuators. We expect this study will provide critical data that could feed into the development of an Earth-based flight test of LNA. The COBRA-Pt framework will provide a modular system by which to study any DEV concept in any atmosphere.

Smith, Brandon↗

Implementation of Combined Feather and Surface-Normal Ice Growth Models in LEWICE/X

Experimental observations have shown that discrete rime ice growths called feathers, which grow in approximately the direction of water droplet impingement, play an important role in the growth of ice on accreting surfaces for some thermodynamic conditions. An improved physical model of ice accretion has been implemented in the LEWICE 2D panel-based ice accretion code maintained by the NASA Lewis Research Center. The LEWICE/X model of ice accretion explicitly simulates regions of feather growth within the framework of the LEWICE model. Water droplets impinging on an accreting surface are withheld from the normal LEWICE mass/energy balance and handled in a separate routine; ice growth resulting from these droplets is performed with enhanced convective heat transfer approximately along droplet impingement directions. An independent underlying ice shape is grown along surface normals using the unmodified LEWICE method. The resulting dual-surface ice shape models roughness-induced feather growth observed in icing wind tunnel tests. Experiments indicate that the exact direction of feather growth is dependent on external conditions. Data is presented to support a linear variation of growth direction with temperature and cloud water content. Test runs of LEWICE/X indicate that the sizes of surface regions containing feathers are influenced by initial roughness element height. This suggests that a previous argument that feather region size is determined by boundary layer transition may be incorrect. Simulation results for two typical test cases give improved shape agreement over unmodified LEWICE.

Matthew T Velazquez↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA PACE Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Pulse Response of the Short-Wave Infrared Detection System of the Ocean Color Instrument for the NASA Pace Mission

The Ocean Color Instrument (OCI) on NASA’s Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission is a hyperspectral Earth imager with a spatial resolution of 1 km x 1 km and a spectral resolution of 5 nm in 2.5 nm steps over 342-887 nm. In addition, OCI provides 7 discrete bands in the 940-2260 nm Short-Wave InfraRed (SWIR) range. The front-end optical imager is a rotating mirror-based system that images the ground scene onto a slit with an instantaneous field of view of 16 km x 1 km. For the SWIR bands, the slit-image is re-imaged onto a 16x1 micro-lens array that effectively acts as the focal plane since each lens element is fiber coupled to wavelength filtered InGaAs and HgCdTe Photo Diodes (PDs). The pulse response of the detection system is critical to OCI SWIR performance. We find that PDs introduce an inherent slow tail in the pulse response due to slow diffusion moving carriers in their n and p regions. We show that this introduces response errors ranging from 1 down to 0.01 % for up to tens of science pixels after the pulse depending on the PD design and materials. It is shown that the response is distinctly different for the InGaAs and HgCdTe PDs. We explain how the front-end design can further increase this error. Finally, we detail the cause of the slow pulse response tail, how to model it, its impact on OCI performance and how it is characterized and corrected to meet OCI requirements.

ocean color↗

Wildfire Segmentation From Remotely Sensed Data Using Quantum-Compatible Conditional Vector Quantized-Variational Autoencoders

Wildfires represent a critical environmental hazard with multifaceted implications for ecosystems, communities, and public health [1]. The escalating frequency and intensity of wildfires globally have intensified the urgency for robust segmentation methodologies to facilitate effective mitigation, response, and recovery strategies [2]. Accurate wildfire segmentation is pivotal for delineating fire boundaries, assessing progression patterns, and prioritizing resource allocation during emergency scenarios. Furthermore, precise segmentation enables stakeholders, including policymakers, environmental scientists, and emergency responders, to formulate evidence-based strategies, thereby minimizing socio-economic disruptions and ecological degradation. Consequently, advancing wildfire segmentation techniques through innovative technological interventions remains a paramount research imperative. Although foundational in wildfire segmentation, traditional deterministic models exhibit inherent limitations that compromise their efficacy in dynamic and uncertain environments. These models often operate on rigid algorithms prioritizing deterministic classifications, thereby overlooking the inherent complexities and uncertainties associated with wildfire behavior and satellite data variability. Such deterministic frameworks tend to produce oversimplified representations that fail to capture the intricate nuances of evolving fire dynamics, spatial heterogeneity, and environmental interactions [1]. Consequently, the deterministic approach’s propensity for uncertainty collapsing [1, 3] hampers the accuracy, reliability, and applicability of segmentation outcomes in real-world scenarios. Contrastingly, stochastic models offer a more nuanced and adaptable framework for wildfire segmentation. By integrating probabilistic elements into the modeling paradigm, stochastic approaches, particularly probabilistic approaches such as variational auto encoders (VAEs) [4], facilitate comprehensive uncertainty assessment, enabling researchers to quantify and incorporate uncertainties into segmentation outcomes effectively. This probabilistic nature empowers stochastic models to encapsulate variability, account for data inconsistencies, and adapt to evolving environmental conditions, enhancing segmentation accuracy, reliability, and robustness. Embracing stochastic methodologies thus catalyzes advancements in wildfire science by fostering a more holistic, adaptive, and resilient segmentation framework. Despite VAEs demonstrating significant promise in various applications, they come with inherent limitations that have garnered attention within the machine learning community. One of the primary drawbacks lies in their reliance on static priors, which essentially assume a fixed distribution for latent variables, thereby limiting the model’s flexibility to capture complex data structures effectively [5]. This static nature leads to suboptimal representations, especially when dealing with complex and high-dimensional data. Additionally, VAEs often struggle with generating sharp and realistic samples, a phenomenon commonly referred to as mode collapse [5, 7, 6]. Furthermore, the optimization process in VAEs, which involves balancing the reconstruction loss and the regularization term, can sometimes be challenging to fine-tune [7]. In recent efforts to address these shortcomings, alternative approaches like Vector Quantized Variational Auto encoders(VQ-VAEs) [7], address the challenges by incorporating discrete latent variables and leveraging techniques that enhance the quality and diversity of generated samples while maintaining efficient training dynamics. VQ-VAEs propose a dynamic prior distribution generation mechanism that diverges from the static priors commonly associated with traditional VAEs. This dynamic approach allows for more adaptive and context-aware latent variable representations, thereby potentially capturing complex data structures more effectively. Unlike autoregressive prior models such as PixelCNN, which, despite their ability to model dependencies across data dimensions, suffer from significant computational inefficiencies and lack flexibility in handling diverse datasets. In our work, we propose to use a generative quantum-compatible approach to help alleviate the shortcomings of autoregressive prior model in VQ-VAEs. Restricted Boltzmann Machines (RBMs) are a viable alternative prior model that can learn prior distributions in a faster and more flexible manner. In this research endeavor, we meticulously curate a state-of-the-art dataset leveraging satellite MODIS data in conjunction with VIIRS fire masks, derived from Fire Radiative Power (FRP), thereby encapsulating diverse wildfire scenarios and environmental contexts. We developed a conditional VQ-VAE architecture with the RBM prior model that is trained in a supervised manner for segmenting wildfire masks. This innovative approach synergistically harnesses deep learning capabilities, enabling the generation of segmentation maps characterized by heightened precision, granularity, and contextual relevance. Furthermore, replacing the autoregressive prior learning method proposed by the original VQ-VAE with a prior density approximation via quantum-compatible RBM facilitates expedited inference processes, augments flexibility in prior sampling, optimizes computational efficiency and establishes a groundbreaking benchmark in wildfire segmentation methodologies.

quantum machine learning↗

Path Planning Algorithms for the Adaptive Sensor Fleet

The Adaptive Sensor Fleet (ASF) is a general purpose fleet management and planning system being developed by NASA in coordination with NOAA. The current mission of ASF is to provide the capability for autonomous cooperative survey and sampling of dynamic oceanographic phenomena such as current systems and algae blooms. Each ASF vessel is a software model that represents a real world platform that carries a variety of sensors. The OASIS platform will provide the first physical vessel, outfitted with the systems and payloads necessary to execute the oceanographic observations described in this paper. The ASF architecture is being designed for extensibility to accommodate heterogenous fleet elements, and is not limited to using the OASIS platform to acquire data. This paper describes the path planning algorithms developed for the acquisition phase of a typical ASF task. Given a polygonal target region to be surveyed, the region is subdivided according to the number of vessels in the fleet. The subdivision algorithm seeks a solution in which all subregions have equal area and minimum mean radius. Once the subregions are defined, a dynamic programming method is used to find a minimum-time path for each vessel from its initial position to its assigned region. This path plan includes the effects of water currents as well as avoidance of known obstacles. A fleet-level planning algorithm then shuffles the individual vessel assignments to find the overall solution which puts all vessels in their assigned regions in the minimum time. This shuffle algorithm may be described as a process of elimination on the sorted list of permutations of a cost matrix. All these path planning algorithms are facilitated by discretizing the region of interest onto a hexagonal tiling.

Stoneking, Eric↗