Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion model”

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,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A Coupled Soil-Atmosphere Model of H2O2 on Mars

The Viking Gas Chromatograph Mass Spectrometer failed to detect organic compounds on Mars, and both the Viking Labeled Release and the Viking Gas Exchange experiments indicated a reactive soil surface. These results have led to the widespread belief that there are oxidants in the martian soil. Since H2O2 is produced by photochemical processes in the atmosphere of Mars, and has been shown in the laboratory to reproduce closely the Viking LR results, it is a likely candidate for a martian soil oxidant. Here, we report on the results of a coupled soil/atmosphere transport model for H202 on Mars. Upon diffusing into the soil, its concentration is determined by the extent to which it is adsorbed and by the rate at which it is catalytically destroyed. An analytical model for calculating the distribution of H202 in the martian atmosphere and soil is developed. The concentration of H202 in the soil is shown to go to zero at a finite depth, a consequence of the nonlinear soil diffusion equation. The model is parameterized in terms of an unknown quantity, the lifetime of H202 against heterogeneous catalytic destruction in the soil. Calculated concentrations are compared with a H202 concentration of 30 nmoles/cu cm, inferred from the Viking Labeled Release experiment. A significant result of this model is that for a wide range of H202 lifetimes (up to 105 years), the extinction depth was found to be less than 3 m. The maximum possible concentration in the top 4 cm is calculated to be approx. 240 nmoles/cu cm, achieved with lifetimes of greater than 1000 years. Concentrations higher than 30 nmoles/cu cm require lifetimes of greater than 4.3 terrestrial years. For a wide range of H202 lifetimes, it was found that the atmospheric concentration is only weakly coupled with soil loss processes. Losses to the soil become significant only when lifetimes are less than a few hours. If there are depths below which H202 is not transported, it is plausible that organic compounds, protected from an oxidizing environment, may still exist. They would have been deposited by meteors, or be the organic remains of past life.

Bullock, Mark A.↗
NASA NTRS1983

The problem of cooling the cold Io torus

Two models are developed for the transport of ions inward from Io's orbit, the first using the assumption that radial diffusion is the dominant transport mechanism while the second uses a combination of diffusive and convective transport. The models include thermal and number density transport, radiation, ionization and pickup of local neutrals, recombination, charge exchange, and Coulomb interactions. It is found that pure diffusive transport can account for the dramatic inward depletion of the torus only by invoking recombination or by postulating a massive increase in the production rate of torus ions sometime prior to the Voyager mission. Recombination is determined to be far too slow to be the cause of the observed density decrease inside of Io. The model combining convection and diffusion is shown to reasonably match the data, but only if a diffusion coefficient 100 times less than that derived from Pioneer observations is used. In addition, it is shown that the Pioneer derived diffusion rate combined with Voyager temperature and density measurements indicate a large nonradiative sink of energy in the inner torus.

Richardson, J. D.↗
DOE OSTI2017

Experimental and Theoretical Analysis of Ink Dispersion Stability for Polymer Electrolyte Fuel Cell Applications

The aggregate size in fuel cell catalyst inks depends on the type of dispersion medium, particle concentration, and addition of stabilizing agents. In this work, ink stability and particle size of carbon black and carbon black/Nafion dispersions in four non-aqueous media, viz., methanol, ethanol, isopropanol and ethyl acetate are studied. Based on visual inspection, isopropanol is found to be the best medium for dispersion of carbon black inks. To rationalize this observation, a semi-empirical model based on diffusion-limited aggregation was developed to evaluate the rate of particle aggregation and predict the ink stability time for each dispersion medium. The proposed model supports the experimental observation by qualitatively predicting the same relationship between carbon stability and the dispersion media. The model also showed that the dielectric constant of the dispersion medium and the particle zeta potential are primarily responsible for the ink stability. Particle size for the different inks was determined by dynamic light scattering with and without dilution. Experimental results show that Nafion is a strong stabilizing agent, increasing the ink stability and decreasing the particle size of carbon aggregates. The beneficial effects of Nafion are independent of its concentration and are observed even at Nafion volume fractions of 10 wt%. The interaction energy is found to be a strong function of the surface potential for the dispersion medium with a higher dielectric constant.

25 ENERGY STORAGE↗
NASA NTRS1993

Near-contact diffusion and compensation in extrinsic photoconductors

Carrier diffusion from heavily doped contact regions has been experimentally monitored with resistivity measurements of variable thickness Si:P extrinsic photoconductors (n(+) -n-n(+)). Modeling of free carrier diffusion at the interface of heavily doped contacts and highly resistive bulk at low temperatures predicts extended regions (5-30 microns) of excess carriers in high purity materials so that, in thin device structures, free carrier diffusion profiles from each contact will overlap and determine the resistivity of the device. In this work, a decrease in resistivity of four orders of magnitude was observed in a 5-micron thick structure compared to a 10-micron thick device. The resistivity of an ohmic structure in the thin limit is strongly dependent on the bulk and near-contact compensation, and resistivity measurements can be used as a sensitive measure of compensation at interfaces or in the tail of implanted layers.

Olsen, C. S.↗
NASA NTRSDate not supplied

Modeling the Uptake of Hydrogen Chloride onto Interior Spacecraft Materials

Hydrogen chloride (HCl) is a major combustion product from the pyrolysis of polyvinyl chloride (PVC) insulated electrical wires, a common spacecraft fire safety concern. Models at two different scales were developed to predict HCl uptake on anodized, chromate conversion coated (Iridite), and bare aluminum surfaces, as well as on Nomex fabric: a macroscopic one-step global surface reaction model where all the active sites are on the exterior surface, and a pore model where the interior active sites deeper into the oxide layer can also be accessed by HCl. Experiments were performed to calibrate kinetic and diffusion constants in the models. A cast acrylic test cell was used to measure the differences between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. For the materials with a thin (< 200 µm) or no oxide layer, the macroscopic surface reaction model adequately predicts the experimental measurements. For the anodized aluminum with a thicker oxide layer, the pore model provided a better match to experimental results. The results will be discussed with respect to the spacecraft fire safety project (Saffire).

modeling↗
DOE OSTI2020

Multilevel Convergence Analysis of Multigrid-Reduction-in-Time

This study presents a multilevel convergence framework for multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid estimates. The framework provides a priori upper bounds on the convergence of MGRIT V- and F-cycles, with different relaxation schemes, by deriving the respective residual and error propagation operators. The residual and error operators are functions of the time-stepping operator, analyzed directly and bounded in the norm, both numerically and analytically. We present various upper bounds of different computational cost and varying sharpness. These upper bounds are complemented by proposing analytic formulae for the approximate convergence factor of V-cycle algorithms that take the number of fine grid time points, the temporal coarsening factors, and the eigenvalues of the time-stepping operator as parameters. The paper concludes with supporting numerical investigations of parabolic (anisotropic diffusion) and hyperbolic (wave equation) model problems. We assess the sharpness of the bounds and the quality of the approximate convergence factors. Observations from these numerical investigations demonstrate the value of the proposed multilevel convergence framework for estimating MGRIT convergence a priori and for the design of a convergent algorithm. We further highlight that observations in the literature are captured by the theory, including that two-level Parareal and multilevel MGRIT with F-relaxation do not yield scalable algorithms and the benefit of a stronger relaxation scheme. An important observation is that with increasing numbers of levels MGRIT convergence deteriorates for the hyperbolic model problem, while constant convergence factors can be achieved for the diffusion equation. The theory also indicates that L-stable Runge--Kutta schemes are more amendable to multilevel parallel-in-time integration with MGRIT than A-stable Runge--Kutta schemes.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1987

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗
NASA NTRS1992

Self-diffusion of magnesium in spinel and in equilibrium melts - Constraints on flash heating of silicates

An isotopic tracer is used to measure Mg self-diffusion in spinel and coexisting melt at bulk chemical equilibrium. The diffusion coefficients were calculated from the measured isotope profiles using a model that includes the complementary diffusion of Mg-24, Mg-25, and Mg-26 in both phases with the constraint that the Mg content of each phase is constant. The activation energy and preexponential factor for Mg self-diffusion in spinel are, respectively, 384 +/- 7 kJ and 74.6 +/- 1.1 sq cm/s. These data indicate Mg diffusion in spinel is much slower than previous estimates. The activation energy for Mg self-diffusion in coexisting melt is 343 +/- 25 kJ and the preexponential factor is 7791.9 +/- 1.3 sq cm/s. These results are used to evaluate cooling rates of plagioclase-olivine inclusions (POIs) in the Allende meteorite. Given a maximum melting temperature for POIs of about 1500 C, these results show that a 1-micron radius spinel would equilibrate isotopically with a melt within about 60 min.

Sheng, Y. J.↗
NASA NTRS1995

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗
NASA NTRS2000

Photochemical Modeling of CH3 Abundances in the Outer Solar System

Recent measurements of methyl radicals (CH3) in the upper atmospheres of Saturn and Neptune by the Infrared Space Observatory (ISO) provide new constraints to photochemical models of hydrocarbon chemistry in the outer solar system. The derived column abundances of CH3 on Saturn above 10 mbar and Neptune above the 0.2 mbar pressure level are (2.5 - 6.0) x 10(exp 13) / sq cm and (0.7 - 2.8) x 10(exp 13) / sq cm, respectively. We use the updated Caltech/Jet Propulsion Laboratory photochemical model, which incorporates hydrocarbon photochemistry, vertical molecular and bulk atmospheric eddy diffusion, and realistic radiative transfer modeling, to study the CH3 abundances in the upper atmosphere of the giant planets and Titan. We identify the key reactions that control the concentrations of CH3 in the model, such as the three-body recombination reaction, CH3 + CH3 + M yields C2H6 + M. We evaluate and extrapolate the three-body rate constant of this reaction to the low-temperature limit (1.8 x 10(exp -16) T(sup -3.75) e(sup -300/T), T < 300 K) and compare methyl radical abundances in five atmospheres: Jupiter, Saturn, Uranus, Neptune, and Titan. The sensitivity of our models to the rate coefficients for the reactions H + CH3 + M yields CH4 + M, H + C2H3 yields C2H2 + H2, (sup 1)CH2 + H2 yields CH3 + H, and H + C2H5 yields 2CH3, the branching ratios of CH4 photolysis, vertical mixing in the five atmospheres, and Lyman alpha photon enhancement at the orbit of Neptune have all been tested. The results of our model CH3 abundances for both Saturn (5.1 x 10(exp 13) / sq cm) and Neptune (2.2 x 10(exp 13) / sq cm) show good agreement with ISO Short Wavelength Spectrometer measurements. Using the same chemical reaction set, our calculations also successfully generate vertical profiles of stable hydrocarbons consistent with Voyager and ground-based measurements in these outer solar system atmospheres. Predictions of CH3 column concentrations (for p <= 0.2 mbar) in the atmospheres of Jupiter (3.3 x 10(exp 13) /sq cm), Uranus (2.5 x 10(exp 12) / sq cm), and Titan (1.9 x 10(exp 15) / sq cm) may be checked by future observations.

Lee, Anthony Y. T.↗
NASA NTRS2012

Quasi 1D Modeling of Mixed Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of the 2-dimensional bifurcated mixed compression supersonic inlet is being developed. The model utilizes computational fluid dynamics for both the supersonic and subsonic diffusers. The oblique shocks are modeled utilizing compressible flow equations. This model also implements variable geometry required to control the normal shock position. The model is flexible and can also be utilized to simulate other mixed compression supersonic inlet designs. The model was validated both in time and in the frequency domain against the legacy LArge Perturbation INlet code, which has been previously verified using test data. This legacy code written in FORTRAN is quite extensive and complex in terms of the amount of software and number of subroutines. Further, the legacy code is not suitable for closed loop feedback controls design, and the simulation environment is not amenable to systems integration. Therefore, a solution is to develop an innovative, more simplified, mixed compression inlet model with the same steady state and dynamic performance as the legacy code that also can be used for controls design. The new nonlinear dynamic model is implemented in MATLAB Simulink. This environment allows easier development of linear models for controls design for shock positioning. The new model is also well suited for integration with a propulsion system model to study inlet/propulsion system performance, and integration with an aero-servo-elastic system model to study integrated vehicle ride quality, vehicle stability, and efficiency.

Kopasakis, George↗
DOE OSTI2020

Multimodel simulation of vertical gas transfer in a temperate lake

In recent decades, several lake models of varying complexity have been developed and incorporated in numerical weather prediction systems and climate models. To foster enhanced forecasting ability and verification, improvement of these lake models remains essential. This especially applies to the limited simulation capabilities of biogeochemical processes in lakes and greenhouse gas exchanges with the atmosphere. Here we present multi-model simulations of physical variables and dissolved gas dynamics in a temperate lake (Harp Lake, Canada). The five models (ALBM, FLake, LAKE, LAKEoneD, MTCR-1) considered within this most recent round of the Lake Model Intercomparison Project (LakeMIP) all captured the seasonal temperature variability well. In contrast, none of the models is able to reproduce the exact dates of ice-cover and ice-off, leading to considerable errors in the simulation of eddy diffusivity around those dates.We then conducted an additional modeling experiment with a diffusing passive tracer to isolate the effect of the eddy diffusivity on gas concentration. Remarkably, sophisticated k -" models do not demonstrate a significant difference in the vertical diffusion of a passive tracer compared to models with much simpler turbulence closures. All models simulate less intensive spring overturn compared to autumn. Reduced mixing in the models consequently leads to the accumulation of the passive tracer distribution in the water column. The lake models with a comprehensive biogeochemical module, such as ALBM and LAKE, predict dissolved oxygen dynamics adequate to the observed data. However, for the surface carbon dioxide concentration the correlation between modeled (ALBM, LAKE) and observed data is weak (~0.3). Overall our results indicate the need to improve the representation of physical and biogeochemical processes in lake models, thereby contributing to enhanced weather prediction and climate projection capabilities.

Guseva, Sofya↗
DOE OSTI2019

Predictions of diffusion rates of large organic molecules in secondary organic aerosols using the Stokes–Einstein and fractional Stokes–Einstein relations

Information on the rate of diffusion of organic molecules within secondary organic aerosol (SOA) is needed to accurately predict the effects of SOA on climate and air quality. Diffusion can be important for predicting the growth, evaporation, and reaction rates of SOA under certain atmospheric conditions. Often, researchers have predicted diffusion rates of organic molecules within SOA using measurements of viscosity and the Stokes–Einstein relation ($D∝1/η$, where D is the diffusion coefficient and η is viscosity). However, the accuracy of this relation for predicting diffusion in SOA remains uncertain. Using rectangular area fluorescence recovery after photobleaching (rFRAP), we determined diffusion coefficients of fluorescent organic molecules over 8 orders in magnitude in proxies of SOA including citric acid, sorbitol,and a sucrose–citric acid mixture. These results were combined with literature data to evaluate the Stokes–Einstein relation for predicting the diffusion of organic molecules in SOA. Although almost all the data agree with the Stokes–Einstein relation within a factor of 10, a fractional Stokes–Einstein relation ($D∝1/η^ξ$) with ξ =0.93 is a better model for predicting the diffusion of organic molecules in the SOA proxies studied. In addition, based on the output from a chemical transport model, the Stokes–Einstein relation can overpredict mixing times of organic molecules within SOA by as much as 1 order of magnitude at an altitude of ~3 km compared to the fractional Stokes–Einstein relation with ξ =0.93. These results also have implications for other areas such as in food sciences and the preservation of biomolecules.

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

Numerical Simulations of Cathodoluminescence Measurements in Thin-Film Solar Cells

Cathodoluminescence (CL) measurements can be applied to assess grain-boundary (GB) and grain-interior (GI) recombination in thin-film solar cell materials and made quantitative if we can develop CL models that account for material and measurement complexities. Recently, we developed a three-dimensional numerical CL model, based in MATLAB, that simulates the GI CL intensity as a function of four parameters: grain size, GI lifetime, and GB and surface recombination velocities. The model assumes that GB electrostatic potentials are screened by the high excesscarrier densities used in the CL measurement such that transport is governed by ambipolar diffusion. Here, we develop models to address directly GB potentials and their effects on these measurements. First, we transfer the MATLAB-based model to COMSOL software, and then introduce shallow donors to the GBs to produce potentials. We also develop a two-dimensional model in COMSOL to simulate CL GB contrast with GB potentials. Simulations indicate that GB potentials can increase or decrease CL intensities relative to the zero-potential case. However, the high electron-beam currents typically applied in CL measurements minimize the impact of GB potentials.

14 SOLAR ENERGY↗
NASA NTRS1996

Flamelet Model Application for Non-Premixed Turbulent Combustion

The current Final Report contains results of the study which was performed in Scientific Research Center 'ECOLEN' (Moscow, Russia). The study concerns the development and verification of non-expensive approach for modeling of supersonic turbulent diffusion flames based on flamelet consideration of the chemistry/turbulence interaction (FL approach). Research work included: development of the approach and CFD tests of the flamelet model for supersonic jet flames; development of the simplified procedure for solution of the flamelet equations based on partial equilibrium chemistry assumption; study of the flame ignition/extinction predictions provided by flamelet model. The performed investigation demonstrated that FL approach allowed to describe satisfactory main features of supersonic H 2/air jet flames. Model demonstrated also high capabilities for reduction of the computational expenses in CFD modeling of the supersonic flames taking into account detailed oxidation chemistry. However, some disadvantages and restrictions of the existing version of approach were found in this study. They were: (1) inaccuracy in predictions of the passive scalar statistics by our turbulence model for one of the considered test cases; and (2) applicability of the available version of the flamelet model to flames without large ignition delay distance only. Based on the results of the performed investigation, we formulated and submitted to the National Aeronautics and Space Administration our Project Proposal for the next step research directed toward further improvement of the FL approach.

Secundov, A.↗
DOE OSTI2018

Laboratory Photoionization Fronts in Nitrogen Gas: A Numerical Feasibility and Parameter Study

Photoionization fronts play a dominant role in many astrophysical situations but remain difficult to achieve in a laboratory experiment. We present the results from a computational parameter study evaluating the feasibility of the photoionization experiment presented in the design paper by Drake et al. in which a photoionization front is generated in a nitrogen medium. The nitrogen gas density and the Planckian radiation temperature of the X-ray source define each simulation. Simulations modeled experiments in which the X-ray flux is generated by a laser-heated gold foil, suitable for experiments using many kJ of laser energy, and experiments in which the flux is generated by a "z-pinch" device, which implodes a cylindrical shell of conducting wires. The models are run using crash, our block-adaptive-mesh code for multimaterial radiation hydrodynamics. The radiative transfer model uses multigroup, flux-limited diffusion with 30 radiation groups. In addition, electron heat conduction is modeled using a single-group, flux-limited diffusion. In the theory, a photoionization front can exist only when the ratios of the electron recombination rate to the photoionization rate and the electron-impact ionization rate to the recombination rate lie in certain ranges. These ratios are computed for several ionization states of nitrogen. Photoionization fronts are found to exist for laser-driven models with moderate nitrogen densities (~10 21 cm –3 ) and radiation temperatures above 90 eV. For "z-pinch"-driven models, lower nitrogen densities are preferred (<10 21 cm –3 ). Here, we conclude that the proposed experiments are likely to generate photoionization fronts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Borophene Synthesis on Au(111)

Borophene (the first two-dimensional (2D) allotrope of boron) is emerging as a groundbreaking system for boron-based chemistry and, more broadly, the field of low-dimensional materials. Exploration of the phase space for growth is critical because borophene is a synthetic 2D material that does not have a bulk layered counterpart and thus cannot be isolated via exfoliation methods. Furthermore, we report synthesis of borophene on Au(111) substrates. Unlike previously studied growth on Ag substrates, boron diffuses into Au at elevated temperatures and segregates to the surface to form borophene islands as the substrate cools. These observations are supported by ab initio modeling of interstitial boron diffusion into the Au lattice. Borophene synthesis also modifies the surface reconstruction of the Au(111) substrate, resulting in a trigonal network that templates growth at low coverage. This initial growth is composed of discrete borophene nanoclusters, whose shape and size are consistent with theoretical predictions. As the concentration of boron increases, nanotemplating breaks down and larger borophene islands are observed. Spectroscopic measurements reveal that borophene grown on Au(111) possesses a metallic electronic structure, suggesting potential applications in 2D plasmonics, superconductivity, interconnects, electrodes, and transparent conductors.

36 MATERIALS SCIENCE↗