Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A perturbation calculation of the dynamical effects of cosmic rays and interstellar gas on the solar wind

The modulation of galactic cosmic rays and other energetic particles by the solar wind produces a gradient in their pressure, which in turn influences the wind dynamics. The basic equations describing this interaction are presented in the 'hydrodynamic' approximation. A perturbation solution of the equations is presented for the case in which both the galactic cosmic ray pressure, and the pressure of the anomalous cosmic ray component accelerated at the wind termination shock, are small compared with the wind ram pressure. Analytical expressions for the deceleration of the wind and the modification of the shock and its location are derived in this case. These effects are estimated to have a relative magnitude of several percent in the solar wind. The interstellar neutral gas which penetrates the heliosphere is ionized, predominantly by photoionization and charge exchange with the wind, and may also have a significant dynamical effect on the wind. The basic equations describing this interaction are also presented, A perturbation solution is presented under the assumption that the mass loading and momentum-loading of the wind by the interstellar pickup ions is small compared with the wind ram pressure. Analytical expressions for the deceleration of the wind, the contribution of the pickup ions to the wind pressure, and the modification of the termination shock location are derived. Again, with the exception of pickup ion pressure which is large compared with solar wind thermal pressure, the effects are estimated to be several percent in relative magnitude in the solar wind.

Lee, M. A.↗
DOE OSTI2023

TaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-vanadium (Tanu-V). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the element V. For Ta, p semi-core states as valence were chosen. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name case-N where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. and Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15-50 (1996) (4) Kresse, G. and Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2025

TaTi_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-titanium (Ta-Ti). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-niobium (Hf-Nb). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Nb. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Calculation of two- and three-dimensional transonic cascade flow field using the Navier-Stokes equations

A Navier-Stokes analysis employing the time-dependent Linearized Block Implicit scheme (LBI) was applied to two-dimensional and three-dimensional transonic turbulent cascade flows. In general, the geometrical configuration of the turbine blade impacts both the grid construction procedure and the implementation of the numerical algorithm. Since modern turbine blades of interest are characterized by very blunt leading edges, rounded trailing edges and high stacking angles, a robust grid construction procedure is required that can accommodate the severe body shape while resolving regions of large flow gradients. A constructive O-type grid generation technique, suitable for cascades with rounded trailing edges, was developed and used to construct the C3X turbine cascade coordinate grid. Two-dimensional calculations were performed employing the Navier-Stokes procedure for the C3X turbine cascade, and the predicted pressure coefficients and heat transfer rates were compared with the experimental data. Three-dimensional Navier-Stokes calculations were also performed.

Weinberg, B. C.↗
DOE OSTI2018

Optimal beam margins in linac-based VMAT stereotactic ablative body radiotherapy: a Pareto front analysis for liver metastases

We explored the Pareto fronts mathematical strategy to determine the optimal block margin and prescription isodose for stereotactic body radiotherapy (SBRT) treatments of liver metastases using the volumetric-modulated arc therapy (VMAT) technique. Three targets (planning target volumes [PTVs] = 20, 55, and 101 cc) were selected. A single fraction dose of 26 Gy was prescribed (prescription dose [PD]). VMAT plans were generated for 3 different beam energies. Pareto fronts based on (1) different multileaf collimator (MLC) block margin around PTV and (2) different prescription isodose lines (IDL) were produced. For each block margin, the greatest IDL fulfilling the criteria (95% of PTV reached 100%) was considered as providing the optimal clinical plan for PTV coverage. Liver D{sub mean}, V7Gy, and V12Gy were used against the PTV coverage to generate the fronts. Gradient indexes (GI and mGI), homogeneity index (HI), and healthy liver irradiation in terms of D{sub mean}, V7Gy, and V12Gy were calculated to compare different plans. In addition, each target was also optimized with a full-inverse planning engine to obtain a direct comparison with anatomy-based treatment planning system (TPS) results. About 900 plans were calculated to generate the fronts. GI and mGI show a U-shaped behavior as a function of beam margin with minimal values obtained with a +1 mm MLC margin. For these plans, the IDL ranges from 74% to 86%. GI and mGI show also a V-shaped behavior with respect to HI index, with minimum values at 1 mm for all metrics, independent of tumor dimensions and beam energy. Full-inversed optimized plans reported worse results with respect to Pareto plans. In conclusion, Pareto fronts provide a rigorous strategy to choose clinical optimal plans in SBRT treatments. We show that a 1-mm MLC block margin provides the best results with regard to healthy liver tissue irradiation and steepness of dose fallout.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Aspect ratio effects on limited scrape-off layer plasma turbulence

The drift-reduced Braginskii model describing turbulence in the tokamak scrape-off layer is written for a general magnetic configuration with a limiter. The equilibrium is then specified for a circular concentric magnetic geometry retaining aspect ratio effects. Simulations are then carried out with the help of the global, flux-driven fluid three-dimensional code GBS [Ricci et al., Plasma Phys. Controlled Fusion 54, 124047 (2012)]. Linearly, both simulations and simplified analytical models reveal a stabilization of ballooning modes. Nonlinearly, flux-driven nonlinear simulations give a pressure characteristic length whose trends are correctly captured by the gradient removal theory [Ricci and Rogers, Phys. Plasmas 20, 010702 (2013)], that assumes the profile flattening from the linear modes as the saturation mechanism. More specifically, the linear stabilization of ballooning modes is reflected by a 15% increase in the steady-state pressure gradient obtained from GBS nonlinear simulations when going from an infinite to a realistic aspect ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Dosimetric study of surface applicators of HDR brachytherapy GammaMed Plus equipment

The cone type surface applicators used in HDR brachytherapy for treatment of small skin lesions are an alternative to be used with both electron beams and orthovoltage X-ray equipment. For a good treatment planning is necessary to know the dose distribution of these applicators, which can be obtained by experimental measurement and Monte Carlo simulation as well. In this study the dose distribution of surface applicators of 3 and 3.5 cm diameter, respectively of HDR brachytherapy GammaMed Plus equipment has been estimated using the Monte Carlo method, MCNP code. The applicators simulated were placed on the surface of a water phantom of 20 × 20 × 20 cm and the dose was calculated at depths from 0 to 3 cm with increments of 0.25 mm. The dose profiles obtained at depth show the expected gradients for surface therapy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2015

Crystal structure optimisation using an auxiliary equation of state

Standard procedures for local crystal-structure optimisation involve numerous energy and force calculations. It is common to calculate an energy–volume curve, fitting an equation of state around the equilibrium cell volume. This is a computationally intensive process, in particular, for low-symmetry crystal structures where each isochoric optimisation involves energy minimisation over many degrees of freedom. Such procedures can be prohibitive for non-local exchange-correlation functionals or other “beyond” density functional theory electronic structure techniques, particularly where analytical gradients are not available. We present a simple approach for efficient optimisation of crystal structures based on a known equation of state. The equilibrium volume can be predicted from one single-point calculation and refined with successive calculations if required. The approach is validated for PbS, PbTe, ZnS, and ZnTe using nine density functionals and applied to the quaternary semiconductor Cu{sub 2}ZnSnS{sub 4} and the magnetic metal-organic framework HKUST-1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

A global wave-driven magnetohydrodynamic solar model with a unified treatment of open and closed magnetic field topologies

We describe, analyze, and validate the recently developed Alfvén Wave Solar Model, a three-dimensional global model starting from the top of the chromosphere and extending into interplanetary space (out to 1-2 AU). This model solves the extended, two-temperature magnetohydrodynamics equations coupled to a wave kinetic equation for low-frequency Alfvén waves. In this picture, heating and acceleration of the plasma are due to wave dissipation and to wave pressure gradients, respectively. The dissipation process is described by a fully developed turbulent cascade of counterpropagating waves. We adopt a unified approach for calculating the wave dissipation in both open and closed magnetic field lines, allowing for a self-consistent treatment in any magnetic topology. Wave dissipation is the only heating mechanism assumed in the model; no geometric heating functions are invoked. Electron heat conduction and radiative cooling are also included. We demonstrate that the large-scale, steady state (in the corotating frame) properties of the solar environment are reproduced, using three adjustable parameters: the Poynting flux of chromospheric Alfvén waves, the perpendicular correlation length of the turbulence, and a pseudoreflection coefficient. We compare model results for Carrington rotation 2063 (2007 November-December) with remote observations in the extreme-ultraviolet and X-ray ranges from the Solar Terrestrial Relations Observatory, Solar and Heliospheric Observatory, and Hinode spacecraft and with in situ measurements by Ulysses. The results are in good agreement with observations. This is the first global simulation that is simultaneously consistent with observations of both the thermal structure of the lower corona and the wind structure beyond Earth's orbit.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Gradient corrections to the exchange-correlation free energy

We develop the first-order gradient correction to the exchange-correlation free energy of the homogeneous electron gas for use in finite-temperature density functional calculations. Based on this, we propose and implement a simple temperature-dependent extension for functionals beyond the local density approximation. These finite-temperature functionals show improvement over zero-temperature functionals, as compared to path-integral Monte Carlo calculations for deuterium equations of state, and perform without computational cost increase compared to zero-temperature functionals and so should be used for finite-temperature calculations. Furthermore, while the present functionals are valid at all temperatures including zero, non-negligible difference with zero-temperature functionals begins at temperatures above 10 000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling a Nb3Sn Cryounit in GPT at UITF

Nb₃Sn is a prospective material for future superconducting RF (SRF) accelerator cavities. The material can achieve higher quality factors, higher temperature operation and potentially higher accelerating gradients (E_{acc} 96 MV/m) compared to conventional niobium. In this work, we performed modeling of the Upgraded Injector Test Facility (UITF) at Jefferson Lab utilizing newly constructed Nb₃Sn cavities. We studied the effects of the buncher cavity and varied the gun voltages from 200-500 keV. We have calibrated and optimized the SRF cavity gradients and phases for the Nb₃Sn five-cell cavities energy gains with the framework of General Particle Tracer (GPT). Our calculations show the beam goes cleanly through the unit. There is full energy gain out of the second SRF cavity but not from the first SRF cavity due to non-relativistic phase shifts.

Pokharel, S.↗
DOE OSTI2013

Summary Statistics for LC_TA03_130719: Micro-CT Data Acquired at LLNL, Specimen 3 of 3

Using x-ray micro computed tomography (MicroCT), we have characterized the linear attenuation coefficients (LAC), μ, of a sample of a dry powder material, lithium carbonate (LC). The specimen was prepared at Lawrence Livermore National Laboratory (LLNL), loaded into a 60mL low density polyethylene (LDPE) bottle. After completed packing, the specimen was scanned following the protocol for MicroCT measurements under Test Plan 79. This particular specimen, LC_TA03_130719, recorded the bulk packing density (mass of sample divided by volume of sample) shown above. Two additional preparations were made and analyzed. We used the computer program IMGREC to reconstruct the CT images. The values of the key parameters used in the x-ray data capture and image reconstruction are given in this report. Additional experimental details may be found in the SOP and a separate document. To characterize the statistical distribution of LAC values in each CT image, we first isolated an ~80% region or segment of volume elements (“voxels”) lying completely within the sample, away from the walls of the container. We then calculated the mean value, standard deviation and entropy for (a) the high and low energy image segments and for (b) their digital gradient images 2 . The statistics of the initial image of LAC values are called “first order statistics;” those of the gradient image, “second order statistics.” See Seetho for details of the analysis used to obtain the numbers reported in this document.

42 ENGINEERING↗
DOE OSTI2013

Summary Statistics for TA_MP03_130620: Micro-CT Data Acquired at LLNL, Specimen 3 of 3

Using x-ray micro computed tomography (MicroCT), we have characterized the linear attenuation coefficients (LAC), μ, of a sample of a dry powder material, tartaric acid (TA). The specimen was prepared at Lawrence Livermore National Laboratory (LLNL), loaded into a 60mL low density polyethylene (LDPE) bottle. After completed packing, the specimen was scanned following the protocol for MicroCT measurements under Test Plan 79. This particular specimen, TA_MP03_130620, recorded the bulk packing density (mass of sample divided by volume of sample) shown above. Two additional preparations were made and analyzed. We used the computer program IMGREC to reconstruct the CT images. The values of the key parameters used in the x-ray data capture and image reconstruction are given in this report. Additional experimental details may be found in the SOP and a separate document. To characterize the statistical distribution of LAC values in each CT image, we first isolated an ~80% region or segment of volume elements (“voxels”) lying completely within the sample, away from the walls of the container. We then calculated the mean value, standard deviation and entropy for (a) the high and low energy image segments and for (b) their digital gradient images 2 . The statistics of the initial image of LAC values are called “first order statistics;” those of the gradient image, “second order statistics.” See Seetho for details of the analysis used to obtain the numbers reported in this document.

42 ENGINEERING↗
DOE OSTI2013

Summary Statistics for TA_MH03_130828: Micro-CT Data Acquired at LLNL, Specimen 3 of 3

Using x-ray micro computed tomography (MicroCT), we have characterized the linear attenuation coefficients (LAC), μ, of a sample of a dry powder material, tartaric acid (TA). The specimen was prepared at Lawrence Livermore National Laboratory (LLNL), loaded into a 60mL low density polyethylene (LDPE) bottle. After completed packing, the specimen was scanned following the protocol for MicroCT measurements under Test Plan 79. This particular specimen, TA_MH03_130828, recorded the bulk packing density (mass of sample divided by volume of sample) shown above. Two additional preparations were made and analyzed. We used the computer program IMGREC to reconstruct the CT images. The values of the key parameters used in the x-ray data capture and image reconstruction are given in this report. Additional experimental details may be found in the SOP and a separate document. To characterize the statistical distribution of LAC values in each CT image, we first isolated an ~80% region or segment of volume elements (“voxels”) lying completely within the sample, away from the walls of the container. We then calculated the mean value, standard deviation and entropy for (a) the high and low energy image segments and for (b) their digital gradient images. The statistics of the initial image of LAC values are called “first order statistics;” those of the gradient image, “second order statistics.” See Seetho for details of the analysis used to obtain the numbers reported in this document.

42 ENGINEERING↗
DOE OSTI2013

Summary Statistics for TA_MH01_130826: Micro-CT Data Acquired at LLNL, Specimen 1 of 3

Using x-ray micro computed tomography (MicroCT), we have characterized the linear attenuation coefficients (LAC), μ, of a sample of a dry powder material, tartaric acid (TA). The specimen was prepared at Lawrence Livermore National Laboratory (LLNL), loaded into a 60mL low density polyethylene (LDPE) bottle. After completed packing, the specimen was scanned following the protocol for MicroCT measurements under Test Plan 79. This particular specimen, TA_MH01_130826, recorded the bulk packing density (mass of sample divided by volume of sample) shown above. Two additional preparations were made and analyzed. We used the computer program IMGREC to reconstruct the CT images. The values of the key parameters used in the x-ray data capture and image reconstruction are given in this report. Additional experimental details may be found in the SOP and a separate document. To characterize the statistical distribution of LAC values in each CT image, we first isolated an ~80% region or segment of volume elements (“voxels”) lying completely within the sample, away from the walls of the container. We then calculated the mean value, standard deviation and entropy for (a) the high and low energy image segments and for (b) their digital gradient images. The statistics of the initial image of LAC values are called “first order statistics;” those of the gradient image, “second order statistics.” See Seetho for details of the analysis used to obtain the numbers reported in this document.

42 ENGINEERING↗
DOE OSTI2013

Summary Statistics for TA_MH02_130827: Micro-CT Data Acquired at LLNL, Specimen 2 of 3

Using x-ray micro computed tomography (MicroCT), we have characterized the linear attenuation coefficients (LAC), μ, of a sample of a dry powder material, tartaric acid (TA). The specimen was prepared at Lawrence Livermore National Laboratory (LLNL), loaded into a 60mL low density polyethylene (LDPE) bottle. After completed packing, the specimen was scanned following the protocol for MicroCT measurements under Test Plan 79. This particular specimen, TA_MH03_130827, recorded the bulk packing density (mass of sample divided by volume of sample) shown above. Two additional preparations were made and analyzed. We used the computer program IMGREC to reconstruct the CT images. The values of the key parameters used in the x-ray data capture and image reconstruction are given in this report. Additional experimental details may be found in the SOP and a separate document. To characterize the statistical distribution of LAC values in each CT image, we first isolated an ~80% region or segment of volume elements (“voxels”) lying completely within the sample, away from the walls of the container. We then calculated the mean value, standard deviation and entropy for (a) the high and low energy image segments and for (b) their digital gradient images. The statistics of the initial image of LAC values are called “first order statistics;” those of the gradient image, “second order statistics.” See Seetho for details of the analysis used to obtain the numbers reported in this document.

42 ENGINEERING↗