Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled system”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tightly-coupled camera/LiDAR integration for point cloud generation from GNSS/INS-assisted UAV mapping systems

Unmanned aerial vehicles (UAVs) equipped with integrated global navigation satellite systems/inertial navigation systems (GNSS/INS) together with cameras and/or LiDAR sensors are being widely used for topographic mapping in a variety of applications such as precision agriculture, coastal monitoring, and archaeological documentation. Integration of image-based and LiDAR point clouds can provide a comprehensive 3D model of the area of interest. For such integration, ensuring a good alignment between data from the different sources is critical. Although many works have been conducted on this topic, there is still a need for a rigorous integration approach that minimizes the discrepancy between camera and LiDAR data caused by inaccurate system calibration parameters and/or trajectory artifacts. This study proposes an automated tightly-coupled camera/LiDAR integration workflow for GNSS/INS-assisted UAV systems. The proposed strategy is conducted in three main steps. First, an image-based point cloud is generated using a LiDAR/GNSS/INS-assisted structure from motion (SfM) strategy. Then, feature correspondences between image-based and LiDAR point clouds are automatically identified. Finally, an integrated-bundle adjustment procedure including image points, LiDAR raw measurements, and GNSS/INS information is conducted to minimize the discrepancy between point clouds from different sensors while estimating system calibration parameters and refining the trajectory information. The proposed SfM strategy and integration framework are evaluated using five datasets. The SfM results show that using LiDAR data can facilitate feature matching and further increase the number of reconstructed 3D points. The experimental results also illustrate that the developed automated camera/LiDAR integration strategy is capable of accurately estimating system calibration parameters to achieve good alignment among camera/LiDAR data from single/multiple systems. Finally, an absolute accuracy in the range of 3–5 cm is achieved for the image/LiDAR point clouds after the integration process.

42 ENGINEERING↗
DOE OSTI2024

Local atomic and magnetic structure of multiferroic (Sr,Ba)⁢(Mn,Ti)⁢O 3

Here, we present a detailed study of the local atomic and magnetic structure of the type-I multiferroic perovskite system (Sr,Ba)(Mn,Ti)O 3 using x-ray and neutron pair distribution function (PDF) analysis, polarized neutron scattering, and muon spin relaxation (μSR) techniques. The atomic PDF analysis reveals widespread nanoscale tetragonal distortions of the crystal structure even in the paraelectric phase with average cubic symmetry, corresponding to incipient ferroelectricity in the local structure. Magnetic PDF analysis, polarized neutron scattering, and μSR likewise confirm the presence of short-range antiferromagnetic correlations in the paramagnetic state, which grow in magnitude as the temperature approaches the magnetic transition. We show that these short-range magnetic correlations coincide with a reduction of the tetragonal (i.e., ferroelectric) distortion in the average structure, suggesting that short-range magnetism can play an important role in magnetoelectric and/or magnetostructural phenomena even without genuine long-range magnetic order. The reduction of the tetragonal distortion scales linearly with the local magnetic order parameter, pointing to spontaneous linear magnetoelectric coupling in this system. These findings provide greater insight into the multiferroic properties of (Sr,Ba)(Mn,Ti)O 3 and demonstrate the importance of investigating the local atomic and magnetic structure to gain a deeper understanding of the intertwined degrees of freedom in multiferroics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Regional environmental controllers influence continental scale soil carbon stocks and future carbon dynamics

Abstract Understanding the influence of environmental factors on soil organic carbon (SOC) is critical for quantifying and reducing the uncertainty in carbon climate feedback projections under changing environmental conditions. We explored the effect of climatic variables, land cover types, topographic attributes, soil types and bedrock geology on SOC stocks of top 1 m depth across conterminous United States (US) ecoregions. Using 4559 soil profile observations and high-resolution data of environmental factors, we identified dominant environmental controllers of SOC stocks in 21 US ecoregions using geographically weighted regression. We used projected climatic data of SSP126 and SSP585 scenarios from GFDL-ESM 4 Earth System Model of Coupled Model Intercomparison Project phase 6 to predict SOC stock changes across continental US between 2030 and 2100. Both baseline and predicted changes in SOC stocks were compared with SOC stocks represented in GFDL-ESM4 projections. Among 56 environmental predictors, we found 12 as dominant controllers across all ecoregions. The adjusted geospatial model with the 12 environmental controllers showed an R 2 of 0.48 in testing dataset. Higher precipitation and lower temperatures were associated with higher levels of SOC stocks in majority of ecoregions. Changes in land cover types (vegetation properties) was important in drier ecosystem as North American deserts, whereas soil types and topography were more important in American prairies. Wetlands of the Everglades was highly sensitive to projected temperature changes. The SOC stocks did not change under SSP126 until 2100, however SOC stocks decreased up to 21% under SSP585. Our results, based on environmental controllers of SOC stocks, help to predict impacts of changing environmental conditions on SOC stocks more reliably and may reduce uncertainties found in both, geospatial and Earth System Models. In addition, the description of different environmental controllers for US ecoregions can help to describe the scope and importance of global and local models.

58 GEOSCIENCES↗
DOE OSTI2023

FORCE Integration with DRAFT and IDAES

Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to produce multiple commodities and improve economic viability under uncertain market conditions. Technical and economic analysis of IES requires modeling of complex processes with software models having enough fidelity to capture real-world dynamics while still being capable of running using reasonable computing resources. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Recent developments for the FORCE toolset have centered on strengthening the flexibility and modularity to link with external models and other available software to enhance IES simulations. Adding versatility to the FORCE toolset improves capability resulting in better techno-economic simulation of nuclear and IES components (both in potential higher fidelity and accuracy to expected performance after deployment). It also helps leverage existing work in the IES field and improve efficiency in national code development. This report focuses on two such endeavors: integration of the Dynamic Reliability Analysis Framework Tool (DRAFT) and the Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES) software packages into FORCE.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Atomic layer etching (ALE) of III-nitrides

Atomic layer etching (ALE) was performed on (Al, In, Ga)N thin films using a cyclic process of alternating Cl 2 gas absorption and Ar + ion bombardment in an inductively coupled plasma etcher system. The etch damage was characterized by comparing photoluminescence of blue single quantum well light-emitting diodes before and after the etch as well as bulk resistivities of etched p-doped layers. It was found that etched surfaces were smooth and highly conformal, retaining the step-terrace features of the as-grown surface, thus realizing ALE. Longer exposures to the dry etching increased the bulk resistivity of etched surfaces layers slightly, with a damaged depth of ∼55 nm. With further optimization and damage recovery, ALE is a promising candidate for controlled etching with atomic accuracy. Finally, it was found that Al 0.1 Ga 0.9 N acts as an etch barrier for the ALE etch, making it a suitable etch to reveal buried V-defects in III-nitride light emitting diodes.

Atomic-layer etching↗
DOE OSTI2026

Phonon modulation of strongly coupled gold tetrahedral plasmonic nanoparticles and a carbocyanine J-aggregate

Coupling exciton and plasmon excitations to form polaritons are of great interest for manipulating energy transfer at the nanoscale via the formation of hybrid light–matter states. In this study, we successfully couple gold tetrahedral nanoparticles with the J-aggregate forming dye 5,5′,6,6′-tetrachloro-1,1′-diethyl-3,3′-di(4-sulfobutyl)-benzimidazolocarbocyanine (TDBC) to form a strongly coupled colloidal polariton system with a Rabi splitting energy of ∼206 meV. These gold tetrahedra exhibit coherent phonon modes upon photoexcitation, which produce transient oscillations of the LSPR energy for isolated tetrahedra. Transient absorption measurements of the polariton system were performed and show how these coherent phonon modes influence the polariton states’ extinction. We found that the oscillation period increases by 0.5 ± 0.14 ps upon surface deposition of TDBC dye, demonstrating LSPR sensitivity to the refractive index environment. Shifts in the plasmon resonance due to the coherent phonon modes transiently alters LSPR alignment with the J-aggregate exciton peak, resulting in shifts of the hybrid polariton states’ oscillator strength.

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

The Madden–Julian Oscillation in the Energy Exascale Earth System Model Version 1

Abstract The present study examines the characteristics of the Madden–Julian Oscillation (MJO) events represented in the Energy Exascale Earth System Model version 1 (E3SMv1), DOE’s new Earth system model. The coupled E3SMv1 realistically simulates the eastward propagation of precipitation and moist static energy (MSE) anomalies associated with the MJO. As in observations, horizontal moisture advection and longwave radiative feedback are found to be the dominant processes in E3SMv1 that lead to the eastward movement and maintenance of the MJO MSE anomalies, respectively. Modulation of the diurnal cycle of precipitation in the Maritime Continent region by the MJO is also well represented in the model despite systematic biases in the magnitude and phase of the precipitation diurnal cycle. On the MJO impact over the midlatitude, E3SMv1 reasonably captures the pattern of the MJO teleconnections across the North Pacific and North America, with improvement in the performance in a high‐resolution version, despite the magnitude being a bit weaker than the observed. Regarding the interannual variability of the MJO, the El Niño‐Southern Oscillation (ENSO) modulation of the zonal extent of MJO’s eastward propagation, as well as associated changes in the mean state moisture gradient in the tropical west Pacific, are well reproduced in the model. However, MJO in E3SMv1 exhibits no sensitivity to the Quasi‐Biennial Oscillation (QBO), with the MJO propagation characteristics being almost identical between easterly QBO and westerly QBO years. Processes that have been suggested as critical to MJO simulation are also examined by utilizing recently developed process‐oriented diagnostics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Exchange‐Interaction‐Like Behavior in Ferroelectric Bilayers

Abstract Interlayer coupling in materials, such as exchange interactions at the interface between an antiferromagnet and a ferromagnet, can produce exotic phenomena not present in the parent materials. While such interfacial coupling in magnetic systems is widely studied, there is considerably less work on analogous electric counterparts (i.e., akin to electric “exchange‐bias‐like” or “exchange‐spring‐like” interactions between two polar materials) despite the likelihood that such effects can also engender new features associated with anisotropic electric dipole alignment. Here, electric analogs of such exchange interactions are reported, and their physical origins are explained for bilayers of in‐plane polarized Pb 1−x Sr x TiO 3 ferroelectrics. Variation of the strontium content and thickness of the layers provides for deterministic control over the switching properties of the bilayer system resulting in phenomena analogous to an exchange‐spring interaction and, leveraging added control of these interactions with an electric field, the ability to realize multistate‐memory function. Such observations not only hold technological promise for ferroelectrics and multiferroics but also extend the similarities between ferromagnetic and ferroelectric materials to include the manifestation of exchange‐interaction‐like phenomena.

Kavle, Pravin↗
DOE OSTI2020

Revealing the Structural Evolution of Green Rust Synthesized in Ionic Liquids by In Situ Molecular Imaging

Abstract Ionic liquids are green solvents that have wide applications in material synthesis, catalysis, and separation. A model switchable ionic liquid (SWIL) consisting of 1,8‐diazabicycloundec‐7‐ene (DBU) and 1‐hexanol with carbon dioxide (CO 2 ) as the trigger gas is chosen to synthesize nanocrystalline green rust. Under anoxic conditions, by adding iron (II) acetate (Fe(C 2 H 3 O 2 ) 2 ) and methanol to the degassed SWIL, a nanoparticulate green rust with carbonate (nano‐GR) is synthesized. More importantly, the molecular structural change of the ionic liquid resulting in green rust crystallines is first characterized using in situ liquid using time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) coupled with the system for analysis at the liquid vacuum interface (SALVI), a vacuum compatible microfluidic reactor in this work. Dimers and ionic clusters consisting ferric and ferrous ions are identified in liquid ToF‐SIMS spectra. Spectral principal component analysis is used to confirm that these ion pairs including dimers and cluster ions are different from nonionic liquids. The results show that liquid ToF‐SIMS can be a useful tool to study complex liquids at the molecular level providing insights into predicative synthesis of nanomaterials using environmentally friendly green solvents and giving unique visualization of the evolving material interface due to nucleation.

Shen, Yanjie↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2022

Decoding the dynamics of poleward shifting climate zones using aqua-planet model simulations

Growing evidence indicates that the atmospheric and oceanic circulation experiences a systematic poleward shift in a warming climate. However, the complexity of the climate system, including the coupling between the ocean and the atmosphere, natural climate variability and land-sea distribution, tends to obfuscate the causal mechanism underlying the circulation shift. Here, using an idealised coupled aqua-planet model, we explore the mechanism of the shifting circulation, by isolating the contributing factors from the direct CO 2 forcing, the indirect ocean surface warming, and the wind-stress feedback from the ocean dynamics. We find that, in contrast to the direct CO 2 forcing, ocean surface warming, in particular an enhanced subtropical ocean warming, plays an important role in driving the circulation shift. This enhanced subtropical ocean warming emerges from the background Ekman convergence of surface anomalous heat in the absence of the ocean dynamical change. Additionally, it expands the tropical warm water zone, causes a poleward shift of the mid-latitude temperature gradient, hence forces a corresponding shift in the atmospheric circulation and the associated wind pattern. The shift in wind, in turn drives a shift in the ocean circulation. Our simulations, despite being idealised, capture the main features of the observed climate changes, for example, the enhanced subtropical ocean warming, poleward shift of the patterns of near-surface wind, sea level pressure, storm tracks, precipitation and large-scale ocean circulation, implying that increase in greenhouse gas concentrations not only raises the temperature, but can also systematically shift the climate zones poleward.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Derivation of K-matrix reaction theory in a discrete basis formalism

The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann–Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present a simpler alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Strain-controlled magnetic ordering in 2D carbon metamaterials

Here, we use ab initio spin-polarized density functional theory to study the magnetic order in a Kagomé-like 2D metamaterial consisting of pristine or substitutionally doped phenalenyl radicals polymerized into a nanoporous, graphene-like structure. In this and in a larger class of related structures, the constituent polyaromatic hydrocarbon molecules can be considered as quantum dots that may carry a net magnetic moment. The structure of this porous system and the coupling between the quantum dots may be changed significantly by applying moderate strain, thus allowing to control the magnetic order and the underlying electronic structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An experimentally validated evaporative phase change heat transfer model for low mass flux applications using R134a in plate heat exchangers

Plate frame heat exchangers are common in liquid-coupled vapor compression systems due to their compactness and ease of maintenance. In some of these systems, the refrigerant within the evaporator can enter the heat exchanger as a subcooled liquid. As a result, the refrigerant passes through three different phases: single-phase liquid, two-phase fluid, and a superheated vapor. Unfortunately, there have been limited prior studies that have identified the best method to predict performance under these circumstances. Furthermore, previous investigations have not evaluated refrigerant evaporation in large industrial sized plate frame heat exchangers under conditions experienced in the present study experimental test facility, especially at mass fluxes below 7 kg m –2 s –1 . In the present investigation, a model was developed to predict the performance of plate and frame heat exchangers when the fluid enters as a subcooled liquid and exits as a superheated vapor. The model used full-sized plate heat exchanger geometry that was discretized into 20 sections to accurately capture local heat transfer and pressure drop effects. The model was validated using an R134a counter flow heat exchanger used in a turbo-compression cooling system test facility that had subcooled liquid entering the evaporator at low mass fluxes. A variety of empirical correlations were evaluated to determine what combination yielded the best predictive capability over the following range of conditions: 5.8 < G < 6.8 kg m –2 s –1 , 2 < q” < 2.8 kW m –2 , and 10 °C < T sat < 15 °C. When the correlations used a fixed control volume length in the heat flux calculation, the most accurate combination resulted in a mean absolute error of 5.5%. Finally, future studies can use the approach described here to optimize heat exchanger size and performance.

42 ENGINEERING↗
DOE OSTI2023

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SMR safety through HTTF modeling and benchmark efforts for code validation for gas-cooled reactor applications

Accurate modeling and simulation tools for thermal-hydraulics calculations are a key element needed to design and license new advanced reactors including Small Modular Reactors (SMR) and Microreactors. Uncertainties in modeling and simulation can have significant safety and economic implications. The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled prismatic-block nuclear reactors. High-quality measurement data is available from the HTTF that is suitable for a thermal-hydraulics code validation benchmark for gas-cooled reactor simulations. Here, this paper summarizes individual HTTF modeling efforts to date for tool validation at Idaho National Laboratory (INL), Argonne National Laboratory (ANL), Oregon State University (OSU) and Canadian Nuclear Laboratories (CNL) using system thermal-hydraulics codes, Computational Fluid Dynamics (CFD) codes and system-CFD code couplings. Also, the paper introduces the ongoing OECD Nuclear Energy Agency (NEA) High Temperature Gas Reactor Thermal-Hydraulics (HTGR T/H) benchmark that allows for better comparisons of results between different international modeling teams. The benchmark provides well defined computational problems that include code-to-code comparisons and comparisons to measured data. These problems provide an avenue for quantifying accuracy and identifying sources of uncertainty in thermal-hydraulics calculations, including in measured thermophysical properties, as part of validation for gas-cooled reactor simulation tools.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2020

Identification of the molecular-weight growth reaction network in counterflow flames of the C 3 H 4 isomers allene and propyne

The reaction networks responsible for aromatics formation in counterflow flames of the C 3 H 4 isomers allene and propyne are identified through a combined experimental and modeling study. Mole fraction profiles of near-atmospheric pressure (933 mbar) diffusion flames fueled by the C 3 H 4 isomers are analyzed by means of a newly assembled, chemically detailed kinetic mechanism. The experiment consists of a counterflow burner system that is coupled to a high-resolution time-of-flight molecular-beam mass spectrometer with single-photon ionization via synchrotron-generated vacuum-ultraviolet photons. Flame-sampled, mass-specific photoionization efficiency curves are used to identify the presence of aliphatically substituted aromatic species in addition to the commonly considered pericondensed ring structures. The new mechanism describes the formation and growth of aromatics through repetitive sequences of radical–radical and radical–molecule reactions that include C 1 6 intermediates. Higher concentrations of aromatic species are observed in the allene flame and the new mechanism captures the observed experimental trends very accurately. The results indicate the importance of the aliphatically substituted aromatics and of ring-enlargement reactions for the growth reactions. According to the model simulations, radical+radical recombination and PAH-radical+molecule reactions play an important role in PAH growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗