Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Image Processing”

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

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

At least 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

A Fast Implementation of the ISOCLUS Algorithm

Unsupervised clustering is a fundamental tool in numerous image processing and remote sensing applications. For example, unsupervised clustering is often used to obtain vegetation maps of an area of interest. This approach is useful when reliable training data are either scarce or expensive, and when relatively little a priori information about the data is available. Unsupervised clustering methods play a significant role in the pursuit of unsupervised classification. One of the most popular and widely used clustering schemes for remote sensing applications is the ISOCLUS algorithm, which is based on the ISODATA method. The algorithm is given a set of n data points (or samples) in d-dimensional space, an integer k indicating the initial number of clusters, and a number of additional parameters. The general goal is to compute a set of cluster centers in d-space. Although there is no specific optimization criterion, the algorithm is similar in spirit to the well known k-means clustering method in which the objective is to minimize the average squared distance of each point to its nearest center, called the average distortion. One significant feature of ISOCLUS over k-means is that clusters may be merged or split, and so the final number of clusters may be different from the number k supplied as part of the input. This algorithm will be described in later in this paper. The ISOCLUS algorithm can run very slowly, particularly on large data sets. Given its wide use in remote sensing, its efficient computation is an important goal. We have developed a fast implementation of the ISOCLUS algorithm. Our improvement is based on a recent acceleration to the k-means algorithm, the filtering algorithm, by Kanungo et al.. They showed that, by storing the data in a kd-tree, it was possible to significantly reduce the running time of k-means. We have adapted this method for the ISOCLUS algorithm. For technical reasons, which are explained later, it is necessary to make a minor modification to the ISOCLUS specification. We provide empirical evidence, on both synthetic and Landsat image data sets, that our algorithm's performance is essentially the same as that of ISOCLUS, but with significantly lower running times. We show that our algorithm runs from 3 to 30 times faster than a straightforward implementation of ISOCLUS. Our adaptation of the filtering algorithm involves the efficient computation of a number of cluster statistics that are needed for ISOCLUS, but not for k-means.

Memarsadeghi, Nargess↗
DOE OSTI2022

Quantitative spatiotemporal mapping of thermal runaway propagation rates in lithium-ion cells using cross-correlated Gabor filtering

Abuse testing of lithium-ion batteries is widely performed in order to develop new safety standards and strategies. However, testing methodologies are not standardised across the research community, especially with failure mechanisms being inherently difficult to reproduce. High-speed X-ray radiography is proven to be a valuable tool to capture events occurring during cell failure, but the observations made remain largely qualitative. We have therefore developed a robust image processing toolbox that can quantify, for the first time, the rate of propagation of battery failure mechanisms revealed by high-speed X-ray radiography. Using Gabor filter, the toolbox selectively tracks the electrode structure at the onset of failure. This facilitated the estimation of the displacement of electrodes undergoing abuse via nail penetration, and also the tracking of objects, such as the nail, as it propagates through a cell. Further, by cross-correlating the Gabor signals, we have produced practical, illustrative spatiotemporal maps of the failure events. From these, we can quantify the propagation rates of electrode displacement prior to the onset of thermal runaway. The highest recorded acceleration (≈514 mm s -2 ) was when a nail penetrated a cell radially (perpendicular to the electrodes) as opposed to axially (parallel to the electrodes). The initiation of thermal runaway was also resolved in combination with electrode displacement, which occurred at a lower acceleration (≈108 mm s -2 ). Our assistive toolbox can also be used to study other types of failure mechanisms, extracting otherwise unattainable kinetic data. Ultimately, this tool can be used to not only validate existing theoretical mechanical models, but also standardise battery failure testing procedures.

25 ENERGY STORAGE↗
NASA NTRS2013

A Combined Very Large Telescope and Gemini Study of the Atmosphere of the Directly Imaged Planet, Beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet Beta Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While Beta Pic b's near-IR colors mimic those of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects-i.e., κ And b and 1RXJ 1609B-match Beta Pic b's JHKsL photometry and its 3.1 micron and 5 micron photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large (approx. 60 micron)dust grains fail to reproduce the Beta Pic b spectrum. However, models incorporating thick clouds similar to those found forHR8799 bcde, but also with small (a fewmicrons) modal particle sizes, yield fits consistent with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles (a = 4 micron), we derive atmosphere parameters of log(g) = 3.8 +/- 0.2 and Teff = 1575-1650 K, an inferred mass of 7+4 −3 MJ, and a luminosity of log(L/L) approx. −3.80 +/- 0.02. The best-estimated planet radius, is approx. equal to 1.65 +/- 0.06 RJ, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if Beta Pic b is younger than is approx. equal to 7 Myr, consistent with a late formation well after its host star's birth approx. 12+8 −4 Myr ago.

overluminous↗
DOE OSTI2015

Effects of pulse duration on magnetostimulation thresholds

Purpose: Medical imaging techniques such as magnetic resonance imaging and magnetic particle imaging (MPI) utilize time-varying magnetic fields that are subject to magnetostimulation limits, which often limit the speed of the imaging process. Various human-subject experiments have studied the amplitude and frequency dependence of these thresholds for gradient or homogeneous magnetic fields. Another contributing factor was shown to be number of cycles in a magnetic pulse, where the thresholds decreased with longer pulses. The latter result was demonstrated on two subjects only, at a single frequency of 1.27 kHz. Hence, whether the observed effect was due to the number of cycles or due to the pulse duration was not specified. In addition, a gradient-type field was utilized; hence, whether the same phenomenon applies to homogeneous magnetic fields remained unknown. Here, the authors investigate the pulse duration dependence of magnetostimulation limits for a 20-fold range of frequencies using homogeneous magnetic fields, such as the ones used for the drive field in MPI. Methods: Magnetostimulation thresholds were measured in the arms of six healthy subjects (age: 27 ± 5 yr). Each experiment comprised testing the thresholds at eight different pulse durations between 2 and 125 ms at a single frequency, which took approximately 30–40 min/subject. A total of 34 experiments were performed at three different frequencies: 1.2, 5.7, and 25.5 kHz. A solenoid coil providing homogeneous magnetic field was used to induce stimulation, and the field amplitude was measured in real time. A pre-emphasis based pulse shaping method was employed to accurately control the pulse durations. Subjects reported stimulation via a mouse click whenever they felt a twitching/tingling sensation. A sigmoid function was fitted to the subject responses to find the threshold at a specific frequency and duration, and the whole procedure was repeated at all relevant frequencies and pulse durations. Results: The magnetostimulation limits decreased with increasing pulse duration (T{sub pulse}). For T{sub pulse} < 18 ms, the thresholds were significantly higher than at the longest pulse durations (p < 0.01, paired Wilcoxon signed-rank test). The normalized magnetostimulation threshold (B{sub Norm}) vs duration curve at all three frequencies agreed almost identically, indicating that the observed effect is independent of the operating frequency. At the shortest pulse duration (T{sub pulse} ≈ 2 ms), the thresholds were approximately 24% higher than at the asymptotes. The thresholds decreased to within 4% of their asymptotic values for T{sub pulse} > 20 ms. These trends were well characterized (R{sup 2} = 0.78) by a stretched exponential function given by B{sub Norm}=1+αe{sup −(T{sub p}{sub u}{sub l}{sub s}{sub e}/β){sup γ}}, where the fitted parameters were α = 0.44, β = 4.32, and γ = 0.60. Conclusions: This work shows for the first time that the magnetostimulation thresholds decrease with increasing pulse duration, and that this effect is independent of the operating frequency. Normalized threshold vs duration trends are almost identical for a 20-fold range of frequencies: the thresholds are significantly higher at short pulse durations and settle to within 4% of their asymptotic values for durations longer than 20 ms. These results emphasize the importance of matching the human-subject experiments to the imaging conditions of a particular setup. Knowing the dependence of the safety limits to all contributing factors is critical for increasing the time-efficiency of imaging systems that utilize time-varying magnetic fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Molecular Simulations of the Chain Length Dependent Adsorption of C7‐C14 n ‐Alkanes in ZIF‐8

Recent experiments show that the diffusivities of C8 and C10 n-alkanes in ZIF-8 are higher than those of C7 and C9, respectively. We investigated this unusual ‘odd-even' effect by simulating the adsorption of C7-C14 n-alkanes in ZIF-8 using hybrid Monte Carlo molecular simulations. The resultant adsorption isotherms, guest-host energies, isosteric heats, and chain length distributions are analyzed for trends among the n-alkanes. ZIF-8 cages filled with n-alkanes are characterized using a combination of image processing and data science to quantify the differences in packing that occur with chain length. In conclusion, results indicate that packing changes drastically from C7 to C8 and from C9 to C10, which is consistent with the diffusion trends previously reported in experiments.

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

Fractal Scaling of Explosively Driven Product Gases

ABSTRACT Characterization of the interface between explosive product gases and ambient air in an explosion is a complicated task due to the turbulent mixing and inherently three‐dimensional expansion of the interface. This study aims to quantify the evolution of the interface as a temporally varying Hausdorff dimension. Two test series were conducted with Composition C‐4 charges with masses of 105 and 880 g. Imaging data were collected from the time of detonation until shock wave detachment using ultra‐high‐speed cameras. Gas cloud profiles were extracted using automated image processing algorithms, and the Hausdorff dimension of these two‐dimensional slices of the gas cloud was then estimated using boxcounting algorithms. When scaled with standard gas dynamic nondimensional scalings, the Hausdorff dimension of all explosive events appears to collapse towards a single curve. The fireball was initially nonfractal and began to develop fractal properties as the shock wave separated from the detonation products. Artificial perturbation of the charged surface had no detectable impact on the evolution of the Hausdorff dimension in the early development of the fireball outside of error, despite visible phenomenological differences in the early development of mixing on the fireball surface.

42 ENGINEERING↗
DOE OSTI2025

Innovative Method for Reliable Measurement of PEM Water Electrolyzer Component Resistances

Understanding the sheet resistance of porous electrodes is essential for improving the performance of polymer electrolyte membrane (PEM) water electrolyzers and related technologies. Despite its importance, existing methods often fail to provide reliable and comprehensive data, especially for porous materials with complex morphologies and non‐uniform thicknesses. This study introduces a robust and straightforward method for determining the sheet resistance of porous electrodes using a novel probe concept based on industrial printed circuit board (PCB) technology. This probe measures resistance across ten distances, ranging from 250 µm to 2500 µm, enabling local mapping of resistance. The study focuses on the sheet resistance of key components in PEM water electrolyzers, including the gas diffusion layer (GDL), porous transport layer (PTL), and catalyst layers deposited on a membrane. Additionally, an image‐processing‐based method is presented to obtain the thickness distribution of the studied catalyst layers, facilitating a detailed analysis of the electrical in‐plane resistivity with thickness variations. Overall, this methodology has the potential to expedite material integration and bridge the gap between electrode engineering and single‐cell testing, thereby advancing the development of PEM water electrolyzers.

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

Root identification in minirhizotron imagery with multiple instance learning

In this study, multiple instance learning (MIL) algorithms to automatically perform root detection and segmentation in minirhizotron imagery using only image-level labels are proposed. Root and soil characteristics vary from location to location, and thus, supervised machine learning approaches that are trained with local data provide the best ability to identify and segment roots in minirhizotron imagery. However, labeling roots for training data (or otherwise) is an extremely tedious and time-consuming task. This paper aims to address this problem by labeling data at the image level (rather than the individual root or root pixel level) and train algorithms to perform individual root pixel level segmentation using MIL strategies. Three MIL methods (multiple instance adaptive cosine coherence estimator, multiple instance support vector machine, multiple instance learning with randomized trees) were applied to root detection and compared to non-MIL approaches. The results show that MIL methods improve root segmentation in challenging minirhizotron imagery and reduce the labeling burden. In our results, multiple instance support vector machine outperformed other methods. The multiple instance adaptive cosine coherence estimator algorithm was a close second with an added advantage that it learned an interpretable root signature which identified the traits used to distinguish roots from soil and did not require parameter selection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Toplogical optimization of structures using Fourier representations

The minimization of compliance subject to a mass constraint is the topology optimization design problem of interest. The goal is to determine the optimal configuration of material within an allowed volume. Our approach builds upon the well-known density method in which the decision variable is the material density in every cell in a mesh. In it’s most basic form the density method consists of three steps: 1) the problem is convexified by replacing the integer material indicator function with a volume fraction, 2) the problem is regularized by filtering the volume fraction field to impose a minimum length scale; 3) the filtered volume fraction is penalized to steer the material distribution toward binary designs. The filtering step is used to yield a mesh-independent solution and to eliminate checkerboard instabilities. In image processing terms this is a low-pass filter, and a consequence is that the decision variables are not independent and a change of basis could significantly reduce the dimension of the nonlinear programming problem. Based on this observation, we represent the volume fraction field with a truncated Fourier representation. Furthermore, this imposes a minimal length scale on the problem, eliminates checkerboard instabilities, and also reduces the number of decision variables by over 100 × (two dimensions) or 1000 × (three dimensions).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Are we at a crossroads or a plateau? Radiomics and machine learning in abdominal oncology imaging

Advances in radiomics and machine learning have driven a technology boom in the automated analysis of radiology images. For the past several years, expectations have been nearly boundless for these new technologies to revolutionize radiology image analysis and interpretation. In this editorial, I compare the expectations with the realities with particular attention to applications in abdominal oncology imaging. I explore whether these technologies will leave us at a crossroads to an exciting future or to a sustained plateau and disillusionment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Body composition changes after left gastric artery embolization in overweight and obese individuals

Purpose: To determine the effects of left gastric artery embolization (LGAE) on computed tomography (CT) body composition change. Materials and methods: Sixteen overweight or obese patients who had abdominal CT scans before and after LGAE for gastric bleeding were retrospectively reviewed. Body composition analysis was performed with semiautomated imaging processing algorithms (MATLAB 13.0, Math Works, MA). Adipose tissue and lean skeletal muscle were measured using threshold attenuation values. Total body fat index (BFI), subcutaneous fat index (SFI), visceral fat index (VFI), intramuscular fat index (IMFI), and skeletal muscle index (SMI) were determined ([tissue area (cm)]{sup 2}/[height (m)]{sup 2}). Excess body weight (EBW) was determined based on the Lorentz formula for ideal body weight. Results: Mean follow-up was 1.5 ± 0.8 months. Following LGAE, patients experienced significantly decreased body weight (p = 0.003), BMI (p = 0.005), EBW (p = 0.003), BFI (p = 0.03), SFI (p = 0.03), and SMI (p < 0.001). Changes in VFI and IMFI did not significantly change (p = 0.13 and p = 0.83, respectively). Conclusions: Patients who underwent LGAE had significant unintended weight loss as a result of decreased body fat and skeletal muscle. Body composition analysis can readily assess the extent of fat loss and identify muscle wasting.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Saline as the Sole Contrast Agent for Successful MRI-guided Epidural Injections

Purpose. To assess the performance of sterile saline solution as the sole contrast agent for percutaneous magnetic resonance imaging (MRI)-guided epidural injections at 1.5 T. Methods. A retrospective analysis of two different techniques of MRI-guided epidural injections was performed with either gadolinium-enhanced saline solution or sterile saline solution for documentation of the epidural location of the needle tip. T1-weighted spoiled gradient echo (FLASH) images or T2-weighted single-shot turbo spin echo (HASTE) images visualized the test injectants. Methods were compared by technical success rate, image quality, table time, and rate of complications. Results. 105 MRI-guided epidural injections (12 of 105 with gadolinium-enhanced saline solution and 93 of 105 with sterile saline solution) were performed successfully and without complications. Visualization of sterile saline solution and gadolinium-enhanced saline solution was sufficient, good, or excellent in all 105 interventions. For either test injectant, quantitative image analysis demonstrated comparable high contrast-to-noise ratios of test injectants to adjacent body substances with reliable statistical significance levels (p < 0.001). The mean table time was 22 {+-} 9 min in the gadolinium-enhanced saline solution group and 22 {+-} 8 min in the saline solution group (p = 0.75). Conclusion. Sterile saline is suitable as the sole contrast agent for successful and safe percutaneous MRI-guided epidural drug delivery at 1.5 T.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

Three-Dimensional Rotational Angiography of the Foot in Critical Limb Ischemia: A New Dimension in Revascularization Strategy

Purpose. To evaluate the additional value of three-dimensional rotational angiography (3DRA) of the foot compared with digital subtraction angiography (DSA) in patients with critical limb ischemia (CLI). Technique. For 3DRA, the C-arm was placed in the propeller position with the foot in an isocentric position. The patient's unaffected foot was positioned in a footrest outside the field of view. For correct timing of 3DRA, the delay from contrast injection in the popliteal artery at the level of knee joint to complete pedal arterial enhancement was assessed using DSA. With this delay, 3DRA was started after injection of 15 ml contrast. Imaging of the 3DRA could directly be reconstructed and visualized.Materials and MethodsPatients undergoing 3DRA of the foot were prospectively registered. DSA and 3DRA images were scored separately for arterial patency and presence of collaterals. Treatment strategies were proposed based on DSA with and without the availability of 3DRA. Results. Eleven patients underwent 3DRA of the foot. One 3DRA was not included because the acquisition was focused on the heel instead of the entire foot. Diagnostic quality of 3DRA was good in all ten patients. 3DRA compared with DSA showed additional patent arteries in six patients, patent plantar arch in three patients, and collaterals between the pedal arteries in five patients. Additional information from 3DRA resulted in a change of treatment strategy in six patients. Conclusion, 3DRA of the foot contains valuable additional real-time information to better guide peripheral vascular interventions in patients with CLI and nonhealing tissue lesions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

A new method of relative radiometric calibration for hyperspectral imaging based on skylight monitor

A new method of relative radiometric calibration based on the skylight monitor is presented for airborne hyperspectral remote sensing to eliminate the negative effect on data resulting from the time-varying skylight during a long flight. In terms of large area operation of mapping or spectral imaging, a static radiation model hardly reveals actual cases as external environment keeps changing. Compared with common considerable factors, fluctuation of skylight is the primary one that may observably affect spectral imaging continuously leading to radiation changes upon flight lines. To get rid of the influence of skylight on spectral imagery and monitoring the variation during the flight, we assemble a set of skylight monitor fixed to the aircraft cabin for the sake of recording the semi-spherical skylight keeping synchronization with hyperspectral imaging system and the position and orientation system. The skylight monitor consists of a cosine corrector on the top receiving the skylight, a stable optical fiber and a visible and near-infrared linear array modular spectrometer. The corresponding matched data of skylight with each frame of images is used to correct related digital number after stripe non-uniformity correction, achieving mean relative standard deviation (coefficient of variation or CV) of under 6% for all 256 bands. The final result shows a satisfactory calibrated image of an appropriate evaluation criteria with contrast of the original image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Developing a system for in vivo imaging of maize roots containing iodinated contrast media in soil using synchrotron XCT and XRF

We sought to develop a novel experimental system which enabled application of iodinated contrast media to in vivo plant roots intact in soil and was compatible with time-resolved synchrotron X-ray computed tomography imaging. The system was developed to overcome issues of low contrast to noise within X-ray computed tomography images of plant roots and soil environments, the latter of which can complicate image processing and result in the loss of anatomical information. To demonstrate the efficacy of the system we employ the novel use of both synchrotron X-ray computed tomography and synchrotron X-ray fluorescence mapping to capture the translocation of the contrast media through root vasculature into the leaves. With the application of contrast media we identify fluid flow in root vasculature and visualise anatomical features, which are otherwise often only observable in ex vivo microscopy, including: the xylem, metaxylem, pith, fibres in aerenchyma and leaf venation. We are also able to observe interactions between aerenchyma cross sectional area and solute transport in the root vasculature with depth. Our novel system was capable of successfully delivering sufficient contrast media into root and leaf tissues such that anatomical features could be visualised and internal fluid transport observed. We propose that our system could be used in future to study internal plant transport mechanisms and parameterise models for fluid flow in plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Enhancing the Range and Reliability of the Spacer Layer Imaging Method

The spacer layer imaging method (SLIM) is widely used to measure the thickness of additive and lubricant films, in lubricant development and evaluation, and for fundamental research into elastohydrodynamic lubrication and tribofilm formation mechanisms. The film thickness measurement, as implemented on several popular tribometers, provides powerful, non-destructive in-situ mapping of film topography with nanometre-scale height sensitivity. However, the results can be highly sensitive to experimental procedure, machine condition, and image analysis, in some cases reporting unphysical film thickness trends. The prevailing image analysis techniques make it challenging to interrogate these errors, often hiding their multivariate nonlinear behaviour from the user by spatial averaging. Herein, several common ‘silent errors’ in the SLIM measurement, including colour matching to incorrect fringe orders, and colour drift due to the optical properties of the system or film itself, are discussed, with examples. A robust suite of novel a priori and a posteriori methods to address these issues, and to improve the accuracy and reliability of the measurement, are also presented, including a novel, computationally inexpensive circle-finding algorithm for automated image processing. In combination, these methods allow reliable mapping of films up to at least 800 nm in thickness, representing a significant milestone for the utility of SLIM applied to elastohydrodynamic contact.

EHL film geometry↗
DOE OSTI2019

Characterization of Roughness, Porosity and Thermal Resistances of Continuous Casting Mold Slag Layers Devitrified and Crystallized in Laboratory

This paper presents an experimental investigation about the variation, with temperature, of surface roughness, porosity, degree of crystallization and microstructure, occurring in solid layers of a commercial mold powder slag used for casting medium carbon steel thin slabs. Characterization is done with samples devitrified and crystallized at 698 K, 823 K, 898 K and 973 K (425 °C, 550 °C, 625 °C and 700 °C) for 10,800 seconds, in the furnace of a high-temperature confocal laser scanning microscope. As well as, with slag disks devitrified non-isothermally in a heat flow apparatus based on the application of Fourier's law in steady state conditions, where the high temperature surface is maintained at selected values, between 927 K and 1241 K (654 °C and 968 °C). Surface roughness and porosity determinations are performed by image analysis of confocal laser scanning and scanning electron micrographs, respectively. For the time considered and temperatures ≥ 823 K (550 °C) the slag crystallizes completely, and variations of roughness and porosity with temperature are associated with changes in size and packing of the primary crystalline-phase grains. Crystallization induced porosity is significant and causes slag layer expansion. Hence crystallization does not lead to gap formation with a neighboring surface (e.g., mold wall) by shrinkage. Instead, surface separation arises from development of roughness, and measurements with the heat flow apparatus confirm that interfacial thermal resistance between a slag layer and a metallic wall can be estimated from the sum of the means of maximum peak to valley height of their surfaces. Furthermore, the results establish that slag layers roughness and porosity increase with temperature and reveal effective thermal conductivities in the range of ~ 0.3 to ~ 0.5 W/m K. Low effective thermal conductivities must be a consequence of porous bands, whose porosity depends on temperature and reaches values as high as ~ 0.20 and ~ 0.50.

36 MATERIALS SCIENCE↗