SEARCH · Search NASA

Results for “Parabolic”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗

DOE OSTI2025

Effect of a Fine-Scale Layered Structure of the Atmosphere on Infrasound Signals from Fragmenting Meteoroids

We investigate the influence of a fine-scale (FS) layered structure in the atmosphere on the propagation of infrasound signals generated by fragmenting meteoroids. Using a pseudo-differential parabolic equation (PPE) approach, we model broadband acoustic signals from point sources at altitudes of 35–100 km. The presence of FS fluctuations in the stratosphere (37–45 km) and the lower thermosphere (100–120 km) modifies ray trajectories, causing multiple arrivals and prolonged signal durations at ground stations. In particular, meteoroids fragmenting at 80–100 km can produce two distinct thermospheric arrivals beyond 150 km range, while meteoroids descending to 50 km or below yield weak, long-lived arrivals within the acoustic shadow zone via antiguiding propagation and diffraction. Comparison with observed infrasound data confirms that FS-layered inhomogeneities can account for multi-arrival “N-waves,” broadening potential interpretations of meteoroid signals. The results also apply to other atmospheric-entry objects, such as sample return capsules, emphasizing how FS structure impacts shock wave propagation. In conclusion, our findings advance understanding of wavefield evolution in a layered atmosphere and have broad relevance for global infrasound monitoring of diverse phenomena (e.g., re-entry capsules, rocket launches, and large-scale explosions).

Aeroacoustics↗

DOE OSTI2024

Effects of uranium metal carbon content on hydriding kinetics and corrosion blister number/area at sub-ambient pressures

Carbon is a common impurity in uranium metal, resulting in a number of uranium–carbon inclusion phases that contribute to an increase in metal defect density as carbon content increases. It is widely held that uranium hydride corrosion preferentially nucleates at these defect sites, and that an increase in carbon content will therefore represents an increase in uranium hydride corrosion sites on the metal surface. We hydrided six uranium sources with differing carbon contents to explore whether this assumption holds in a sub-ambient (~ 0.1 atm hydrogen), sealed environment, and report the resulting reaction kinetics and uranium hydride blister benchmarking data. We find that carbon content is not strongly correlated with reaction kinetics terms or the resulting hydride blister number and area, but that there is a tight relationship between corrosion blister number/area and kinetics as is expected. Further, we find that there is a strong trend of decreasing variance in the blister number, blister area, and induction time as carbon content increases (higher carbon content results in more reproducible blister populations). Additionally, we find a narrow band of uranium metal consumption at the end of the parabolic phase of reaction progress (beginning of linear growth phase) of 0.098 ± 0.011 w/w%, a fact that may be useful in assaying hydrogen corrosion of uranium metal within sealed environments generally.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2021

Lessons Learned in Employing Data Analytics to Predict Oxidation Kinetics and Spallation Behavior of High-Temperature NiCr-Based Alloys

Machine learning (ML) can offer many advantages in predicting material properties over traditional materials development methods based solely on limited experimental investigations or physical-based simulations with the capability to reduce development cost, risk, and time. However, so far, limited efforts have been made to predict alloy oxidation kinetics and spallation behavior via ML due to the lack of consistently measured and sufficient experimental data and the inherent complexity in oxidation behavior of multicomponent high-temperature alloys. A previous study reported the ability of ML to predict oxidation kinetics of NiCr-based alloys as a function of alloy composition and operating conditions. Here, the performance of a ML model in predicting rate constants and spallation probability was evaluated in light of the roles of the data distribution of the experimental dataset (data analytics), the alloy composition, the exposure environment and the chosen oxidation approach to extracting kinetic values from the measured mass changes (but using either a simple parabolic law or a statistical cyclic oxidation model). Potential strategies to improve the predictions and enhance the extrapolative capability of the previously trained model will be discussed.

36 MATERIALS SCIENCE↗

DOE OSTI2023

High-Temperature Corrosion of Chromia-Forming Ni-Based Alloys in CO 2 Containing Impurities

Future technologies require structural materials resistant to environmental degradation in high-temperature CO 2 -rich environments. In this work, we exposed several commercially available Ni-based alloys (230, 263, 282, 617, 625, and 740H) to atmospheric pressure gases intended to simulate the compositions expected in future direct-fired supercritical CO 2 power cycles. The alloys were exposed to 95% CO 2 + 4% H 2 O + 1% O 2 and the same gas containing 0.1% SO 2 at temperatures of 600, 650, 700, 750, and 800 °C for 2500 h. With minor exceptions, chromia scales formed on all alloys at all temperatures in the SO 2 -free gas, yielding parabolic growth rates that followed an Arrhenius temperature dependence. Behavior in the SO 2 -containing gas was more complex. Generally, the alloys performed well at temperatures of 650, 750, and 800 °C. While some alloys further performed relatively well across the whole temperature range, several of the alloys experienced chromia failure resulting in non-protective duplex oxide scales and high oxidation rates, at temperatures of 600 and 700 °C. Deviation from protective behavior was associated with internal sulfide formation and, additionally for the case of 600 °C, external sulfate formation. Extensive carburization accompanied growth of the non-protective duplex oxide scales, which made it more difficult for the alloy to recover after initiation of the sulfur-induced accelerated corrosion process. The thermodynamic and kinetic factors influencing the accelerated corrosion in the presence of sulfur are discussed. The results suggest that caution is required when assessing compatibility of Ni-based alloys for CO 2 -based systems when sulfur-based impurities are expected.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Microstructural, Oxidation, and Mechanical Behavior of NbTi-Based Refractory Alloys with 5 to 10 Pct Co, Cr, and Ni Additions

NbTi-based refractory alloys with additions of Co, Cr, and Ni represent an interesting medium-entropy alloy system with potential for protective oxide film formation, high strength, and ductility. This study investigates the microstructural evolution, oxidation behavior, and mechanical properties of NbTi-based alloys containing 5 to 10 at. pct Co, Cr, and Ni. CALPHAD predictions suggest that this composition range can be heat treated to obtain a predominantly body-centered cubic matrix phase. Mechanical properties, including microhardness, yield strength, maximum strength, and specific strength are evaluated through isothermal compression tests conducted between room temperature and 800 °C. The oxidation kinetics of these alloys are assessed through discontinuous oxidation tests. Parabolic oxidation kinetics were observed for NbTi–10Ni and NbTi–5Co, while linear oxidation kinetics were found for NbTi–10Cr and NbTi–10(CoCrNi). Microstructures and oxide layers are characterized using X-ray diffraction, electron backscatter diffraction, energy-dispersive X-ray spectroscopy, and scanning electron microscopy. All alloys exhibit significant mechanical softening between room temperature and 800 °C, with elastic-perfectly plastic flow observed at 800 °C. The addition of 10 pct Cr to NbTi resulted in two BCC phases up to 1050 °C, conflicting with CALPHAD predictions of a single-phase solid solution at this temperature, and resulting in higher flow stress at 800 °C. NbTi–10(CoCrNi) exhibited the lowest flow stress at 800 °C despite having more ‘cocktail effect’ potential and insufficient molar fractions of Co, Cr, or Ni to form a desirable protective oxide film.

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 OSTI2021

Generalized multifractality at spin quantum Hall transition

Generalized multifractality characterizes scaling of eigenstate observables at Anderson-localization critical points. We explore generalized multifractality in 2D systems, with the main focus on the spin quantum Hall (SQH) transition in superconductors of symmetry class C. Relations and differences with the conventional integer quantum Hall (IQH) transition are also studied. Using the field-theoretical formalism of non-linear sigma-model, we derive the pure-scaling operators representing generalizing multifractality and then “translate” them to the language of eigenstate observables. Performing numerical simulations on network models for SQH and IQH transitions, we confirm the analytical predictions for scaling observables and determine the corresponding exponents. Remarkably, the generalized-multifractality exponents at the SQH critical point strongly violate the generalized parabolicity of the spectrum, which implies violation of the local conformal invariance at this critical point.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

A transmissive concentrator photovoltaic module with cells directly cooled by silicone oil for solar cogeneration systems

Hybrid concentrator photovoltaic-thermal systems can cogenerate electricity and heat by beam-splitting incoming concentrated light onto photovoltaic cells and a thermal receiver to increase total conversion efficiency and potentially reduce system cost. To demonstrate this, we have designed and prototyped a transmissive spectrum-splitting concentrator photovoltaic module that maximizes solar energy conversion by utilizing the entire solar spectrum. Visible light is collected using infrared-transmissive triple-junction photovoltaic cells to achieve an in-band module efficiency of 43.3% for light of wavelength λ < 873 nm, while 44.2% of out-of-band light with λ > 873 nm is transmitted through for collection by a thermal receiver. During testing on a dual-axis tracked parabolic concentrator dish at up to 166 suns, cell temperatures were maintained at 119 °C or below via a novel active cooling method. This cooling system strictly flows silicone oil directly across both sides of the cells, without inhibiting optical transmission, as verified through experimentation and simulation. The module was validated outdoors for 572 sun·hrs, and achieved a maximum thermal receiver temperature of 180 °C. 86.1% of incident solar power is collected at 166 suns average concentration collectively among the electrical, cooling, and thermal receiver subsystems. The remaining 13.9% is lost to mirror reflectivity, dish shadowing, receiver reflection, and thermal losses. The ability to directly cool the cells with an inert silicone oil offers the potential for reduced system cost relative to previous transmissive hybrid concentrator photovoltaic-thermal systems, including microfluidic-cooled designs. This solar cogeneration capability is valuable in a wide range of commercial and industrial applications.

14 SOLAR ENERGY↗

DOE OSTI2025

Materials property changes in ETU-10 graphite due to neutron irradiation at elevated temperatures

Graphite grade ETU-10, from IBIDEN Co., Ltd. Has been irradiated in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). The irradiation program was developed to provide a preliminary study the irradiation-induced property changes to the dimensions/volume, elastic properties, strength, electrical resistivity, coefficient of thermal expansion, and thermal diffusivity/conductivity over a range of temperatures and neutron exposures that may be relevant for future nuclear reactors. The irradiation envelope covers a range of irradiation temperatures (300°C–900 °C) and fluences (up to 40 × 10 25 n/m 2 [E > 0.1 MeV] or ~30 dpa) that would be relevant for advanced nuclear reactors. Further, the dimensional change was observed to be anisotropic for an isotropic graphite, the specimen dimensions, volume, Young's modulus, shear modulus, and strength all displayed a parabolic fluence dependence, the electrical resistivity had a rapid rise followed by a decrease and a later increase, at high fluence the mean coefficient of thermal expansions was similar for all irradiation temperatures, and thermal conductivity rapidly decreased followed by a continued loss.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge

In this work, we present the results of two-dimensional modeling of the effects of non-equilibrium excitation and electrode geometry on H 2 /air ignition in a nanosecond plasma discharge. A multiscale adaptive reduced chemistry solver for plasma assisted combustion (MARCS-PAC) based on PASSKEy discharge modeling package and compressible multi-component reactive flow solver ASURF+ is developed and validated. This model is applied to simulate the impact of non-equilibrium plasma excitation and electrode geometry and heat loss on the dynamics of the discharge from streamer to spark and ignition kernel development in a H 2 /air mixture with a pair of cylindrical electrodes. The results show that the plasmagenerated species (N 2 (A), N 2 (B), N 2 (a'), N 2 (C), O( 1 D), O and H) in the spark and afterglow significantly accelerate the ignition kernel development. The increase of discharge voltage at the same total discharge energy promotes the non-equilibrium active species production. It is found that the production of electronically excited species at higher reduced electric field strength is more efficient in enhancing ignition in comparison to the vibrational excitation and heating. Moreover, the 2D simulation clearly reveals that the electric field and active species distribution are highly non-uniform. The streamers are initiated at the sharp outer edges of the negative and positive electrodes by a strong electric field while the electric field is much weaker at the centerline of the electrodes. Furthermore, the simulations reveal that the ignition enhancement is sensitive to the variation of electrode shape, diameter, and gap size due to the changes of electric field distribution and location of streamer formation. A cylindrical electrode produces a larger discharge volume and ignition kernel than the parabolic and spherical electrodes, when the discharge is localized near the axis of the gap. It is found that there is a non-monotonic dependence of ignition kernel size on the electrode diameter and inter-electrode distance. The increase of electrode diameter and gap size above the optimal conditions leads to the reduction of ignition kernel volume, due to the decrease of active species concentration and gas temperature. At a larger electrode surface area and electrode diameter as well as smaller electrode gap size, the heat loss to electrode plays a greater role in reducing the ignition kernel size and slowing ignition kernel development. This work provides insights and guidance to understand the kinetic enhancement of non-equilibrium plasma and the effects of electrode geometries on ignition for the optimization ignitors in advanced engines.

42 ENGINEERING↗

DOE OSTI2025

Direct Discontinuous Galerkin methods for the reacting multi-component flow equations

The Direct Discontinuous Galerkin (DDG (Liu and Yan, 2008)) method and a counterpart with Interface Correction (DDGIC (Danis and Yan, 2022)) are extended to compute diffusion terms that arise when solving the compressible multi-component flow equations in thermochemical nonequilibrium. Thermodynamic properties, transport properties, chemical reaction rates, and energy exchange terms are computed using Mutation++ (Scoggins et al., 2020). The DG method is applied on unstructured grids, where the accuracy and convergence rates can be sensitive to the numerical method chosen for parabolic terms. A method for determining the homogeneity tensor of the flow equations required for DDGIC is shown. The convergence properties of the DDG methods are studied and compared to the Interior Penalty (IP) method. A number of numerical experiments are conducted to assess the accuracy and performance of the method. The numerical results and convergence studies indicate that DDG and DDGIC provide accurate solutions and perform well for general flows in thermochemical nonequilibrium.

Diffusion↗

DOE OSTI2021

Water vapor oxidation of SiC layer in surrogate TRISO fuel particles

Under accidental conditions for high temperature gas-cooled reactors (HTGR), the SiC layer in tri-structural-isotropic (TRISO) fuel particles can be exposed to water vapor. In this study, oxidation behaviors of surrogate TRISO fuel particles were investigated in a He-20 vol% water vapor mixed atmosphere at temperatures up to 1600 °C. The growth of the crystalline oxide passivation layer with cracks and pores followed a parabolic law with time, where the maximum was 2.3 μm at 1600 °C. The oxide layer and the SiC surface under the oxide became flattened with increasing temperature, as a function of the silica viscosity and the diffusion path of water vapor. Volatilization of the oxide layer was analyzed using a mechanistic model that an inert gas in oxidizing atmospheres could influence the magnitude of volatilization. The fracture load and strength of the oxidized and thinned SiC layer were numerically estimated to decrease from 2.27 to 1.69 N and from 317 to 299 MPa, respectively, with the SiC thickness decrease from 35 to 32 μm. This prediction indicates that the oxidized SiC layer should retain fission products. Additionally, the mechanical integrity of each layer in the TRISO fuel particle after oxidation was evaluated. The results in this work provide important data for the safety analysis of accidental scenarios in HTGRs.

36 MATERIALS SCIENCE↗

DOE OSTI2020

High temperature oxidation behaviors of bulk SiC with low partial pressures of air and water vapor in argon

Oxidation behaviors of bulk SiC at 1000–1400 °C in Ar-20 vol% H 2 O-10 vol% air were investigated. The thickness of the passivation oxide layer showed a parabolic increase with time. Dense and uniform layers formed below 1200 °C, whereas pores and cracks were observed at higher temperatures. An analytical model was proposed to estimate the oxide layer growth and SiC mass loss rate. The SiC layer should be effective in retaining the fission products after oxidation by excessive water and air. This work is the first detailed study of SiC layer oxidation during loss of coolant accidents (LOCA) for HTGRs.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effect of exposure temperature on the corrosion behavior of a FeNiCrCuAl high entropy alloy in supercritical water

The corrosion behavior of a FeNiCrCuAl high entropy alloy exposed to supercritical water at 380–650 °C and 25 MPa was investigated. The weight gain exhibited a nearly parabolic kinetics. At 380 °C, Al 2 O 3 precipitated on Fe 3 O 4 particles in FeCr-rich regions and intensified as the temperature increased to 650 °C, eventually resulting in an Al 2 O 3 particle. A mixed Fe 3 O 4 and NiAl 2 O 4 outer layer developed on NiAl-rich regions at 380 °C, whereas an Al 2 O 3 layer formed at 650 °C. Cu 2 O formed along the boundaries of NiAl-rich regions, regardless of exposure temperatures. In conclusion, the corrosion mechanism at different temperatures was elucidated.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗

DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

SEM petrography of dispersed organic matter in black shales: A review

Organic matter (OM)-hosted pores are important constituents of the pore system of black shales and play a crucial role in determining their methane adsorption capacity and porosity. OM-hosted pores are generally observed and described with scanning electron microscope (SEM) on Ar ion-milled surfaces. However, SEM imaging is not able to reliably distinguish OM types and relate the observed pores to specific macerals. Partly because of this inability to relate organic pores to macerals, the evolution of organic porosity during thermal maturation remains poorly understood. In this paper, we review the petrographic characteristics of dispersed organic matter (DOM) in black shales under the SEM. Organic petrographic classification of DOM developed for reflected-light microscopy is so far the most practical method when describing DOM in black shales under the SEM because this classification has information on the origin of DOM. Therefore, correlative microscopy (combination of reflected-light and electron microscopy) is the most effective method to identify both OM types and OM-hosted pores. This method, however, is not readily available to most researchers. Although identifying OM on the basis of SEM observations is a challenging task, it is achievable provided there is a good understanding of the studied shales, especially their thermal maturity and original OM composition. Therefore, the overall objective of this paper is to review petrographic characteristics of DOM in black shales under the SEM to provide some guidelines for identifying DOM from SEM observations. We also review factors that control the formation and preservation of OM-hosted pores. OM-hosted pores consist of primary and secondary organic pores. Primary organic pores are pores inherited from the biological structure of the original OM. Secondary organic pores develop during hydrocarbon generation and expulsion from oil-prone OM and are hosted by solid bitumen or pyrobitumen. The development of secondary organic pores is controlled by thermal maturity and OM type, and their preservation is subject to thermal maturity, OM content, and mineralogical composition. The presented view of the evolution of micropore and mesopore characteristics of OM with thermal maturity is based on data from the literature. The specific surface area and pore volume of OM in black shales follow a parabolic pattern with increasing thermal maturity (quantified via vitrinite reflectance, R o ). The initial increase reflects development of OM-hosted pores, and the subsequent decrease is due to denser stacking of aromatic units in the macromolecular structure of OM, with maximum values (specific surface area ~ 300 m 2 /g and pore volume ~ 0.3 cm 3 /g) reached at R o values in the 2.5-3.5% range. The contribution of OM-hosted pores to the pore characteristics of black shales depends on OM content, OM type, and thermal maturity.

03 NATURAL GAS↗