Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tomographic reconstruction”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Laser-driven high-resolution MeV x-ray tomography

The need for high-resolution MeV x-ray tomography to observe the three-dimensional structure of dense, large-sized objects is rapidly increasing for the non-destructive evaluation of critical additively manufactured parts, national security, and other applications. We report a demonstration of high-resolution MeV computed tomography of a dense, large object with a laser-driven x-ray source. A record detector-limited MeV radiograph resolution of < 200 µm as determined with the Bennett approximation of the point spread function was achieved by irradiating millimeter-thick tungsten targets with 300 TW femtosecond laser pulses at a 0.5 Hz repetition rate. A tungsten alloy step wedge spectrometer indicates that the peak of the x-ray emission is between 1 and 2 MeV, with an endpoint energy of 19 MeV. To illustrate the radiographic imaging capability of the system, a tomographic reconstruction of a nickel superalloy turbine blade ( maximum $ρr$ = 139 g/cm 2 ) with sub-millimeter resolution was performed using 2160 individual radiographs. The small x-ray source size opens the prospect of extremely high-resolution tomographs of large, dense objects. This laser-driven approach has major advantages for non-destructive evaluation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

4D Imaging of ZnO-Coated Nanoporous Al 2 O 3 Aerogels by Chemically Sensitive Ptychographic Tomography: Implications for Designer Catalysts

The 3D chemical structure (4D spectromicroscopy) of nanoporous Al 2 O 3 aerogels coated with ZnO by atomic layer deposition (ALD) was studied by multienergy scanning transmission X-ray microscopy. These materials are representative of a class of designer catalysts in which the nanoporous support is prepared separately from the active catalyst material, which is subsequently introduced by ALD, thereby allowing independent optimization of the morphology, chemistry, and spatial distribution of the support and catalyst. The samples studied were prepared by Ga ion and Xe plasma focused ion beam (FIB) milling as well as drop casting from water suspension. Zn L and Al K edge spectra of six samples with three different ZnO loadings were measured to investigate how loading and different sample preparation methods affect the 3D distribution of the ZnO and Al 2 O 3 . Scanning transmission X-ray microscopy (STXM) and ptychographic imaging at two energies each at the Zn L 3 and Al K edge were measured. The ptychography data were analyzed by using the SHARP reconstruction code to generate quantitative 2D chemical maps of the Al 2 O 3 and the ZnO. The STXM and ptychography maps were then measured at a sequence of tilt angles, covering up to 160° of rotation. The 3D structure of the ZnO and Al 2 O 3 was derived from the tilt series data by tomographic reconstruction using a compressed sensing algorithm. A twodimensional spatial resolution (half-period) of 6 nm, measured by Fourier ring correlation, and a 3D spatial resolution (half-period) of 9 nm, measured by Fourier shell correlation, were achieved when using the COSMIC beamline at the Advanced Light Source (ALS). The results show that for all of the ZnO loadings studied there is nonuniform coverage of the ZnO on the Al 2 O 3 aerogel framework. In addition, we found that both FIB methods create sample artifacts, although the distortion was less with Xe plasma than Ga ion FIB.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Evaluation of Electron Tomography Capabilities for Shale Imaging

Despite the advantageous resolution of electron tomography (ET), reconstruction of three-dimensional (3D) images from multiple two-dimensional (2D) projections presents several challenges, including small signal-to-noise ratios, and a limited projection range. This study evaluates the capabilities of ET for thin sections of shale, a complex nanoporous medium. A numerical phantom with 1.24 nm pixel size is constructed based on the tomographic reconstruction of a Barnett shale. A dataset of 2D projection images is numerically generated from the 3D phantom and studied over a range of conditions. First, common reconstruction techniques are used to reconstruct the shale structure. The reconstruction uncertainty is quantified by comparing overall values of storage and transport metrics, as well as the misclassification of pore voxels compared to the phantom. We then select the most robust reconstruction technique and we vary the acquisition conditions to quantify the effect of artifacts. We find a strong agreement for large pores over the different acquisition workflows, while a wider variability exists for nanometer-scale features. Furthermore, the limited projection range and reconstruction are identified as the main experimental bottlenecks, thereby suggesting that sample thinning, advanced holders, and advanced reconstruction algorithms offer opportunities for improvement.

3D scanning transmission electron microscopy↗
NASA NTRS1993

Tomographic decoding algorithm for a nonoverlapping redundant array

A tomographic reconstruction algorithm is developed for the nonoverlapping redundant array X-ray imaging system whereby the background contributions from out-of-focus planes can be eliminated. The algorithm makes use of two constraints derived from the physical characteristics of the nonoverlapping redundant array system in tandem with the correlation decoding process. It is simple, direct, and noniterative. Tomographic images of computer-generated planar and three-dimensional objects are provided to illustrate the effectiveness of the algorithm.

Yin, Lo I.↗
NASA NTRS2002

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗
NASA NTRS2003

Hierarchies of Models: Toward Understanding Planetary Nebulae

Stars like our sun (initial masses between 0.8 to 8 solar masses) end their lives as swollen red giants surrounded by cool extended atmospheres. The nuclear reactions in their cores create carbon, nitrogen and oxygen, which are transported by convection to the outer envelope of the stellar atmosphere. As the star finally collapses to become a white dwarf, this envelope is expelled from the star to form a planetary nebula (PN) rich in organic molecules. The physics, dynamics, and chemistry of these nebulae are poorly understood and have implications not only for our understanding of the stellar life cycle but also for organic astrochemistry and the creation of prebiotic molecules in interstellar space. We are working toward generating three-dimensional models of planetary nebulae (PNe), which include the size, orientation, shape, expansion rate and mass distribution of the nebula. Such a reconstruction of a PN is a challenging problem for several reasons. First, the data consist of images obtained over time from the Hubble Space Telescope (HST) and spectra obtained from Kitt Peak National Observatory (KPNO) and Cerro Tololo Inter-American Observatory (CTIO). These images are of course taken from a single viewpoint in space, which amounts to a very challenging tomographic reconstruction. Second, the fact that we have two disparate and orthogonal data types requires that we utilize a method that allows these data to be used together to obtain a solution. To address these first two challenges we employ Bayesian model estimation using a parameterized physical model that incorporates much prior information about the known physics of the PN. In our previous works we have found that the forward problem of the comprehensive model is extremely time consuming. To address this challenge, we explore the use of a set of hierarchical models, which allow us to estimate increasingly more detailed sets of model parameters. These hierarchical models of increasing complexity are akin to scientific theories of increasing sophistication, with each new model/theory being a refinement of a previous one by either incorporating additional prior information or by introducing a new set of parameters to model an entirely new phenomenon. We apply these models to both a simulated and a real ellipsoidal PN to initially estimate the position, angular size, and orientation of the nebula as a two-dimensional object and use these estimates to later examine its three-dimensional properties. The efficiency/accuracy tradeoffs of the techniques are studied to determine the advantages and disadvantages of employing a set of hierarchical models over a single comprehensive model.

Knuth, Kevin H.↗
DOE OSTI2022

Semi-Automated, Object-Based Tomography of Dislocation Structures

The characterization of the three-dimensional arrangement of dislocations is important for many analyses in materials science. Dislocation tomography in transmission electron microscopy is conventionally accomplished through intensity-based reconstruction algorithms. Although such methods work successfully, a disadvantage is that they require many images to be collected over a large tilt range. Here, we present an alternative, semi-automated object-based approach that reduces the data collection requirements by drawing on the prior knowledge that dislocations are line objects. Our approach consists of three steps: (1) initial extraction of dislocation line objects from the individual frames, (2) alignment and matching of these objects across the frames in the tilt series, and (3) tomographic reconstruction to determine the full three-dimensional configuration of the dislocations. Drawing on innovations in graph theory, we employ a node-line segment representation for the dislocation lines and a novel arc-length mapping scheme to relate the dislocations to each other across the images in the tilt series. We demonstrate the method for a dataset collected from a dislocation network imaged by diffraction-contrast scanning transmission electron microscopy. Based on these results and a detailed uncertainty analysis for the algorithm, we discuss opportunities for optimizing data collection and further automating the method.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Focused‐beam X‐ray fluorescence and diffraction microtomographies for mineralogical and chemical characterization of unsectioned extraterrestrial samples

Abstract This study describes the application of new synchrotron X‐ray fluorescence (XRF) and diffraction (XRD) microtomographies for the 3‐D visualization of chemical and mineralogical variations in unsectioned extraterrestrial samples. These improved methods have been applied to three compositionally diverse chondritic meteorite samples that were between 300 and 400 μm in diameter, including samples prepared from fragments of the CR2 chondrite LaPaz Icefield (LAP) 02342, H5 chondrite MacAlpine Hills (MAC) 88203, and the CM2 chondrite Murchison. The synchrotron‐based XRF and XRD tomographies used are focused‐beam techniques that measure the intensities of fluorescent and diffracted X‐rays in a sample simultaneously during irradiation by a high‐energy microfocused incident X‐ray beam. Measured sinograms of the emitted and diffracted intensities were then tomographically reconstructed to generate 2‐D slices of XRF and XRD intensity through the sample, with reconstructed pixel resolution of 1–2 μm, defined by the resolution of the focused incident X‐ray beam. For sample LAP 02342, primary mineral phases that were visualized in reconstructed slices using these techniques included isolated grains of α‐Fe, orthopyroxene, and olivine. For our sample of MAC 88203, XRF/XRD tomography allowed visualization of forsteritic olivine as a primary mineral phase, a vitrified fusion crust at the sample surface, identification of localized Cr‐rich spinels at spatial resolutions of several micrometers, and imaging of a plagioclase‐rich glassy matrix. In the sample of Murchison, major identifiable phases include clinoenstatite‐ and olivine‐rich chondrules, variable serpentine matrix minerals and small Cr‐rich spinels. Most notable in the tomographic analysis of Murchison is the ability to quantitatively distinguish and visualize the complex mixture of serpentine‐group minerals and associated tochilinite–cronstedtite intergrowths. These methods provide new opportunities for spatially resolved characterization of sample texture, mineralogy, crystal structure, and chemical state in unsectioned samples. This provides researchers an ability to characterize such samples internally with minimal disruption of sample micro‐structures and chemistry, possibly without the need for sample extraction from some types of sampling and capture media.

Geochemistry & Geophysics↗
NASA NTRS2003

Three-Dimensional Morphological Analysis of ALH84001 Magnetite Using Electron Tomography

We report here the crystal morphologies of MV-1 and ALH84001 magnetites as calculated by back-projection using electron tomography. In the present study, we used a 300 keV TEM with a field emission gun (Tecnai F-30 from FEI Inc.), equipped with a 2048 x 2048 pixel CCD camera from Gatan Inc. to image magnetite crystals over tilt ranges of approx. +/- 72 deg in 2 deg tilt intervals. The images were aligned for back-projection, either manually, or through the use of fiducial 5 nm Au spheres affixed to the specimen prior to microscopy. Three-dimensional (3-D) reconstructions were computed using weighted back-projection of the tilted views. The tomograms were viewed and analyzed as a series of slices 1.0 nm thick, taken parallel to the specimen-supporting grid, using the IMOD software package. The shape of each magnetite crystal was determined by defining the external contour of a given magnetite in each slice and assembling a stack of these contours in 3-D. To aid in visualization, the stacked contour array was reduced to an optimal mesh by Delaunay triangulation. The surface normal to each of the triangles in the mesh was calculated and the triangle faces colored according to the orientation of that surface normal relative to the principal crystallographic axis of magnetite. Green surfaces correspond to {111} orientations, blue surfaces to {100} orientations, and red surfaces to {110} orientations. Triangles whose surface normal did not correspond to one of the principal axes were colored gray. Within the experimental and numerical uncertainties of the deconvolution, the tomographic reconstruction of both MV-1 and ALH84001 magnetites are equivalent and correspond to a truncated hexa-octahedral morphology.

Thomas-Keprta, Kathie L.↗
NASA NTRS1996

Structural relations between collagen and mineral in bone as determined by high voltage electron microscopic tomography

Aspects of the ultrastructural interaction between collagen and mineral crystals in embryonic chick bone have been examined by the novel technique of high voltage electron microscopic tomography to obtain three-dimensional information concerning extracellular calcification in this tissue. Newly mineralizing osteoid along periosteal surfaces of mid-diaphyseal regions from normal chick tibiae was embedded, cut into 0.25 microns thick sections, and documented at 1.0 MV in the Albany AEI-EM7 high voltage electron microscope. The areas of the tissue studied contained electron dense mineral crystals associated with collagen fibrils, some marked by crystals disposed along their cylindrically shaped lengths. Tomographic reconstructions of one site with two mineralizing fibrils were computed from a 5 degrees tilt series of micrographs over a +/- 60 degrees range. Reconstructions showed that the mineral crystals were platelets of irregular shape. Their sizes were variable, measured here up to 80 x 30 x 8 nm in length, width, and thickness, respectively. The longest crystal dimension, corresponding to the c-axis crystallographically, was generally parallel to the collagen fibril long axis. Individual crystals were oriented parallel to one another in each fibril examined. They were also parallel in the neighboring but apparently spatially separate fibrils. Crystals were periodically (approximately 67 nm repeat distance) arranged along the fibrils and their location appeared to correspond to collagen hole and overlap zones defined by geometrical imaging techniques. The crystals appeared to be continuously distributed along a fibril, their size and number increasing in a tapered fashion from a relatively narrow tip containing smaller and infrequent crystals to wider regions having more densely packed and larger crystals. Defined for the first time by direct visual 3D imaging, these data describe the size, shape, location, orientation, and development of early crystals in normal bone collagen. The results suggest that platelet-shaped crystals are arranged in channels or grooves which are formed by collagen hole zones in register and that crystal sizes may exceed the dimensions of hole zones. Such data agree with those from mineral-matrix interaction in normally calcifying avian tendon obtained by similar high voltage tomographic means, but in addition they indicate a possible gradual and continuous deposition of crystals in collagen of bone unlike tendon and imply that individual collagen fibrils in local regions of osteoid are organized such that they all may be aligned in a coherent manner.

NASA Discipline Number 40-40↗
NASA NTRS2000

Coronal Reconstruction using LASCO and UVCS Observations

The main goal of the research described in the original proposal was to develop methods to quantify coronal and inner-heliospheric velocity fields of the 'quiet' solar wind. For this we planned to use several sources of observations:(1) SOHO/UVCS velocity information in the range 1.5-3 Solar Radii obtained from Doppler dimming observations; (2) projected solar wind velocities (into the plane of the sky) obtained from SOHO/ LASCO images (1.1-30 Solar Radii), primarily derived from two-dimensional correlation tracking techniques; (3) Interplanetary scintillation observations of the heliospheric (> 26 Solar Radii) solar wind velocity from the Solar-Terrestrial Environment Laboratory (STELab) in Nagoya, Japan; (4) Ecliptic in situ observations: data for the ecliptic solar wind are available from the MIT and Los Alamos plasma experiments on the Earth-orbiting IMP-8 spacecraft, from the Cellas instrument on SOHO near the LI Lagrange point, and from the WIND spacecraft; (5) Out-of-ecliptic in situ observations: these data are available primarily from the Los Alamos SWOOPS instrument on Ulysses, which passed over the solar north pole in August 1995, about one year prior to the Whole Sun Month period. Where ever possible we planned to use the first Whole Sun Month as the main time period for the analysis, since we expected that for this period it would be easiest to obtain adequate coverage over the extended period of time required to analyze 'quiet' solar wind patterns. Beyond the observations mentioned above (primarily SOHO data) we extended our selection of data to several events identified in the Yohkoh/SXT data base which directly promised to provide us with clues about the connection between the slow solar wind observed by IPS in the inner heliosphere and their sources in the low corona, in particular active regions. We also obtained valuable results using SWOOPS ill situ observations from the pole-to-pole passage of Ulysses in a comparison with solar wind velocities derived from a tomographic reconstruction of Nagoya IPS observations.

Hick, P. P.↗
DOE OSTI2022

Grazing-Incidence Texture Tomography and Diffuse Reflectivity Tomography of an Organic Semiconductor Device Array

The use of grazing-incidence scattering methods for the characterization of 2D patterned organic thin films is limited due to the elongated 1D footprint of the X-ray beam on the sample. However, this characteristic feature can be turned into an advantage, when combined with tomographic reconstruction. In this pilot study we show, how the use of a chosen texture reflection and a diffuse reflectivity signal can each provide 2D images of the deposits, simultaneously revealing the organic film's crystal orientation and the location of the metal electrodes in a field-effect transistor structure from a single sequence of diffraction images.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Laboratory-Based Micro-X-ray Computed Tomography of Energy Materials at Idaho National Laboratory

Abstract The Idaho National Laboratory (INL) has implemented laboratory-based micro-X-ray computed tomography in a laboratory equipped for the examination of highly radioactive samples. This capability provides nondestructive three-dimensional volumetric information on samples to inform subsequent traditional destructive examinations as well as real-world inputs for high-fidelity scientific modeling. Samples can be imaged with spatial resolutions ranging from several hundred nm/voxel up to ~ 100 µm/voxel. The best usable spatial resolution achieved to date is 384 nm/voxel with this instrument, while the highest radiological dose rate of a sample imaged is ~ 60 R/h β/γ on contact. Advanced data analysis, including custom tomographic reconstruction and segmentation methods, have also been developed to support this capability. In addition to traditional digital X-ray radiography and tomography, this instrument is also able to visualize in situ tensile and compression testing as well as perform diffraction contrast tomography. This work describes the X-ray computed tomography post-irradiation examination capabilities at INL, as well as detailing a variety of applications this instrument has examined.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Demonstration and analysis of volumetric additive manufacturing via sub-orbital spaceflight testing

Computed Axial Lithography (CAL) represents a significant advancement in the emerging field of Volumetric Additive Manufacturing (VAM). CAL addresses key limitations of traditional photopolymer additive manufacturing technologies, by eliminating the need for layering and support structures. Unlike conventional methods, CAL prints components by illuminating all points within a desired geometry simultaneously, using tomographic reconstruction to form the object in a single step. This unique approach eliminates the relative motion between the object and the precursor material, enabling faster printing speeds and reducing the waste associated with support structures. However, CAL parts require post-processing steps before they can be utilized. CAL's core attributes make it particularly suited for In-Space Manufacturing (ISM), due to its fast fabrication times, wide breadth of materials it can use, and minimized footprint. CAL has been successfully demonstrated in microgravity during parabolic flight experiments. However to fully validate and understand CAL's behaviour in microgravity, all manufacturing and post-processing steps must be integrated. In June 2024, we conducted SpaceCAL Mission 3, testing this entire workflow on a suborbital flight aboard Virgin Galactic's SpaceShipTwo. During ~140 s of microgravity, the system autonomously manufactured and post-processed four parts using PEGDA700 resin. Post-flight analysis showed that 2/4 parts were recognisable, while others were distorted due to bubble formation from residual water droplets, off-axis optical aberrations, and non-uniform solvent rinsing. Despite these limitations, this study represents the first integrated CAL workflow in space, providing an initial experimental demonstration and analysis for closed-loop in-space manufacturing.

Waddell, Taylor [Department of Mechanical Engineer↗
DOE OSTI2024

From irregular to regular eutectic growth in the Al-Al 3 Ni system: In situ observations during directional solidification

Here we investigate the irregular eutectic growth dynamics of the Al-Al 3 Ni alloy, in which one of the solid phases (Al 3 Ni) grows faceted from the liquid. Leveraging in situ optical microscopy and synchrotron transmission x-ray microscopy, we address the question of the degree of coupling between Al and Al 3 Ni at the growth front and that of the shape of the microstructures left behind in the bulk solid during directional solidification. Real-time optical observations bring evidence for a morphological transition from a eutectic-grain dependent, irregular eutectic growth at low solidification velocity V (typically 1 μm s -1 ), to a weakly anisotropic, regular growth at higher V ( reaching 10 μm s -1 ). Unprecedented x-ray nano-imaging of the solid–liquid interface, and 3D characterization of the growth patterns, were made possible by a new DS setup at Brookhaven National Laboratory’s NSLS-II. At low V, the leading tips of partly faceted Al 3 Ni crystals are observed to grow not far ahead of the Al growth front. Correlating in situ images and postmortem 3D tomographic reconstructions reveals that the presence of faceted and non-faceted regions of Al 3 Ni crystals in the solid is a direct consequence of coupling and decoupling during DS, respectively. Upon increasing V, the lead distance of Al 3 Ni vanishes, and the shape of Al 3 Ni ceases to be governed by faceted growth. These observations cast light onto the basic mechanisms (faceted growth, diffusive coupling, and the dynamics of trijunctions) governing a faceted to rod-like transition upon increasing V in the Al 3 Ni system, with broad implications to a large class of irregular eutectics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulation toolkit for digital material characterization of large image-based microstructures

In this paper, an efficient image-based simulation toolkit for material characterization is presented, which is scalable to work from personal computers to workstations. The effective thermal conductivity, elasticity, and permeability are evaluated employing a computational homogenization framework based on the Finite Element Method (FEM). Two complementary open-source packages are presented: one developed in Python, which can convert digital images into voxel meshes (pyTomoviewer); the other developed in Julia, that can run numerical simulations to compute effective material properties (chpack). Also, a CUDA C version of chpack is provided (chfem_gpu). They were designed to deal with large multi-phase models, so strategies were devised to minimize their memory footprint, while avoiding a high toll on execution time. The voxel-based approach significantly simplifies the FEM meshes and allows efficient matrix-free implementations. In that sense, to handle large linear systems of equations, the element-by-element (EBE) technique is adopted, in conjunction with a low-memory implementation of the Preconditioned Conjugate Gradient (PCG) method. Finally, the code was thoroughly tested on an artificial geometry made of a square array of cylinders, for which analytical solutions exist, as well as on a real micro-tomographic reconstruction of FiberForm TM , a carbon preform commonly used in thermal protection systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Spray process of multi-component gasoline surrogate fuel under ECN Spray G conditions

As modern gasoline direct injection (GDI) engines utilize sophisticated injection strategies, a detailed understanding of the air-fuel mixing process is crucial to further improvements in engine emission and fuel economy. In this study, a comprehensive evaluation of the spray process of single-component iso-octane (IC8) and multi-component gasoline surrogate E00 (36 % n-pentane, 46 % iso-octane, and 18 % n-undecane, by volume) fuels was conducted using an Engine Combustion Network (ECN) Spray G injector. High-speed extinction, schlieren, and microscopy imaging campaigns were carried out under engine-like ambient conditions in a spray vessel. Experimental results including liquid/vapor penetration, local liquid volume fraction, droplet size, and projected liquid film on the nozzle tip were compared under ECN G1 (573 K, 3.5 kg/m 3 ), G2 (333 K, 0.5 kg/m 3 ), and G3 (333 K, 1.01 kg/m 3 ) conditions. In addition to the experiments, preferential evaporation process of the E00 fuel was elucidated by Large–Eddy Simulations (LES). The three-dimensional liquid volume fraction measurement enabled by the computed tomographic reconstruction showed substantial plume collapse for E00 under the G2 and G3 conditions having wider plume growth and plume-to-plume interaction due to the fuel high vapor pressure. Here, the CFD simulation of E00 showed an inhomogeneity in the way fuel components vaporized, with more volatile components carried downstream in the spray after the end of injection. The high vapor pressure of E00 also results in ~4 μm smaller average droplet diameter than IC8, reflecting a higher rate of initial vaporization even though the final boiling point temperature is higher. Consistent with high vapor pressure, E00 had a wider plume cone angle and enhanced interaction with the wall to cover the entire surface of the nozzle tip in a film. However, the liquid fuel underwent faster evaporation, so the final projected tip wetting area was smaller than the IC8 under the flash-boiling condition.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Effects of detailed geometry and real fluid thermodynamics on Spray G atomization

We present recent results toward the quantification of spray characteristics at engine conditions for an eight-hole counter-bored (stepped) GDI injector-Spray G in the ECN denomination. This computational study is characterized by two novel features: the detailed description of a real injector's internal surfaces via tomographic reconstruction; and a general equation of state that represents the thermodynamic properties of homogeneous liquid-vapor mixtures. The combined level-set moment-of-fluid approach, coupled to an embedded boundary formulation for moving solid walls, makes it possible to seamlessly connect the injector's internal flow to the spray. The Large Eddy Simulation (LES) discussed here presents evidence of partial hydraulic flipping and, during the closing transient, string cavitation. Results are validated by measurements of spray density profiles and droplet size distribution.

42 ENGINEERING↗