Computational Fluid Dynamics and Magnetohydrodynamics Modeling of Liquid Metal Flows for Fusion Blanket Applications
Explore the source record for details and available documents.
Magidanga, Vailet [ORNL] (ORCID:0009000410659658)↗
SEARCH · Engineering Papers
Results for “Liquid dynamics”
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.
Structural studies on choline-carboxylate bio-ionic liquids by x-ray scattering and molecular dynamics
We report a X-ray diffraction and molecular dynamics study on three choline-based bio-ionic liquids, choline formate, [Ch] [For], choline propanoate, [Ch][Pro], and choline butanoate, [Ch][But]. For the first time, this class of ionic liquids has been investigated by X-ray diffraction. Experimental and theoretical structure factors have been compared for each term of the series. Local structural organization has been obtained from ab initio calculations through static models of isolated ion pairs and dynamic simulations of small portions of liquids through twelve, ten, and nine ion pairs for [Ch][For], [Ch][Pro], and [Ch][But], respectively. All the theoretical models indicate that cations and anions are connected by strong hydrogen bonding and form stable ion pairs in the liquid that are reminiscent of the static ab initio ion pairs. Different structural aspects may affect the radial distribution function, like the local structure of ion pairs and the conformation of choline. When small portions of liquids have been simulated by dynamic quantum chemical methods, some key structural features of the X-ray radial distribution function were well reproduced whereas the classical force fields here applied did not entirely reproduce all the observed structural features.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A mathematical model for the study of slosh effects on rendezvous dynamics.
Mathematical model for studying liquid slosh effects on rendezvous dynamics by incorporating sloshing dynamics into vehicle motion equation
MATHEMATICAL MODEL↗
Manipulating molecular order in nematic liquid crystal capillary bridges via surfactant adsorption: guiding principles from dissipative particle dynamics simulations
The ability of liquid crystals (LCs) to change orientational order is used in applications, ranging from sensors to displays. The aim of this work is to computationally investigate how surfactant adsorption on cylindrical LC bridges can be used to control such orientational order. Building from classical fundamental lessons, understanding the ordering of mesogens along a preferred axis with the help of molecular modelling contributes to investigations of systems that could be a platform for LC-based sensing applications. The coarse-grained dissipative particle dynamics (DPD) simulation method is implemented here, because it allows us to quantify the effect of molecular features on the properties of meso-scopic systems containing LC bridges, an aqueous solvent, and surfactants at various concentrations. Three surfactant types are modelled with short, medium, and long tail lengths, respectively. All surfactants adsorb at the LC–water interface. It is found that the length of the surfactant hydrophobic tail determines the effectiveness by which the LC order is affected. Short tails are not as effective as long ones. Surfactants with long tails affect the LC order, but, in agreement with experiments, predominantly only within a short distance from the LC–water interface. For these surfactants, the surface density at the LC–water interface is an important knob that can be used to control the order of the LCs. As the effective LC–surfactant interactions change, so does the distribution of the surfactants at the interface. Consistent with theoretical expectations, the results presented here elucidate the effect of molecular features on the anchoring mechanism between surfactants and mesogens within cylindrical bridges dispersed in aqueous systems and could be helpful for designing novel surface-active compounds in the development of advanced sensing devices based on LCs.
74 ATOMIC AND MOLECULAR PHYSICS↗
Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores
Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations
We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Dynamics and statics of nonaxisymmetric and symmetric liquid bridges
This program of theoretical and experimental ground-based research focuses on the understanding of the dynamics and stability limits of nonaxisymmetric and symmetric liquid bridges. There are three basic objectives: First, to determine the stability limits of nonaxisymmetric liquid bridges held between non-coaxial parallel disks, Second, to examine the dynamics of nonaxisymmetric bridges and nonaxisymmetric oscillations of initially axisymmetric bridges. The third objective is to experimentally investigate the vibration sensitivity of liquid bridges under terrestrial and low gravity conditions. Some of these experiments will require a low gravity environment and the ground-based research will culminate in a definitive flight experiment.
Alexander, J. Iwan D.↗
Time Resolved Visualization of A Liquid Jet in an RDE Using MHz Rate Diesel PLIF
The characterization of the dynamic response of liquid jets to transient detonation wave passage is critical for optimization and modeling of liquid fueled rotating detonation combustors. In this work, a rotating detonation combustor (RDC) is operated on hydrogen and air to sustain stable detonation waves that acts as a detonation driver and interact in a one-way coupled manner with a single liquid fuel jet that propagates into the combustion chamber with cycle periods of ~ 250 μs. Diesel is utilized as a realistic fuel surrogate with higher aromatic compounds to enable fluorescence excitation, using the 355 nm third harmonic output of a burst-mode Nd:YAG laser, imaged at rates up to 1 MHz. By optimizing the technique to accommodate orders of magnitude variations in the fuel density throughout the injection process, the PLIF data enable quantitative measurements including the refill time, the relative recovery between liquid and gaseous jets and jet trajectory across various momentum flux ratio. As the passage of the detonation wave imparts significant changes in the momentum flux ratio, the qualitative liquid break-up process and spatial distribution varies significantly in time. As the injection system recovers ~ 70% of the cycle period and return to a quasi-steady position and allow comparisons with theoretical jet trajectories. These data, enabled by ultra-high-speed PLIF imaging, represent some of the first detailed measurements for quantifying the dynamic response and recovery of liquid jets exposed to periodic detonations in an operating RDC.
Propulsion↗
Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation
Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Insights into the Reactivity of Brookite TiO 2 Nanorods in Liquid Water from Ab Initio Molecular Dynamics Simulations
Brookite TiO 2 , a rare natural polymorph of TiO 2 , has been reported to be an excellent photocatalyst for the production of hydrogen from water and aqueous alcohol solutions, especially when it is reduced and synthesized in the form of nanorods. Here, we investigate the reactivity of stoichiometric and reduced brookite nanorods in liquid water using ab initio molecular dynamics and hybrid density functional theory calculations. Our simulations show a much higher water dissociation fraction on reduced nanorods than on stoichiometric ones, with an accumulation of the resulting bridging hydroxyls (O br H) and terminal hydroxyls (Ti –OH) on different facets of the nanorod. ObrH groups accumulate preferentially on low-energy (210) facets, where they are stabilized by adjacent reduced Ti (Ti 3+ ) sites, while Ti –OH groups prefer to form at the four-fold coordinated Ti atoms on high-energy (010) facets. This hydroxylation pattern also favors the spatial localization of excited holes on the (010) facets. This coupling between water-induced surface chemistry and charge separation underpins the enhanced photocatalytic activity of brookite nanorods, providing useful information for the design of more efficient TiO2-based nanostructures for solar-driven hydrogen evolution.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Liquid-based gating mechanism with tunable multiphase selectivity and antifouling behaviour
Living organisms make extensive use of micro- and nanometre-sized pores as gatekeepers for controlling the movement of fluids, vapours and solids between complex environments. The ability of such pores to coordinate multiphase transport, in a highly selective and subtly triggered fashion and without clogging, has inspired interest in synthetic gated pores for applications ranging from fluid processing to 3D printing and lab-on-chip systems(1-10). But although specific gating and transport behaviours have been realized by precisely tailoring pore surface chemistries and pore geometries(6,11-17), a single system capable of controlling complex, selective multiphase transport has remained a distant prospect, and fouling is nearly inevitable(11,12). Here we introduce a gating mechanism that uses a capillary-stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state. Theoretical modelling and experiments demonstrate that for each transport substance, the gating threshold-the pressure needed to open the pores-can be rationally tuned over a wide pressure range. This enables us to realize in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow us to dynamically modulate gas-liquid sorting in a microfluidic flow and to separate a three-phase air-water-oil mixture, with the liquid lining ensuring sustained antifouling behaviour. Because the liquid gating strategy enables efficient long-term operation and can be applied to a variety of pore structures and membrane materials, and to micro- as well as macroscale fluid systems, we expect it to prove useful in a wide range of applications.
Hou, X↗
Confining Liquids inside Carbon Nanotubes: Accelerated Molecular Dynamics with Spliced, Soft-Core Potentials and Simulated Annealing
Understanding emergent phenomena of fluids under physical confinement requires the development of advanced tools for rapid and accurate simulation of their physiochemical properties. Simulating liquid molecules commensurate in size with the nanoscale enclosures that confine them is a key challenge. In this work, we demonstrate an accelerated molecular dynamics simulation technique that combines soft-core potentials (SCP) and simulated annealing (SA) to analyze confined liquids. This integrated SCP/SA method relies on a new spliced soft-core potential (SSCP), which enables tunable accuracy with respect to the target hard-core potential (HCP). SCP/SA enables the packing of enclosures with bulk material in a controlled, thermodynamically consistent manner. The enhanced SSCP accuracy is a critical feature of SCP/SA, enabling a smooth transition between the SCP and the HCP at a desired SCP hardness. We applied SCP/SA to the problem of filling a carbon nanotube (CNT) in periodic boundary conditions with a popular ionic liquid (IL), 1- butyl-3-methylimidazolium hexafluorophosphate $[\text{BMIM}^+][\text{PF}^–_6]$. We performed a series of triplicate simulations on systems with varying CNT diameter and charge to demonstrate SCP/SA’s versatility. Beyond this IL/CNT system, the SCP/SA simulation framework has a broad range of potential applications, not limited to nanoscale enclosures and interfaces, including both solid-state and biological systems.
36 MATERIALS SCIENCE↗
Effects of polarizability on the structural and thermodynamics properties of [C{sub n}mim][Gly] ionic liquids (n = 1–4) using EEM/MM molecular dynamic simulations
An extended electronegativity equalization method/molecular mechanics (EEM/MM) model for ionic liquids is used to investigate the structures and properties of 1-alkyl-3-methylimidazolium glycine ionic liquids [C{sub n}mim][Gly] (n = 1–4) with alkyl substituents of different lengths. The EEM/MM model describes the electrostatic interactions of atoms and their changes in different ambient environments. This property is the most outstanding characteristic of the model. EEM parameters (i.e., valence electronegativities and valence hardness parameters) are calibrated using linear regression and least-squares methods, which can accurately predict the gas-phase properties of [C{sub n}mim]{sup +}, [Gly]{sup −}, and [C{sub n}mim][Gly] ion pairs. We utilize the EEM/MM force field to systematically investigate the effects of polarizability on the accuracy of [C{sub n}mim][Gly] properties predicted through the molecular dynamic simulations. EEM/MM explicitly describes the atom-based polarizability of [C{sub n}mim][Gly]; thus, the densities, enthalpies of vaporization, self-diffusion coefficients, and conductivities of the [C{sub n}mim][Gly] are consistent with the experimental values. The calculated radial distribution functions provide a mechanistic understanding of the effects of polarizability on ionic aggregations in amino acid ionic liquids. The effects of alkyl chain length on the diffusion coefficient and conductivity are also discussed.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Phonon dynamics in the generalized Kitaev spin liquid
Candidate materials for the Kitaev spin liquid generically have residual interactions beyond the Kitaev coupling. It therefore becomes necessary to understand how signatures of the quantum spin liquid, e.g., fractionalization of the spin excitations, are affected by the presence of these interactions. Recently it was shown that phonon dynamics is an indirect but effective probe to study fractionalized excitations in the Kitaev spin liquid. Ultrasound experiments can measure sound attenuation, which should show characteristic temperature and angular dependence of the sound-attenuation coefficient if the scattering of phonons happens predominantly on Majorana fermions. So far the computation of the sound attenuation was only done in the pure spin-phonon coupled Kitaev model, without taking into account residual interactions. In order to understand experimental signatures, here we present a mean-field study of the sound attenuation in the generalized J-K_Γ model, which is relevant to candidate materials. Further, our findings show that as long as the system is in the spin liquid phase, characteristic features of the sound attenuation remain observable even in the presence of residual interactions.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Coarsening Dynamics and Marangoni Effects in Thin Liquid Crystal Bubbles in Microgravity
The Observation and Analysis of Smectic Islands in Space (OASIS) flight hardware was successfully launched on SpaceX-6 on April 15, 2015 and was operated in the Microgravity Science Glovebox (MSG) on board the International Space Station (ISS). The OASIS project comprises a series of experiments that probe the interfacial and hydrodynamic behavior of spherical-bubble freely suspended liquid crystal (FSLC) membranes in space. These are the thinnest known stable condensed phase structures, making them ideal for studies of two-dimensional (2D) coarsening dynamics and thermocapillary phenomena in microgravity. The OASIS experimental investigation was carried out using four different smectic A and C liquid crystal materials in four separate sample chambers housed inside the MSG. In this report, we present the behavior of collective dynamics on 2D bubble surface, including the equilibrium spatial organization and interaction of islands in electric fields and temperature gradients, and the diffusion and coalescence-driven coarsening dynamics of island emulsions in microgravity. We have observed spontaneous bubble thickening behavior caused by gradients between the bubble-blowing needle and ambient air temperatures. A uniform, thicker band forms during coarsening as a result of non-uniform heating by the LED illumination panels. These are proposed to be a result of Marangoni convection on the bubble surface.
microgravity↗
A Method for Obtaining Liquid–Solid Adsorption Rates from Molecular Dynamics Simulations: Applied to Methanol on Pt(111) in H 2 O
Adsorption is an important step in heterogeneous catalysis as it predetermines how many reactant molecules can participate in a surface reaction per unit time. While the rate of adsorption processes is well studied in gas–solid adsorption in both theory and experiment, such rates are still not well studied for liquid–solid adsorption. This is partly because the ever-changing configurations of liquid-phase solvent molecules impede the ability to study a molecule approaching a surface from a liquid phase by either experiment or theory. In this work, we develop a method using molecular dynamics (MD) simulations to study the rate of adsorption in liquid–solid adsorption processes. Specifically, we use MD to model the diffusion of a methanol molecule in aqueous solvent and its adsorption to a Pt(111) surface. We find that by approximating the solute motion as following the same displacement rates as a random walk model, the adsorbed and non-adsorbed states of the methanol molecule near the Pt(111) surface can be discerned and quantified. In particular, this methodology enables extracting a sticking coefficient and a macroscopically relatable adsorption rate. This method can be applied to arbitrary types of reactants and surfaces, as well as different liquid environments, thus providing a general tool for predicting quantitative adsorption rates of liquid–solid adsorption systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structural ordering at solid-liquid interfaces in Al-Sm system: A molecular-dynamics study
The structural ordering at solid-liquid interfaces far from equilibrium is studied with molecular dynamics simulations for the Al-Sm system. Using the van-Hove self-correlation function as the criterion to identify attachment/detachment events that occur at the interface, we are able to determine the time-dependent interface position, and characterize the detailed interfacial structure ordering surrounding the attached atoms. For the interface between an undercooled Al90Sm10 liquid and a metastable cubic structure, the solid induces the crystalline order of the cubic phase in the liquid layers, promoting the continuous growth of the crystal phase. When the same liquid is put in contact with f.c.c. Al, Sm from the liquid can still attach to the solid interface despite its insolubility in the Al lattice. Non-f.c.c. order is revealed surrounding the attached Sm atoms. Lastly, we show that the local structure ordering at interface is highly correlated to solid packing and liquid ordering.
36 MATERIALS SCIENCE↗
Dynamical Structure Factor of the Three-Dimensional Quantum Spin Liquid Candidate NaCaNi 2 F 7
We study the dynamical structure factor of the spin-1 pyrochlore material NaCaN 2 F 7 , which is well described by a weakly perturbed nearest-neighbour Heisenberg Hamiltonian, Our three approaches—molecular dynamics simulations, stochastic dynamical theory, and linear spin wave theory—reproduce remarkably well the momentum dependence of the experimental inelastic neutron scattering intensity as well as its energy dependence with the exception of the lowest energies. Here, we discuss two surprising aspects and their implications for quantum spin liquids in general: the complete lack of sharp quasiparticle excitations in momentum space and the success of the linear spin wave theory in a regime where it would be expected to fail for several reasons.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗