Search NASA⌕ Search

SEARCH · Search NASA

Results for “post processing”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Biochar–polymer composites for 3D printing: a review

Biochar, a bio-based co-product of biofuel production via thermochemical conversion, holds potential as a filler for polymer composites to reduce costs, improve thermomechanical properties, and aid in environmental remediation. 3D-printed biochar composites have received growing interest over the past few years but have experienced difficulties such as poor layer adhesion and nozzle clogging. Currently, no literature review examines 3D-printed biochar composites and related biochar properties in-depth. This work summarizes and discusses recent studies on 3D-printed polymer and biochar composites and examines their mechanical, thermal, and additional properties that result from each study. Technical challenges in printability, such as nozzle clogging from particle size and biochar aggregation, are also discussed. Furthermore, this work discusses the variability of biochar properties resulting from the pyrolysis conditions and feedstock choice in relation to potential 3D printing outcomes. In particular, several studies reported that high lignin feedstocks could be candidates for 3D printing. The post-processing approaches of the biochar via physical and chemical methods are also introduced. Ball milling appears to hold the most promise for physical treatments due to its tunability of particle size, surface area, and functional groups, while chemical treatments with acids or alkalis are used to tailor biochar porosity and wettability. Overall, it was determined that future research needs to be done relating biochar production and post-processing methods to resulting 3D printing parameters as the number of studies is limited.

Day, Rachel [Auburn University, Auburn, AL]↗
DOE OSTI2024

ThunderBoltz API

The ThunderBoltz application programming interface (API) code is written in Python and is comprised of a set of tools to facilitate compilation of the ThunderBoltz code, as well as fast assembly and formatting of input files for the ThunderBoltz code, post-processing tools of ThunderBoltz results, plotting tools, and runs/schedules the ThunderBoltz code executable for calculations. The ThunderBoltz API code is utilized for importing and manipulating input cross section sets, input conditions, and any other simulation settings made available within the ThunderBoltz input deck via user-defined settings or via automatic generation. The API comes with a set of plotting capabilities of input cross sections, results from ThunderBoltz, and post-processed results carried out with the API.

Park, Ryan↗
DOE OSTI2015

FEASIBILITY OF HELIOSPHERIC IMAGING FROM NEAR EARTH

Imaging solar wind structures via Thomson scattered sunlight has proved important to understanding the inner heliosphere. The principal challenge of heliospheric imaging is background subtraction: typical solar wind features are fainter than the zodiacal light and starfield by 2-3 orders of magnitude. Careful post-processing is required to separate the solar wind signal from the static background. Remnant background, and not photon noise, is the dominant noise source in current STEREO data. We demonstrate that 10× shorter exposure times would not strongly affect the noise level in these data. Further, we demonstrate that current processing techniques are sufficient to separate not only the existing background of the STEREO images but also diffuse variable backgrounds such as are expected to be seen from low Earth orbit. We report on a hare-and-hounds style study, demonstrating blind signal extraction from STEREO/HI-2 data that have been degraded by the addition of large-scale, time-dependent artifacts to simulate viewing through airglow or high-altitude aurora. We demonstrate removal of these effects via image processing, with little degradation compared to the original. Even with as few as three highly degraded source images over 48 hr, it is possible to detect and track large coronal mass ejections more than 40° from the Sun. This implies that neither the high altitude aurora discovered by Coriolis/SMEI, nor airglow effects seen from low Earth orbit, are impediments to a hypothetical next-generation heliospheric imager in low Earth orbit; and also that post-processing is as important to heliospheric image qualitiy as are optical contamination effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

HydraGNN_Predictive_GFM_2024 - Ensemble of predictive graph foundation models for ground state atomistic materials modeling

We provide the ensemble of fifteen pre-trained graph foundation models (GFMs) for atomistic materials modeling applications. Each one of the fifteen GFMs has been trained on five open-source datasets that (once aggregated) amount to over 154 million atomistic structures, which cover over two-thirds of the natural elements of the periodic table and that comprises a broad set of organic and inorganic compounds. This vast set of atomistic structures comprises ground state configurations that are dynamically stable (i.e., equilibrated structures with atomic forces approximately close to zero values) as well as dynamically unstable structures (i.e., non-equilibrium structures with non-negligible non-zero values of atomic forces). The ensemble of datasets aggregated does NOT include excited states. The datasets have been curated to remove atomistic structures with spectral norm of the force tensor above 100 eV/angstrom. Moreover, a linear term of the energy was computed for each dataset using a linear regression model that uses the chemical concentration of each natural element as regressor. The linear term predicted by the linear regression model has been subtracted from each original energy value to perform a re-alignment of the energy values across different electronic structures approximation theories performed to generate the diverse multi-source, multi-fidelity datasets. The folder "ADIOS_files" contains the set of pre-processed datasets in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used for the development and training of GFMs in this work. The "ADIOS_files" directory contains 6 sub-directories named as follows: - ANI1x-v3.bp - MPTrj-v3.bp - OC2020-20M-v3.bp - OC2020-v3.bp - OC2022-v3.bp - qm7x-v3.bp Each sub-directory contains the pre-processed datasets converted in Adaptable I/O System (ADIOS) format (https://www.exascaleproject.org/research-project/adios/) that have been used to the development, training, and performance testing of the ensemble go predictive graph foundation models. Each GFM was developed using HydraGNN (https://github.com/ORNL/HydraGNN) as underlying graph neural network (GNN) architecture. The multi-task learning (MTL) capability of HydraGNN was used to simultaneously train the GFMs on labeled values for direct predictions of energy (a total system property of an atomistic structure that measures the chemical stability) and atomic forces (an atomic level property of an atomistic structure that measures the dynamical stability). The hyper parameters of the GFM have been tuned using scalable hyperparameter optimization (HPO) algorithms implemented in the software DeepHyper (https://github.com/deephyper/deephyper). The pre-training of each HPO trial was performed using distributed data parallelism (DDP) to scale the training across 128 compute nodes of the exascale OLCF supercomputer Frontier. Each HPO trial was trained only for 10 epochs and an early stopping was performed to avoid wasting significant computational resources on GNN architectures that were clearly underperforming. For each HPO trial, the 'omnistat' tool developed by (AMD Research - Advanced Micro Device) was used to measure the total energy consumption in kWh. The ensemble of GFMs was obtained by selecting the fifteen best performing HPO trials. Four models have been selected for their clear advantage in accuracy, and these are the GFMs with IDs 229, 156, 147, 260. Additional eleven models have been selected based on judicious balance between accuracy and energy consumption needed for training, and these are the GFMs with IDs 165, 78, 137, 1, 175, 171, 181, 67, 179, 167, 351. Each selected GFM of the ensemble was continued to cumulate a total of at most 30 epochs. In some cases, the total number of epochs actually performed was les than 30 due to two combined factors: (1) the size of the GFM (i.e., the number of model parameters to train) and (2) the total wall-clock time for which the computational resources could be allocated on OLCF-Frontier. The "Ensemble_of_models" directory contains 15 sub-directories named as follows: - gfm_0.229 - gfm_0.156 - gfm_0.147 - gfm_0.260 - gfm_0.165 - gfm_0.78 - gfm_0.137 - gfm_0.1 - gfm_0.175 - gfm_0.171 - gfm_0.181 - gfm_0.67 - gfm_0.179 - gfm_0.167 - gfm_0.351 Each one of these sub-directories refers to one of the fifteen HPO trials that have been selected to continue the pre-training with at most 30 epochs. With each sub-directory associated with a specific HPO trial, the following files can be found: - config.json: file for argument parsing to develop and train an HydraGNN architecture - gfm_0.ID_epoch_N.pk: file with model parameters for HPO ID trial after N epochs of training The ensemble of fifteen GFM architectures was used for (1) ensemble averaging to stabilize the predictions of energy and atomic forces after pre-training for post-processing analysis and (2) ensemble uncertainty quantification (UQ). The code used to develop, pre-train, and load the pre-trained models for post-processing analysis is available on the ORNL-GitHub at the following link: https://github.com/ORNL/HydraGNN/tree/Predictive_GFM_2024

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of thermal treatments on electrochemical behavior of binder jetted 17-4 PH stainless steel

Binder jetted 17–4 PH stainless steel was post-processed to relate heat treatment, microstructure, and corrosion in 3.5 wt% NaCl. Specimens were sintered at 1380 or 1400 °C, solution-annealed at 1055 °C for 1 h, and aged at 482 °C for 1 h. Here, as-sintered parts showed α′-martensitic matrix with a δ-ferrite network and Cu-rich precipitates in ferrite; inclusions (MnS, NbC) promoted localized attack. solutionizing redistributed elements and reduced ferrite, while aging generated coherent Cu nano-precipitates. Corrosion resistance was highly sensitive to post-processing in which aged specimens exhibited the lowest corrosion current density and formed a thicker, stable Cr 2 O 3 -rich passive film, whereas sintered specimens degraded most. Pitting potential depended on sintering temperature and microstructure, with 1400 °C sintering yielding more positive pitting potentials and the best overall performance after aging. Although the 1400-solutionized condition showed a relatively noble pitting response versus 1380-solutionized, it displayed unstable corrosion kinetics attributed to an imperfect passive film linked to higher NbC density. XPS depth profiles corroborated these trends, showing thicker, more continuous Cr 2 O 3 in aged states and discontinuous/thinning oxides in less resistant conditions. Practically, high-temperature sintering (∼1400 °C) followed by solutionizing and aging is recommended, with further gains expected from reducing NbC/MnS populations and porosity via powder and process control.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Sci-Thur PM - Colourful Interactions: Highlights 04: A Fast Quantitative MRI Acquisition and Processing Pipeline for Radiation Treatment Planning and Simulation

MRI-only Radiation Treatment Planning (RTP) is becoming increasingly popular because of a simplified work-flow, and less inconvenience to the patient who avoids multiple scans. The advantages of MRI-based RTP over traditional CT-based RTP lie in its superior soft-tissue contrast, and absence of ionizing radiation dose. The lack of electron-density information in MRI can be addressed by automatic tissue classification. To distinguish bone from air, which both appear dark in MRI, an ultra-short echo time (UTE) pulse sequence may be used. Quantitative MRI parametric maps can provide improved tissue segmentation/classification and better sensitivity in monitoring disease progression and treatment outcome than standard weighted images. Superior tumor contrast can be achieved on pure T{sub 1} images compared to conventional T{sub 1}-weighted images acquired in the same scan duration and voxel resolution. In this study, we have developed a robust and fast quantitative MRI acquisition and post-processing work-flow that integrates these latest advances into the MRI-based RTP of brain lesions. Using 3D multi-echo FLASH images at two different optimized flip angles (both acquired in under 9 min, and 1mm isotropic resolution), parametric maps of T{sub 1}, proton-density (M{sub 0}), and T{sub 2}{sup *} are obtained with high contrast-to-noise ratio, and negligible geometrical distortions, water-fat shifts and susceptibility effects. An additional 3D UTE MRI dataset is acquired (in under 4 min) and post-processed to classify tissues for dose simulation. The pipeline was tested on four healthy volunteers and a clinical trial on brain cancer patients is underway.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

The Fifth International Workshop on Ice Nucleation phase 2 (FIN-02): laboratory intercomparison of ice nucleation measurements

Abstract. The second phase of the Fifth International Ice Nucleation Workshop (FIN-02)involved the gathering of a large number of researchers at the KarlsruheInstitute of Technology's Aerosol Interactions and Dynamics of the Atmosphere(AIDA) facility to promote characterization and understanding of icenucleation measurements made by a variety of methods used worldwide.Compared to the previous workshop in 2007, participation was doubled,reflecting a vibrant research area. Experimental methods involved sampling ofaerosol particles by direct processing ice nucleation measuring systems fromthe same volume of air in separate experiments using different ice nucleatingparticle (INP) types, and collections of aerosol particle samples ontofilters or into liquid for sharing amongst measurement techniques thatpost-process these samples. In this manner, any errors introduced bydifferences in generation methods when samples are shared across laboratorieswere mitigated. Furthermore, as much as possible, aerosol particle sizedistribution was controlled so that the size limitations of different methodswere minimized. The results presented here use data from the workshop toassess the comparability of immersion freezing measurement methods activatingINPs in bulk suspensions, methods that activate INPs in condensation and/orimmersion freezing modes as single particles on a substrate, continuous flowdiffusion chambers (CFDCs) directly sampling and processing particles wellabove water saturation to maximize immersion and subsequent freezing ofaerosol particles, and expansion cloud chamber simulations in which liquidcloud droplets were first activated on aerosol particles prior to freezing.The AIDA expansion chamber measurements are expected to be the closestrepresentation to INP activation in atmospheric cloud parcels in thesecomparisons, due to exposing particles freely to adiabatic cooling. The different particle types used as INPs included the minerals illite NX andpotassium feldspar (K-feldspar), two natural soil dusts representative of arable sandy loam(Argentina) and highly erodible sandy dryland (Tunisia) soils, respectively,and a bacterial INP (Snomax ® ). Consideredtogether, the agreement among post-processed immersion freezing measurementsof the numbers and fractions of particles active at different temperaturesfollowing bulk collection of particles into liquid was excellent, withpossible temperature uncertainties inferred to be a key factor in determiningINP uncertainties. Collection onto filters for rinsing versus directly intoliquid in impingers made little difference. For methods that activatedcollected single particles on a substrate at a controlled humidity at orabove water saturation, agreement with immersion freezing methods was good inmost cases, but was biased low in a few others for reasons that have not beenresolved, but could relate to water vapor competition effects. AmongstCFDC-style instruments, various factors requiring (variable) highersupersaturations to achieve equivalent immersion freezing activation dominatethe uncertainty between these measurements, and for comparison with bulkimmersion freezing methods. When operated above water saturation to includeassessment of immersion freezing, CFDC measurements often measured at orabove the upper bound of immersion freezing device measurements, but oftenunderestimated INP concentration in comparison to an immersion freezingmethod that first activates all particles into liquid droplets prior tocooling (the PIMCA-PINC device, or Portable Immersion Mode Cooling chAmber–Portable Ice Nucleation Chamber), and typically slightly underestimated INPnumber concentrations in comparison to cloud parcel expansions in the AIDAchamber; this can be largely mitigated when it is possible to raise therelative humidity to sufficiently high values in the CFDCs, although this isnot always possible operationally. Correspondence of measurements of INPs among direct sampling andpost-processing systems varied depending on the INP type. Agreement was bestfor Snomax ® particles in the temperature regimecolder than −10 ∘ C, where their ice nucleation activity is nearlymaximized and changes very little with temperature. At temperatures warmer than −10 ∘ C, Snomax ® INP measurements (allvia freezing of suspensions) demonstrated discrepancies consistent withprevious reports of the instability of its protein aggregates that appear tomake it less suitable as a calibration INP at these temperatures. ForArgentinian soil dust particles, there was excellent agreement across allmeasurement methods; measures ranged within 1 order of magnitude for INPnumber concentrations, active fractions and calculated active site densitiesover a 25 to 30 ∘ C range and 5 to 8 orders of correspondingmagnitude change in number concentrations. This was also the case for alltemperatures warmer than −25 ∘ C in Tunisian dust experiments. Incontrast, discrepancies in measurements of INP concentrations or active sitedensities that exceeded 2 orders of magnitude across a broad range of temperaturemeasurements found at temperatures warmer than −25 ∘ C in a previous study werereplicated for illite NX. Discrepancies also exceeded 2 orders of magnitude attemperatures of −20 to −25 ∘ C for potassium feldspar (K-feldspar), but these coincidedwith the range of temperatures at which INP concentrations increase rapidly atapproximately an order of magnitude per 2 ∘ C cooling forK-feldspar. These few discrepancies did not outweigh the overall positive outcomes of theworkshop activity, nor the future utility of this data set or future similarefforts for resolving remaining measurement issues. Measurements of the samematerials were repeatable over the time of the workshop and demonstratedstrong consistency with prior studies, as reflected by agreement of databroadly with parameterizations of different specific or general (e.g., soildust) aerosol types. The divergent measurements of the INP activity of illiteNX by direct versus post-processing methods were not repeated for otherparticle types, and the Snomax ® datademonstrated that, at least for a biological INP type, there is no expectedmeasurement bias between bulk collection and direct immediately processedfreezing methods to as warm as −10 ∘ C. Since particle size rangeswere limited for this workshop, it can be expected that for atmosphericpopulations of INPs, measurement discrepancies will appear due to thedifferent capabilities of methods for sampling the full aerosol sizedistribution, or due to limitations on achieving sufficient watersupersaturations to fully capture immersion freezing in direct processinginstruments. Overall, this workshop presents an improved picture of presentcapabilities for measuring INPs than in past workshops, and providesdirection toward addressing remaining measurement issues.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Scalable Hybrid Learning Techniques for Scientific Data Compression

Data compression is becoming critical for storing scientific data because many scientific applications need to store large amounts of data and post process this data for scientific discovery. Unlike image and video compression algorithms that limit errors to primary data (PD), scientists require compression techniques that accurately preserve derived quantities of interest (QoIs). Here, this article presents a physics-informed compression technique implemented as an end-to-end, scalable, GPU-based pipeline for data compression that addresses this requirement. Our hybrid compression technique combines machine learning techniques and standard compression methods. Specifically, we combine an autoencoder, an error-bounded lossy compressor to provide guarantees on raw data error, and a constraint satisfaction post-processing step to preserve the QoIs within a minimal error (generally less than floating point error). The effectiveness of the data compression pipeline is demonstrated by compressing nuclear fusion simulation data generated by a large-scale fusion code, XGC, which produces hundreds of terabytes of data in a single day. Our approach works within the ADIOS framework and results in compression by a factor of more than 150 while requiring only a few percent of the computational resources necessary for generating the data, making the overall approach highly effective for practical scenarios.

ITER↗
DOE OSTI2018

Correlation of Microstructure to Creep Response of Hot Isostatically Pressed and Aged Electron Beam Melted Inconel 718

In this paper, creep rupture samples were fabricated by additive manufacturing (AM) via electron beam melting (EBM) to study the effects of grain morphology (equiaxed/columnar grains) and loading direction (longitudinal/transverse) with respect to build direction on creep deformation at 923 K (650 °C) with applied stresses of 580 and 600 MPa. The observed minimum creep rates and creep rupture lives of EBM Inconel 718 after post-processing by hot isostatic pressing (HIP) were found to be comparable to wrought material. The material with equiaxed grains exhibited low creep strain (2 pct) and short creep lifetimes (800 hours), whereas longer times (approximately 4500 hours) and high creep strain (up to 23 pct) were observed for material with columnar grains. The high stress exponent (n > 14) reflected the resistance to dislocation motion by γ” particles during creep. A precipitate-free zone (PFZ) was observed around the grain boundary δ phase. Creep damage occurred as voids and cracks in the PFZ. Finally, optimal post-processing of EBM Inconel 718 material should be explored to prevent δ phase embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Constrained Versions of DEDICOM for Use in Unsupervised Part-Of-Speech Tagging

This reports describes extensions of DEDICOM (DEcomposition into DIrectional COMponents) data models [3] that incorporate bound and linear constraints. The main purpose of these extensions is to investigate the use of improved data models for unsupervised part-of-speech tagging, as described by Chew et al. [2]. In that work, a single domain, two-way DEDICOM model was computed on a matrix of bigram fre- quencies of tokens in a corpus and used to identify parts-of-speech as an unsupervised approach to that problem. An open problem identi ed in that work was the com- putation of a DEDICOM model that more closely resembled the matrices used in a Hidden Markov Model (HMM), speci cally through post-processing of the DEDICOM factor matrices. The work reported here consists of the description of several models that aim to provide a direct solution to that problem and a way to t those models. The approach taken here is to incorporate the model requirements as bound and lin- ear constrains into the DEDICOM model directly and solve the data tting problem as a constrained optimization problem. This is in contrast to the typical approaches in the literature, where the DEDICOM model is t using unconstrained optimization approaches, and model requirements are satis ed as a post-processing step.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Experimental characterization of hot-electron emission and shock dynamics in the context of the shock ignition approach to inertial confinement fusion

We report on planar target experiments conducted on the OMEGA-EP laser facility performed in the context of the shock ignition (SI) approach to inertial confinement fusion. The experiment aimed at characterizing the propagation of strong shock in matter and the generation of hot electrons (HEs), with laser parameters relevant to SI (1-ns UV laser beams with I ~10 16 W/cm 2 ). Time-resolved radiographs of the propagating shock front were performed in order to study the hydrodynamic evolution. The hot-electron source was characterized in terms of Maxwellian temperature, T h , and laser to hot-electron energy conversion efficiency η using data from different x-ray spectrometers. The post-processing of these data gives a range of the possible values for T h and η [i.e., T h [keV] ϵ (20, 50) and η ϵ (2%, 13%)]. These values are used as input in hydrodynamic simulations to reproduce the results obtained in radiographs, thus constraining the range for the HE measurements. According to this procedure, we found that the laser converts ~10% ± 4% of energy into hot electrons with T h = 27 ± 8 keV. The paper shows how the coupling of different diagnostics and numerical tools is required to sufficiently constrain the problem, solving the large ambiguity coming from the post-processing of spectrometers data. The effect of the hot electrons on the shock dynamics is then discussed, showing an increase in the pressure around the shock front. Furthermore, the low temperature found in this experiment without pre-compression laser pulses could be advantageous for the SI scheme, but the high conversion efficiency may lead to an increase in the shell adiabat, with detrimental effects on the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gravitation and Mesh Adaption

The Gravitation and Mesh Adaptation (GaMA) toolbox consists of a set of Matlab classes for modeling the environment around asteroids and comets. A variety of gravitational models and supporting algorithms from are consolidated alongside a custom meshing utility tailored for the application. This allows the user to work within a single streamlined environment to import and manipulate surface definitions, probe the dynamical environment, integrate trajectories, and post-process results. For other applications, modified surface meshes can be exported. Gravity Models: 1) Werner's analytic polyhedron model 2) Mascon model 3) Gottlieb's spherical harmonic model 4) Approximate polyhedron models 5) Curvilinear surface models 6) Custom composite models Meshing Features: 1) Array-based half-edge data structure 2) Supports curvilinear surface definitions up to degree 4 3) Ray-tracing 4) Coarsening 5) Feature-based refinement 6) Projection 7) Smoothing 8) Mesh quality validity tests 9) Mesh repair Additional Features: 1) Solar radiation pressure model 2) Distant 3rd body model 3) Collision detection 4) Trajectory integration and post-processing 5) Visualization of surface fields

Pearl, Jason↗
DOE OSTI2016

Characterizing the nano-structure and defect structure of nano-scaled non-ferrous structural alloys

The presence and interaction of nanotwins, geometrically necessary dislocations, and grain boundaries play a key role in the mechanical properties of nanostructured crystalline materials. Therefore, it is vital to determine the orientation, width and distance of nanotwins, the angle and axis of grain boundary misorientations as well as the type and the distributions of dislocations in an automatic and statistically meaningful fashion in a relatively large area. In this paper, such details are provided using a transmission electron microscope-based orientation microscopy technique called ASTAR™/precession electron diffraction. The remarkable spatial resolution of this technique (~ 2 nm) enables highly detailed characterization of nanotwins, grain boundaries and the configuration of dislocations. This orientation microscopy technique provides the raw data required for the determination of these parameters. The procedures to post-process the ASTAR™/PED datasets in order to obtain the important (and currently largely hidden) details of nanotwins as well as quantifications of dislocation density distributions are described in this study. - Highlights: • EBSD cannot characterize defects such as dislocations, grain boundaries and nanotwins in severely deformed metals. • TEM based orientation microscopy technique called ASTAR™/PED was used to resolve the problem. • Locations and orientations of nanotwins, dislocation density distribution and grain boundary characters can be resolved. • This work provides the bases for further studies on the interactions between dislocations, grain boundaries and nanotwins. • The computation part is explained sufficiently which helps the readers to post process their own data.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-FG-BRB-03: Basis Vector Model Based Method for Proton Stopping Power Estimation From Experimental Dual Energy CT Data

Purpose: This work aims at reducing the uncertainty in proton stopping power (SP) estimation by a novel combination of a linear, separable basis vector model (BVM) for stopping power calculation (Med Phys 43:600) and a statistical, model-based dual-energy CT (DECT) image reconstruction algorithm (TMI 35:685). The method was applied to experimental data. Methods: BVM assumes the photon attenuation coefficients, electron densities, and mean excitation energies (I-values) of unknown materials can be approximated by a combination of the corresponding quantities of two reference materials. The DECT projection data for a phantom with 5 different known materials was collected on a Philips Brilliance scanner using two scans at 90 kVp and 140 kVp. The line integral alternating minimization (LIAM) algorithm was used to recover the two BVM coefficient images using the measured source spectra. The proton stopping powers are then estimated from the Bethe-Bloch equation using electron densities and I-values derived from the BVM coefficients. The proton stopping powers and proton ranges for the phantom materials estimated via our BVM based DECT method are compared to ICRU reference values and a post-processing DECT analysis (Yang PMB 55:1343) applied to vendorreconstructed images using the Torikoshi parametric fit model (tPFM). Results: For the phantom materials, the average stopping power estimations for 175 MeV protons derived from our method are within 1% of the ICRU reference values (except for Teflon with a 1.48% error), with an average standard deviation of 0.46% over pixels. The resultant proton ranges agree with the reference values within 2 mm. Conclusion: Our principled DECT iterative reconstruction algorithm, incorporating optimal beam hardening and scatter corrections, in conjunction with a simple linear BVM model, achieves more accurate and robust proton stopping power maps than the post-processing, nonlinear tPFM based DECT analysis applied to conventional reconstructions of low and high energy scans. Funding Support: NIH R01CA 75371; NCI grant R01 CA 149305.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Fast decay of classification error in variational quantum circuits

Variational quantum circuits (VQCs) have shown great potential in near-term applications. However, the discriminative power of a VQC, in connection to its circuit architecture and depth, is not understood. To unleash the genuine discriminative power of a VQC, we propose a VQC system with the optimal classical post-processing—maximum-likelihood estimation on measuring all VQC output qubits. Via extensive numerical simulations, we find that the error of VQC quantum data classification typically decays exponentially with the circuit depth, when the VQC architecture is extensive—the number of gates does not shrink with the circuit depth. This fast error suppression ends at the saturation towards the ultimate Helstrom limit of quantum state discrimination. On the other hand, non-extensive VQCs such as quantum convolutional neural networks are sub-optimal and fail to achieve the Helstrom limit, demonstrating a trade-off between ansatz complexity and classification performance in general. To achieve the best performance for a given VQC, the optimal classical post-processing is crucial even for a binary classification problem. To simplify VQCs for near-term implementations, we find that utilizing the symmetry of the input properly can improve the performance, while oversimplification can lead to degradation.

Zhang, Bingzhi↗
DOE OSTI2022

Multi-scale characterization of supersolidus liquid phase sintered H13 tool steel manufactured via binder jet additive manufacturing

Additive manufacturing (AM) of H13 tool steel by binder jet 3D printing (BJ3DP) followed by pressureless supersolidus liquid phase sintering (SLPS) provides a low-cost alternative manufacturing method for components with intricate geometric features. However, the microstructure-mechanical property relationships for BJ3DP-SLPS produced H13 tool steel are not well understood, which makes it challenging to develop printing and post-processing methods that maximize part performance. In this work, we leverage atom probe tomography and transmission electron microscopy along with thermodynamic calculations to rationalize the microstructure-mechanical property relationships in as-sintered BJ3DP H13 tool steel. Here, we report for the first time, the presence of a continuous eutectic film-like carbide in H13 along with the more commonly observed cuboidal MX carbides in the prior liquid channels of the microstructure. Further, atom probe tomography revealed the interconnected nature of the MX carbides that appear to be discrete in two-dimensional micrographs. These continuous eutectic carbides and interconnected MX carbides result in brittle failure of the material. Characterization of these microstructural features will be critical in developing appropriate post-processing heat treatments for the improved mechanical performance of BJ3DP H13.

36 MATERIALS SCIENCE↗
DOE OSTI2023

AmeriFlux US-GL1 Stannard Rock

This is the AmeriFlux version of the carbon flux data for the site US-GL1 Stannard Rock. Site Description - Stannard Rock is located 39 km from the nearest shore (Keweenaw Peninsula) in Lake Superior, 44 miles NNE of Marquette, Michigan, and 24 miles ESE of Manitou Island. The site is located on the historic Stannard Rock Lighthouse, which was completed in 1882. Eddy covariance instrumentation was installed in 2008 by a network of scientists from both US and Canada, eventually to be called the Great Lakes Evaporation Network (GLEN). The intent of GLEN has been to provide observations of over-lake meteorology and evaporation, improve forecasting of Great Lakes water levels, and support a wide variety of stakeholders, including the National Weather Service (NWS), Environment and Climate Change Canada, National Oceanic and Atmospheric Administration, U.S. Coast Guard, recreational boaters and commercial shipping, emergency management officials, and the Great Lakes research community. The eddy covariance station along with other ancillary meteorological instrumentation is located at an approximate elevation of 39.2 meters above mean lake water level. Meteorological data from the lighthouse are sent to the National Data Buoy Center, where they can be viewed in real-time at http://www.ndbc.noaa.gov/station_page.php?station=stdm4 Uncorrected half-hour fluxes were computed directly on the logger using a 30-minute block averaging period as high-frequency data were not available. The half hour flux measurements downloaded from the datalogger were post-processed with the following filters and corrections. Latent and sensible heat and carbon dioxide fluxes were corrected with 2-D coordinate rotation. Although the primary objective of data collection at Stannard Rock was to quantify the evaporative flux, additional preliminary measurements of the carbon dioxide concentration and flux from the LI-7500 are also included in this dataset but it is advised to use the carbon data with caution. Carbon data reported in this dataset includes turbulent fluxes of CO2 with no storage correction (FC, µmol m-2 s-1) and CO2 density in mole fraction of wet air (CO2), which was originally output on the datalogger as average CO2 density (mg m-2 s-1) and converted into µmol mol-1 using air temperature and pressure in post-processing. Webb, Pearman, and Leuning terms were applied to account for density fluctuations for water vapor and CO2. Sonic path length, high-frequency attenuation and sensor separation were accounted for according to Horst and Massman, and the oxygen absorption correction for the KH2O hygrometer was also applied. Latent and sensible heat fluxes were assumed to be unrealistic above an absolute value of 1000 W m-2 and were removed. Both carbon flux (FC) and carbon dioxide mole fraction in wet air (CO2) and were assumed to be unrealistic above 1000 µmol m-2 s-1 and 1000 µmol mol-1 respectively. Spikes in latent and sensible heat and carbon fluxes and densities (often due to periods of precipitation) were identified by computing the mean and standard deviation over a moving, overlapping window of 336 half-hours (7 days), similar to Shao et al., and were removed when the flux was more than 1.5 standard deviations from the moving window’s mean. While Vickers and Mahrt use a threshold of 3.5 standard deviations from the mean, a conservative value of 1.5 was chosen due to the noisy nature of over-lake data at this particular site. This process was repeated twice for latent and sensible heat, and carbon dioxide fluxes and densities and therefore it is possible that some real, realistic data was filtered out in this process. No detrending was performed. As per AmeriFlux standards, no friction velocity (USTAR, m s-1) filtering was performed.

Spence, Chris↗
DOE OSTI2025

AmeriFlux FLUXNET-1F US-GL1 Stannard Rock

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-GL1 Stannard Rock. This is the FLUXNET version of the carbon flux data for the site US-GL1 Stannard Rock produced by applying the standard ONEFlux (1F) software. Site Description - Stannard Rock is located 39 km from the nearest shore (Keweenaw Peninsula) in Lake Superior, 44 miles NNE of Marquette, Michigan, and 24 miles ESE of Manitou Island. The site is located on the historic Stannard Rock Lighthouse, which was completed in 1882. Eddy covariance instrumentation was installed in 2008 by a network of scientists from both US and Canada, eventually to be called the Great Lakes Evaporation Network (GLEN). The intent of GLEN has been to provide observations of over-lake meteorology and evaporation, improve forecasting of Great Lakes water levels, and support a wide variety of stakeholders, including the National Weather Service (NWS), Environment and Climate Change Canada, National Oceanic and Atmospheric Administration, U.S. Coast Guard, recreational boaters and commercial shipping, emergency management officials, and the Great Lakes research community. The eddy covariance station along with other ancillary meteorological instrumentation is located at an approximate elevation of 39.2 meters above mean lake water level. Meteorological data from the lighthouse are sent to the National Data Buoy Center, where they can be viewed in real-time at http://www.ndbc.noaa.gov/station_page.php?station=stdm4 Uncorrected half-hour fluxes were computed directly on the logger using a 30-minute block averaging period as high-frequency data were not available. The half hour flux measurements downloaded from the datalogger were post-processed with the following filters and corrections. Latent and sensible heat and carbon dioxide fluxes were corrected with 2-D coordinate rotation. Although the primary objective of data collection at Stannard Rock was to quantify the evaporative flux, additional preliminary measurements of the carbon dioxide concentration and flux from the LI-7500 are also included in this dataset but it is advised to use the carbon data with caution. Carbon data reported in this dataset includes turbulent fluxes of CO2 with no storage correction (FC, µmol m-2 s-1) and CO2 density in mole fraction of wet air (CO2), which was originally output on the datalogger as average CO2 density (mg m-2 s-1) and converted into µmol mol-1 using air temperature and pressure in post-processing. Webb, Pearman, and Leuning terms were applied to account for density fluctuations for water vapor and CO2. Sonic path length, high-frequency attenuation and sensor separation were accounted for according to Horst and Massman, and the oxygen absorption correction for the KH2O hygrometer was also applied. Latent and sensible heat fluxes were assumed to be unrealistic above an absolute value of 1000 W m-2 and were removed. Both carbon flux (FC) and carbon dioxide mole fraction in wet air (CO2) and were assumed to be unrealistic above 1000 µmol m-2 s-1 and 1000 µmol mol-1 respectively. Spikes in latent and sensible heat and carbon fluxes and densities (often due to periods of precipitation) were identified by computing the mean and standard deviation over a moving, overlapping window of 336 half-hours (7 days), similar to Shao et al., and were removed when the flux was more than 1.5 standard deviations from the moving window’s mean. While Vickers and Mahrt use a threshold of 3.5 standard deviations from the mean, a conservative value of 1.5 was chosen due to the noisy nature of over-lake data at this particular site. This process was repeated twice for latent and sensible heat, and carbon dioxide fluxes and densities and therefore it is possible that some real, realistic data was filtered out in this process. No detrending was performed. As per AmeriFlux standards, no friction velocity (USTAR, m s-1) filtering was performed.

Spence, Chris [Environment and Climate Change Cana↗