Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fourier imaging”

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 649 records · Page 36

FLIT: A Generic Fortran Library based on Interfaces and Templates

This Fortran code consists of multiple modules with a focus on simplifying array operations, image processing, and numerical computation especially for computational geophysics applications. We intend to use this code to demonstrate the application and usefulness of Fortran interface and templates for generic programming, especially for computational geophysics and seismology applications. The code has several notable features. Firstly, it is based on a modularized structure, where each module contains multiple functions but with a focus of functionality. Secondly, it heavily uses interfaces and templates for improving the genericness and convenience of the resulting code, where a same function interface can enclose a group of functions that perform the same functionality but with inputs/output variables of different data types. Thirdly, it includes a variety of generic functions with an emphasis on array operations, such as rotation, flipping, cropping, padding, fast Fourier transform, Gaussian blurring, interpolation, and so on. We name this package FLIP – a generic Fortran Library based on Interfaces and Templates.

Gao, Kai↗

Water chemistry in flume channel and hyporheic zone (i.e., porewater) associated with: “Rethinking Aerobic Respiration in the Hyporheic Zone Under Variation in Carbon and Nitrogen Stoichiometry”

Dissolved oxygen (DO), total organic carbon (TOC), total nitrogen (TN), molecular data for organic matter, and biochemical reactions for surface water and porewater (i.e., hyporheic zone) collected from a water recirculating flume located at the University of Texas, Austin. The flume contained real river water from Lower Colorado River(Austin, TX) and clean sand. Hyporheic exchange in the flume was induced through The study aims to understand relationships between aerobic metabolism of organic matter and molecular characteristics of organic matter, such as thermodynamic signature and nitrogen content, through the extent of the hyporheic zone at 10 cm- resolution, and through time. During the experiment, organic matter (dry leaves) was added to the flume and removed after 24 hours. The water samples were collected before the addition of leaves, at the time of removal of leaves, and at hour 72. The water samples were analyzed using ultrahigh resolution Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) and total organic carbon (TOC) and total nitrogen (TN) analysis. Dissolved oxygen content throughout the surface water and the hyporheic zone of the flume was measured with a large planar optode. This data package is associated with the publication ’ Rethinking Aerobic Respiration in the Hyporheic Zone Under Variation in Carbon and Nitrogen Stoichiometry’ published in Environmental Science and Technology (Turețcaia et al., 2023 https://doi.org/10.1021/acs.est.3c04765). The dataset is comprised of five folders (1) Diss_O2_pic, (2) input_files (3) output_files; (4) python_code; and (5) R_code . Diss_O2_pic contains siximages of dissolved oxygen distribution in a bedform at hours 0, 24, and 72 of the experiment conducted in a large recirculation flume. Images are in separate R and G channels (i.e., RGB). The input_files contains (1) a csv file with FTICR peaks identified within each sample, (2) a csv file with molecular information pertinent to FTICR data with Gibbs free energy calculations adjusted for environmental temperature, (3) a csv file containing concentrations of non-purgeable organic carbon measured throughout the experiment , (4) a csv file containing concentrations of total nitrogen measured throughout the experiment, (5) a csv file containing total biochemical reactions (i.e., transformations) identified in the dataset, (6) a csv containing transformation profiles, and (7) a csv file containing transformations with formulas, and (8) a jpg file with schematic representation of locations for sample collection. The output_files contains (1) and xlsx file containing percent biochemical reactions containing nitrogen identified across all 39 sample, (2) a csv file of merged FTICR data and molecular information files, (3) a csv files containing average Gibbs free energy within sampling domains and at each sampling location, (4) a csv file with average concentrations of dissolved oxygen across sampling locations at hour 0, (5) a csv file with average concentrations of dissolved oxygen across sampling locations at hour 24, (6) a csv file with average concentrations of dissolved oxygen across sampling locations at hour 72, (7) a csv file with percent chemical classes identified across sampling locations at hour 0, (8) a csv file with percent chemical classes identified across sampling locations at hour 24, (9) a csv file with percent chemical classes identified across sampling locations at hour 72, and (10) a csv file containing percent nitrogen containing biochemical reactions identified across sampling locations at hours 0, 24, and 72. The python_code contains seven ipynb files which are Jupyter Notebooks used for data analysis and figures generation. The R_code contains 3 R files with R code used for data analysis and figures generation. This data package contains the processed data used in the associated manuscript. This data has not been previously published.

54 ENVIRONMENTAL SCIENCES↗

NASA Tech Briefs, February 2011

Topics covered include: Multi-Segment Radius Measurement Using an Absolute Distance Meter Through a Null Assembly; Fiber-Optic Magnetic-Field-Strength Measurement System for Lightning Detection; Photocatalytic Active Radiation Measurements and Use; Computer Generated Hologram System for Wavefront Measurement System Calibration; Non-Contact Thermal Properties Measurement with Low-Power Laser and IR Camera System; SpaceCube 2.0: An Advanced Hybrid Onboard Data Processor; CMOS Imager Has Better Cross-Talk and Full-Well Performance; High-Performance Wireless Telemetry; Telemetry-Based Ranging; JWST Wavefront Control Toolbox; Java Image I/O for VICAR, PDS, and ISIS; X-Band Acquisition Aid Software; Antimicrobial-Coated Granules for Disinfecting Water; Range 7 Scanner Integration with PaR Robot Scanning System; Methods of Antimicrobial Coating of Diverse Materials; High-Operating-Temperature Barrier Infrared Detector with Tailorable Cutoff Wavelength; A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane; Thermally Conductive Tape Based on Carbon Nanotube Arrays; Two Catalysts for Selective Oxidation of Contaminant Gases; Nanoscale Metal Oxide Semiconductors for Gas Sensing; Lightweight, Ultra-High-Temperature, CMC-Lined Carbon/Carbon Structures; Sample Acquisition and Handling System from a Remote Platform; Improved Rare-Earth Emitter Hollow Cathode; High-Temperature Smart Structures for Engine Noise Reduction and Performance Enhancement; Cryogenic Scan Mechanism for Fourier Transform Spectrometer; Piezoelectric Rotary Tube Motor; Thermoelectric Energy Conversion Technology for High-Altitude Airships; Combustor Computations for CO2-Neutral Aviation; Use of Dynamic Distortion to Predict and Alleviate Loss of Control; Cycle Time Reduction in Trapped Mercury Ion Atomic Frequency Standards; and A (201)Hg+ Comagnetometer for (199)Hg+ Trapped Ion Space Atomic Clocks.

Source record↗

Applications of a Forward-Looking Interferometer for the On-board Detection of Aviation Weather Hazards

The Forward-Looking Interferometer (FLI) is a new instrument concept for obtaining measurements of potential weather hazards to alert flight crews. The FLI concept is based on high-resolution Infrared (IR) Fourier Transform Spectrometry (FTS) technologies that have been developed for satellite remote sensing, and which have also been applied to the detection of aerosols and gases for other purposes. It is being evaluated for multiple hazards including clear air turbulence (CAT), volcanic ash, wake vortices, low slant range visibility, dry wind shear, and icing, during all phases of flight. Previous sensitivity and characterization studies addressed the phenomenology that supports detection and mitigation by the FLI. Techniques for determining the range, and hence warning time, were demonstrated for several of the hazards, and a table of research instrument parameters was developed for investigating all of the hazards discussed above. This work supports the feasibility of detecting multiple hazards with an FLI multi-hazard airborne sensor, and for producing enhanced IR images in reduced visibility conditions; however, further research must be performed to develop a means to estimate the intensities of the hazards posed to an aircraft and to develop robust algorithms to relate sensor measurables to hazard levels. In addition, validation tests need to be performed with a prototype system.

West, Leanne↗

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗

Lightning-channel morphology by return-stroke radiation field waveforms

Simultaneous video and wideband electric field recordings of 32 cloud-to-ground lightning flashes in Florida were analyzed to show the formation of new channels to ground can be detected by examination of the return-stroke radiation fields alone. The return-stroke E and dE/dt waveforms were subjectively classified according to their fine structure. Then the video images were examined field by field to identify each waveform with a visible channel to ground. Fifty-five correlated waveforms and channel images were obtained. Of these, all 34 first-stroke waveforms (multiple jagged E peaks, noisy dE/dt), 8 of which were not radiated by the chronologically first stroke in the flash, came from new channels to ground (not previously seen on video). All 18 subsequent-stroke waveforms (smoothly rounded E and quiet dE/dt after initial peak) were radiated by old channels (illuminated by a previous stroke). Two double-ground waveforms (two distinct first-return-stroke pulses separated by tens of microseconds or less) coincided with video fields showing two new channels. One `anomalous-stroke' waveform (beginning like a first stroke and ending like a subsequent) was produced by a new channel segment to ground branching off an old channel. This waveform classification depends on the presence or absence of high-frequency fine structure. Fourier analysis shows that first-stroke waveforms contain about 18 dB more spectral power in the frequency interval from 500 kHz to at least 7 MHz than subsequent-stroke waveforms for at least 13 microseconds after the main peak.

Willett, J. C.↗

Design of coupled mace filters for optical pattern recognition using practical spatial light modulators

Spatial light modulators (SLMs) are being used in correlation-based optical pattern recognition systems to implement the Fourier domain filters. Currently available SLMs have certain limitations with respect to the realizability of these filters. Therefore, it is necessary to incorporate the SLM constraints in the design of the filters. The design of a SLM-constrained minimum average correlation energy (SLM-MACE) filter using the simulated annealing-based optimization technique was investigated. The SLM-MACE filter was synthesized for three different types of constraints. The performance of the filter was evaluated in terms of its recognition (discrimination) capabilities using computer simulations. The correlation plane characteristics of the SLM-MACE filter were found to be reasonably good. The SLM-MACE filter yielded far better results than the analytical MACE filter implemented on practical SLMs using the constrained magnitude technique. Further, the filter performance was evaluated in the presence of noise in the input test images. This work demonstrated the need to include the SLM constraints in the filter design. Finally, a method is suggested to reduce the computation time required for the synthesis of the SLM-MACE filter.

Rajan, P. K.↗

Relationship between Horizontal Wind Velocity and Normalized Surface Cross Section Using Data from the HIWRAP Dual-Frequency Airborne Radar

The High-Altitude Imaging Wind and Rain Airborne Profiler (HIWRAP) dual-frequency conically scanning airborne radar provides estimates of the range-profiled mean Doppler and backscattered power from the precipitation and surface. A velocity–azimuth display analysis yields near-surface estimates of the mean horizontal wind vector υh in cases in which precipitation is present throughout the scan. From the surface return, the normalized radar cross section (NRCS) is obtained, which, by a method previously described, can be corrected for path attenuation. Comparisons between υh and the attenuation-corrected NRCS are used to derive transfer functions that provide estimates of the wind vector from the NRCS data under both rain and rain-free conditions. A reasonably robust transfer function is found by using the mean NRCS (⟨NRCS⟩) over the scan along with a filtering of the data based on a Fourier series analysis of υh and the NRCS. The approach gives good correlation coefficients between υh and ⟨NRCS⟩ at Ku band at incidence angles of 30° and 40°. The correlation degrades if the Ka-band data are used rather than the Ku band.

R Meneghini↗

Hydrogen leak detection in the Space Shuttle

This study focuses on a helium gas jet flowing into room air. Measurements of helium concentration and velocity in the jet-air mixture are reported. The objective is to learn about jet characteristics so that dynamically similar hydrogen leaks may be located in the Space Shuttle. The hazardous gas detection system (HGDS) in the mobile launch pad uses mass spectrometers to monitor the shuttle environment for leaks. The mass spectrometers are fed by long sample tubes which draw gas from the payload bay, mid body, aft engine compartment and external tank. The overall purpose of this study is to improve the HGDS especially in its potential for locating hydrogen leaks. A rapid-response leak detection experiment was designed, built, and tested, following on the work done in this program last summer. The apparatus included a Perkin Elmer MGA-1200 mass spectrometer and air velocity transducer, both monitored by a Macintosh IIFX computer using LabVIEW software. A jet of helium flowing into the lab air simulated a gas leak. Steady helium or hydrogen-nitrogen jets were logged for concentration and velocity, and the power spectral density of each was computed. Last year, large eddies and vortices were visually seen with Schlieren imaging, and they were detected in the time plots of the various instruments. The response time of the MGA-1200 was found in the range of 0.05 to 0.1 sec. Pulsed concentration waves were clearly detected at 25 cycles per sec by spectral analysis of MGA data. No peaks were detected in the power spectrum, so in the present study, 10 Hz bandwidth-averaged power levels were examined at regular frequency intervals. The practical consequences of last year's study are as follows: sampling frequency should be increased above the present rate of 1 sample per second so that transients could be observed and analyzed with frequency response methods. Many more experiments and conditions were observed in this second summer, including the effects of orifice diameter, jet velocity, sample tube design, radial effects, vertical flow, and low hydrogen concentrations (1 percent). A frequent observation was that the power spectrum, calculated from the Fourier transform of concentration fluctuations, gives a separate piece of information from concentration. Many of the tests suggest that power is high where mixing occurs at the helium-air interface. This fact is apparently independent of the concentration level, which could be high or low, but depends on the sample location relative to the jet (leak) origin, whereas high concentration may be due to a strong leak far away or a small leak close to the sample tube. If the power is low for any concentration level, this would signify helium is arriving at the sample tube by diffusion, not chaotic mixing caused by the jet interaction with air. The practical result is to propose a modification of the HGDL mass spectrometer data sampling and software so that sampling rates could be capable of observing at least 25 Hz fluctuations.

Barile, Ronald G↗

A dual high-energy radiography platform with 15 μ m resolution at the National Ignition Facility

To study matter at extreme densities and pressures, we need mega laser facilities such as the National Ignition Facility as well as creative methods to make observations during timescales of a billionth of a second. To facilitate this, we developed a platform and diagnostic to characterize a new point-projection radiography configuration using two micro-wires irradiated by a short pulse laser system that provides a large field of view with up to 3.6 ns separation between images. In this work, we used tungsten-carbide solid spheres as reference objects and inferred characteristics of the back-lighter source using a forward-fitting algorithm. The resolution of the system is inferred to be 15 μm (using 12.5 μm diameter wires). The bremsstrahlung temperature of the source is 70–300 keV, depending on laser energy and coupling efficiency. By adding the images recorded on multiple stacked image plates, the signal-to-noise of the system is nearly doubled. The imaging characterization technique described here can be adapted to most point-projection platforms where the resolution, spectral contrast, and signal-to-noise are important.

47 OTHER INSTRUMENTATION↗

The effects of free surfaces on deformation twinning in HCP metals

Deformation twinning is a predominant mode of plastic deformation for hexagonal close packed metals, like Mg and Ti. The heterogenous microstructure and the local stresses associated with twinning play a key role in their mechanical response and fracture. Surface analyses, like electron microscopy, are frequently employed to spatially map microstructure and micromechanical fields in order to study twinning behavior. However, these measurements are inherently influenced by the vicinity of the free surface. In this work, an elasto-visco-plastic fast-Fourier-transform (EVP-FFT) polycrystal modeling approach is employed to investigate the effects of free surfaces on twin development before and after loading. We compare calculated micromechanical fields on free surfaces with those calculated inside the bulk and, in some cases, experimental surface measurements. The results indicate that the creation of free surfaces can promote twin propagation and growth and can influence twin morphology by causing a twin lamella to become larger, more blunted and irregular. The structure along the twin boundaries are also affected, due to the higher driving stresses that extend prismatic-basal and basal-prismatic facets. Furthermore, free surfaces invoke different slip activities in the twin and the surrounding parent crystal by enhancing basal, prismatic and pyramidal slip in some localized regions, while reducing slip in others. We demonstrate that the simulated free-surface effects lead to better qualitative and quantitative agreement with experimental measurements from scanning electron microscopy and digital image correlation.

36 MATERIALS SCIENCE↗

Response Surface Methodology As a New Approach for Finding Optimal MALDI Matrix Spraying Parameters for Mass Spectrometry Imaging

Automated spraying devices have become ubiquitous in laboratories employing matrixassisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in part because they permit control of a number of matrix application parameters that can easily be reproduced for intra- and interlaboratory studies. Determining the optimal parameters for MALDI matrix application, such as temperature, flow rate, spraying velocity, number of spraying cycles, and solvent composition for matrix application, is critical for obtaining high-quality MALDI-MSI data. However, there are no established approaches for optimizing these multiple parameters simultaneously. Instead optimization is performed iteratively (i.e., one parameter at a time), which is time consuming and can lead to overall non-optimal settings. In this report, we demonstrate the use a novel experimental design and the response surface methodology to optimize five parameters of MALDI matrix application using a robotic sprayer. Thirtytwo combinations of MALDI matrix spraying conditions were tested, which allowed us to elucidate relationships between each of the application parameters as determined by MALDI-MS (specifically, using a 15 Tesla Fourier transform ion cyclotron resonance mass spectrometer). As such, we were able to determine the optimal automated spraying parameters that minimized signal de-localization and enabled high MALDI sensitivity. We envision this optimization strategy can be utilized for other MALDI-MSI applications, tissue types, and matrix application approaches.

60 APPLIED LIFE SCIENCES↗

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗

NASA Tech Briefs, April 2011

Topics covered include: Amperometric Solid Electrolyte Oxygen Microsensors with Easy Batch Fabrication; Two-Axis Direct Fluid Shear Stress Sensor for Aerodynamic Applications; Target Assembly to Check Boresight Alignment of Active Sensors; Virtual Sensor Test Instrumentation; Evaluation of the Reflection Coefficient of Microstrip Elements for Reflectarray Antennas; Miniaturized Ka-Band Dual-Channel Radar; Continuous-Integration Laser Energy Lidar Monitor; Miniaturized Airborne Imaging Central Server System; Radiation-Tolerant, SpaceWire-Compatible Switching Fabric; Small Microprocessor for ASIC or FPGA Implementation; Source-Coupled, N-Channel, JFET-Based Digital Logic Gate Structure Using Resistive Level Shifters; High-Voltage-Input Level Translator Using Standard CMOS; Monitoring Digital Closed-Loop Feedback Systems; MASCOT - MATLAB Stability and Control Toolbox; MIRO Continuum Calibration for Asteroid Mode; GOATS Image Projection Component; Coded Modulation in C and MATLAB; Low-Dead-Volume Inlet for Vacuum Chamber; Thermal Control Method for High-Current Wire Bundles by Injecting a Thermally Conductive Filler; Method for Selective Cleaning of Mold Release from Composite Honeycomb Surfaces; Infrared-Bolometer Arrays with Reflective Backshorts; Commercialization of LARC (trade mark) -SI Polyimide Technology; Novel Low-Density Ablators Containing Hyperbranched Poly(azomethine)s; Carbon Nanotubes on Titanium Substrates for Stray Light Suppression; Monolithic, High-Speed Fiber-Optic Switching Array for Lidar; Grid-Tied Photovoltaic Power System; Spectroelectrochemical Instrument Measures TOC; A Miniaturized Video System for Monitoring Drosophila Behavior; Hydrofocusing Bioreactor Produces Anti-Cancer Alkaloids; Creep Measurement Video Extensometer; Radius of Curvature Measurement of Large Optics Using Interferometry and Laser Tracker n-B-pi-p Superlattice Infrared Detector; Safe Onboard Guidance and Control Under Probabilistic Uncertainty; General Tool for Evaluating High-Contrast Coronagraphic Telescope Performance Error Budgets; Hidden Statistics of Schroedinger Equation; Optimal Padding for the Two-Dimensional Fast Fourier Transform; Spatial Query for Planetary Data; Higher Order Mode Coupling in Feed Waveguide of a Planar Slot Array Antenna; Evolutionary Computational Methods for Identifying Emergent Behavior in Autonomous Systems; Sampling Theorem in Terms of the Bandwidth and Sampling Interval; Meteoroid/Orbital Debris Shield Engineering Development Practice and Procedure; Self-Balancing, Optical-Center-Pivot, Fast-Steering Mirror; Wireless Orbiter Hang-Angle Inclinometer System; and Internal Electrostatic Discharge Monitor - IESDM.

Source record↗

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

High-resolution cryo-EM reconstructions in the presence of substantial aberrations

Here, an analysis is performed of how uncorrected antisymmetric aberrations, such as coma and trefoil, affect cryo-EM single-particle reconstruction (SPR) results, and an analytical formula quantifying information loss owing to their presence is inferred that explains why Fourier-shell coefficient-based statistics may report significantly overestimated resolution if these aberrations are not fully corrected. The analysis is validated with reference-based aberration refinement for two cryo-EM SPR data sets acquired with a 200 kV microscope in the presence of coma exceeding 40 µm, and 2.3 and 2.7 Å reconstructions for 144 and 173 kDa particles, respectively, were obtained. The results provide a description of an efficient approach for assessing information loss in cryo-EM SPR data acquired in the presence of higher order aberrations, and address inconsistent guidelines regarding the level of aberrations that is acceptable in cryo-EM SPR experiments.

3D reconstruction and image processing↗

Accurate real space iterative reconstruction (RESIRE) algorithm for tomography

Tomography has made a revolutionary impact on the physical, biological and medical sciences. The mathematical foundation of tomography is to reconstruct a three-dimensional (3D) object from a set of two-dimensional (2D) projections. As the number of projections that can be measured from a sample is usually limited by the tolerable radiation dose and/or the geometric constraint on the tilt range, a main challenge in tomography is to achieve the best possible 3D reconstruction from a limited number of projections with noise. Over the years, a number of tomographic reconstruction methods have been developed including direct inversion, real-space, and Fourier-based iterative algorithms. Here, we report the development of a real-space iterative reconstruction (RESIRE) algorithm for accurate tomographic reconstruction. RESIRE iterates between the update of a reconstructed 3D object and the measured projections using a forward and back projection step. The forward projection step is implemented by the Fourier slice theorem or the Radon transform, and the back projection step by a linear transformation. Our numerical and experimental results demonstrate that RESIRE performs more accurate 3D reconstructions than other existing tomographic algorithms, when there are a limited number of projections with noise. Furthermore, RESIRE can be used to reconstruct the 3D structure of extended objects as demonstrated by the determination of the 3D atomic structure of an amorphous Ta thin film. We expect that RESIRE can be widely employed in the tomography applications in different fields. Finally, to make the method accessible to the general user community, the MATLAB source code of RESIRE and all the simulated and experimental data are available at https://zenodo.org/record/7273314.

97 MATHEMATICS AND COMPUTING↗