Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface roughness”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Timescales and Processes of Methane Hydrate Formation and Breakdown, With Application to Geologic Systems

Abstract Gas hydrate is an ice‐like form of water and low molecular weight gas stable at temperatures of roughly −10°C to 25°C and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one sixth of Earth's methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean‐atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean‐atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth's near‐surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some field studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

58 GEOSCIENCES↗
DOE OSTI2022

Interlayer spacing in pillared and grafted MCM-22 type silicas: density functional theory analysis versus experiment

Pillaring of synthetic layered crystalline silicates and aluminosilicates provides a strategy to enhance their adsorption and separation performance, and can facilitate the understanding of such behavior in more complex natural clays. In this study, we perform the first-principles density functional theory calculations for the pillaring of the pure silica polymorph of an MCM-22 type molecular sieve. Starting with a precursor material MCM-22P with fully hydroxylated layers, a pillaring agent, (EtO) 3 SiR, can react with hydroxyl groups (–OH) on adjacent internal surfaces, 2(–OH) + (EtO) 3 SiR + H 2 O → (–O) 2 SiOHR + 3EtOH, to form a pillar bridging these surfaces, or with a single hydroxyl, –OH + (EtO) 3 SiR + 2H 2 O → (–O)Si(OH) 2 R + 3EtOH, grafting to one surface. For computational efficiency, we replace the experimental organic ligand, R, by a methyl group. We find that the interlayer spacing in MCM-22 is reduced by 2.66 Å relative to weakly bound layers in the precursor MCM-22P. Including (–O) 2 SiR bridges for 50% (100%) of the hydroxyl sites in MCM-22P increases the interlayer spacing relative to MCM-22 by 2.52 Å (2.46 Å). For comparison, we also analyze the system where all –OH groups in MCM-22P are replaced by non-bridging grafted (–O)Si(OH) 2 R which results in a smaller interlayer spacing expansion of 2.17 Å relative to MCM-22. Our results for the interlayer spacing in the pillared materials are compatible with experimental observations for a similar MCM-22 type material with low Al content (Si : Al = 51 : 1) of an expansion relative to MCM-22 of roughly 2.8 Å and 2.5 Å from our x-ray diffraction and scanning transmission electron microscopy analyses, respectively. The latter analysis reveals significant variation in individual layer spacings.

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

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Brittle-to-Ductile Transition Temperature Decreases at the Microscale in Muscovite Mica

Characterizing the structural and chemical transitions of muscovite mica across various scales and thermal profiles is essential for ensuring the reliability of high-voltage power insulation, aerospace thermal barriers, and nanoscale semiconductor substrates under extreme operational stress. Here, we investigate notch fracture toughness $K^{n}_{Ic}$ of muscovite mica as a function of feature size w and temperature T with microscale shear testing. $K^{n}_{Ic}$ decreased and then stayed constant as w increased, likely from fabrication-induced changes to the microstructure and roughness. Moreover, $K^{n}_{Ic}$ initially decreased and then increased as temperature T increased, with the reversal in trend indicative of a brittle-to-ductile transition. Complementary Raman spectroscopy and molecular dynamics showed that the observed trends were due to an increase in bond lengths and a decrease in critical stresses/strains. In comparing the microscale results to macroscale counterparts, it was surmised that the decrease in brittle-to-ductile transition temperature was due to an increase in the dehydroxylation rate.

DelRio, Frank W. [Sandia National Laboratories (SN↗
DOE OSTI2024

Shallow Rate-Redox Potential Scaling in Aqueous Molecular Oxygen Reduction Electrocatalysis Across a Family of Iron Macrocycles

Rate-overpotential scaling relationships have been employed widely to understand trends in oxygen reduction reaction (ORR) electrocatalysis by dissolved metal macrocycles in organic electrolytes. Similar scaling relationships remain unknown for surface-adsorbed ORR electrocatalysts in the acidic aqueous environments germane to proton-exchange membrane (PEM) fuel cells. Herein, we examine ORR catalysis in aqueous perchloric acid media for a structurally diverse array of iron macrocycle complexes adsorbed on Vulcan carbon black. The macrocycles encompass Fe– N 4 , Fe–N 2 N' 2 and Fe–N x C 4-x motifs bearing pyrrolic, pyridinic, and N-heterocyclic carbene (NHC) moieties in the primary ligation sphere, giving rise to a 670 mV range in Fe(III/II) redox potentials, E Fe(III/II) . Experimental Tafel data in the micropolarization regime were extrapolated to the E Fe(III/II) to furnish estimated per-site-normalized current density (j per-site ) values that span ~4.6 orders of magnitude across the family of compounds. Despite the structural diversity of this family of compounds, extrapolated j per-site values correlate with the Fe(III/II) redox potentials in a roughly log-linear fashion with a shallow scaling factor of approximately 145 mV/decade. Further, these findings highlight that negative shifts in E Fe(III/II) lead to diminishing returns in catalytic rate promotion and suggest that changes to the primary ligating environment in a macrocycle are insufficient to break fundamental rate-potential scaling relationships in aqueous ORR catalysis. Together these studies motivate the further development of higher-potential iron complexes that employ motifs beyond the equatorial ligation plane to enhance ORR catalysis.

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

Saturation and decay of an induced scattering instability driven by static field asymmetries in a non-neutral plasma

Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Functional variability in specific root respiration translates to autotrophic differences in soil respiration in a temperate deciduous forest

CO 2 release from forest soils (R s ) is a prominent flux in the global carbon cycle. Rs is derived from roots (autotrophic respiration, R a ) and microbial (heterotrophic) respiration and is highly dynamic, as it depends on edaphic and environmental conditions as well as root functional traits and microbial community composition. It is unclear how root functional traits affect root and microbial respiration rates; however, their consideration may help parse out the relative contributions of root and microbial respiration to R s . At a temperate forest site, root systems of 3–4 functional root orders and their surrounding surface soil were carefully excavated and placed into custom trays designed to repeatedly measure R s in situ on eight temperate tree species that varied in their root functional strategies and mycorrhizal affinity. R s was measured bi-weekly to monthly for nearly one year using a custom chamber attached to a gas exchange system. R s varied over time, ranging from 0.3 to 12 µmol m -2 s -1 . Comparable root systems of the same species were excised from the soil and specific root respiration rates (R r ) were measured. Rr ranged from 2.5 to 9.0 nmol g -1 s -1 and was negatively correlated with root tissue density and positively related to root tissue nitrogen concentration. Using R r to estimate R a , we estimate that R a accounts for <10%, on average 2–3%, of R s for individual root systems (averaging 1.2 g dry biomass) housed in surrounding soil (average 1.3 kg dry mass) in situ; thus, Ra was roughly 20 times greater than Rh per unit mass. The contribution of R a peaked in the fall and coincided with leaf senescence of the forest canopy. A soil-sterilizing experimental treatment designed to help isolate R a in situ reduced bacterial biomass and shifted fungal community composition, but there was no reduction in Rs of the in-situ root-soil tray systems. The relative R a to R s ratio increased with root functional strategies characterized by greater specific root length and tip abundance, but also to greater root tissue density. The ratio of R a to R s also increased with warmer soil temperatures and decreased slightly with increasing soil moisture. We discuss how incorporating root functional traits as modulators of the autotrophic contribution to R s could be considered when modeling total soil CO 2 efflux from forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Superconductivity and magnetic and transport properties of single-crystalline CaK ( Fe 1 – x Cr x ) 4 As 4

Members of the CaK (Fe 1–x Cr x ) 4 As 4 series have been synthesized by high-temperature solution growth in single-crystalline form and characterized by x-ray diffraction, elemental analysis, and magnetic and transport measurements. The effects of Cr substitution on the superconducting and magnetic ground states of CaKFe 4 As 4 (T c = 35 K) have been studied. These measurements show that the superconducting transition temperature decreases monotonically and is finally suppressed below 1.8 K as x is increased from 0 to 0.038. For x-values greater than 0.012, signatures of a magnetic transition can be detected in magnetic measurements with the associated features in the transport measurements becoming detectable for x ≥ 0.038. The magnetic transition temperature increases in a roughly linear manner as Cr substitution increases. A temperature-composition (T – x) phase diagram is constructed, revealing a half-dome of superconductivity with the magnetic transition temperature, T*, appearing near 22 K for x ~ 0.017 and rising slowly up to 60 K for x ~ 0.077. The T – x phase diagrams for CaK (Fe 1–x T x ) 4 As 4 for T = Cr and Mn are essentially the same despite the nominally different band filling; this is in marked contrast to T = Co and Ni series for which the T – x diagrams scale by a factor of 2, consistent with the different changes in band filling Co and Ni would produce when replacing Fe. Superconductivity of CaK (Fe 1–x Cr x ) 4 As 4 is also studied as a function of magnetic field. A clear change in $H^{'}_{c2}$ (T)/T c , where $H^{'}_{c2}$ (T) is dH c2 (T)/dT, at x ~ 0.012 is observed and probably is related to a change of the Fermi surface due to magnetic order. Further, coherence length and the London penetration depths are also calculated based on H c1 and H c2 data. Both of them as a function of x show changes near x = 0.012, again consistent with Fermi surface changes associated with the magnetic ordering seen for higher x-values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Modeling DOE Standard Canister Configurations with Updated Surface Chemistry

Road-ready and final disposition packaging configurations for the advanced test reactor (ATR) fuel currently specifies storage within helium backfilled DOE sealed standard canisters. The aluminum cladding of the ATR fuel contains an oxyhydroxide layer of boehmite/bayerite that generates hydrogen when subjected to irradiation. Understanding the effect of this hydrogen buildup over time to important for long term storage considerations. Previous modeling efforts have built a coupled CFD-chemical model to simulate the temperature gas phase concentrations within the DOE sealed standard canisters. These models have coupled the temperature conditions to both the gas phase radiolysis chemistry and the surface chemistry associated with the oxyhydroxide layer. Continued experimental work has identified trends for aluminum surrogate samples with oxyhydroxide layers for tests done at higher dose rates. At low initial doses a fast generation rate of hydrogen occurs which starts to roll over to a lower generation rate as the total dose applied in increased. A previous study created a small-scale chemical model was built to replicate a mini-canister surrogate system at SRNL as well as for the smaller capsule tests performed at INL. The previous iteration of the model utilized step function for its G-values for the hydrogen generation over a 200-year period. This model makes an update to the surface chemistry to account for theorized surface chemistry reactions allowing for oxygen to remain bounded to the oxyhydroxide layer. As additional experiments continue, the kinetic fits may be adjusted to adapt to new data. Three canister configurations are modeled – the base 18-inch, 15-foot DOR standard canister with 30 ATR fuel elements, an 18-inch, 10-foot DOE standard canister with 32 ATR fuel elements and a 24-inch, 10-foot DOE standard canister with 40 ATR fuel elements. In previous reports, the primary sensitivity of the canister conditions was identified as the decay heat of the fuel and the dried condition of the fuel. Only these parameters are studied in the report. For the nominal case with dried fuel, the hydrogen generation is only about 2%, so it is even less than the flammability limit if it were exposed to oxygen. For the densely packed 18-inch case the total hydrogen concentration is only around 4% with the nominal decay heat. For the densely packed 24-inch case the total hydrogen concentration is only 2.5%, roughly 20% higher than the original packing design, and still under flammability conditions if exposed to oxygen.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Tropical extreme droughts drive long-term increase in atmospheric CO 2 growth rate variability

The terrestrial carbon sink slows the accumulation of carbon dioxide (CO 2 ) in the atmosphere by absorbing roughly 30% of anthropogenic CO 2 emissions, but varies greatly from year to year. The resulting variations in the atmospheric CO 2 growth rate (CGR) have been related to tropical temperature and water availability. The apparent sensitivity of CGR to tropical temperature ( ${{{{{{\rm{\gamma }}}}}}}_{{{{{{\rm{CGR}}}}}}}^{{{{{{\rm{T}}}}}}}$ ) has changed markedly over the past six decades, however, the drivers of the observation to date remains unidentified. Here, we use atmospheric observations, multiple global vegetation models and machine learning products to analyze the cause of the sensitivity change. We found that a threefold increase in ${{{{{{\rm{\gamma }}}}}}}_{{{{{{\rm{CGR}}}}}}}^{{{{{{\rm{T}}}}}}}$ emerged due to the long-term changes in the magnitude of CGR variability (i.e., indicated by one standard deviation of CGR; STD CGR ), which increased 34.7% from 1960-1979 to 1985-2004 and subsequently decreased 14.4% in 1997-2016. We found a close relationship (r 2 = 0.75, p < 0.01) between STD CGR and the tropical vegetated area (23°S – 23°N) affected by extreme droughts, which influenced 6-9% of the tropical vegetated surface. A 1% increase in the tropical area affected by extreme droughts led to about 0.14 Pg C yr -1 increase in STD CGR . The historical changes in STD CGR were dominated by extreme drought-affected areas in tropical Africa and Asia, and semi-arid ecosystems. The outsized influence of extreme droughts over a small fraction of vegetated surface amplified the interannual variability in CGR and explained the observed long-term dynamics of ${{{{{{\rm{\gamma }}}}}}}_{{{{{{\rm{CGR}}}}}}}^{{{{{{\rm{T}}}}}}}$ .

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2021

The De-Icing Comparison Experiment (D-ICE): a study of broadband radiometric measurements under icing conditions in the Arctic

Abstract. Surface-based measurements of broadband shortwave (solar) and longwave (infrared) radiative fluxes using thermopile radiometers are made regularly around the globe for scientific and operational environmental monitoring. The occurrence of ice on sensor windows in cold environments – whether snow, rime, or frost – is a common problem that is difficult to prevent as well as difficult to correct in post-processing. The Baseline Surface Radiation Network (BSRN) community recognizes radiometer icing as a major outstanding measurement uncertainty. Towards constraining this uncertainty, the De-Icing Comparison Experiment (D-ICE) was carried out at the NOAA Atmospheric Baseline Observatory in Utqiaġvik (formerly Barrow), Alaska, from August 2017 to July 2018. The purpose of D-ICE was to evaluate existing ventilation and heating technologies developed to mitigate radiometer icing. D-ICE consisted of 20 pyranometers and 5 pyrgeometers operating in various ventilator housings alongside operational systems that are part of NOAA's Barrow BSRN station and the US Department of Energy Atmospheric Radiation Measurement (ARM) program North Slope of Alaska and Oliktok Point observatories. To detect icing, radiometers were monitored continuously using cameras, with a total of more than 1 million images of radiometer domes archived. Ventilator and ventilator–heater performance overall was skillful with the average of the systems mitigating ice formation 77 % (many >90 %) of the time during which icing conditions were present. Ventilators without heating elements were also effective and capable of providing heat through roughly equal contributions of waste energy from the ventilator fan and adiabatic heating downstream of the fan. This provided ∼0.6 ∘C of warming, enough to subsaturate the air up to a relative humidity (with respect to ice) of ∼105 %. Because the mitigation technologies performed well, a near complete record of verified ice-free radiometric fluxes was assembled for the duration of the campaign. This well-characterized data set is suitable for model evaluation, in particular for the Year of Polar Prediction (YOPP) first Special Observing Period (SOP1). We used the data set to calculate short- and long-term biases in iced sensors, finding that biases can be up to +60 W m−2 (longwave) and −211 to +188 W m−2 (shortwave). However, because of the frequency of icing, mitigation of ice by ventilators, cloud conditions, and the timing of icing relative to available sunlight, the biases in the monthly means were generally less than the aggregate uncertainty attributed to other conventional sources in both the shortwave and longwave.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2023

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

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

A Modified Vegetation Photosynthesis and Respiration Model (VPRM) for the Eastern USA and Canada, Evaluated With Comparison to Atmospheric Observations and Other Biospheric Models

Atmospheric CO 2 measurements from a dense surface network can help to evaluate terrestrial biosphere model (TBM) simulations of Net Ecosystem Exchange (NEE) with two key benefits. First, gridded CO 2 flux estimates can be evaluated over regional scales, not possible using flux tower observations at discrete locations for model evaluation. Second, TBM ability to explain atmospheric CO 2 fluctuations due to the biosphere can be directly tested, an important objective for anthropogenic emissions monitoring using atmospheric observations. Here, we customize the Vegetation Photosynthesis and Respiration Model (VPRM) for an eastern North American domain with strong biological activity upwind of urban areas. Parameters are optimized using flux tower observations from a historical database with sites in (and near) the domain. In addition, the respiration model (originally a linear function of temperature) is modified to account for impacts of changing foliage, non-linear temperature, and water stress. Flux estimates from VPRM, the Carnegie-Ames-Stanford Approach (CASA) model and the Simple Biosphere Model v4 (SiB4), are convolved with footprints from atmospheric transport models for evaluation with CO 2 observations at 21 towers in the domain, with roughly half of the towers used here for the first time. Results show that the new respiration model in VPRM helps to correct a growing season sink bias in the atmosphere associated with underestimated summertime respiration using the original model with annual parameters. The new VPRM also better explains fine-scale atmospheric CO 2 variability compared to other TBMs, due to higher resolution diagnostic phenology, the new respiration model, domain-specific parameters, and high-quality input data sets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.

54 ENVIRONMENTAL SCIENCES↗