Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID MECHANICS”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Influence of magnetic angle on the E × B drift in a magnetic presheath

In this work, we provide experimental measurements of the E X B flow above a non-magnetic, absorbing boundary in a magnetized plasma. Measurements are taken as a function of the angle between the background magnetic field and the boundary normal. The measurements are compared to computational predictions for oblique magnetic fields based on Chodura's model. Ion flow measurements were obtained with laser induced fluorescence and the presheath potential structure was measured with an emissive probe. The ions were found to accelerate to nearly 30% of the sound speed parallel to the boundary at oblique angles of the magnetic field (~80°) and the ion speed at the electrostatic sheath edge was found to decrease with angle. The edge of the magnetic presheath has been experimentally determined based on the appearance of the E X B drift and is found to be independent of the angle of the magnetic field. We also demonstrate that laser induced fluorescence is suitable as a non-perturbative diagnostic to measure the electric field in a magnetized sheath.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Multiphysics simulation of microscale copper printing by confined electrodeposition using a nozzle array

3D printing of metals at the microscale and nanoscale is crucial to produce high-density interconnects and intricate structures in electronic devices. Conventional melting and sintering processes are not suitable for these scales due to a reliance on individual metal particles in the size range of tens of micrometers. Confined electrodeposition (CED) is an established alternative to conventional metal 3D printing processes in which an electrolyte is used to selectively induce deposition of the metal on the printing surface. However, commercialization and efficiency of this process have been limited due to a reliance on sub-micrometer nozzles to achieve desirable deposition rates and single nozzle to achieve uniformity of printed structures. Here, we address these challenges by computationally analyzing an array of microscale nozzles. We tailor the convection within the electrolyte to alter both deposition rate and geometric uniformity of the printed structures. The results show that for large nozzles the evaporation alone is not sufficient to obtain high deposition rates, yet an external pressure can be used to increase deposition and alter uniformity (thickness) of printed structures. Our results can be used to design and analyze new experiments toward parallel multi-nozzle deposition using CED toward high-throughput metal printing.

3D printing↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗
DOE OSTI2025

Planar shock-induced bubble collapse and jetting in water captured via X-ray phase contrast imaging

Shock wave–bubble interactions in water manifest rich dynamics driven by a combination of strong pressure and density mismatches. They have a wide variety of applications, including the injection of pharmaceuticals, and through scaling, enable the exploration of various aspects of high-energy-density systems such as inertial confinement fusion. In this work, the interaction between a micrometric nitrogen bubble and a planar shock wave, characterized by a Mach number of M = 1.24 and a peak pressure of p max = 0.57 , is experimentally recorded using ultra-high-speed x-ray phase contrast imaging. Highly resolved radiographs provide access to all phase discontinuities along the beam path, offering quantities such as the time-varying bubble size, the speed of a jet produced during the bubble collapse, and the time evolution of the shock wave front, which are critical benchmark data for numerical scheme validation. This study addresses the lack of well-characterized, repeatable, and high spatiotemporal resolution experiments at negative Atwood numbers by providing shock–bubble visualization and corresponding numerical simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Transport at high ${\beta_p}$ and development of candidate steady state scenarios for ITER

On DIII-D, the high β p scenario has an internal transport barrier (ITB), β N ~β p ~3,q 95 ~10, and very high normalized confinement H 98,y2 ~1.6. Recently, plasmas starting with these conditions have been dynamically driven to q 95 ~6 and β p ~2, where we find the ITB and high performance persist for five energy confinement times. These conditions are projected to meet the ITER steady-state goal of Q = 5. The ITB is maintained at lower β p with a strong reverse shear, consistent with predictions that negative central shear can lower the β p threshold for the ITB. There are two observed confinement states in the high β p scenario: H-mode confinement state with a high edge pedestal, and an enhanced confinement state with a low pedestal and an ITB. It has been observed in a scan of external resonant magnetic perturbation amplitude that when there are no large type-I ELMs, there is no transition to enhanced confinement. This is consistent with the proposed mechanism for ITB formation being a type-I ELM. Quasilinear gyro-Landau fluid predictive modeling of ITER suggests that only a modest reverse shear is required to achieve the ITB formation necessary for Q=5 when electromagnetic physics including the kinetic ballooning mode (KBM) is incorporated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗
DOE OSTI2022

Opportunities and Challenges in Passive Thermal-Fluid and Energy Systems

Abstract This article focuses on passive systems that are used in energy and thermal-fluid applications. These passive systems do not have moving parts and are reliable and cost-effective. Fluid motion in these passive devices could be driven by capillary force, gravity, osmotic pressure, and/or concentration gradient. The fundamental mechanisms and limitations of transport phenomena for passive systems are highlighted, followed by their applications in heat pipes, fuel cells, thermal energy storage, and desalination systems. The capabilities of the passive systems are limited by the balance between the driving force and transport resistance. Based on the fundamental understanding of fluid flow and phase change in passive systems, this study proposes associated transport phenomena and quantitative criteria to determine the maximum heat transfer rate, the transport distance, and minimum pore size of wick structures (when relevant) in these passive devices. This article concludes with the discussion of challenges and future opportunities of passive systems.

Engineering↗
DOE OSTI2022

Behavior of Tidally Triggered Earthquakes Depends on Fluid Conditions

Seismicity at The Geysers geothermal field and in Oklahoma is heavily influenced by industrial activities related to energy production, though the mechanism in which earthquakes are induced or triggered is different. At The Geysers, much of the seismicity is linked to thermoelastic stresses caused by injecting cold water into hot rocks, whereas in Oklahoma, the seismicity is linked to a reduction of confining stress on faults due to increasing pore pressure resulting from wastewater injections. Here, in this study, we show that these contrasting conditions are also evident in tidally triggered earthquakes. At The Geysers, earthquakes preferentially occur during maximum extensional stress, which does not occur at the same time as maximum shear stress on optimally oriented faults in the regional stress field. In Oklahoma, earthquakes preferentially occur during maximum shear stress on optimally oriented faults rather than maximum extensional stress. The magnitude of tidal extensional stress is naturally much greater than tidal shear stress. However, in a fluid-saturated environment, pore pressure responds to changes in volume, which can counteract or reduce the effect of the applied stress. The difference in behavior at these two sites is indicative of the level of coupling between applied stress and pore pressure, corresponding to unsaturated conditions at The Geysers and high pore pressure in Oklahoma.

58 GEOSCIENCES↗
DOE OSTI2022

Effect of interfacial vibrational coupling on surface wettability and water transport

Here, we report that the atomic-scale vibrational coupling at the solid-fluid interface can substantially alter the interfacial properties such as wettability and fluid slip. The wettability of water droplets on substrates subjected to various vibrational frequencies is studied using molecular dynamics simulation. The contact angle increases (i.e., becomes more hydrophobic) when the oscillation frequency of the substrate matches the intermolecular bending frequency of liquid water. We investigate the underlying mechanism by examining the dynamics of water molecules at the interface and find that the temporal contact between the solid and fluid is shorter when the frequencies match, resulting in weak solid-fluid adsorption. We further report that the vibrational match at the interface reduces wall-fluid friction and enhances water transport through the nanopore. Our findings demonstrate the importance of the atomic-scale vibrational coupling at the solid-fluid interface on the physicochemical behavior of nanodevices and biological nanochannels.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor – dynamic fluid-body interaction modeling

To regulate primary coolant flow in the Advanced Test Reactor (ATR), a butterfly valve was installed between the primary coolant pumps and the reactor core. If the mechanical connection between the valve's disk and its shaft ever fails, the disk will rotate freely. Rapid disk rotation may induce pressure surges that could damage upstream pipes. In the present work, the rotational trajectory and pressure evolution during a disk free-rotation scenario were analyzed via the dynamic fluid-body interaction (DFBI) approach in STAR-CCM+, with the movement of a solid being driven by the forces and moment/torque imposed by its surrounding fluid. Assuming a large initial opening angle, the disk accelerates slowly, but swiftly passes the closed position. As a result of the sudden valve closure, a pressure surge occurs in the upstream pipes, exceeding their maximum allowable pressure. Furthermore, the disk does not stabilize at the closed position but continually oscillates due to the unsteady nature of the coolant flow. Because of the significant and continuous water hammer effect, a fix to the butterfly valve is being implemented to prevent rapid valve closure due to failure at the valve's disk/shaft mechanical connection.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A level-set immersed boundary method for reactive transport in complex topologies with moving interfaces

A simulation framework based on the level-set and the immersed boundary methods (LS-IBM) has been developed for reactive transport problems in porous media involving a moving solid-fluid interface. The interface movement due to surface reactions is tracked by the level-set method, while the immersed boundary method captures the momentum and mass transport at the interface. The proposed method is capable of accurately modeling transport near evolving boundaries in Cartesian grids. The framework formulation guarantees second order accuracy in space. Since the interface velocity is only defined at the moving boundary, an interface velocity propagation method is also proposed. The method can be applied to other moving interface problems of the “Stefan” type. Here, we validate the proposed LS-IBM both for flow and transport close to an immersed object with reactive boundaries as well as for crystal growth. Lastly, the proposed method provides a powerful tool to model more realistic problems involving moving reactive interfaces in complex domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Using defects as a ‘fossil record’ to help interpret complex processes during additive manufacturing: as applied to raster-scanned electron beam powder bed additively manufactured Ti–6Al–4V

Abstract Defects in parts produced by additive manufacturing, instead of simply being perceived as deleterious, can act as important sources of information associated with the complex physical processes that occur during materials deposition and subsequent thermal cycles. Indeed, they act as materials-state ‘fossil’ records of the dynamic AM process. The approach of using defects as epoch-like records of prior history has been developed while studying additively manufactured Ti–6Al–4V and has given new insights into processes that may otherwise remain either obscured or unquantified. Analogous to ‘epochs,’ the evolution of these defects often is characterized by physics that span across a temporal length scale. To demonstrate this approach, a broad range of analyses including optical and electron microscopy, X-ray computed tomography, energy-dispersive spectroscopy, and electron backscatter diffraction have been used to characterize a raster-scanned electron beam Ti–6Al–4V sample. These analysis techniques provide key characteristics of defects such as their morphology, location within the part, complex compositional fields interacting with the defects, and structures on the free surfaces of defects. Observed defects have been classified as banding, spherical porosity, and lack of fusion. Banding is directly related to preferential evaporation of Al, which has an influence on mechanical properties. Lack-of-fusion defects can be used to understand columnar grain growth, fluid flow of melt pools, humping, and spattering events. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2020

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗