Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial transferability”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Comparison of theory with the experimental characterization of the spatial frequency response of interferometers using a binary pseudo-random array sample

Experimental evaluations of the surface height response of an interference microscope using a binary pseudo-random array test sample are compared with a theory based on a Fourier optics model. Measurements of key instrument characteristics, including the illumination, imaging, and obscuring apertures of three different Mirau objectives, support the theoretical calculations. Agreement between experimental and theoretical modeling confirms the predictability of the spatial frequency response for the purpose of specification and optimization of instrument configuration for specific metrology tasks. The results also provide confidence in methods of compensating for the decrease in instrument response with spatial frequency.

Calibration↗
DOE OSTI2026

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗
DOE OSTI2020

1-D imaging of rotation-vibration non-equilibrium from pure rotational ultrafast coherent anti-Stokes Raman scattering

We present one-dimensional (1-D) imaging of rotation-vibration non-equilibrium measured by two-beam pure rotational hybrid femtosecond/picosecond coherent anti-Stokes Raman scattering (fs/ps CARS). Simultaneous measurements of the spatial distribution of molecular rotation-vibration non-equilibrium are critical for understanding molecular energy transfer in low temperature plasmas and hypersonic flows. However, non-equilibrium CARS thermometry until now was limited to point measurements. The red shift of rotational energy levels by vibrational excitation was used to determine the rotational and vibrational temperatures from 1-D images of the pure rotational spectrum. Vibrational temperatures up to 5500 K were detected in a CH4/N2 nanosecond-pulsed pin-to-pin plasma within 2 mm near the cathode. As a result, this approach enables study of non-equilibrium systems with 40 µm spatial resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

2D perovskite-based high spatial resolution X-ray detectors

Abstract X-ray radiography is the most widely used imaging technique with applications encompassing medical and industrial imaging, homeland security, and materials research. Although a significant amount of research and development has gone into improving the spatial resolution of the current state-of-the-art indirect X-ray detectors, it is still limited by the detector thickness and microcolumnar structure quality. This paper demonstrates high spatial resolution X-ray imaging with solution-processable two-dimensional hybrid perovskite single-crystal scintillators grown inside microcapillary channels as small as 20 µm. These highly scalable non-hygroscopic detectors demonstrate excellent spatial resolution similar to the direct X-ray detectors. X-ray imaging results of a camera constructed using this scintillator show Modulation Transfer Function values significantly better than the current state-of-the-art X-ray detectors. These structured detectors open up a new era of low-cost large-area ultrahigh spatial resolution high frame rate X-ray imaging with numerous applications.

2D perovskite↗
DOE OSTI2024

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

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

A Holistic Algorithmic Approach to Improving Accuracy, Robustness, and Computational Efficiency for Atmospheric Dynamics

Atmospheric weather and climate models must perform simulations very quickly to be useful. Therefore, modelers have traditionally focused on reducing computations as much as possible. However, in our new era of increasingly compute-capable hardware, data movement is now the prohibiting expense. This study examines the computational benefits of a new algorithmic approach to modeling atmospheric dynamics on scales relevant to weather and climate simulation. Rather than minimizing computations, this new approach considers the larger problem more holistically, including spatial accuracy, temporal accuracy, robustness (i.e., oscillations), on-node efficiency, and internode data transfers together at once. Numerical experiments demonstrate how computations can be strategically increased to simultaneously address each of these constraints while reducing data movement to adapt to modern accelerated hardware. The new algorithm can achieve at times up to 80% peak floating point throughput in single precision on the Nvidia Tesla V100 GPU, where the traditional approach is shown to only achieve single-digit floating point efficiency. Further, the new algorithm is twice as fast as a standard Runge--Kutta time integrator, and high-order accuracy with Weighted Essentially Non-Oscillatory (WENO) limiting came at less than 30% additional runtime cost on a GPU, thus increasing the accuracy per degree of freedom.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Supramolecular Control of Redox-Mediated Interactions for Ion-Selective Electroseparations

This report summarizes efforts from the Su Group at the University of Illinois at Urbana-Champaign in researching and developing materials and processes for the recovery of rare earth elements and other metal elements that are critical to the United States economy and national security, under DOE Award No. DE-SC0021409. This report highlights efforts in the recovery of: • Transition metal elements (including arsenic, chromium, molybdenum, rhenium, vanadium): Redox-active metallopolymers are able to capture these elements in their oxyanion form. Direct control of target ion capture/release through electrochemical-mediated switching of metallocene binding sites reduces chemical input and lowers generation of waste in element recovery processes. We study the effects of metallopolymer structure on targeted ion binding through a combination of experimental screening of ion selectivity from multi-component mixtures, computational efforts which elucidate the charge-transfer origins of metallocene binding and provide insight on how binding site structure influences selectivity, and state-of-the-art neutron reflectivity of metallocene polymer films to understand the interplay of polymer-ion (charge-transfer), polymer-solvent (polymer hydrophilicity/hydrophobicity), and ion-solvent (ion hydration) interactions with spatial resolution of polymer films. • Rare earth elements (including yttrium, neodymium, europium, gadolinium, dysprosium, cerium): Here we employ redox-copolymers for recovery of rare earth elements by electrochemically regenerated ion-exchange. Rare earth elements require an alternative route to capture/release than directly with metallopolymers since rare earth elements typically exist in cationic form in water. To address this, copolymerization of traditional ion-exchange moieties with metallocene moieties allowed for capture of rare earth elements with the ion-exchange moieties and electrochemical-mediated release of rare earth elements through oxidation of metallocene moieties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Lessons Learned from Ecosystem-Scale Experimental Field Studies (Workshop Report)

Efforts to understand and predict ecosystem responses to environmental change require long-term, large-scale, spatially representative experiments and observations that capture natural variability, test predictive models, and generate transferable knowledge. Such studies are indispensable for unraveling the complexities of terrestrial ecosystems and their responses to disturbances and evolving environmental conditions, while generating the data necessary for developing mechanistic models and predictive tools that inform decision-making processes. Having a rich history of designing and executing large-scale ecosystem experiments, the U.S. Department of Energy’s Environmental System Science program convened a workshop in January 2025 that brought together leaders in the field to distill critical lessons from decades of experience in large-scale experiments. The workshop aimed to (1) provide an ecosystem experiment primer for best practices, thus ensuring a high scientific return on investment for funding agencies, and (2) offer a robust framework for the design and management of future research initiatives. This report synthesizes insights and experiences from workshop participants and is structured to capture the entire research life cycle, from goal setting and design to operations, adaptive management, team dynamics, collaborations, and the often overlooked aspect of decommissioning. By synthesizing decision-making and lessons learned across diverse research approaches, the report aims to provide a template of essential factors to consider when designing successful long-term, large-scale ecosystem experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Modeling the co-assembly of binary nanoparticles

Abstract In this work, we present a binary assembly model that can predict the co-assembly structure and spatial frequency spectra of monodispersed nanoparticles with two different particle sizes. The approach relies on an iterative algorithm based on geometric constraints, which can simulate the assembly patterns of particles with two distinct diameters, size distributions, and at various mixture ratios on a planar surface. The two-dimensional spatial-frequency spectra of the modeled assembles can be analyzed using fast Fourier transform analysis to examine their frequency content. The simulated co-assembly structures and spectra are compared with assembled nanoparticles fabricated using transfer coating method are in qualitative agreement with the experimental results. The co-assembly model can also be used to predict the peak spatial frequency and the full-width at half-maximum bandwidth, which can lead to the design of the structure spectra by selection of different monodispersed particles. This work can find applications in fabrication of non-periodic nanostructures for functional surfaces, light extraction structures, and broadband nanophotonics.

Mohanty, Saurav (ORCID:0009000302592855)↗
DOE OSTI2022

Impacts of Sub‐Grid Topographic Representations on Surface Energy Balance and Boundary Conditions in the E3SM Land Model: A Case Study in Sierra Nevada

Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Direct Optoelectronic Imaging of 2D Semiconductor–3D Metal Buried Interfaces

The semiconductor-metal junction is one of the most critical factors for high performance electronic devices. In two-dimensional (2D) semiconductor devices, minimizing voltage drop at this junction is of particularly challenging and important. Despite numerous studies on contact resistance in 2D semiconductors, the exact nature of the buried interface under a 3D metal remains elusive and unclear. Herein, we report the direct observation of 2D semiconductor-metal buried interface using recently developed metal-assisted transfer technique to expose the buried interface which is then directly investigated using scanning probe techniques. We characterize the spatially varying electronic and optical properties of this buried interface with < 20 nm resolution. To be specific, potential, conductance and photoluminescence at the buried metal/MoS 2 interface are correlated under variety of metal deposition conditions as well as type of metal contacts. We observe that direct evaporation of Au on MoS 2 induces large strain of ~5% in the MoS 2 which coupled with charge transfer leads to degenerate doping of the MoS 2 underneath the contact. These factors lead to improvement of contact resistance to record values of 138 kΩ μm measured using local conductance probes. This approach was adopted to investigate MoS 2 -In/Au alloy interface, which records lowest contact resistance of 63 kΩ μm. Overall, our results highlight that the MoS 2 /Metal interface is very sensitive to method of contact making and provides a universal strategy to characterize buried contact interfaces involving 2D semiconductors.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗
DOE OSTI2020

Photochemical Upconversion in Water Using Cu(I) MLCT Excited States: Role of Energy Shuttling at the Micellar/Water Interface

Photochemical upconversion (UC) through triplet–triplet annihilation (TTA), which employs a visible absorbing triplet photosensitizer and an annihilator, is a process that generates a high energy photon from two lower energy photons. TTA-UC has been largely developed in pure organic solvents and solid-state polymeric constructs while featuring near exclusive use of rare and expensive metals within the photosensitizer. In this current investigation, we demonstrate that TTA-UC from the long lifetime earth-abundant photosensitizer [Cu(dsbtmp) 2 ](PF) 6 (dsbtmp = 2,9-di(sec-butyl)-3,4,7,8-tetramethyl-1,10-phenanthroline), abbreviated as Cu-PS, functions in water through encapsulation within a cationic-based assembly. Cetyltrimethylammonium bromide was the surfactant of choice as it electrostatically binds the negatively charged water-soluble 10-phenylanthracene-9-carboxylate (PAC) acceptor/annihilator and ultimately facilitates energy transfer across the interface. Efficient triplet–triplet energy transfer (TTET) from Cu-PS to the PAC acceptor was achieved in this aqueous assembly. Unfortunately, the hindered mobility of the PAC moieties ultimately hampered the annihilation process, and this was reflected in attenuated TTA rates and efficiencies. The combined experimental data illustrated that the water-soluble PAC acceptor was able to vectorially deliver the excited-state energy stored in Cu-PS across the interface into the bulk aqueous solution by engaging in excited-state electron transfer with methyl viologen acceptors. Lastly, these results are important for remotely operating photoredox reactions in water while rendering a photosensitizer spatially isolated in the hydrophobic core of a micelle.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nanoscale imaging of quantum dot dimers using time-resolved super-resolution microscopy combined with scanning electron microscopy

Abstract Time-resolved super-resolution microscopy was used in conjunction with scanning electron microscopy to image individual colloidal CdSe/CdS semiconductor quantum dots (QD) and QD dimers. The photoluminescence (PL) lifetimes, intensities, and structural parameters were acquired with nanometer scale spatial resolution and sub-nanosecond time resolution. The combination of these two techniques was more powerful than either alone, enabling us to resolve the PL properties of individual QDs within QD dimers as they blinked on and off, measure interparticle distances, and identify QDs that may be participating in energy transfer. The localization precision of our optical imaging technique was ∼3 nm, low enough that the emission from individual QDs within the dimers could be spatially resolved. While the majority of QDs within dimers acted as independent emitters, at least one pair of QDs in our study exhibited lifetime and intensity behaviors consistent with resonance energy transfer from a shorter lifetime and lower intensity donor QD to a longer lifetime and higher intensity acceptor QD. For this case, we demonstrate how the combined super-resolution optical imaging and scanning electron microscopy data can be used to characterize the energy transfer rate.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Electron-mediated entanglement of two distant macroscopic ferromagnets within a nonequilibrium spintronic device

Using the nascent concept of quantum spin-transfer torque, we demonstrate that a current pulse can be harnessed to entangle quantum localized spins of two spatially separated ferromagnets (FMs) which are initially unentangled. The envisaged setup is composed of a spin-polarizer (FM p ) and a spin-analyzer (FM a ) FM layers separated by a normal metal (NM) spacer. The injection of a current pulse into the device leads to a time-dependent superposition of many-body states characterized by a high degree of entanglement between the spin degrees of freedom of the two distant FM layers. The nonequilibrium dynamics are due to the transfer of spin angular momentum from itinerant electrons to the localized spins via a quantum spin-torque mechanism that remains active even for collinear but antiparallel arrangements of the FM p and FM a magnetizations (a situation in which the conventional spin torque is absent). We quantify the mixed-state entanglement generated between the FM layers by tracking the time evolution of the full density matrix and analyzing the build-up of the mutual logarithmic negativity over time. The effect of decoherence and dissipation in the FM layers due to coupling to bosonic baths at finite temperature, the use of multielectron current pulses, and the dependence on the number of spins are also considered in an effort to ascertain the robustness of our predictions under realistic conditions. Lastly, we propose a “current-pump–x-ray-probe” scheme, utilizing ultrafast x-ray spectroscopy, that can witness nonequilibrium and transient entanglement of the FM layers by extracting its time-dependent quantum Fisher information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A multiscale packed-bed reactor model for sustainable ethylene production via chemical looping oxidative coupling of methane

The rising global warming concerns and shale gas discovery have prompted research in the direction of greenhouse gas (GHG), such as methane, reduction and conversion. Oxidative coupling of methane (OCM) offers a pathway to low carbon-intense valorization of methane while producing ethylene, a chemical regarded as central to the petrochemical industry. Even after decades of OCM discovery, researchers keep understanding the process and underlying chemical reactions in a pursuit to achieve industrial viability for OCM. Here, in general, OCM suffers from low C 2 selectivity, yield and reactor temperature runaways due to highly exothermic nature of its reactions. Computational Fluid Dynamics (CFD) tools help analyze spatial gradients within the reactor to deeply understand the diffusion of species, mass and heat transfer phenomena. Furthermore, challenges associated with scaling up such as hot spot formation and parametric sensitivity can be addressed without having to expend on costly experiments. The current paper presents a multiscale packed-bed reactor CFD model coupled with a chemical kinetic model for the chemical looping OCM. The CFD model includes two scales i.e., macroscale for catalyst bed and microscale for individual pellets. Moreover, a chemical kinetic model based on 10 gas-phase reactions is integrated with the CFD model. An additional surface reaction for the formation of gas-phase oxygen from catalyst surface is added to account for the absence of feed oxygen. The model is calibrated against experimental results. The calibrated model captures trends in CH 4 conversion, C 2 selectivity and C 2 yield within a ± 4.35 % range across a temperature range of 700-900 °C. Moreover, model fidelity is evaluated by varying key computational parameters such as mesh resolution and time step size. The model is also verified by varying the inlet methane concentration and the gas hourly space velocity (GHSV) and comparing the results with literature. A sensitivity analysis and scale-up of the current model is undergoing.

Chemical looping↗
DOE OSTI2022

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Estimation of extreme temperatures in direct solar methane pyrolysis within a porous medium

Porous media have wide application in renewable energy conversion processes, such as solar-thermal fuels production and decarbonization. Heat transport mechanisms within porous media can be highly complex, particularly under extreme conditions encountered in concentrated solar thermal reactors in which direct measurement of temperature is challenging. Here, we implement and report an inverse heat conduction model to estimate the temperature distribution throughout a porous substrate domain in a direct solar methane pyrolysis process. By solving a two-dimensional heat transfer problem and applying an inverse optimization algorithm, we estimate the quasi-steady state spatial temperature distribution in a fibrous porous carbon substrate. The results are validated indirectly by experimentally measured graphite deposition and a simplified reaction kinetic model.

finite difference method↗