Numerical modeling of fiber orientation in multi-layer, isothermal material-extrusion big area additive manufacturing
Not Available
36 MATERIALS SCIENCE↗
SEARCH · Engineering Papers
Results for “ISOTHERMAL LAYER”
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.
Global fluid turbulence simulations in the scrape-off layer of a stellarator island divertor
Isothermal fluid turbulence simulations have been performed in the edge and scrape-off layer (SOL) of an analytic stellarator configuration with an island divertor, thereby providing numerical insight into edge turbulence in regions around islands in a stellarator. The steady-state transport follows the a curvature drive that is inverse to the major radius ( $1/R$ ) toward the outboard side, but large fluctuations are present throughout the island divertor region, with the average wavelength of similar size to the island width. The system exhibits a prominent $m=2$ , $n=5$ mode, where m is the poloidal mode number and n is the toroidal mode number, although other modes are present. The amplitude and radial extent of the density fluctuations are similar throughout the edge and SOL, but can decrease near island O-points. The fluctuations exhibit a predominantly positive skewness on the outboard midplane, indicating blob-like perturbations for the transport into the outer SOL. It is determined that a point on the separatrix is generally more correlated with regions outside of the SOL than a nearby reference point which does not lie on the separatrix.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Growth of ZrO 2 films on mesoporous silica sieve via atomic layer deposition
The atomic layer deposition of ZrO 2 thin films over the mesopores of a mesoporous silica sieve, SBA-15 (Santa Barbara Amorphous-15), using the tetrakis(dimethylamino)zirconium(IV) (TDMAZ) as ZrO 2 precursor was tested and characterized using N 2 adsorption-desorption isotherms, taking advantage of the well-defined shape and size distribution of the pores in the SBA-15. Three samples were prepared, using 2, 4, and 6 deposition cycles to control the thickness of the films. It was determined that, as the average size of the pores decreased (with the increasing number of cycles), their size distribution remained narrow, indicating a homogeneous distribution of the ZrO 2 throughout the surfaces of the SBA-15 pores. This conclusion was confirmed by transmission electron microscopy and X-ray photoelectron spectroscopy. Further, the deposition rate appears to slow down after 4 cycles, an observation that we explain by different chemisorption properties of TDMAZ molecules over SBA-15 and ZrO 2 surface and by relating the kinetic diameter of TDMAZ with the pore diameter of SBA-15 after 4 deposition cycles: mass transport limitations have become significant at this point. They may be affected by complex factors like electronic effects caused by exposed ZrO 2 surface and the multi-directional adsorption of precursor molecules in the pore walls.
36 MATERIALS SCIENCE↗
CO Oxidation Catalyzed by Au Dispersed on SBA-15 Modified with TiO2 Films Grown via Atomic Layer Deposition (ALD)
It has been established that gold, when in nanoparticle (NP) form and in contact with reducible oxides, can promote oxidation reactions under mild conditions. Here, we report results from our exploration of the catalytic oxidation of carbon monoxide using catalysts where Au NPs were combined with thin titanium oxide films deposited on SBA-15 using atomic layer deposition (ALD). Both orders of deposition, with TiO2 added either before or after Au dispersion, were tested for two titania film thicknesses amounting to about half and full TiO2 monolayers. The resulting catalysts were characterized using various techniques, mainly electron microscopy and N2 adsorption–desorption isotherms, and the kinetics of the oxidation of CO with O2 were followed using infrared absorption spectroscopy. A synergy between the Au and TiO2 phases as it relates to the bonding and conversion of CO was identified, the tuning of which could be controlled by varying the synthetic parameters. The ALD of TiO2 films proved to be an effective way to maximize the Au-TiO2 interface sites, and with that help with the activation of molecular oxygen.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
From Langmuir isotherm to Brunauer–Emmett–Teller isotherm
Abstract Developed for multilayer adsorption, the Brunauer–Emmett–Teller ( BET) isotherm considers the adsorption of the first layer as an equilibrium chemical reaction between adsorbate molecules and adsorption sites and the adsorption of the second and subsequent layers as a condensation–evaporation process. Following the recent development of an activity‐based formulation for the Langmuir isotherm for monolayer adsorption, we present an activity‐based formulation for the BET isotherm in which species concentrations are replaced with species activities. Capturing the adsorbent surface heterogeneity for the adsorption of the first layer, the resulting thermodynamic BET isotherm is shown to accurately represent pure component adsorption isotherms over the relative pressure range of zero to unity or prior to the onset of capillary condensation. The thermodynamic BET isotherm should facilitate accurate estimation of monolayer adsorption capacity and the corresponding adsorbent surface area.
Vyawahare, Pradeep↗
Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)
Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.
36 MATERIALS SCIENCE↗
Generalized Brunauer–Emmett–Teller Isotherm for Mixed-Gas Multilayer Adsorption Equilibria
Here, this work presents a rigorous thermodynamic framework to calculate mixed-gas multilayer adsorption equilibria from single-component isotherms and vapor–liquid equilibria. Named the generalized Brunauer–Emmett–Teller (gBET) isotherm, the newly formulated isotherm considers adsorbent surface heterogeneity, competitive adsorption on the monolayer, condensation–evaporation on the subsequent layers, and adsorbed phase nonideality for the monolayer and the subsequent layers. The monolayer adsorbed phase nonideality is tracked by using an area-based adsorption nonrandom two-liquid activity coefficient model. The adsorbed phase composition and corresponding nonideality in the subsequent layers are calculated at the dew point condition of the mixed gas with either an equation of state or an activity coefficient model for the vapor–liquid equilibria. The proposed model is validated with six single and three binary multilayer adsorption equilibrium systems, and the model results are compared against those from the classical BET isotherm for single-component adsorption and ideal adsorbed solution theory for mixed-gas adsorption equilibria.
09 BIOMASS FUELS↗
Effect of layer bending on montmorillonite hydration and structure from molecular simulation
Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm 2 system that included dehydrated and hydrated regions in the same interlayer.
58 GEOSCIENCES↗
Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis
The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.
36 MATERIALS SCIENCE↗
High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)
This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.
54 ENVIRONMENTAL SCIENCES↗
On the pressure and temperature dependence of adsorption densities and other thermodynamic properties in gas shales
Black shale source rocks have become a key natural gas resource in the US and China. Unlike in conventional hydrocarbon reservoirs, a significant fraction of the gas in shales is adsorbed onto the inordinate specific surface areas (SSA) of micropores inside organic matter with densities exceeding those of bulk gas by an order of magnitude. Yet, the distribution of adsorbed versus bulk gas remains poorly understood. Experiments measure excess adsorption isotherms, which are corrected to absolute adsorption amounts by a factor that involves the density, ρ abs , of the adsorption layer(s). Constant values of ρ abs are widely used in the literature but, we argue, are inconsistent with the pressure- and temperature-dependence of adsorption behavior. In this work, we propose a more reliable approach that assumes a constant volume of each adsorption layer that can be determined from the SSA of the substrate. We measure nine high-pressure methane excess adsorption isotherms up to 15 MPa for three shale samples from China and Germany at temperatures of 65, 75, and 95 °C, and low-pressure carbon dioxide and nitrogen isotherms to estimate the SSA. We also investigate another 24 isotherms at 35, 50, and 60° C or 8 Chinese shale samples and 3 isotherms at 60, 90, and 120° C for an Argentinian shale sample from the literature. A new variation of the multi-layer Ono-Kondo lattice model is introduced to derive the temperature and pressure dependent adsorption layer densities, which are subsequently used to obtain the absolute adsorption isotherms of interest. Finally, this practical methodology has the potential to significantly improve our estimates of shale gas-in-place resources, and adsorption modeling more broadly.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Nonclassical Behavior in Competitive Ion Adsorption at a Charged Solid–Water Interface
Ion adsorption at solid-water interfaces is commonly described by interactions between specific surface sites and adsorbed ions in classical models. However, energetic contributions from non-site-specific ion-ion interactions have not been explored systematically. Here, we report non-classical behaviors observed during competitive adsorption between Sr 2+ and Na + /Rb + at the negatively-charged muscovite mica (001)-water interface, revealing apparent controls of adsorbed ion speciation over the interfacial reactivity. In the absence of competing cations, Sr 2+ adsorbs in approximately equivalent proportions of inner-sphere and outer-sphere complexes, whereas it adsorbs predominantly as an outer-sphere complex in the presence of Na + /Rb + . This transformation of adsorbed Sr 2+ speciation significantly decreases its adsorption strength, as indicated by ~15-fold shift in the Sr 2+ adsorption edge concentration, compared to that calculated from a classical Langmuir isotherm model developed based on site-specific interactions. Finally, these observations highlight the importance of non-site-specific interactions in controlling the energetics of chemical reactions at the charged interface.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Unraveling the core of fuel cell performance: engineering the ionomer/catalyst interface
The biggest obstacle to the widespread implementation of polymer electrolyte membrane fuel cells (PEMFCs) is their cost, primarily due to the use of platinum catalysts. The high intrinsic catalyst activity exhibited on a rotating disk electrode (RDE) is rarely realized in a membrane electrode assembly (MEA), which is a long-standing challenge for PEMFCs and a cause of low catalyst utilization. To translate the high RDE performance of a catalyst into a MEA, the design of an ideal ionomer/catalyst interface is proposed: a thin, conformal ionomer film covers the maximum surface of a Pt nanoparticle and thus simultaneously maximizes catalyst utilization, (i.e., high mass activity and electrochemically active surface area) and O 2 diffusion rate (i.e., high current density performance) without compromising proton conduction. Building such an interface is a long-standing challenge due to the lack of interaction between the ionomer and catalyst particles, resulting in large ionomer agglomerates and inhomogeneous ionomer coverage over the catalyst nanoparticle, with consequent poor fuel cell performance. In this work, this ionomer/catalyst interface has been engineered, utilizing the electrostatic attraction between positively charged catalyst and negatively charged ionomer particles in a catalyst ink and preserved in a solid catalyst layer. As a result, this interface leads to previously unachieved proton exchange membrane fuel cell performance in terms of both catalyst utilization (75% vs. 45%) and peak/rated power density (i.e., 1.430/0.930 W cm -2 , H 2 /air, cathode Pt loading: 0.1 mgPt cm -2 ) for pure Pt catalysts, even better than those of Pt alloy catalysts. This work demonstrates the formation of an interface in the liquid phase (using ultra-small-angle X-ray scattering in combination with cryo-TEM, isothermal–titration–calorimetry) and the preservation of the interface in the solid catalyst layer (using TEM) and estimates the effective coverage and thickness of the ionomer film (using limiting current density, RDE and fuel cell performance).
25 ENERGY STORAGE↗
Single-Atom Doping at the Molecule–Metal Interface: How Rh Affects Surface Explosion Kinetics
Controlling chemical reactions at interfaces is central to many fields, including atmospheric, environmental, and biological chemistry as well as catalysis and corrosion. We performed a fundamental study of how catalytically active Rh atoms placed at and above a molecule-metal interface influence the surface reaction kinetics of tartaric acid (TA) decomposition. Specifically, TA decomposition on Cu(100) exhibits autocatalytic decomposition kinetics involving a slow initiation step followed by fast decomposition of the molecular layer. These so-called “surface explosions” are extremely sensitive to initial conditions, making it fundamentally interesting to study how the placement of reactive atoms affects the reaction rate. Temperature-programmed reaction (TPR) experiments reveal that Rh atoms embedded in the Cu(100) surface beneath TA (TA/Rh/Cu) or atop the TA layer (Rh/TA/Cu) reduce the decomposition temperature and modify the reaction rate constants. Isothermal TPR analysis reveals that Rh enhances both the initiation rate constant, k i , and the explosion rate constant, k e , facilitating explosive decomposition at lower temperatures. This result provides evidence that both the initiation and explosion steps occur at the metal-molecule interface and are accelerated by the presence of Rh at this interface. This study illustrates how small amounts of reactive elements may be used to control nonlinear kinetic processes at interfaces.
desorption↗
Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions
With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.
58 GEOSCIENCES↗
Water Uptake by Silica Nanopores: Impacts of Surface Hydrophilicity and Pore Size
Water uptake at 298 K by two mesoporous silicas with different pore sizes was studied using volumetric vapor sorption. Here, through variation of sample pretreatment temperature and time, the number of surface hydroxyl groups was varied, leading to marked changes in the water sorption behavior. With increasing surface hydroxyl density, the adsorption of water to the pore surfaces increases, leading to formation of a thicker and denser water sorption layer. The BET model was used to measure surface hydroxyl density from low pressure parts of water adsorption isotherms. Adsorption of water molecules scales with the number density of surface hydroxyl groups, and monolayer coverage is reached at reduced pressure of ca. 0.3. The onset of adsorption pore condensation of water shifts to lower reduced pressure with increasing surface hydroxyl density, indicating growing thickness of adsorption layers. However, the water pore condensation step of the desorption branch shifts to smaller reduced pressure with increasing surface hydroxylation, reducing the adsorption-desorption hysteresis width. An adsorbed phase model is utilized to calculate the distribution of pore water molecules between adsorbed and pore condensation phases, and characterize the density and thickness of the water sorption phase. In the smaller pores, the adsorbed phase forms a sparse monolayer, while approximately a bilayer is formed in the wider pores. The previously reported silica nanopore underfilling by water is confirmed and rationalized by a significantly reduced water density in the adsorbed phase with respect to bulk, while the pore core is filled with water of approximately bulk density.
36 MATERIALS SCIENCE↗
Effect of resistivity on poloidal asymmetries in electron density in the edge and scrape off layer in 3D full-F gyro-fluid simulations with FELTOR
In this contribution, we analyze the poloidal asymmetries in electron density in the edge and scrape-off layer (SOL) in a COMPASS-sized, diverted tokamak with the 3D full-F, isothermal, electromagnetic, and gyro-fluid model FELTOR. The study is performed for different simulations that span over 2 orders of magnitude in resistivity. The poloidal asymmetries of density are evaluated in the edge and the SOL relative to the outer midplane (OMP), where the highest densities are usually found. In the closed magnetic surfaces, the relative poloidal asymmetry with respect to the OMP is not larger than 20% independently of the plasma resistivity. For the open field lines in the SOL, the relative density asymmetry can range from 55% for the highest resistivity to around 40% for the lowest. The lowest densities inside the separatrix are found between the inner midplane and the top of the magnetic configuration, away from the X-point. In the SOL it is usually close to the X-point in the high field side. The observations in the closed field lines are consistent with ballooning transport but in the SOL oppose experimental evidence for high-density plasmas. In conclusion, this indicates the necessity for more complex physics to reproduce the experimental observations in the SOL for higher density plasmas, such as neutral-plasma dynamics and realistic divertor conditions.
3D full-f gyrofluid↗
Molecular Structure of Adsorbed Water Phases in Silica Nanopores
The adsorption of water vapor in silica nanopores with different pore morphologies and surface hydrophilicities was studied to quantify the densities and thicknesses of the water sorption layers and deduce their molecular structures. Water adsorption to surface hydroxyls is described by a multilayer sorption model. At low pressure, the water adsorption isotherms are largely independent of pore size and the adsorbed amounts scale with the surface hydroxyl density. Adsorbed-phase densities corresponding to the adsorption of two water molecules per surface hydroxyl group are found in the first adsorbed water layer for a wide range of surface hydroxyl densities. The densities and layer thickness values found in narrow pores indicate that patchy adsorbed layers form if not enough water molecules exist for a full layer, which coexist with dry pore surface regions. This behavior indicates cooperative adsorption effects, i.e., a preference for the formation of hydrogen bonds between water molecules bound to surface hydroxyls. In narrow pores, pore condensation limits further growth of the sorption layer, and in larger pores and at the planar quartz surface, a second adsorbed water layer is formed, which can hold up to approximately 4 additional water molecules per surface hydroxyl group. The water molecules in these thicker adsorption layers arrange such that the sorption layer density is similar to the bulk water density. Pore confinement limitations on the sorption layers are observed in pores with radii of as large as 8 nm. Molecular dynamics modeling reveals two preferential orientations for water molecules adsorbing to surface hydroxyl groups and suggests an intralayer structuring in the adsorbed monolayer. Furthermore, adsorbed water molecules in the sorption layer are bonded to the surface hydroxyl group via the donation and acceptance of hydrogen bonds.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗