Engineering Papers⌕ Search

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.

At least 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Stabilized Electrode/Electrolyte Interphase by a Saturated Ionic Liquid Electrolyte for High-Voltage NMC532/Si-Graphite Cells

Nonaqueous electrolyte has become one of the technical barriers in enabling Li-ion battery comprising of a high voltage cathode and high capacity anode. Here, we demonstrate a saturated piperidinum bis(fluorosulfonyl)imide ionic liquid (IL) with a LiFSI salt not only supports the redox reaction on the cathode at high voltages, but also shows exceptional kinetic stability on the lithiated anode as evidenced by its improved cycling performance in a NMC532/Si-graphite full cells cycled between 4.6 and 3.0 V. On the basis of the spectroscopic/microscopic analysis and molecular dynamics (MD) simulations, the superior performance of the cells is attributed to the formation of solid-electrolyte-interphase on both electrode as well as unique solvation structure where a deadlocked coordination network is established at the saturated state, which prevents transition metal dissolution into the electrolyte via a solvation process.

25 ENERGY STORAGE↗
DOE OSTI2020

A review of the application of 2D isotropic-anisotropic correlation NMR spectroscopy in structural studies of chalcogenide glasses

The application of solid-state high-resolution NMR spectroscopy in the structural investigation of chalcogenide glasses in Ge/As/P/Si-X (X = S, Se,Te) systems has remained challenging even for the spin-1/2 nuclides ( 29 Si, 31 P, 77 Se, 125 Te), owing to their low natural abundance (except for 31 P), slow spin-lattice relaxation rate and large CSA and chemical shift distribution induced line broadening effects. However, most of these deleterious effects can be successfully overcome in two-dimensional (2D) isotropic-anisotropic correlation NMR experiments, especially when performed at high magnetic field and in conjunction with the Car-Purcell-Meiboom-Gill (CPMG) echo train acquisition. Here, we present a short introduction to the basic principles of such experiments and review their applications over the last decade in deciphering various short- and intermediate- range structural characteristics of chalcogenide glasses in S-Se, Se-Te, Ge-Se, As-Se and Si-Se systems as well as in investigating the molecular dynamics in a P-Se supercooled liquid. We anticipate possible future applications of these 2D isotropic-anisotropic correlation NMR experiments, particularly in conjunction with density functional theory-based calculations of NMR chemical shift tensor parameters and additional signal enhancement schemes, in addressing complex structural correlations and distributions in chalcogenide glasses.

2D NMR↗
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2020

Reconfigurable structure and tunable transport in synchronized active spinner materials

Ensembles of actuated colloids are excellent model systems to explore emergent out-of-equilibrium structures, complex collective dynamics, and design rules for the next generation materials. Here, we demonstrate that ferromagnetic microparticles suspended at an air-water interface and energized by an external rotating magnetic field spontaneously form dynamic ensembles of synchronized spinners in a certain range of the excitation field parameters. Each spinner generates strong hydrodynamic flows, and collective interactions of the multiple spinners promote a formation of dynamic lattices. On the basis of experiments and simulations, we reveal structural transitions from liquid to nearly crystalline states in this novel active spinner material and demonstrate that dynamic spinner lattices are reconfigurable, capable of self-healing behavior and that the transport of embedded inert cargo particles can be remotely tuned by the parameters of the external excitation field. Our findings provide insights into the behavior of active spinner materials with reconfigurable structural order and tunable functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large-Scale Atomistic Simulations: Investigating Free Expansion

The experiment investigates free expansion of a supercritical fluid into a two-phase liquid-vapor coexistence region. A huge molecular dynamics simulation (6 billion Lennard-Jones atoms) was run on 5760 GPUs (33% of LLNL Sierra) using LAMMPS/Kokkos software. This improved visualization workflow and started preliminary simulations of aluminum using SNAP machine learning potential.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The role of defects and solid/liquid interfacial interactions on controlling anisotropic growth of novel two-dimensional materials

Atomic-scale understanding about how interfacial properties are modified in the presence of defects and how they are correlated to the properties of nucleation and growth is critical for controlled materials growth on solid/liquid interfaces. The goal of this proposal is to establish preliminary data on the critical role of surface defects and chemical compositions of solvent on the structure, dynamics, reactivity, and free energetics at graphene/ionic liquid (IL) interfaces and their impact on two-dimensional (2D) growth of transition metal oxides on the surface. This study aims to develop generalized physical/chemical principles for directional growth at solid/liquid interfaces to achieve atomic-level control of two-dimensional materials synthesis. The proposed study of the effects of surface defects on the interfacial structure and reactivity will improve our understanding of the critical roles of heterogeneity, disorder, and solvents on 2D materials synthesis and lead to better control of directional growth of materials at solid/liquid interfaces.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Tribology experiment

A two-dimensional concept for Spacelab rack 7 was developed to study the interaction of liquid lubricants and surfaces under static and dynamic conditions in a low-gravity environment fluid wetting and spreading experiments of a journal bearing experiments, and means to accurately measure and record the low-gravity environment during experimentation are planned. The wetting and spreading process of selected commercial lubricants on representative surface are to the observes in a near-zero gravity environment.

Wall, W. A.↗
NASA NTRS1988

New theoretical results for the Lehmann effect in cholesteric liquid crystals

The Lehmann effect arising in a cholesteric liquid crystal drop when a temperature gradient is applied parallel to its helical axis is investigated theoretically using a local approach. A pseudoscalar quantity is introduced to allow for cross couplings which are absent in nematic liquid crystals, and the statics and dissipative dynamics are analyzed in detail. It is shown that the Lehmann effect is purely dynamic for the case of an external electric field and purely static for an external density gradient, but includes both dynamic and static coupling contributions for the cases of external temperature or concentration gradients.

Brand, Helmut R.↗
NASA NTRS1990

A dynamic technique for measuring surface tension at high temperatures in microgravity environment

The feasibility of a dynamic technique for measuring surface tension of liquid metals at high temperatures in a microgravity environment was demonstrated. The basic method involves heating a tubular specimen resistively from ambient temperature through its melting point in about 1 sec by passing an electrical current pulse through it, while simultaneously recording the pertinent experimental quantities. Static equilibrium for the molten specimen is achieved in a microgravity environment by splitting the current after it passes through the specimen tube and returning a fraction along the tube axis, and the remaining fraction outside the specimen. Adjustments to the current split enable a balance between the magnetic and surface tension forces acting on the specimen. Values for surface tension are determined from measurements of the equilibrium dimensions of the molten specimen tube, and the magnitudes of the currents. Rapid melting experiments, performed during microgravity simulations with the NASA KC-135 aircraft, yield a value for the surface tension of copper at its melting point which is in agreement with literature data. Measurements of surface tension of a refractory metal (tantalum) are underway.

Miiller, A. P.↗
NASA NTRS1989

Design considerations for a micro-g superfluid helium fluid acquisition system

The general description, the operation, and the design of a superfluid helium (SFHe) fluid acquisition system (FAS) for use under microgravity conditions is presented. For the type of FAS considered here, where fine-mesh woven screens are used to retain flowing SFHe within a gallery arm (flow) channel, those forces which determine the flow dynamics are the micro-g accelerations, liquid surface tension, and tensile strength and cumulative pressure drops along a flow path that begins at the bulk liquid and ends at the entrance to a pump. For this case, the dimensionless number, N(T) is written as the ratio between the pressure drop across the screen and the surface tension forces at the screen for low fluid velocities. Static Bond number measurements have bene taken for SFHe using 325 x 2300 twilled Dutch screen and have indicated a screen pore hydraulic radius of 0.00031 cm.

Lee, J. M.↗
NASA NTRS1992

A numerical model for atomization-spray coupling in liquid rocket thrust chambers

The physical process of atomization is an important consideration in the stable operation of liquid rocket engines. Many spray combustion computational fluid dynamics (CFD) codes do not include an atomization sub-model but assume arbitrary drop size distributions, drop initial locations, and velocities. A method of coupling an atomization model with the spray model in a REFLEQS CFD code is presented. This method is based on a jet-embedding technique in which the equations governing the liquid jet core are solved separately using the surrounding gas phase conditions. The droplet initial conditions are calculated using a stability analysis appropriate for the atomization regime of liquid jet break-up. This novel coupling model is used to analyze the SSME fuel preburner single injector flow. Results of the diffusion flame characteristics in a single injection element will be presented. The effect of relative velocity, mixture ratio and droplet initial conditions will be shown. The predictions of present atomization model is compared with that of the widely used CICM correlation. The results are also compared with the predictions of volume-of-fluid method.

M. G. Giridharan↗
NASA NTRS1991

Longitudinal and transverse modes of slosh wave excitation in rotating dewar associated with gravity jitters

Characteristics of slosh waves based on the dynamical behavior of oscillations at the liquid-vapor interface have been investigated. Twelve case studies of slosh wave excitation due to various frequencies of gravity jitters under different rotating speeds of the propellant tank and different levels of background gravity environment have been simulated. The study shows that slosh waves excited inside the spacecraft propellant tank are characterized by the lowest frequency of the waves initiated, frequencies of the gravity jitters imposed on the propellant system, the levels of background gravity environment, and dewar rotating speeds. Conditions for suppression and amplification of the slosh waves are discussed.

Hung, R. J.↗
NASA NTRS1995

Effect of baffle on slosh reaction forces in rotating liquid helium subjected to a lateral impulse in microgravity

Sloshing dynamics within a partially filled rotating Dewar of superfluid He II are investigated in response to a lateral impulse. The study investigates several factors, including how the rotating bubble of superfluid He II reacts to the impulse in microgravity, how the amplitudes of slosh reaction forces act on the Dewar with various rotating speeds, how the frequencies of the sloshing modes excited differ in terms of differences in rotating speeds, and how the sloshing dynamics differ with and without a baffle. The numerical computation of sloshing dynamics is based on the noninertial frame spacecraft-bound coordinates. Results of the simulations are illustrated.

Hung, R. J.↗
NASA NTRS2000

Towards Flange-to-Flange Turbopump Simulations for Liquid Rocket Engines

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbopump geometry through numerical simulation will be of significant value toward design. This will help to improve safety of future space missions. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the MPI and MLP versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology. INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbopump simulations, moving boundary capability and an efficient time-accurate integration methods are build in the flow solver. To handle the geometric complexity and moving boundary problems, overset grid scheme is incorporated with the solver that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. The performance of the two time integration schemes for time-accurate computations is investigated. For an unsteady flow which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiterations procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in our computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive. The current geometry for the LOX boost turbopump has various rotating and stationary components, such as inducer, stators, kicker, hydrolic turbine, where the flow is extremely unsteady. Figure 1 shows the geometry and computed surface pressure of the inducer. The inducer and the hydrolic turbine rotate in different rotational speed.

Kiris, Cetin↗
NASA NTRS1999

Drop Ejection From an Oscillating Rod

The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.

Wilkes, E. D.↗
NASA NTRS2000

Supercomputing Aspects for Simulating Incompressible Flow

The primary objective of this research is to support the design of liquid rocket systems for the Advanced Space Transportation System. Since the space launch systems in the near future are likely to rely on liquid rocket engines, increasing the efficiency and reliability of the engine components is an important task. One of the major problems in the liquid rocket engine is to understand fluid dynamics of fuel and oxidizer flows from the fuel tank to plume. Understanding the flow through the entire turbo-pump geometry through numerical simulation will be of significant value toward design. One of the milestones of this effort is to develop, apply and demonstrate the capability and accuracy of 3D CFD methods as efficient design analysis tools on high performance computer platforms. The development of the Message Passage Interface (MPI) and Multi Level Parallel (MLP) versions of the INS3D code is currently underway. The serial version of INS3D code is a multidimensional incompressible Navier-Stokes solver based on overset grid technology, INS3D-MPI is based on the explicit massage-passing interface across processors and is primarily suited for distributed memory systems. INS3D-MLP is based on multi-level parallel method and is suitable for distributed-shared memory systems. For the entire turbo-pump simulations, moving boundary capability and efficient time-accurate integration methods are built in the flow solver, To handle the geometric complexity and moving boundary problems, an overset grid scheme is incorporated with the solver so that new connectivity data will be obtained at each time step. The Chimera overlapped grid scheme allows subdomains move relative to each other, and provides a great flexibility when the boundary movement creates large displacements. Two numerical procedures, one based on artificial compressibility method and the other pressure projection method, are outlined for obtaining time-accurate solutions of the incompressible Navier-Stokes equations. The performance of the two methods is compared by obtaining unsteady solutions for the evolution of twin vortices behind a flat plate. Calculated results are compared with experimental and other numerical results. For an unsteady flow, which requires small physical time step, the pressure projection method was found to be computationally efficient since it does not require any subiteration procedure. It was observed that the artificial compressibility method requires a fast convergence scheme at each physical time step in order to satisfy the incompressibility condition. This was obtained by using a GMRES-ILU(0) solver in present computations. When a line-relaxation scheme was used, the time accuracy was degraded and time-accurate computations became very expensive.

Kwak, Dochan↗
NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗