Studies of the flow, bonding and damping characteristics of a squeeze film under dynamic conditions Final report, 15 Jun. 1968 - 31 Jun. 1970
Flow, bonding and damping characteristics of liquid squeeze film under dynamic conditions
SEARCH · Engineering Papers
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.
Flow, bonding and damping characteristics of liquid squeeze film under dynamic conditions
We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.
A physics-based mathematical model for simulating heat and moisture transport on the skin and within clothing was developed to improve the accuracy of human thermal modeling predictions. In a previous study, a set of sweating hotplate experiments was specifically created to obtain test data for model validation. The clothing moisture model was subsequently coupled with a human thermal model and tested by simulating human experiments described in the open literature, which served to verify the model’s ability to predict whole body thermo-physiologic state in work-rest cycles during which sweating was known to occur. In this study, a human subject test was specifically designed to validate the model’s ability to predict local and whole body sweat evaporation and accumulation on the surface of the body. The human subject test incorporated twelve clothed male subjects of moderate physical fitness, ages 18-30, which were prescribed a work-rest cycle in a climate-controlled chamber (25 °C, 50% RH). Two activity levels were considered during the active stage (4 MET vs. 6.5 MET), and two sets of ending ambient conditions were considered for the final rest stage (25 °C vs. 15 °C, 50% RH). Model predictions were compared to the measured local skin temperatures, core temperatures, and whole body sweat losses. In general, the model was able to reproduce the trends observed in the human subject trial, but more importantly, provided insight into how dynamic changes in liquid content can impact the evolution of local and whole body thermo-physiologic state. This study demonstrates the ability of an advanced moisture model to capture distinct physics attributable to moisture accumulation, allowing for increased accuracy in predicting human thermal state and for interpreting the results of measurements derived from human subject tests.
The Rectilinear Control System can be used to simulate a launch vehicle during liftoff. Several control schemes have been developed that can control different dynamic models of the rectilinear plant. A robust control algorithm was developed that can control a pendulum to maintain an inverted position. A fluid slosh tank will be attached to the pendulum in order to test robustness in the presence of unknown slosh characteristics. The rectilinear plant consists of a DC motor and three carts mounted in series. Each cart s weight can be adjusted with brass masses and the carts can be coupled with springs. The pendulum is mounted on the first cart and an adjustable air damper can be attached to the third cart if desired. Each cart and the pendulum have a quadrature encoder to determine position. Full state feedback was implemented in order to develop the control algorithm along with a state estimator to determine the velocity states of the system. A MATLAB program was used to convert the state space matrices from continuous time to discrete time. This program also used a desired phase margin and damping ratio to determine the feedback gain matrix that would be used in the LabVIEW program. This experiment will allow engineers to gain a better understanding of liquid propellant slosh dynamics, therefore enabling them to develop more robust control algorithms for launch vehicle systems
The promise of high energy density lithium–sulfur batteries with long cycle life is currently tempered by the rapid degradation of lithium-metal anodes with cycling. An in-depth understanding of the dynamical behavior in liquid electrolytes, including the mechanisms underlying depletion of lithium inventory and evolution of lithium interphases, is crucial to make Li–S batteries a reality. We use here an anode-free full cell configuration, pairing a Li 2 S cathode with a bare nickel current collector with no lithium metal on it, to quantitatively estimate the lithium inventory loss per cycle. Lithium inventory loss is shown to be the main factor limiting the overall cyclability of Li–S batteries. Time-of-flight secondary ion mass spectrometry measurements on the deposited lithium reveal the presence of substantial metallic lithium even after most of the active lithium inventory has been depleted. The trapped metallic lithium is rendered electrochemically inactive by the growth of a resistive electrolyte decomposition interphase on the lithium surface. The bulk of the deposited lithium is shown to be composed of various fully reduced interphasial components, including several hydrogen-containing species that show a substantial reduction in intensity with cycling. This indicates considerable gas evolution and is also correlated with the loss of lithium inventory. The use of an anode-free full cell configuration provides a framework for accurate assessment of the dynamics of lithium inventory depletion and characterization of the accompanying interphasial evolution with cycling. The insights gained will prove invaluable to the development of strategies for extending the cycle life of energy-dense Li–S batteries.
This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.
Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.
Tests of cryogenic chilldown in a vertical stainless-steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters a room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFDbased modeling framework employing the mixture multiphase model was developed to simulate this chilldown process and was applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a range of test conditions are simulated to assess the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and considered in the CFD framework. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. In addition, simulated distributions of the vapor volume fraction are presented, but the current lack of experimental volume fraction data prevents comparisons.
Tests of cryogenic chilldown in a vertical stainless steel tube are simulated using computational fluid dynamics (CFD). Subcooled liquid nitrogen enters the room-temperature tube at one end, resulting in film boiling and eventual quenching of the tube walls. A CFD-based modeling framework employing the mixture multiphase model has been developed to simulate this chilldown process and applied to validate it against experimental data. Extensive sensitivity studies have been performed to create a consistent modeling framework, and a two test conditions are simulated to test the modeling capability. The importance of the rewetting temperature and its relation to the foam limit temperature is discussed and is considered in the CFD model. The presented results demonstrate a qualitatively consistent representation of the film boiling and quenching regime. However, the quantitative accuracy of the simulations varies with the test conditions. Simulated vapor volume fraction distributions are also presented, although the current lack of experimental volume fraction data does not permit a comparison.
CHESS, chopper spectrometer examining small samples, is a planned direct geometry neutron chopper spectrometer designed to detect and analyze weak signals intrinsic to small cross sections (e.g., small mass, small magnetic moments, or neutron absorbing materials) in powders, liquids, and crystals. CHESS is optimized to enable transformative investigations of quantum materials, spin liquids, thermoelectrics, battery materials, and liquids. The broad dynamic range of the instrument is also well suited to study relaxation processes and excitations in soft and biological matter. The 15 Hz repetition rate of the Second Target Station at the Spallation Neutron Source enables the use of multiple incident energies within a single source pulse, greatly expanding the information gained in a single measurement. Furthermore, the high flux grants an enhanced capability for polarization analysis. This enables the separation of nuclear from magnetic scattering or coherent from incoherent scattering in hydrogenous materials over a large range of energy and momentum transfer. This paper presents optimizations and technical solutions to address the key requirements envisioned in the science case and the anticipated uses of this instrument.
We report a study that elucidates key physicochemical phenomena that underlie the dynamics of adsorbate-induced ordering transitions in micrometer-thick films of nematic liquid crystals (LCs) supported on chemically functionalized surfaces. By flowing a gas containing dimethyl methylphosphonate (DMMP) over a film of nematic 4-cyano-4'-pentylbiphenyl (5CB) supported on a surface decorated with aluminum perchlorate, we quantify the dynamics of the optical response (intensity of transmitted light) of the LC film as a function of key experimental parameters, including linear flow rate of the gas over the surface of the film and the concentration of the DMMP within the gas. Building from prior infrared spectroscopic analysis revealing that the optical response of the LC is due to an ordering transition triggered by competitive ligand-exchange reactions at the chemically functionalized surface (i.e., between the phosphonate and nitrile groups of the DMMP and LC, respectively, with the aluminum perchlorate-decorated surface), we interpret the measurements reported herein within the framework of a simple model of mass transport. The model reveals that the rate-limiting process underlying the response of the LC is the mass transport of the DMMP across the concentration boundary layer formed on the vapor side of the LC–vapor interface, and not the mass transport of DMMP across the LC film nor the intrinsic ligand exchange kinetics. The analysis also reveals that the response of the LC is triggered once the concentration of DMMP within the LC film rises to approximately 0.25 mM. By flowing a stream of gas free of DMMP over a LC film pre-exposed to DMMP, we established that the relaxation of the LC film to the initial state is also largely controlled by mass transport of the DMMP across the vapor side of the LC–vapor interface. Analysis of the relaxation of the system leads to an estimate of the sensitivity of the supported LC film to the concentration of DMMP in the gas phase (approximately 150 ppb for the non-optimized system reported in this paper). These results and others reported in this paper identify key physical phenomena that dictate the dynamics of adsorbate-induced ordering transitions in LC systems, and provide important guidance to the design of LCs as chemo-responsive materials for use in sensors.
We consider the genesis and dynamics of interfacial instability in vertical gas-liquid flows, using as a model the two-dimensional channel flow of a thin falling film sheared by counter-current gas. The methodology is linear stability theory (Orr-Sommerfeld analysis) together with direct numerical simulation of the two-phase flow in the case of nonlinear disturbances. We investigate the influence of two main flow parameters on the interfacial dynamics, namely the film thickness and pressure drop applied to drive the gas stream. To make contact with existing studies in the literature, the effect of various density contrasts is also examined. Energy budget analyses based on the Orr-Sommerfeld theory reveal various coexisting unstable modes (interfacial, shear, internal) in the case of high density contrasts, which results in mode coalescence and mode competition, but only one dynamically relevant unstable interfacial mode for low density contrast. A study of absolute and convective instability for low density contrast shows that the system is absolutely unstable for all but two narrow regions of the investigated parameter space. Direct numerical simulations of the same system (low density contrast) show that linear theory holds up remarkably well upon the onset of large-amplitude waves as well as the existence of weakly nonlinear waves. For high density contrasts, corresponding more closely to an air-water-type system, linear stability theory is also successful at determining the most-dominant features in the interfacial wave dynamics at early-to-intermediate times. Nevertheless, the short waves selected by the linear theory undergo secondary instability and the wave train is no longer regular but rather exhibits chaotic motion. The same linear stability theory predicts when the direction of travel of the waves changes — from downwards to upwards. We outline the practical implications of this change in terms of loading and flooding. The change in direction of the wave propagation is represented graphically in terms of a flow map based on the liquid and gas flow rates and the prediction carries over to the nonlinear regime with only a small deviation.
We performed long-time all-atom molecular dynamics simulations of cationic polymerized ionic liquids with eight mobile counterions, systematically varying size and shape to probe their influence on the decoupling of conductivity from polymer segmental dynamics. We demonstrated rigorous identification of the dilatometric glass-transition temperature ( T g ) for polymerized ionic liquids using an all-atom force field. Polymer segmental relaxation rates are presumed to be consistent for different materials at the same glass-transition-normalized temperature ( T g / T ), allowing us to extract a relative order of decoupling by examining conductivity at the same T g / T . Size, or ionic volume, cannot fully explain decoupling trends, but within certain geometric and chemical-specific classes, small ions generally show a higher degree of decoupling. This size effect is not universal and appears to be overcome when structural results reveal substantial coordination delocalization. Here, we also reveal a universal inverse correlation between ion-association structural relaxation time and absolute conductivity for these polymerized ionic liquids, supporting the ion-hopping interpretation of ion mobility in polymerized ionic liquids.
Calculations of charged systems in periodic boundary conditions (PBC) are problematic because there are spurious interactions between the charges in different periodic images that can affect the physical picture. In addition, the intuitive limit of Coulomb interactions decaying to zero as the interacting charges are placed at infinite separation no longer applies, and for example total energies become undefined. Leveraging subsystem density functional theory (also known as density embedding) we define an impurity model that embeds a finite neutral or charged subsystem within an extended (infinite) surrounding subsystem. The combination of the impurity model and a consistent choice of the Coulomb reference provides us with an algorithm for evaluating the ionization potential (IP) in extended systems. Additionally, we demonstrate our approach in a pilot calculation of the IP of liquid water, based on a configuration from a prior ab initio molecular dynamics (AIMD) simulation of liquid water (Genova et al., J. Chem. Phys. 2016, 144, 234105). The calculations with the impurity model capture the broadening on the ionization energies introduced by the interactions between the water molecules. Furthermore, the calculated average IP value (10.5 eV) compare favorably to experiments (9.9-10.06 eV) and very recent simulations based on the GW approximation (10.55 eV), while at the same time outperforming density embedding calculations carried out with a naïve handling of the electrostatic interactions (about 7 eV).
Here, the morphology of DNA is strongly influenced by its surrounding environment, including factors such as pH, salt type and valency, and the presence of polymers. Inorganic salts are known to reduce the DNA chain length through mechanisms like electrostatic screening and ion bridging. In contrast, ionic liquids, a new class of organic salts, have previously been found to increase the DNA chain length, indicating a distinct mode of interaction between the ionic liquid and DNA chains. This study utilizes self-assembled DNA-AuNPs as a model system to examine changes in the DNA chain morphology and the nanoscale interaction mechanisms in an ionic liquid environment. The DNA chain lengths are measured in solution using X-ray scattering measurements at varying concentrations of two imidazolium ([$BMIM$] acetate and [$EMIM$] acetate) based ionic liquids. Additionally, Molecular Dynamics (MD) simulations are performed mimicking the experimental system. Our results suggest an interplay of electrostatic and groove-binding interactions governing the DNA chain morphology, which depends on IL concentration and the composition of the DNA chains. It has been found that for DNA chains with majority ssDNA, electrostatic interaction dominate, however with increasing composition of double strands, the DNA chains exhibit compaction due to a non-electrostatic hydrophobic groove-binding mechanism.
Not provided.
Materials are usually classified as solids or liquids, based on their structural stability, dynamic response, and rheological properties. Structured liquid, a new state of matter, has attracted much attention in recent years. Different with either solid or liquid, structured liquid combines the desirable characteristics of fluids with the structural stability of a solid, showing a myriad of potential applications in encapsulation, biphasic reactors, and programmable liquid constructs. Here, a brief review is given, by introducing a new strategy to structure liquids based on the formation, assembly, and jamming of nanoparticle surfactants (NPSs) at liquid-liquid interfaces. Furthermore, the interfacial packing of the NPSs can be effectively manipulated using external triggers, endowing the structured liquids with adaptiveness and responsiveness to changes in their external environment.
Highlights: • Chiral transformations with account for chiral anomaly were applied to the CDW model. • Two-fluid hydrodynamics was constructed of the order parameter and the normal liquid. • Topologically nontrivial dynamics appears under applied fields or charge injection. • Spontaneous spacetime vorticity under the driving current was modeled numerically. • Periodic sequence of amplitude nodes with 2π phase increments was demonstrated. - Abstract: Many recent experiments addressed manifestations of electronic crystals, particularly the charge density waves, in nano-junctions, under electric field effect, at high magnetic fields, together with real space visualizations by STM and micro X-ray diffraction. This activity returns the interest to stationary or transient states with static and dynamic topologically nontrivial configurations: electronic vortices as dislocations, instantons as phase slip centers, and ensembles of microscopic solitons. Describing and modeling these states and processes calls for an efficient phenomenological theory which should take into account the degenerate order parameter, various kinds of normal carriers and the electric field. Here we notice that the commonly employed time-depend Ginzburg–Landau approach suffers with violation of the charge conservation law resulting in unphysical generation of particles which is particularly strong for nucleating or moving electronic vortices. We present a consistent theory which exploits the chiral transformations taking into account the principle contribution of the fermionic chiral anomaly to the effective action. The resulting equations clarify partitions of charges, currents and rigidity among subsystems of the condensate and normal carriers. On this basis we perform the numerical modeling of a spontaneously generated coherent sequence of phase slips – the spacetime vortices – serving for the conversion among the injected normal current and the collective one.