Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “image simulation”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Monte Carlo simulation of a novel water-equivalent electronic portal imaging device using plastic scintillating fibers

Purpose: Most electronic portal imaging devices (EPIDs) developed so far use a thin Cu plate/phosphor screen to convert x-ray energies into light photons, while maintaining a high spatial resolution. This results in a low x-ray absorption and thus a low quantum efficiency (QE) of approximately 2-4% for megavoltage (MV) x-rays. A significant increase of QE is desirable for applications such as MV cone-beam computed tomography (MV-CBCT). Furthermore, the Cu plate/phosphor screen contains high atomic number (high-Z) materials, resulting in an undesirable over-response to low energy x-rays (due to photoelectric effect) as well as high energy x-rays (due to pair production) when used for dosimetric verification. Our goal is to develop a new MV x-ray detector that has a high QE and uses low-Z materials to overcome the obstacles faced by current MV x-ray imaging technologies. Methods: A new high QE and low-Z EPID is proposed. It consists of a matrix of plastic scintillating fibers embedded in a water-equivalent medium and coupled to an optically sensitive 2D active matrix flat panel imager (AMFPI) for image readout. It differs from the previous approach that uses segmented crystalline scintillators made of higher density and higher atomic number materials to detect MV x-rays. The plastic scintillating fibers are focused toward the x-ray source to avoid image blurring due to oblique incidence of off-axis x-rays. When MV x-rays interact with the scintillating fibers in the detector, scintillation light will be produced. The light photons produced in a fiber core and emitted within the acceptance angle of the fiber will be guided toward the AMFPI by total internal reflection. A Monte Carlo simulation has been used to investigate imaging and dosimetric characteristics of the proposed detector under irradiation of MV x-rays. Results: Properties, such as detection efficiency, modulation transfer function, detective quantum efficiency (DQE), energy dependence of detector response, and water-equivalence of dose response have been investigated. It has been found that the zero frequency DQE of the proposed detector can be up to 37% at 6 MV. The detector, also, is water-equivalent with a relatively uniform response to different energy x-rays as compared to current EPIDs. Conclusions: The results of our simulations show that, using plastic scintillating fibers, it is possible to construct a water-equivalent EPID that has a better energy response and a higher detection efficiency than current flat panel based EPIDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

POLARIZATION STRUCTURE IN THE CORE OF 1803+784: A SIGNATURE OF RECOLLIMATION SHOCKS?

The extragalactic radio source 1803+784 is a BL Lac object that shows rapid variability and strong linear polarization. Very long baseline interferometry observations reveal that the core possesses a distinctive structure in polarized intensity with two maxima on axis and two minima symmetrically placed on either side. The approximately radial pattern of electric field polarization rods is reminiscent of the results obtained earlier by Cawthorne for conical shocks, but, individually, these do not reproduce the main features of the polarized intensity images. In numerical simulations and experiments, these shocks occur in pairs and help to stabilize jets as they adjust to changes in environment. Here, the polarization resulting from such structures is investigated using an approximate, analytical approach, by making some simple assumptions about the nature of the flow between two such shock waves. For fairly small viewing angles, it is found that a reasonable representation of the core polarization of 1803+784 can be obtained. The similarity between the observed and model polarizations supports the view that the core structure in 1803+784 represents a recollimation shock, and that such shock waves may be responsible for the first disturbance and hence brightening of the quiescent flow in astrophysical jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

TH-CD-206-12: Image-Based Motion Estimation for Plaque Visualization in Coronary Computed Tomography Angiography

Purpose: Visualization and quantification of coronary artery calcification and atherosclerotic plaque benefits from coronary artery motion (CAM) artifact elimination. This work applies a rigid linear motion model to a Volume of Interest (VoI) for estimating motion estimation and compensation of image degradation in Coronary Computed Tomography Angiography (CCTA). Methods: In both simulation and testbench experiments, translational CAM was generated by displacement of the imaging object (i.e. simulated coronary artery and explanted human heart) by ∼8 mm, approximating the motion of a main coronary branch. Rotation was assumed to be negligible. A motion degraded region containing a calcification was selected as the VoI. Local residual motion was assumed to be rigid and linear over the acquisition window, simulating motion observed during diastasis. The (negative) magnitude of the image gradient of the reconstructed VoI was chosen as the motion estimation objective and was minimized with Covariance Matrix Adaptation Evolution Strategy (CMAES). Results: Reconstruction incorporated the estimated CAM yielded signification recovery of fine calcification structures as well as reduced motion artifacts within the selected local region. The compensated reconstruction was further evaluated using two image similarity metrics, the structural similarity index (SSIM) and Root Mean Square Error (RMSE). At the calcification site, the compensated data achieved a 3% increase in SSIM and a 91.2% decrease in RMSE in comparison with the uncompensated reconstruction. Conclusion: Results demonstrate the feasibility of our image-based motion estimation method exploiting a local rigid linear model for CAM compensation. The method shows promising preliminary results for the application of such estimation in CCTA. Further work will involve motion estimation of complex motion corrupted patient data acquired from clinical CT scanner.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Invariant discovery of features across multiple length scales: Applications in microscopy and autonomous materials characterization

Physical imaging is a foundational characterization method in areas from condensed matter physics and chemistry to astronomy and spans length scales from atomic to universe. Images encapsulate crucial data regarding atomic bonding, materials microstructures, and dynamic phenomena such as microstructural evolution and turbulence, among other phenomena. The challenge lies in effectively extracting and interpreting this information. Variational Autoencoders (VAEs) have emerged as powerful tools for identifying the underlying factors of variation in image data, providing a systematic approach to distilling meaningful patterns from complex data sets. However, a significant hurdle in their application is the definition and selection of appropriate descriptors reflecting local structures. Here, we introduce the scale-invariant VAE approach (SI-VAE) based on the progressive training of the VAE with the descriptors sampled at different length scales. The SI-VAE allows the discovery of the length scale-dependent factors of variation in the system. Here, we illustrate this approach using the ferroelectric domain images and generalize it to the movies of the electron-beam induced phenomena in graphene and topography evolution across combinatorial libraries. This approach can further be used to initialize the decision making in automated experiments including structure–property discovery and can be applied across a broad range of imaging methods. This approach is universal and can be applied to any spatially resolved data including both experimental imaging studies and simulations, and can be particularly useful for exploration of phenomena such as turbulence and scale-invariant transformation fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Forming an Edge Dislocation Image at Anomalous X-ray Transmission

The features of forming section and projection topography images of an edge dislocation normal to the crystal surface are considered for the case of X-ray anomalous transmission. Experimental images are analyzed using numerical simulation of the diffraction experiment. A new mechanism for image formation of defects arranged near the output crystal surface has been proposed. Differences between topography images of dislocations and rosettes of local misorientations of reflecting planes have been established.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Towards an improved understanding of the influence of subpixel vegetation structure on pixel-level spectra: a simulation approach

The planned NASA Hyperspectral Infrared Imager (HyspIRI) mission, equipped with an imaging spectrometer that has the capability of monitoring ecosystems globally, will provide an unprecedented opportunity to address scientific challenges related to ecosystem function and change. However, uncertainty remains around the impact of subpixel vegetation structure, in combination with the point spread function, on pixel-level imaging spectroscopy data. We estimated structural parameters, e.g., leaf area index (LAI), canopy cover, and tree location, from HyspIRI spectral data, with the goal of assessing how subpixel variation in these parameters impact pixel-level imaging spectroscopy data. The fine-scale variability of real vegetation structure makes this a challenging endeavor. Therefore, we utilized a simulation-based approach to counter the time-consuming and often destructive sampling needs of vegetation structural analysis and to simultaneously generate synthetic HyspIRI data pre-launch. Three virtual scenes were constructed, corresponding to the actual vegetation structure in the National Ecological Observatory Network’s (NEON) Pacific Southwest Domain (Fresno, CA). These included an oak savanna, a dense coniferous forest, and a conifer-manzanita-mixed forest. Simulated spectroscopy data for these scenes were then generated using the Digital Imaging and Remote Sensing Image Generation (DIRSIG) model. Simulations first were used to verify the physical model, virtual scene geometrical information, and simulation parameters. This was followed by simulations of HyspIRI data, where within-pixel structural variability was introduced, e.g., by iteratively changing per-pixel canopy cover and tree placement, tree clustering, leaf area index (LAI), etc., between simulation runs for the virtual scenes. Finally, narrow-band vegetation indices (VIs) were extracted from the data in an attempt to describe the variability of the subpixel structural parameters; this was done in order to assess VI robustness to changes in structural “levels”, as well as placement of trees/canopies within the instrument’s instantaneous field-of-view (IFOV). Our ultimate goal is not only to better understand how such subpixel variability influence imaging spectroscopy outputs, but also to better estimate vegetation structural parameters using spectra. We constructed regression models for LAI (R2 = 0.92) and canopy cover (R2 = 0.97) with narrow-band VIs via this simulation approach. Our models ultimately are intended to improve the HyspIRI mission’s ability to monitor global vegetation structure.

Van Aardt, Jan↗
NASA NTRSDate not supplied

Simulations of OSIRIS-REx Atmospheric Entry for Comparisons with Data from the Imaging Campaign

Results from pre-flight flow and radiation simulations for the OSIRIS-REx Sample Return Capsule along a candidate entry trajectory are presented. These results are from pre-flight simulations in support of the observation campaign. Details of the methodology for computing temporal variation of spectral irradiance on an instrument onboard an aerial observation platform are presented along with a sampling of results. An attempt is made to use the methodology presented to assess the level of elemental sodium in the shock layer. The final paper will be based on the actual flight trajectory and will include comparisons of predictions with calibrated spectrometric data.

Spectrometric Imaging↗
DOE OSTI2015

Compressed Streak Imaging Beam Diagnostic

We are planning a new beam diagnostic based on Compressed Ultrafast Photography (CUP)1.A foil inserted in the beam path is used to generate a continuous optical image of the beam, which is the same technique used on the existing DARHT II beam imaging diagnostic. This existing diagnostic compresses the beam image with anamorphic lenses into four views that are streaked simultaneously on two streak cameras. In our new beam diagnostic design the anamorphic lenses are replaced with regular lenses and the beam image is stamped with a pseudo-random mask pattern. The full image is imaged onto a streak camera with its entrance slit expanded. Modifying the CUP technique, which uses a single streaked image, our design splits the image into four rotated copies and puts all four images on two streak cameras. A data cube of multiple image frames is reconstructed from the streak camera data through the use of the TwiST algorithm2 combined with total variation denoising3. The additional rotated images improve spatial resolution and reduce noise and image artifacts compared to using a single streaked image. In reconstructions of simulated data, fine detail in the beam profile can be seen and there is a remarkable absence of image artifacts compared to the existing DARHT II beam imaging diagnostic. Performance was evaluated using a pseudo-MTF derived from a simulated wave pattern test object.

Marks, Daniel G.↗
DOE OSTI2015

Mapping Patterns of Ipsilateral Supraclavicular Nodal Metastases in Breast Cancer: Rethinking the Clinical Target Volume for High-risk Patients

Purpose: To map the location of metastatic supraclavicular (SCV) lymph nodes (LNMs) in breast cancer patients with SCV node involvement and determine whether and where the radiation therapy clinical target volume (CTV) of this region could be modified in high-risk subsets. Methods and Materials: Fifty-five patients with metastatic SCV LNMs were eligible for geographic mapping and atlas coverage analysis. All LNMs and their epicenters were registered proportionally by referencing the surrounding landmarks onto simulation computed tomography images of a standard patient. CTVs based on selected SCV atlases, including the one by the Radiation Therapy Oncology Group (RTOG) were contoured. A modified SCV CTV was tried and shown to have better involved-node coverage and thus theoretically improved prophylaxis in this setting. Results: A total of 50 (91%) and 45 (81.8%) patients had LNMs in the medial and lateral SCV subregions, respectively. Also, 36 patients (65.5%) had LNMs located at the junction of the jugular-subclavian veins. All nodes were covered in only 25.5% to 41.8% of patients by different atlases. The RTOG atlas covered all nodes in 25.5% of patients. Stratified by the nodes in all the patients as a whole, 49.2% to 81.3% were covered, and the RTOG atlas covered 62.6%. The lateral and posterior borders were the most overlooked locations. Modification by extending the borders to natural anatomic barriers allowed the new CTV to cover all the nodes in 81.8% of patients and encompass 96.1% of all the nodes. Conclusions: According to the distribution of SCV LNMs, the extent of existing atlases might not be adequate for potential metastatic sites in certain groups of patients. The extension of the lateral and posterior CTV borders in high-risk or recurrent patients might be a reasonable approach for increasing coverage. However, additional data in more homogeneous populations with localized disease are needed before routine application.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Star formation rate and extinction in faint z ∼ 4 Lyman break galaxies

We present a statistical detection of 1.5 GHz radio continuum emission from a sample of faint z ∼ 4 Lyman break galaxies (LBGs). To constrain their extinction and intrinsic star formation rate (SFR), we combine the latest ultradeep Very Large Array 1.5 GHz radio image and the Hubble Space Telescope Advanced Camera for Surveys (ACS) optical images in the GOODS-N. We select a large sample of 1771 z ∼ 4 LBGs from the ACS catalog using B {sub F435W}-dropout color criteria. Our LBG samples have I {sub F775W} ∼ 25-28 (AB), ∼0-3 mag fainter than M{sub UV}{sup ⋆} at z ∼ 4. In our stacked radio images, we find the LBGs to be point-like under our 2'' angular resolution. We measure their mean 1.5 GHz flux by stacking the measurements on the individual objects. We achieve a statistical detection of S {sub 1.5} {sub GHz} = 0.210 ± 0.075 μJy at ∼3σ for the first time on such a faint LBG population at z ∼ 4. The measurement takes into account the effects of source size and blending of multiple objects. The detection is visually confirmed by stacking the radio images of the LBGs, and the uncertainty is quantified with Monte Carlo simulations on the radio image. The stacked radio flux corresponds to an obscured SFR of 16.0 ± 5.7 M {sub ☉} yr{sup –1}, and implies a rest-frame UV extinction correction factor of 3.8. This extinction correction is in excellent agreement with that derived from the observed UV continuum spectral slope, using the local calibration of Meurer et al. This result supports the use of the local calibration on high-redshift LBGs to derive the extinction correction and SFR, and also disfavors a steep reddening curve such as that of the Small Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy

Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.

electron inelastic mean free path↗
DOE OSTI2019

An asymmetric hybrid cryptosystem using equal modulus and random decomposition in hybrid transform domain

In this paper, an asymmetric hybrid cryptosystem with coherent superposition, equal modulus and random decomposition in hybrid transform domain is proposed. To further strengthen the security of the cryptosystem, a hyperchaotic system is used as a pixel-swapping procedure. The hybrid transform is generated by utilizing fractional Fourier transform of various orders and Walsh transform. The hyperchaotic framework’s parameters and starting conditions alongside the fractional orders of the fractional Fourier transform extend the key-space and consequently give extra strength to the proposed cryptosystem. The designed cryptosystem has an extended key-space to avoid any brute-force attack and is non-linear in nature. The scheme is validated on gray-scale images. Computer based simulations have been performed to verify the validity and the performance of the proposed cryptosystem against different types of attacks. Results demonstrate that the proposed cryptosystem not only offers higher protection against noise attacks but is also invulnerable to special attack.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Multiphase Model of Semisolid Slurry Generation and Isothermal Holding During Cooling Slope Rheoprocessing of A356 Al Alloy

In the present paper, we present an experimentally validated 3D multiphase and multiscale solidification model to understand the transport processes involved during slurry generation with a cooling slope. In this process, superheated liquid alloy is poured at the top of the cooling slope and allowed to flow along the slope under the influence of gravity. As the melt flows down the slope, it progressively loses its superheat, starts solidifying at the melt/slope interface with formation of solid crystals, and eventually exits the slope as semisolid slurry. In the present simulation, the three phases considered are the parent melt as the primary phase, and the solid grains and air as secondary phases. The air phase forms a definable air/liquid melt interface as the free surface. After exiting the slope, the slurry fills an isothermal holding bath maintained at the slope exit temperature, which promotes further globularization of microstructure. The outcomes of the present model include prediction of volume fractions of the three different phases considered, grain evolution, grain growth, size, sphericity and distribution of solid grains, temperature field, velocity field, macrosegregation and microsegregation. In addition, the model is found to be capable of making predictions of morphological evolution of primary grains at the onset of isothermal coarsening. The results obtained from the present simulations are validated by performing quantitative image analysis of micrographs of the rapidly oil-quenched semisolid slurry samples, collected from strategic locations along the slope and from the isothermal slurry holding bath.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Potts Model parameter study of particle size, Monte Carlo temperature, and “Particle-Assisted Abnormal Grain Growth”

A Potts Model was proposed to account for temperature and particle heterogeneity-dependent Zener Pinning. This was accomplished by relating the random fluctuations of grain boundary position allowed by the Potts Model switching probability at finite simulation temperature to experimentally observed grain growth stagnation behavior. We assume that these fluctuations arise from random fluctuations in the thermodynamic and kinetic properties of the grain boundaries and/or interfaces as they change with temperature. As an application of this model, the grain growth kinetics of U-10 wt. % Mo nuclear fuels were simulated with the input of microstructural images during different heat treatment processes. Simulated average grain growth behavior is in good agreement with experiments.

Frazier, William E.↗
DOE OSTI2024

Uncertainty quantification and propagation in lithium-ion battery electrodes using bayesian convolutional neural networks

The complex nature of manufacturing processes stipulates electrodes to possess high variability with increased heterogeneity during production. X-ray computed tomography imaging has proved to be critical in visualizing the complicated stochastic particle distribution of as-manufactured electrodes in lithium-ion batteries. However, accurate prediction of their electrochemical performance necessitates precise evaluation of kinetic and transport properties from real electrodes. Image segmentation that characterizes voxels to particle/pore phase is often meticulous and fraught with subjectivity owing to a myriad of unconstrained choices and filter algorithms. Here we utilize a Bayesian convolutional neural network to tackle segmentation subjectivity and quantify its pertinent uncertainties. Otsu inter-variance and Blind/Referenceless Imaging Spatial Quality Evaluator are used to assess the relative image quality of grayscale tomograms, thus evaluating the uncertainty in the derived microstructural attributes. We analyze how image uncertainty is correlated with the uncertainties and magnitude of kinetic and transport properties of an electrode, further identifying pathways of uncertainty propagation within microstructural attributes. The coupled effect of spatial heterogeneity and microstructural anisotropy on the uncertainty quantification of transport parameters is also understood. This work demonstrates a novel methodology to extract microstructural descriptors from real electrode images through quantification of associated uncertainties and discerning the relative strength of their propagation, thus facilitating feedback to manufacturing processes from accurate image based electrochemical simulations.

25 ENERGY STORAGE↗
DOE OSTI2015

Modeling Pancreatic Tumor Motion Using 4-Dimensional Computed Tomography and Surrogate Markers

Purpose: To assess intrafractional positional variations of pancreatic tumors using 4-dimensional computed tomography (4D-CT), their impact on gross tumor volume (GTV) coverage, the reliability of biliary stent, fiducial seeds, and the real-time position management (RPM) external marker as tumor surrogates for setup of respiratory gated treatment, and to build a correlative model of tumor motion. Methods and Materials: We analyzed the respiration-correlated 4D-CT images acquired during simulation of 36 patients with either a biliary stent (n=16) or implanted fiducials (n=20) who were treated with RPM respiratory gated intensity modulated radiation therapy for locally advanced pancreatic cancer. Respiratory displacement relative to end-exhalation was measured for the GTV, the biliary stent, or fiducial seeds, and the RPM marker. The results were compared between the full respiratory cycle and the gating interval. Linear mixed model was used to assess the correlation of GTV motion with the potential surrogate markers. Results: The average ± SD GTV excursions were 0.3 ± 0.2 cm in the left-right direction, 0.6 ± 0.3 cm in the anterior-posterior direction, and 1.3 ± 0.7 cm in the superior-inferior direction. Gating around end-exhalation reduced GTV motion by 46% to 60%. D95% was at least the prescribed 56 Gy in 76% of patients. GTV displacement was associated with the RPM marker, the biliary stent, and the fiducial seeds. The correlation was better with fiducial seeds and with biliary stent. Conclusions: Respiratory gating reduced the margin necessary for radiation therapy for pancreatic tumors. GTV motion was well correlated with biliary stent or fiducial seed displacements, validating their use as surrogates for daily assessment of GTV position during treatment. A patient-specific internal target volume based on 4D-CT is recommended both for gated and not-gated treatment; otherwise, our model can be used to predict the degree of GTV motion.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2017

Toward Personalized Dosimetry with {sup 32}P Microparticle Therapy for Advanced Pancreatic Cancer

Purpose: To develop a Monte Carlo model for patient-specific dosimetry of {sup 32}P microparticle localized internal radionuclide therapy for advanced pancreatic cancer. Methods and Materials: Spherical tumor geometries and a pancreatic phantom were modeled, as well as different 3-dimensional non-uniform clinical pancreatic geometries based on patient-specific ultrasound images. The dosimetry simulations modeled the dose distribution due to the energy spectrum of emitted beta particles. Results: The average dose for small (3-cm diameter) and large (6-cm diameter) spherical tumors was 111 Gy (for 7.6 MBq administered activity) and 128 Gy (for 58 MBq), respectively. For the clinical 3-dimensional geometries, on the basis of patient data, the mean doses delivered to the tumor were calculated to be in the range 102 to 113 Gy, with negligible dose to the pancreas for the smallest tumor volumes. The calculated dose distributions are highly non-uniform. For the largest tumor studied, the pancreas received approximately 6% of the tumor dose (5.7 Gy). Importantly, we found that because the smallest tumor studied exhibited the most dynamic changes in volume in response to the treatment, the dose to tumor and pancreas is significantly underestimated if a static tumor volume is assumed. Conclusions: These results demonstrate the dosimetry of {sup 32}P microparticle localized internal radionuclide therapy for pancreatic cancer and the possibility of developing personalized treatment strategies. The results also highlight the importance of considering the effects of non-uniform dose distributions and dynamic change of tumor mass during treatment on the dosimetry of the tumor and critical organs.

62 RADIOLOGY AND NUCLEAR MEDICINE↗