Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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 1,135 records · Page 63

Mathematical Modeling of an Oscillating Droplet

Oscillating droplets are of interest in a number of disciplines. A practical application is the oscillating drop method, which is a technique for measuring surface tension and viscosity of liquid metals. It is especially suited to undercooled and highly reactive metals, because it is performed by electromagnetic levitation. The natural oscillation frequency of the droplets is related to the surface tension of the material, and the decay of oscillations is related to its viscosity. The fluid flow inside the droplet must be laminar in order for this technique to yield good results. Because no experimental method has yet been developed to visualize flow in electromagnetically-levitated oscillating metal droplets, mathematical modeling is required to determine whether or not turbulence occurs. Three mathematical models of the flow: (1) assuming laminar conditions, (2) using the k-epsilon turbulence model, and (3) using the RNG turbulence model, respectively, are compared and contrasted to determine the physical characteristics of the flow. It is concluded that the RNG model is the best suited for describing this problem. The goal of the presented work was to characterize internal flow in an oscillating droplet of liquid metal, and to verify the accuracy of the characterization by comparing calculated surface tension and viscosity.

Berry, S.↗

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↗

3D structure of anisotropic flow in small collision systems at energies available at the BNL Relativistic Heavy Ion Collider

Here, we present (3 + 1)-dimensional [(3 + 1)D] dynamical simulations of asymmetric nuclear collisions at the BNL Relativistic Heavy Ion Collider (RHIC). Employing a dynamical initial state model coupled to (3 + 1)D viscous relativistic hydrodynamics, we explore the rapidity dependence of anisotropic flow in the RHIC small system scan at 200 GeV center-of-mass energy. We calibrate parameters to describe central 3 He + Au collisions and make extrapolations to d + Au and p + Au collisions. Our calculations demonstrate that approximately 50% of the v 3 (p T ) difference between the measurements by the STAR and PHENIX Collaborations can be explained by the use of reference flow vectors from different rapidity regions. This emphasizes the importance of longitudinal flow decorrelation for anisotropic flow measurements in asymmetric nuclear collisions, and the need for (3 + 1)D simulations. We also present results for the beam energy dependence of particle spectra and anisotropic flow in d + Au collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Implications of Upwells as Hydrodynamic Jets in a Pulse Jet Mixed System

This report evaluates the physics of the upwell flow in pulse jet mixed systems in the Hanford Tank Waste Treatment and Immobilization Plant (WTP). Although the initial downward flow and radial flow from pulse jet mixers (PJMs) has been analyzed in some detail, the upwells have received considerably less attention despite having significant implications for vessel mixing. Do the upwells behave like jets? How do the upwells scale? When will the central upwell breakthrough? What proportion of the vessel is blended by the upwells themselves? Indeed, how the physics of the central upwell is affected by multiple PJMs (e.g., six in the proposed mixing vessels), non-Newtonian rheology, and significant multicomponent solids loadings remain unexplored. The central upwell must satisfy several criteria to be considered a free jet. First, it must travel for several diameters in a nearly constant direction. Second, its velocity must decay with the inverse of elevation. Third, it should have an approximately Gaussian profile. Fourth, the influence of surface or body forces must be negligible. A combination of historical data in a 12.75 ft test vessel, newly analyzed data from the 8 ft test vessel, and conservation of momentum arguments derived specifically for PJM operating conditions demonstrate that the central upwell satisfies these criteria where vigorous breakthrough is achieved. An essential feature of scaling from one vessel to the next is the requirement that the underlying physics does not change adversely. One may have confidence in scaling if (1) correlations and formulas capture the relevant physics; (2) the underlying physics does not change from the conditions under which it was developed to the conditions of interest; (3) all factors relevant to scaling have been incorporated, including flow, material, and geometric considerations; and (4) the uncertainty in the relationships is sufficiently narrow to meet required specifications. Although the central upwell satisfies these criteria when vigorous breakthrough is achieved, not all available data follow the free jet profile for the central upwell, particularly at lower nozzle velocities. Alternative flow regimes are considered and new models for cloud height, “cavern height,” and the rate of jet penetration (jet celerity) are benchmarked against data to anchor scaling analyses. This analytical modeling effort to provide a technical basis for scaling PJM mixed vessels has significant implications for vessel mixing, because jet physics underlies “cavern” height, cloud height, and the volume of mixing considerations. A new four-parameter cloud height model compares favorably to experimental results. This model is predictive of breakthrough in 8 ft vessel tests with the two-part simulant. Analysis of the upwell in the presence of yield stresses finds evidence of expanding turbulent jets, confined turbulent jets, and confined laminar flows. For each, the critical elevation at which jet momentum depletes is predicted, which compare favorably to experimental cavern height data. Partially coupled momentum and energy balances suggest that these are limiting cases of a gradual transition from a turbulent expanding flow to a confined laminar flow. This analysis of the central upwell alone lays essential groundwork for complete analysis of mode three mixing (i.e., breakthrough with slow peripheral mixing). Consideration of jet celerity shows that the rate of jet penetration is a governing consideration in breakthrough to the surface. Estimates of the volume of mixing are presented. This analysis shows that flow along the vessel wall is sluggish such that the central upwell governs the volume of mixing. This analysis of the central upwell alone lays essential groundwork for complete analysis of mode three mixing and estimates of hydrogen release rates from first principles.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Synchrotron X-ray based particle image velocimetry to measure multiphase streamflow and densitometry

Synchrotron-based radiation techniques are increasingly used in the field of flow condition measurements. As a complicated multiphase flow, cavitation has been investigated in a wide variety of industrial fields, to prevent damages on spillways at large dams or intensify the bacteria eradication rate during water treatments. To control the detrimental effects of cavitating flows, the dynamics of and mechanisms affecting cavitation development shall be visualized and identified in real-time. Due to the large reflection and scattering effects of multiphase flow, standard visible-light imaging technologies cannot penetrate vapor bubbles and provide extremely limited information. In this study, the synchrotron source emitted high-flux X-ray pulses are used to capture flow motion and visualize internal structures with sufficiently high spatial and temporal resolutions. The image processing procedures on raw x-ray images are developed using 2D Fourier transform (FT) and wavelet transform (WT) to implement the contrast enhancement and de-noising. Through the edge detection algorithm, tracking either seeded particles or phase interfaces inside the opaque multiphase flow can be achieved to perform the particle image velocimetry. Further, with decomposing particle-only and particle-eliminated subplots from the X-ray phase-contrast image, the time-resolved velocity and void fraction fields are obtained simultaneously, paving a way for further flow condition analysis such as densitometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Saturation scalings of toroidal ion temperature gradient turbulence

The emerging understanding of instability-driven plasma-turbulence saturation in terms of energy transfer to stable modes in the same scale range as the instability is employed to derive a saturation theory for the toroidal branch of ion temperature gradient turbulence that provides the scaling of turbulence and zonal flow levels for all physical parameters. The theory is based on the eigenmode decomposition of a nonlinear fluid model, which is subjected to a statistical closure and simplified via an ordering expansion consistent with zonal-flow catalyzed energy transfer from the unstable mode to the stable mode at large scale. Solution of the closed energy balance equations yields a turbulence level that is proportional to the ratio of the zonal flow damping rate and the inverse of the triplet correlation time of the zonal-flow catalyzed wavenumber triplet interaction. The zonal flow energy is proportional to the ratio of the growth rate and the inverse triplet correlation time. The saturation scalings are applied to the ion heat flux, showing that it has a factor proportional to the quasilinear heat flux and a correction factor that includes the inverse of the triplet correlation time and a reduction due to the stable mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic phases, clustering, and chain formation for driven disk systems in the presence of quenched disorder

In this paper, we numerically examine the dynamic phases and pattern formation of two-dimensional monodisperse repulsive disks driven over random quenched disorder. We show that there is a series of distinct dynamic regimes as a function of increasing drive, including a clogged or pile-up phase near depinning, a homogeneous disordered flow state, and a dynamically phase separated regime consisting of high-density crystalline regions surrounded by a low density of disordered disks. At the highest drives the disks arrange into one-dimensional moving chains. The phase separated regime has parallels with the phase separation observed in active matter systems, but arises from a distinct mechanism consisting of the combination of nonequilibrium fluctuations with density-dependent mobility. We discuss the pronounced differences between this system and previous studies of driven particles with longer-range repulsive interactions moving over random substrates, such as superconducting vortices or electron crystals, where dynamical phase separation and distinct one-dimensional moving chains are not observed. Finally, our results should be generic to a broad class of systems in which the particle-particle interactions are short ranged, such as sterically interacting colloids or Yukawa particles with strong screening driven over random pinning arrays, superconducting vortices in the limit of small penetration depths, or quasi-two-dimensional granular matter flowing over rough landscapes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessing the Role of Hydrodynamics in Enhancing Height-Above-the-Nearest-Drainage Derived Synthetic Rating Curves: A Comparative Study in the Wu River Basin, Taiwan

The conventional approach to generating synthetic rating curves (SRC) using the Height-Above-the-Nearest-Drainage (HAND) method typically relies on the assumption of uniform flow, such as Manning's equation, to establish stage-discharge ratings. The zero-physics application of the uniform flow equation is insufficient for capturing detailed hydraulic features (e.g., backwater effect) and neglects the hydraulic effects from adjacent channels. This lack of hydrodynamic computation can impact the accuracy and effectiveness of riverine flood risk estimation and management. To reduce this foreseeable error, we introduce the HAND-hd workflow, which integrates sophisticated hydrodynamic computations in the production of HAND-based SRC with hydrodynamic features (SRC hd ). The results indicate that SRC hd demonstrates consistent agreement with both gauge observations and benchmark solutions. Additionally, the comparative analysis suggests that SRC hd provides notable improvements in stage-discharge ratings over conventional HAND-based SRCs, particularly in channels with mild bed gradients, where it reduces water stage prediction errors and percent biases. In steeper channel segments, SRC hd maintains comparable accuracy to conventional methods. The comprehensive evaluation in this study emphasizes the potential discrepancies and inaccuracies associated with the adoption of the uniform flow assumption in the conventional HAND-SRCs and addresses the necessity of including hydrodynamic physics in the application of HAND-based SRC (e.g., inundation map) in channels with mild gradients.

54 ENVIRONMENTAL SCIENCES↗

Method for Comparing Porous Carbon Electrode Performance in Redox Flow Batteries

In this paper, we present a method for comparing porous electrodes for use in redox flow batteries. By comparing physical and electrochemical properties of six carbon cloth electrodes (CCEs) including hydraulic permeability, Ohmic resistance, and mass transport overpotential, we show how optimal engineering conditions such as compression and electrolyte flow rates can be derived. These results are combined and applied to a quinone-based redox flow battery, leading to a peak power density around 0.4 W cm –2 for several different CCEs. This method can be applied broadly to other electrode architectures and flow battery chemistries, and provides a comprehensive framework for their comparison.

25 ENERGY STORAGE↗

Current Sheet Magnetic Model for the Solar Corona

A new magnetic model is developed and compared with previous models and the observed solar corona. An attempt is made to more accurately compute the three-dimensional currents flowing in the solar corona. Physical reasons are given that require most of the large scale currents flowing in the solar corona to lie near thin sheets. The current sheets are not constrained into any particular geometry or symmetry as in the previous models of Altschuler and Newkirk (1969) and Schatten et al. (1969). A comparison with the axisymmetric, isothermal MHD solution of Pneuman and Kopp (1970) suggests that the model is able to simulate to high accuracy an isothermal corona. A comparison of the model with the May 30, 1965, solar eclipse and the November 12, 1966, solar eclipse shows the model is capable of computing many features including the polar plume orientations as well as radial and nonradial streamers in the solar corona.

Schatten, K. H.↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid. inverted annular. dispersed and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered.

Antar, B. N.↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum, and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid, inverted annular, dispersed, and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered. The model described above was solved numerically through a mixed finite difference scheme with forward time marching. The inverted annular regime was resolved using a semi-implicit finite differencing while the dispersed regime was solved explicitly. Also, a staggered mesh was used in which the velocity was resolved at mesh boundaries while all other field variables were resolved at the mesh centroids. Different mesh sizes were used depending on the region of solution. A coarse mesh was used in the dispersed flow region while a much finer mesh was used in both the inverted annular flow region and the tube walls.

B. N. Antar↗

Wettability control on multiphase flow in patterned microfluidics

Significance The simultaneous flow of multiple fluid phases through a porous solid occurs in many natural and industrial processes—for example, rainwater infiltrates into soil by displacing air, and carbon dioxide is stored in deep saline aquifers by displacing brine. It has been known for decades that wetting—the affinity of the solid to one of the fluids—can have a strong impact on the flow, but the microscale physics and macroscopic consequences remain poorly understood. Here, we study this in detail by systematically varying the wetting properties of a microfluidic porous medium. Our high-resolution images reveal the fundamental control of wetting on multiphase flow, elucidate the inherently 3D pore-scale mechanisms, and help explain the striking macroscopic displacement patterns that emerge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tetranuclear Polypyridylruthenium(II) Complexes as Selective Nucleic Acid Stains for Flow Cytometric Analysis of Monocytic and Epithelial Lung Carcinoma Large Extracellular Vesicles

Selective staining of extracellular vesicles (EVs) is a major challenge for diagnostic and therapeutic applications. Herein, the EV labeling properties of a new class of tetranuclear polypyridylruthenium(II) complexes, Rubb7-TNL and Rubb7-TL, as phosphorescent stains are described. These new stains have many advantages over standard stains to detect and characterize EVs, including: high specificity for EV staining versus cell staining; high phosphorescence yields; photostability; and a lack of leaching from EVs until incorporation with target cells. As an example of their utility, large EVs released from control (basal) or lipopolysaccharide (LPS)-stimulated THP-1 monocytic leukemia cells were studied as a model of immune system EVs released during bacterial infection. Key findings from EV staining combined with flow cytometry were as follows: (i) LPS-stimulated THP-1 cells generated significantly larger and more numerous large EVs, as compared with those from unstimulated cells; (ii) EVs retained native EV physical properties after staining; and (iii) the new stains selectively differentiated intact large EVs from artificial liposomes, which are models of cell membrane fragments or other lipid-containing debris, as well as distinguished two distinct subpopulations of monocytic EVs within the same experiment, as a result of biochemical differences between unstimulated and LPS-stimulated monocytes. Comparatively, the staining patterns of A549 epithelial lung carcinoma-derived EVs closely resembled those of THP-1 cell line-derived EVs, which highlighted similarities in their selective staining despite their distinct cellular origins. This is consistent with the hypothesis that these new phosphorescent stains target RNA within the EVs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A kinetic theory treatment of heat transfer in plane Poiseuille flow with uniform pressure

Plane compressible Poiseuille flow with uniform pressure (Couette flow with stationary boundaries) is revisited where the Lees two-steam method with the Enskog equation of change is applied. Single particle velocity distribution functions are chosen, which preserve the essential physical features of this flow with arbitrary but uniform plate temperatures and gas pressure. Lower moments are shown to lead to expressions for the parameter functions, molecular number densities, and temperatures which are entirely in agreement with those obtained in the analysis of Lees for compressible plane Couette flow in the limit of low Mach number and vanishing mean gas velocity. Important simplifications result, which are helpful in gaining insight into the power of kinetic theory in fluid mechanics. The temperature distribution, heat flux, as well as density, are completely determined for the whole range of Knudson numbers from free molecular flow to the continuum regime, when the pressure level is specified.

Bahrami, Parviz A.↗

An Assessment of Artificial Compressibility and Pressure Projection Methods for Incompressible Flow Simulations

Performance of the two commonly used numerical procedures, one based on artificial compressibility method and the other pressure projection method, are compared. These formulations are selected primarily because they are designed for three-dimensional applications. The computational procedures are compared by obtaining steady state solutions of a wake vortex and unsteady solutions of a curved duct flow. For steady computations, artificial compressibility was very efficient in terms of computing time and robustness. For an unsteady flow which requires small physical time step, pressure projection method was found to be computationally more efficient than an artificial compressibility method. This comparison is intended to give some basis for selecting a method or a flow solution code for large three-dimensional applications where computing resources become a critical issue.

Kwak, Dochan↗