Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compositional simulation”

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,639 records · Page 91

Chemical Composition and PCT Data for the Initial Set of Hanford Enhanced Waste Loading Glasses

In this report, the Savannah River National Laboratory provides chemical analyses and Product Consistency Test results for 20 simulated high level waste glasses fabricated by the Pacific Northwest National Laboratory. The results of these analyses will be used as part of efforts to revise or extend the validation ranges of the current Hanford Waste Treatment and Immobilization Plant glass property models to cover a broader span of waste compositions. The measured chemical composition data are reported and compared with the targeted values for each component for each glass. Two components of the study glasses, fluorine and silver, were not measured since each of these species would have required the use of an additional preparation method and their measured values were likely to be near or below analytical detection limits. Some of the glasses were difficult to prepare for chemical analysis. A sodium peroxide fusion dissolution method was successful in completely dissolving the glasses. Components present in the glasses in minor concentrations can be difficult to measure using this dissolution method due to dilution requirements. The use of a lithium metaborate preparation method for the minor components (planned for use since it is typically successful in digesting Defense Waste Processing Facility HLW glasses) resulted in an unacceptable amount of undissolved solids remaining in the sample solutions. An acid dissolution method was used instead, which provided more thorough dissolution of the glasses, although a small amount of undissolved material remained for some of the study glasses. The undissolved material was analyzed to determine those components of the glasses that did not fully dissolve. These components (e.g., calcium and chromium) were present in sufficient quantities to be reported from the measurements resulting from the sodium peroxide fusion preparation method, which did not leave undissolved material. Overall, the analyses resulted in sums of oxides that ranged from about 98 to 101.5 wt % for the study glasses, indicating excellent recovery of all the components in the chemical composition analyses. Comparisons of the targeted and measured chemical compositions indicated that, in general, the measured values for the glasses met the targeted concentrations. Exceptions were Cr 2 O 3 , MgO, and P 2 O 5 . The measured values for Cr 2 O 3 were somewhat low when compared to the targeted values for all of the study glasses targeting Cr 2 O 3 concentrations above 0.5 wt %. Many of the measured MgO and P 2 O 5 values were below the targeted values for those glasses that contained these components. Two of the study glasses exhibited differences from the targeted compositions that may indicate a batching error. Glasses EWG-HAI-Centroid-2 and EWG-OL-1672 had measured values for Al 2 O 3 and SiO 2 that were lower than the targeted values, and measured values for B 2 O 3 that were higher than the targeted values. Glass EWG-HAI-Centroid-2 also had a measured value for Fe 2 O 3 that was lower than the targeted value. A review of the PCT data, including standards and blanks, revealed no issues with the performance of the tests. The PCT results were normalized to both the targeted and measured compositions of the study glasses. Comparisons of the normalized PCT results for both the quenched and Canister Centerline Cooled versions of the study glasses are made with the Environmental Assessment benchmark glass for reference.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Simulation of integrated pollutant removal (IPR) water-treatment system using ASPEN Plus

Capturing CO2 from fossil fuel combustion provides an opportunity for tapping a significant water source which can be used as service water for a capture-ready power plant and its peripherals. Researchers at the National Energy Technology Laboratory (NETL) have patented a process—Integrated Pollutant Removal (IPR®)—that uses off-the-shelf technology to produce a sequestration ready CO2 stream from an oxy-combustion power plant. Water condensed from oxy-combustion flue gas via the IPR system has been analyzed for composition and an approach for its treatment—for in-process reuse and for release—has been outlined. A computer simulation model in ASPEN Plus has been developed to simulate water treatment of flue gas derived wastewater from IPR systems. At the field installation, water condensed in the IPR process contains fly ash particles, sodium (largely from spray-tower buffering) and sulfur species as well as heavy metals, cations, and anions. An IPR wastewater treatment system was modeled using unit operations such as equalization, coagulation and flocculation, reverse osmosis, lime softening, crystallization, and pH correction. According to the model results, 70% (by mass) of the inlet stream can be treated as pure water, the other 20% yields as saleable products such as gypsum (CaSO4) and salt (NaCl) and the remaining portion is the waste. More than 99% of fly ash particles are removed in the coagulation and flocculation unit and these solids can be used as filler materials in various applications with further treatment. Results discussed relate to a slipstream IPR installation and are verified experimentally in the coagulation/flocculation step.

20 FOSSIL-FUELED POWER PLANTS↗

WE-EF-BRA-02: A Monte Carlo Study of Macroscopic and Microscopic Dose Descriptors for Kilovoltage Cellular Dosimetry

Purpose: To investigate how doses to cellular (microscopic) targets depend on cell morphology, and how cellular doses relate to doses to bulk tissues and water for 20 to 370 keV photon sources using Monte Carlo (MC) simulations. Methods: Simulation geometries involve cell clusters, single cells, and single nuclear cavities embedded in various healthy and cancerous bulk tissue phantoms. A variety of nucleus and cytoplasm elemental compositions are investigated. Cell and nucleus radii range from 5 to 10 microns and 2 to 9 microns, respectively. Doses to water and bulk tissue cavities are compared to nucleus and cytoplasm doses. Results: Variations in cell dose with simulation geometry are most pronounced for lower energy sources. Nuclear doses are sensitive to the surrounding geometry: the nuclear dose in a multicell model differs from the dose to a cavity of nuclear medium in an otherwise homogeneous bulk tissue phantom by more than 7% at 20 keV. Nuclear doses vary with cell size by up to 20% at 20 keV, with 10% differences persisting up to 90 keV. Bulk tissue and water cavity doses differ from cellular doses by up to 16%. MC results are compared to cavity theory predictions; large and small cavity theories qualitatively predict nuclear doses for energies below and above 50 keV, respectively. Burlin’s (1969) intermediate cavity theory best predicts MC results with an average discrepancy of 4%. Conclusion: Cellular doses vary as a function of source energy, subcellular compartment size, elemental composition, and tissue morphology. Neither water nor bulk tissue is an appropriate surrogate for subcellular targets in radiation dosimetry. The influence of microscopic inhomogeneities in the surrounding environment on the nuclear dose and the importance of the nucleus as a target for radiation-induced cell death emphasizes the potential importance of cellular dosimetry for understanding radiation effects. Funded by the Natural Sciences and Engineering Research Council of Canada (NSERC), the Canada Research Chairs Program (CRC), and the Ontario Ministry of Training, Colleges and Universities.

60 APPLIED LIFE SCIENCES↗

LATE-TIME EVOLUTION OF COMPOSITE SUPERNOVA REMNANTS: DEEP CHANDRA OBSERVATIONS AND HYDRODYNAMICAL MODELING OF A CRUSHED PULSAR WIND NEBULA IN SNR G327.1-1.1

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology: a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for a mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock (RS), whichcan occur as a result of a density gradient in the ambient medium and/or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a ∼17,000-year-old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar’s motion. We also show that the RS/PWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to γ-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs.

79 ASTRONOMY AND ASTROPHYSICS↗

THE EFFECT OF METALLICITY-DEPENDENT T-τ RELATIONS ON CALIBRATED STELLAR MODELS

Mixing length theory is the predominant treatment of convection in stellar models today. Usually described by a single free parameter, α, the common practice is to calibrate it using the properties of the Sun, and apply it to all other stellar models as well. Asteroseismic data from Kepler and CoRoT provide precise properties of other stars which can be used to determine α as well, and a recent study of stars in the Kepler field of view found α to vary with metallicity. Interpreting α obtained from calibrated stellar models, however, is complicated by the fact that the value for α depends on the surface boundary condition of the stellar model, or T-τ relation. Calibrated models that use typical T-τ relations, which are static and insensitive to chemical composition, do not include the complete effect of metallicity on α. We use three-dimensional radiation-hydrodynamic simulations to extract metallicity-dependent T-τ relations and use them in calibrated stellar models. We find the previously reported α-metallicity trend to be robust, and not significantly affected by the surface boundary condition of the stellar models.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying the effect of tow architecture variability on the performance of biaxially braided composite tubes

Targeted experiments are performed with stereoscopic digital image correlation to quantify the effect of irregular braid structure on the local surface mechanical response of carbon-fiber reinforced epoxy composite tubes. Virtual models comparing both nonuniform and uniform tow structure of the outer braid of each specimen are integrated into numerical simulations similar to experiments. Experimental and numerical results agree well. Statistical tests are used to support correlations between as-manufactured braid variation ranges to local increases in stress up to 18% compared to an ideal, uniform structure in the elastic regime. Results support improved manufacturing and modeling efforts of defect sensitive braided composites.

36 MATERIALS SCIENCE↗

Atmospheric Escape From TOI-700 d: Venus versus Earth Analogs

The recent discovery of an Earth-sized planet (TOI-700 d) in the habitable zone of an early-type M-dwarf by the Transiting Exoplanet Survey Satellite constitutes an important advance. In this Letter, we assess the feasibility of this planet to retain an atmosphere – one of the chief ingredients for surface habitability – over long timescales by employing state-of-the-art magnetohydrodynamic models to simulate the stellar wind and the associated rates of atmospheric escape. We take two major factors into consideration, namely, the planetary atmospheric composition and magnetic field. In all cases, we determine that the atmospheric ion escape rates are potentially a few orders of magnitude higher than the inner Solar system planets, but TOI-700 d is nevertheless capable of retaining a 1 bar atmosphere over gigayear timescales for certain regions of the parameter space. The simulations show that the unmagnetized TOI-700 d with a 1 bar Earth-like atmosphere could be stripped away rather quickly (< 1 gigayear), while the unmagnetized TOI-700 d with a 1 bar CO 2 -dominated atmosphere could persist for many billions of years; we find that the magnetized Earth-like case falls in between these two scenarios. We also discuss the prospects for detecting radio emission of the planet (thereby constraining its magnetic field) and discerning the presence of an atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffusion of Selective Penetrants in Interfacially Modified Block Copolymers from Molecular Dynamics Simulations

In order to show the influence of the interface on structure and dynamics of microphase separated polymer systems, we study interfacially modified AB block copolymers with small molecule penetrants. The polymers have a random midblock or tapered midblock whose composition varies from pure A to pure B (or from pure B to pure A for an inverse taper) between two pure blocks of A and B. We perform simple coarse-grained molecular dynamics simulations of symmetric polymers that form lamellae. With normal tapering, both polymer and penetrant diffusion parallel to the lamellae increases as taper length increases. Inverse tapered polymers exist in different conformational states (e.g. stretched versus folded back and forth across the interface) with different dynamic behavior, leading to non-monotonic trends in their diffusion. However, the local mixing of monomers (rather than polymer conformation) appears to be the most important factor in determining penetrant diffusion.

36 MATERIALS SCIENCE↗

Resilient high-temperature reverse osmosis desalination membranes

Conventional thin-film composite (TFC) reverse osmosis (RO) membranes experience irreversible performance loss at high temperatures, restricting their use in industries with high-temperature streams, including oil and gas, pharmaceuticals, electronics, power generation, food production, and hybrid desalination plants. However, the mechanisms driving the performance decline of TFC membranes at high temperatures remain poorly understood. Herein, we combine controlled experiments, molecular dynamics simulations, and micromechanical modeling to elucidate TFC failure mechanisms and to evaluate thermally resilient thin-film cross-linked (TFX) composite membrane. Upon exposure to elevated temperatures (>60°C), salt rejection of TFC dropped from ~99 to <90%, with irreversible structural damage in the polysulfone layer, confirmed by scanning electron microscopy. In contrast, the TFX membrane maintained ~99% salt rejection and showed no signs of physical degradation up to 80°C. Our combined analyses revealed that TFC membrane failure arises from irreversible pore expansion in the thermoplastic polysulfone support, leading to polyamide film rupture and delamination. TFX membranes resist thermal deformation, enabling ultrahigh-temperature RO desalination and water reuse.

Wu, Jishan [Rice Univ., Houston, TX (United States↗

Low-power, ultrafast, and dynamic all-optical tunable plasmonic analog to electromagnetically induced transparency in two resonators side-coupled with a waveguide system

We theoretically and numerically investigate a low-power, ultrafast, and dynamic all-optical tunable plasmonic analog to electromagnetically induced transparency (EIT) in two nanodisk resonators side-coupled to a metal-insulator-metal plasmonic waveguide system. The optical Kerr effect is enhanced by the slow light effect of the plasmonic EIT-like effect and the plasmonic waveguide based on graphene-Ag composite material structures with giant effective Kerr nonlinear coefficient. The optical Kerr effect modulation method is applied to improve tuning rate with response time of subpicoseconds or even femtoseconds. With dynamically tuning the propagation phase of the plasmonic waveguide, π-phase shift of the transmission spectrum in the plasmonic EIT-like system is achieved under excitation of a pump light with an intensity as low as 5.85 MW/cm{sup 2}. The group delay is controlled between 0.09 and 0.4 ps. All observed schemes are analyzed rigorously through finite-difference time-domain simulations and coupled-mode formalism. Results show a new direction toward the low power consumption and ultrafast responses of integration plasmonic photonic devices and all-optical dynamical storage of light devices in optical communication and quantum information processing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Smart Adhesive Joint with High-Definition Fiber-Optic Sensing for Automotive Applications

Structural health monitoring of fiber-reinforced composite-based joints for automotive applications during their manufacturing and on-demand assessment for its durability in working environments is critically needed. High-definition fiber-optic sensing is an effective method to measure internal strain/stress development using minimally invasive continuous sensors. The sensing fiber diameters are in the same order of magnitude when compared to reinforcement (glass, basalt, or carbon fibers) used in polymer composites. They also offer a unique ability to monitor the evolution of residual stresses after repeated thermal exposure with varying temperatures for automotive components/joints during painting using an electrophoretic painting process. In this paper, a high-definition fiber-optic sensor utilizing Rayleigh scattering is embedded within an adhesive joint between a carbon fiber-reinforced thermoset composite panel and an aluminum panel to measure spatially resolved strain development, residual strain, and thermal expansion properties during the electrophoretic paint process-simulated conditions. The strain measured by the continuous fiber-optic sensor was compared with an alternate technique using thermal digital image correlation. The fiber-optic sensor was able to identify the spatial variation of residual strains for a discontinuous carbon fiber-reinforced composite with varying local fiber orientations and resin content.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Geochemical and Geomechanical Effects on Wellbore Cement Fractures

Experimental studies were conducted using batch reactors, X-ray microtomograpy (XMT), and computational fluid dynamics (CFD) simulation to determine changes in cement fracture surfaces, fluid flow pathways, and permeability with geochemical and geomechanical processes. Composite Portland cement-basalt caprock core with artificial fractures was prepared and reacted with CO2-saturated groundwater at 50°C and 10 MPa for 3 to 3.5 months under static conditions to understand the geochemical and geomechanical effects on the integrity of wellbores containing defects. Cement-basalt interface samples were subjected to mechanical stress at 2.7 MPa before the CO2 reaction. XMT provided three-dimensional (3-D) visualization of the opening and interconnection of cement fractures due to mechanical stress. After the CO2 reaction, XMT images revealed that calcium carbonate precipitation occurred extensively within the fractures in the cement matrix, but only partially along fractures located at the cement-basalt interface. The permeability calculated based on CFD simulation was in agreement with the experimentally measured permeability. The experimental results imply that the wellbore cement with fractures is likely to be healed during exposure to CO2-saturated groundwater under static conditions, whereas fractures along the cement-caprock interface are still likely to remain vulnerable to the leakage of CO2. CFD simulation for the flow of different fluids (CO2-saturated brine and supercritical CO2) using a pressure difference of 20 kPa and 200 kPa along ~2 cm-long cement fractures showed that a pressure gradient increase resulted in an increase of CO2 fluids flux by a factor of only ~3-9 because the friction of CO2 fluids on cement fracture surfaces increased with higher flow rate as well. At the same pressure gradient, the simulated flow rate was higher for supercritical CO2 than CO2-saturated brine by a factor of only ~2-3, because the viscosity of supercritical CO2 is much lower than that of CO2-saturated brine. The study suggests that in deep geological reservoirs the geochemical and geomechanical processes have coupled effects on the wellbore cement fracture evolution and fluid flow along the fracture surfaces.

36 MATERIALS SCIENCE↗

How well are hazards associated with derechos reproduced in regional climate simulations?

Abstract. A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km2, caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Shepherd, Tristan (ORCID:0000000186276419)↗

Statistical effects related to low numbers of reacting molecules analyzed for a reversible association reaction A + B = C in ideally dispersed systems: An apparent violation of the law of mass action

Theoretical analysis and computer simulations (Monte Carlo and numerical integration of differential equations) show that the statistical effect of a small number of reacting molecules depends on a way the molecules are distributed among the small volume nano-reactors (droplets in this study). A simple reversible association A + B = C was chosen as a model reaction, enabling to observe both thermodynamic (apparent equilibrium constant) and kinetic effects of a small number of reactant molecules. When substrates are distributed uniformly among droplets, all containing the same equal number of substrate molecules, the apparent equilibrium constant of the association is higher than the chemical one (observed in a macroscopic—large volume system). The average rate of the association, being initially independent of the numbers of molecules, becomes (at higher conversions) higher than that in a macroscopic system: the lower the number of substrate molecules in a droplet, the higher is the rate. This results in the correspondingly higher apparent equilibrium constant. A quite opposite behavior is observed when reactant molecules are distributed randomly among droplets: the apparent association rate and equilibrium constants are lower than those observed in large volume systems, being the lower, the lower is the average number of reacting molecules in a droplet. The random distribution of reactant molecules corresponds to ideal (equal sizes of droplets) dispersing of a reaction mixture. Our simulations have shown that when the equilibrated large volume system is dispersed, the resulting droplet system is already at equilibrium and no changes of proportions of droplets differing in reactant compositions can be observed upon prolongation of the reaction time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Devising a framework to Couple Electromagnetic Geophysical Measurements with Reservoir Simulations to Monitor CO2 Plume Movement

Monitoring CO2 plume movement and extent is critical to any CO2 utilization and storage (CCUS) project but can present significant challenges. Downhole measurements (such as cross-well seismic) can give the resolution needed to track the CO2 plume. That resolution comes at a significant expense in equipment and lost production time, and the resolution is only in the vicinity of that well. Surface seismic is less expensive but has decreasing resolution with depth and may not be sensitive to changes in the CO2 plume. Conversely, electrical conductivity has been shown to be correlated with changes in phase saturation associated with CCUS operations. A novel approach for subsurface monitoring that can be applied with minimal expense and operation disruption is charged wellbore casing controlled-source electromagnetics (CWC-CSEM). This technique injects electrical current into two legacy wells. Surface measurements of the induced electric and magnetic fields are obtained with distributed data loggers. Time-lapse measurements are sensitive to changes in CO2 plume movement, but inversion requires the application of prior knowledge from prior characterization efforts. We propose a methodology where simulations are used to inform CWC-CSEM inversion while simultaneously being validated. As such, uncertainty is reduced, and more robust flow models can be devised. A three-year study is underway at the Bell Creek Field in southeast Montana (DOE-NETL Project DE-FE0028320) to evaluate the effectiveness of CWC-CSEM on an active CO2-EOR operation. Determining the relationship between electrical resistivity and phase saturations, and coupling that to a reservoir model is the thrust of this project. Initially, the relative permeability/capillary pressure relationship and the compositional oil model were identified as model parameters with high uncertainty. To elucidate the impact of these uncertain parameters, simulations were run for ten years with a single five-spot pattern. Final water saturation distributions were converted into resistivity using Archie’s Law and available fluid sampling data to provide an exploration of the time lapse electromagnetic anomaly which can be anticipated using CWC-CSEM methods. This technique has promise as an inexpensive and accurate way to monitor CO2 injection in a CCUS site.

Moodie, Nathan↗

Black carbon absorption at the global scale is affected by particle-scale diversity in composition

Atmospheric black carbon (BC) exerts a strong, but uncertain, warming effect on the climate. BC that is coated with non-absorbing material absorbs more strongly than the same amount of BC in an uncoated particle, but the magnitude of this absorption enhancement (E abs ) is not well constrained. Modelling studies and laboratory measurements have found stronger absorption enhancement than has been observed in the atmosphere. Here, using a particle-resolved aerosol model to simulate diverse BC populations, we show that absorption is overestimated by as much as a factor of two if diversity is neglected and population-averaged composition is assumed across all BC-containing particles. If, instead, composition diversity is resolved, we find E abs =1-1.5 at low relative humidity, consistent with ambient observations. This study offers not only an explanation for the discrepancy between modelled and observed absorption enhancement, but also demonstrates how particle-scale simulations can be used to develop relationships for global-scale models.

Science & Technology - Other Topics↗

Perspective: Size selected clusters for catalysis and electrochemistry

We report that size-selected clusters containing a handful of atoms may possess noble catalytic properties different from nano-sized or bulk catalysts. Size- and composition-selected clusters can also serve as models of the catalytic active site, where an addition or removal of a single atom can have a dramatic effect on their activity and selectivity. In this Perspective, we provide an overview of studies performed under both ultra-high vacuum and realistic reaction conditions aimed at the interrogation, characterization and understanding of the performance of supported size-selected clusters in heterogeneous and electrochemical reactions, which address the effects of cluster size, cluster composition, cluster-support interactions and reaction conditions, the key parameters for the understanding and control of catalyst functionality. Computational modelling based on density functional theory sampling of local minima and energy barriers or ab initio Molecular Dynamics simulations is an integral part of this research by providing fundamental understanding of the catalytic processes at the atomic level, as well as by predicting new materials compositions which can be validated in experiments. Lastly, we discuss approaches which aim at the scale up of the production of well-defined clusters for use in real world applications.

36 MATERIALS SCIENCE↗

A classical molecular dynamics simulation method for the formation of “dry” gels from boro-aluminosilicate glass structures

In contact with water, glass transforms into amorphous and porous structures called gels. A simulation method based on classical molecular dynamics is proposed here to mimic “dry” gels forming from initial oxide glass structures. Six glass compositions were investigated. Two behaviours were evidenced depending on the initial glass composition, and in particular on the quantity of elements removed. If a large quantity of soluble elements (B, Na) was removed, it induced an increase in the average pore size within the gels, and the time needed to stabilise the gel structure increased because more local atomic rearrangements occurred. The gel network displayed a higher proportion of Si-Q4 at the expense of Si-Q3 and a lower average ring size compared to the glass network, irrespective of the glass composition. Surface effects were also highlighted in the dry gels, such as the presence of 3-coordinated Al and a decrease in the average angle Si-O-Si and Al-O-Al. Our findings will be compared to both wet gels and experimental data in further studies, to help find the best procedure to simulate such structures.

36 MATERIALS SCIENCE↗