SEARCH · Engineering Papers

Results for “Dynamic Procedures”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Real-time reconstruction and control of pedestal-top electron density using RMP and gas puff at KSTAR

We report the experimental results of controlling the pedestal-top electron density by applying resonant magnetic perturbation (RMP) with in-vessel control coils and the main gas puff in the 2024-2025 KSTAR experimental campaign. The density is reconstructed using a parameterized $ψ_N$ grid and five channels of line-averaged density measured by the two-colored interferometer (TCI). The reconstruction procedure is accelerated by deploying a multi-layer perceptron to run in approximately 120 µ s and is sufficiently fast for real-time control. A proportional-integral controller was adopted, with the controller gains estimated from the system identification procedure. The experimental results demonstrate that the developed controller can follow a dynamic target while exclusively using both actuators. The absolute percentage errors between the electron density at $ψ_N$ = 0.89 and the target were approximately 1.5% median and a 2.5% average, respectively. The developed controller can even lower the density by using the pump-out mechanism under RMP, and it can follow a more dynamic range of density targets than a single actuator controller. The developed controller will enable experimental scenario exploration within a shot by dynamically setting the density target or maintaining a constant electron density within a discharge.

EFIT↗

DOE OSTI2024

Inverter Model Validation and Calibration Using Phasor Measurement Unit Data

As the penetration of inverter-based renewable energy resources increases in the power grid, especially at the distribution and microgrid levels, the need to accurately represent them in planning studies increases as well. However, due to the lack of well-established standard procedures, and vendor reluctance towards the detailed sharing of proprietary models, automated dynamic model validation and parameter calibration tools for inverter based resources (IBRs) remain scarce. This work presents a model validation and parameter calibration platform for representing IBRs with generic phasor-domain models. Phasor measurements of power system events are used for continuous validation using the data playback method, and model parameters are re-calibrated if a significant mismatch between measurements and model response is observed. Unique features of the proposed platform include- (a) an iterative Bayesian optimization approach towards parameter calibration to address a possible mismatch between the structures of generic models implemented in simulation softwares and actual commercial inverters, (b) error metrics designed to account for a possible mismatch between the time resolution of simulation and measurements, and (c) analysis of the measurement-simulation mismatch to provide guidance to engineering personnel regarding model shortcomings. The performance of the platform has been illustrated using both simulated data and field measurements to validate/calibrate inverter models in GridLAB-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Representation and Impact of Water Head on Power System Planning and Operation

Representing water head information in power system model files, can provide a more realistic model of the system and thereby inform operation and planning personnel in the decision-making process. This article describes a procedure for modifying the power system model files (steady-state and dynamic) to represent water head information. Additionally, the impact of representing the water head on power system reliability studies including contingency analysis, cascading failure analysis and dynamic frequency response analysis has been investigated, using the modified power system models. This paper considers the detailed Western Electricity Coordination Council model during summer and winter conditions as the test system for the impact analysis. Results show that under reduced water head: 1) the number of critical voltage and branch flow violations increases; 2) chances of cascading failure and island formation increases; and 3) frequency nadir decreases as compared to those of the base cases where the water head information is not represented.

13 - HYDRO ENERGY↗

DOE OSTI2024

Measurement of the muon anomalous precession frequency $\omega_a$ in the Fermilab $g-2$ experiment

The Fermilab Muon $g-2$ Experiment was designed to measure the muon's anomalous magnetic moment $a_{\mu} = (g-2)/2$ to 140 parts per billion. The value of $a_{\mu}$ is proportional to the difference frequency $\omega_a=\omega_s-\omega_c$ between the muon's cyclotron frequency and spin precession frequency in the uniform magnetic field of the $g-2$ storage ring. The frequency $\omega_a$ is extracted from the time distribution of the mu-decay positrons recorded by 24 electromagnetic calorimeters positioned around the inner circumference of the storage ring. We will discuss the various approaches to the frequency determination including the reconstruction, fitting of time distributions, and procedures for handling the effects of gain changes, positron pileup and beam dynamics. We also discuss the data consistency checks and the strategy for the averaging of $\omega_a$ across the different analyses.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Markov Chain Monte Carlo Predictions of Neutron-rich Lanthanide Properties as a Probe of r -process Dynamics

Lanthanide element signatures are key to understanding many astrophysical observables, from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis that enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the r-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multidimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of r-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the r-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. Here, we find our mass predictions given outflows that undergo an extended (n,γ)⇌(γ,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

FRET-based dynamic structural biology: Challenges, perspectives and an appeal for open-science practices

Single-molecule FRET (smFRET) has become a mainstream technique for studying biomolecular structural dynamics. The rapid and wide adoption of smFRET experiments by an ever-increasing number of groups has generated significant progress in sample preparation, measurement procedures, data analysis, algorithms and documentation. Several labs that employ smFRET approaches have joined forces to inform the smFRET community about streamlining how to perform experiments and analyze results for obtaining quantitative information on biomolecular structure and dynamics. The recent efforts include blind tests to assess the accuracy and the precision of smFRET experiments among different labs using various procedures. These multi-lab studies have led to the development of smFRET procedures and documentation, which are important when submitting entries into the archiving system for integrative structure models, PDB-Dev. This position paper describes the current ‘state of the art’ from different perspectives, points to unresolved methodological issues for quantitative structural studies, provides a set of ‘soft recommendations’ about which an emerging consensus exists, and lists openly available resources for newcomers and seasoned practitioners. To make further progress, we strongly encourage ‘open science’ practices.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Multi-channel direct detection of light dark matter: theoretical framework

We present a unified theoretical framework for computing spin-independent direct detection rates via various channels relevant for sub-GeV dark matter — nuclear re- coils, electron transitions and single phonon excitations. Despite the very different physics involved, in each case the rate factorizes into the particle-level matrix element squared, and an integral over a target material- and channel-specific dynamic structure factor. We show how the dynamic structure factor can be derived in all three cases following the same procedure, and extend previous results in the literature in several aspects. For electron transitions, we incorporate directional dependence and point out anisotropic target materials with strong daily modulation in the scattering rate. For single phonon excitations, we present a new derivation of the rate formula from first principles for generic spin-independent couplings, and include the first calculation of phonon excitation through electron couplings. We also discuss the interplay between single phonon excitations and nuclear recoils, and clarify the role of Umklapp processes, which can dominate the single phonon production rate for dark matter heavier than an MeV. Our results highlight the complementarity between various search channels in probing different kinematic regimes of dark matter scattering, and provide a common reference to connect dark matter theories with ongoing and future direct detection experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Neutralizing topological obstructions to bubbles of nothing

Theories with compact extra dimensions can exhibit a vacuum instability known as a bubble of nothing. These decay modes can be obstructed if the internal manifold is stabilized by fluxes, or if it carries Wilson lines for background gauge fields, or if the instanton is incompatible with the spin structure. In each of these cases the decay can proceed by adding dynamical charged membranes or gauge fields. We give a general, bottom-up procedure for constructing approximate bubble of nothing solutions in models with internal spheres stabilized by flux and study the influence of the brane tension on the tunneling exponent, finding two branches of solutions that merge at a minimal superextremal value of the tension. In the case of Wilson operators and incompatible fermions, the relevant bubble is shown to be the Euclidean Reissner-Nordstrom black hole, and the ordinary decay exponent is modified by 1/g 2 effects. We examine the Dirac operator on this background and comment on the relevance for models of supergravity with gauged R-symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Approximation of refrigerant thermophysical properties using neural networks to speed up transient thermofluid simulations

Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.

Ma, Jiacheng↗

DOE OSTI2020

Broadening of the drumhead-mode spectrum due to in-plane thermal fluctuations of two-dimensional trapped ion crystals in a Penning trap

Two-dimensional crystals of ions stored in Penning traps are a leading platform for quantum simulation and sensing experiments. For small amplitudes, the out-of-plane motion of such crystals can be described by a discrete set of normal modes called the drumhead modes, which can be used to implement a range of quantum information protocols. However, experimental observations of crystals with Doppler-cooled and even nearground- state-cooled drumhead modes reveal an unresolved drumhead-mode spectrum. In this work, we establish in-plane thermal fluctuations in ion positions as a major contributor to the broadening of the drumhead-mode spectrum. In the process, we demonstrate how the confining magnetic field leads to unconventional in-plane normal modes, whose average potential and kinetic energies are not equal. This property, in turn, has implications for the sampling procedure required to choose the in-plane initial conditions for molecular-dynamics simulations. Further, for current operating conditions of the NIST Penning trap, our study suggests that the two-dimensional crystals produced in this trap undergo in-plane potential-energy fluctuations of the order of 10mK. Our study therefore motivates the need for designing improved techniques to cool the in-plane degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Rancor-HUNTER: Using a Simulator Engine for Realistic Human Performance Modeling of Nuclear Power Operations

The Human Unimodel for Nuclear Technology to Enhance Reliability (HUNTER) is a software system to simulate human performance in support of human reliability analysis (HRA) in nuclear power plants. This paper summarizes recent work to integrate HUNTER with a plant simulator, namely the Rancor Microworld Simulator. Rancor is an offshoot of earlier work at Idaho National Laboratory (INL) to support plant modernization. The graphical software tools used to mimic digital human-system interface upgrades at INL’s Human Systems Simulation Laboratory were linked to the Rancor Microworld Simulator, an INL-developed simplified plant model. HUNTER becomes a “virtual operator” coupled to the Rancor simulator, thereby allowing a tight coupling between a digital human twin and a digital twin of the plant. Rancor-HUNTER may be run through Monte Carlo iterations across a dynamic range of performance shaping factors, thereby producing distributions of human performance in terms of procedure paths, errors instantiations, and task durations. This paper overviews the various unique features of Rancor-HUNTER and presents an example run of Rancor-HUNTER for a startup scenario.

99 - GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Unraveling Adsorbate-Induced Structural Evolution of Iron Carbide Nanoparticles

Iron carbide (Fe x C y ) nanoparticles (NPs) are promising candidates for replacing platinum group metals in industrial applications, such as high-temperature Fischer–Tropsch synthesis. However, due to their amorphous nature, characterization of the active sites has been challenging experimentally and computationally. Here, using a combined density functional theory (DFT), neural network interatomic potential-assisted global optimization, and ensemble learning study, we evaluate dynamic surface changes associated with syngas (H and CO) interactions. For this purpose, we have developed a general procedure that we use to model an experimentally relevant 270-atom Fe 182 C 88 NP using the neural network-assisted stochastic surface walk global optimization algorithm (SSW-NN). Once generated, the Fe 182 C 88 NP active sites and particle morphology are thoroughly characterized before the effects of syngas adsorbate interactions are explored by using DFT and molecular dynamics simulations. Lastly, we explore correlations between geometric and electronic features of the active sites and the adsorption of H (H ads ), using a regularized random forest machine learning algorithm. In doing so, we identified the Fe–C coordination number and p orbital occupancy as the most important descriptors affecting H ads . Furthermore, using a combined ML and quantum chemistry approach, our work demonstrates a general and efficient procedure for generating and probing complex surface phenomena on binary nanoparticles.

Adsorption↗

DOE OSTI2021

A hierarchical gray-box dynamic modeling methodology for direct-expansion cooling systems to support control stability analysis; Méthodologie de modélisation dynamique hiérarchique de type boîte grise pour les systèmes de refroidissement à détente directe afin de soutenir l'analyse de stabilité de la commande

In this paper, a gray-box dynamic modeling approach for direct-expansion cooling systems is presented. The overall approach incorporates a multi-stage training procedure that consists of 1) identification of component sub-models from quasi-steady-state performance data, 2) system model integration with estimation of refrigerant charge and 3) fine tuning of thermal capacitances of the evaporator and condenser to capture the system dynamic responses. Compared to traditional physics-based models, the proposed modeling approach has advantages including reduced engineering efforts in the model development phase, improved computational efficiency and enhanced prediction accuracy. The modeling method was validated using a 3-ton variable-speed heat pump and proved to be capable of accurately predicting the system transient behaviors over a wide range of operating conditions. The established dynamic model was then applied for control stability analysis, with a specific goal of determining a proper control execution time step. The case study results showed that the stable control execution time step could change significantly, from 3 sec to 19 sec, as the operating conditions and control settings vary, and a proper selection of the execution time step is critical to ensure stable and reliable operations.

42 ENGINEERING↗

DOE OSTI2021

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab-initio multiple cloning (AIMC) approach based on the multicongurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to non-adiabatic dynamics with a particular strength in the proper description of decoherence effects due to branching of vibrational wavepackets at a level crossing. In this work, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-eld theory to determine the equations of motion for classical trajectories. The multi-dimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based non-adiabatic coupling approach within the AIMC method, all observable are calculated on-the-y in the non-adiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level non-adiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Evolution of Charge-Lattice Dynamics across the Kuramoto Synchronization Phase Diagram of Quantum Tunneling Polarons in Cuprate Superconductors

Because of its sensitivity to the instantaneous structure factor, S(Q,t = 0), Extended X-ray Absorption Fine Structure (EXAFS) is a powerful tool for probing the dynamic structure of condensed matter systems in which the charge and lattice dynamics are coupled. When applied to hole-doped cuprate superconductors, EXAFS has revealed the presence of internal quantum tunneling polarons (IQTPs). An IQTP arises in EXAFS as a two-site distribution for certain Cu–O pairs, which is also duplicated in inelastic scattering but not observed in standard diffraction measurements. The Cu–Sr pair distribution has been found to be highly anharmonic and strongly correlated to both the IQTPs and to superconductivity, as, for example, in YSr 2 Cu 2.75 Mo 0.25 O 7.54 (T c =84 K). In order to describe such nontrivial, anharmonic charge-lattice dynamics, we have proposed a model Hamiltonian for a prototype six-atom cluster, in which two Cu-apical-O IQTPs are charge-transfer bridged through Cu atoms by an O atom in the CuO 2 plane and are anharmonically coupled via a Sr atom. By applying an exact diagonalization procedure to this cluster, we have verified that our model indeed produces an intricate interplay between charge and lattice dynamics. Then, by using the Kuramoto model for the synchronization of coupled quantum oscillators, we have found a first-order phase transition for the IQTPs into a synchronized, phase-locked phase. Most importantly, we have shown that this transition results specifically from the anharmonicity. Finally, we have provided a phase diagram showing the onset of the phase-locking of IQTPs as a function of the charge-lattice and anharmonic couplings in our model. We have found that the charge, initially confined to the apical oxygens, is partially pumped into the CuO 2 plane in the synchronized phase, which suggests a possible connection between the synchronized dynamic structure and high-temperature superconductivity (HTSC) in doped cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Inference-Engine v0.1.0

Given a pre-trained neural network, Inference-Engine performs maps network inputs to outputs by executing the forward pass through the provided network. Although the predominant programming language for machine-learning is Python, most high-performance computing (HPC) applications are written in Fortran, C, or C++. Inference-Engine aims to support HPC programs and is written in Fortran, a language with a large feature set supporting interoperability with C. This software exposes concurrency in a portable way by using standard language features that some modern Fortran compilers can exploit with various optimizations, including offloading computation to a Graphics Processing Unit (GPU). In particular, this software makes extensive use of Fortran's "do concurrent" parallel loop construct, implicitly parallel array statements, and pure procedures that can be invoked inside "do concurrent" blocks. Inference-Engine also supports dynamic choice of inference methods at runtime. Two current options include one method that uses Fortran's "dot_product" intrinsic function inside "do concurrent" blocks and another method that instead uses Fortran' "matmul" array intrinsic function. We plan to investigate automatic compiler offloading of "do concurrent" calculations to GPUs and compile-time substitution of optimized libraries such as the Basic Linear Algebra Library (BLAS) for "matmul" invocations. We also envision the potential for the choice of which method to use could happen at program launch based on in situ performance measurements on any given platform.

Rouson, Damian↗

DOE OSTI2022

Multiscale water drop contact angles at selected silica surfaces

In this paper, multiscale advancing contact angles for glycerol/water drops at silica surfaces are reported for millidrops, submicron-drops, and nanodrops. Selected silica surfaces were muscovite, silicon, and talc. The contact angles for millidrops (1–2 mm) were determined by the traditional sessile drop technique. For submicron-drops (0.1–1.0 μm), a hollow tip Atomic Force Microscope (AFM) procedure was used. The contact angles for nanodrops (~7 nm) were examined from Molecular Dynamics (MD) simulation. The results were compared to evaluate the effect of drop size on the contact angle. In the case of the hydrophobic talc surface, the 75° advancing contact angle did not vary significantly with drop size. For the hydrophilic muscovite surface, the water drop wet the surface and an advancing contact angle of about 10° was found for the millidrops and submicron-drops. However, for the MD simulated nanodrops, attachment and spreading of the ~7 nm drop created a 2D film of molecular dimensions, the contact angle of which was difficult to define and varied from 0° to 17°. Perhaps of equal interest from the MD simulation results was that the spreading of the glycerol/water nanodrop at the muscovite surface resulted in crystallographic directional transport of water molecules to the extremities of the 2D film. Such separation and segregation left the center of the film with an increased concentration of glycerol. Based on these results, the line tension, which has been found in other investigations to account for contact angle decrease with a decrease in drop size, does not seem to be a significant factor in this study.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗