Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW VELOCITY - DISTRIBUTION”

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 685 records · Page 38

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗

Turbulence generation by shock interaction with a highly nonuniform medium

An initially planar shock wave propagating into a medium of nonuniform density will be perturbed, leading to the generation of postshock velocity perturbations. Using numerical simulations we study this phenomenon in the case of highly nonuniform density (order-unity normalized variance, $\sigma_\rho/\bar{\rho}\sim1$) and strong shocks (shock Mach numbers $\bar{M}_s\gtrsim10$). This leads to a highly disrupted shock and a turbulent postshock flow. We simulate this interaction for a range of shock drives and initial density configurations meant to mimic those which might be presently achieved in experiments. Theoretical considerations lead to scaling relations, which are found to reasonably predict the postshock turbulence properties. The turbulent velocity dispersion and turbulent Mach number are found to depend on the preshock density dispersion and shock speed in a manner consistent with the linear Richtymer-Meshkov instability prediction. We also show a dependence of the turbulence generation on the scale of density perturbations. The postshock pressure and density, which can be substantially reduced relative to the unperturbed case, are found to be reasonably predicted by a simplified analysis that treats the extended shock transition region as a single normal shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterizing skyrmion flow phases with principal component analysis

Principal component analysis (PCA) is a powerful method that can identify patterns in large, complex data sets by constructing low-dimensional order parameters from higher-dimensional feature vectors. There are increasing efforts to use space-and-time-dependent PCA to detect transitions in nonequilibrium systems that are difficult to characterize with equilibrium methods. Here, we demonstrate that feature vectors incorporating the position and velocity information of driven skyrmions moving through random disorder permit PCA to resolve different types of disordered skyrmion motion as a function of driving force and the ratio of the Magnus force to the dissipation. Since the Magnus force creates gyroscopic motion and a finite Hall angle, skyrmions can exhibit a greater range of flow phases than what is observed in overdamped driven systems with quenched disorder. We show that in addition to identifying previously known skyrmion flow phases, PCA detects several additional phases, including different types of channel flow, moving fluids, and partially ordered states. Guided by the PCA analysis, we further characterize the disordered flow phases to elucidate the different microscopic dynamics and show that the changes in the PCA-derived order parameters can be connected to features in bulk transport measures, including the transverse and longitudinal velocity-force curves, differential conductivity, topological defect density, and changes in the skyrmion Hall angle as a function of drive. We discuss how asymmetric feature vectors can be used to improve the resolution of the PCA analysis, and how this technique can be extended to find disordered phases in other nonequilibrium systems with time-dependent dynamics.

36 MATERIALS SCIENCE↗

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Driven probe particle dynamics in a bubble and pattern forming system

We numerically examine the dynamics of a probe particle driven at a constant force through an assembly of particles with competing long-range repulsion and short-range attraction that forms a bubble or stripe state. In the bubble regime, we identify several distinct types of motion, including an elastic or pinned regime where the probe particle remains inside a bubble and drags all other bubbles with it. There is also a plastic bubble phase where the bubble in which the probe particle is trapped is able to move past the adjacent bubbles. At larger drives, there is a breakthrough regime where the probe particle jumps from bubble to bubble and, in some cases, can induce correlated rotations or plastic rearrangements of the particles within the bubbles. At the highest drives, the probe particle moves sufficiently rapidly that the background particles undergo only small distortions. The distinctive dynamic flow states and the transitions between them are accompanied by signatures in the effective drag on the driven particle, jumps in the velocity–force curves, and changes in the time-dependent velocity fluctuations. We map the dynamic phase diagram for this system for varied interaction lengths, bubble sizes, and densities.

36 MATERIALS SCIENCE↗

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗

Influence of conical wire array geometry on plasma flow and temperature profiles of radiatively cooled jets

The influence of conical wire array geometry on the formation and dynamics of pulsed-power driven plasma jets is investigated. In the experiments, the jet becomes isolated from the inflows as it passes through an aperture, allowing the study of its intrinsic evolution for different conical angles. Here, our results show that, regardless of the array opening angle, the jets are supersonic, highly collisional, and exhibit an exponential axial density decay with a characteristic scale length of 𝐿 𝑛 ≈ 3 mm, significantly shorter than the overall length of the jet. In contrast, axial velocity systematically increases with larger array opening angles. The near invariance of the density profile is attributed to a compensating mechanism between geometric divergence and axial acceleration, consistent with an asymptotically steady flow regime. Additionally, temperature measurements reveal ion-electron thermal decoupling near the base of the jet, with equilibration downstream. Analysis of radiative cooling and collisional energy exchange timescales indicates that the plasma evolves in a radiatively stable regime and that the two-temperature profile is fully described by the competition between both methods. These findings provide valuable information into the internal structure and evolution of pulsed-power plasma jets in regimes dominated by geometric shaping and radiative cooling, demonstrating experimental control over flow acceleration and collimation relevant to both laboratory and astrophysical contexts.

Physics - Plasma physics↗

Numerical simulation of frost formation and heat transfer on fin-and-tube heat exchangers in turbulent cross-flow

Frost formation in fin-and-tube heat exchangers in turbulent cross-flow presents significant challenges in industrial refrigeration applications, affecting heat transfer efficiency and operational reliability. The purpose of this work is to investigate frost deposition and growth on a staggered bank of a fin-and-tube freezer coil under turbulent forced convection conditions. The focus here is on investigating conditions that closely replicate real-world scenarios in large walk-in industrial freezers. Using a direct numerical simulation approach, we examine the flow dynamics and thermal behaviour in the presence of frost, considering turbulent regimes characterized by a Reynolds number in the range 1050 ≤ R e D , avg ≤ 4800 , with the characteristic length being the outer diameter of the tube and the velocity being the bulk fluid velocity between the plates (fins). Computational fluid dynamics simulations are employed to resolve the interactions between turbulent airflow and the frost layer. Our approach incorporates a modified immersed boundary method and a slow-time acceleration technique to address the complex dynamic interface between the continuously evolving frost layer and the flowing air stream. Our findings indicate that frost forms more on the sides of the finned surfaces (plates) and less on the tubes themselves. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗

A comprehensive review of nonlinear oscillators in hydrokinetic energy harnessing using flow-induced vibrations

A comprehensive review of hydrokinetic energy converters based on alternating lift technology (ALT) is provided. Emphasis is on nonlinear oscillators based on Flow Induced Vibration (FIV) or Oscillation (FIO). Due to strong coupling in Fluid-Structure Interaction (FSI), and in order to maximize the hydrokinetic harnessed energy, design of nonlinear oscillators and analysis by model tests or computational fluid dynamics dominates this area. Research confirmed that the nonlinear oscillator can harvest energy from a stochastic excitation modeled by a generic wide spectrum, and overcome the most severe oscillator limitations: specifically, the need for continuous frequency tuning due to the narrow bandwidth response, and low efficiency outside the narrow bandwidth oscillator response. This review covers the following aspects of nonlinear oscillators in ALT converters: (1) Geometric changes in oscillator cross-section; e.g., circular, square, rectangular, or trilateral shapes. (2) Passive turbulence control of FIV/FIO. (3) Position based nonlinear stiffness. (4) Multi-cylinder synergistic FIV/FIO. (5) Mechanically linked oscillators. (6) Velocity-based, nonlinear, adaptive harnessing damping.

42 ENGINEERING↗

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↗

Time-dependent probability density function analysis of H-mode transitions

Abstract The first application of time-dependent probability density function (PDF) analysis to the L-H transition in fusion plasmas is presented. PDFs are constructed using Doppler Backscattering data of perpendicular fluctuation velocity, , and turbulence from the edge region of the DIII-D tokamak. These raw time-series data are sliced into millisecond-long sliding time-windows to create PDFs. During the transition, the PDFs develop strong right tails, indicative of turbulence-suppressing localised flows in the plasma edge; such features and other subtle behaviours are explored using novel information geometry techniques. This letter examines the applicability of these techniques to predict L-H transitions and investigate predator-prey self-regulation theories between turbulence and perpendicular velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fabrication and Preliminary Demonstration of Microwave Resonant Cavity Transducer Performance

We are investigating a microwave cavity-based transducer for in-core high-temperature fluid flow sensing in molten salt cooled reactors (MSCR) and sodium fast reactors (SFR). This sensor is a hollow metallic cylindrical cavity, which can be fabricated from stainless steel, and as such is expected to be resilient to radiation, high temperature and corrosive environment of MSCR and SFR. The principle of sensing consists of making one wall of the cylindrical cavity flexible enough so that dynamic pressure, which is proportional to fluid velocity, will cause membrane deflection. Membrane deflection causes cavity volume change, which leads to a shift in the resonant frequency. We have developed an initial design for proof-of-principle testing of the flow sensor performance in microwave K-band. A cylindrical resonator prototype was fabricated from brass for initial tests in water. The cavity size is matched to the flange of a standard WR-42 waveguide. Microwave field is coupled into the resonant cavity through a subwavelength-size aperture. A test article was developed, consisting of a piping Tee with bulkhead WR-42 microwave waveguide installed in leak-proof design. In the test article, the cylindrical cavity is positioned in the center of the pipe. A microwave waveguide circulator was installed in the setup to suppress microwave reflections at the cavity entrance. Preliminary spectral characterization of cavity spectral response was performed with microwave VNA. Applying mechanical pressure to cavity membrane showed a measurable shift in the microwave resonant frequency. We also investigate mechanical integrity of the flowmeter’s membrane through computer simulations. By calculating the stress on the plate due to deflection and, comparing the stress to the material ultimate tensile strength and yield strength, it can be estimated if the plate will fail. The stress on the plate was calculated with an analytic closed form solution model, and with COMSOL Structural Mechanics Module which does not involve any approximations. Both the analytic model and COMSOL model showed that maximum stresses on the plate, which are at the radial boundary of the plate, are three orders of magnitude smaller than the yield strength and ultimate tensile strength. This indicates that the sensor is at a low risk of mechanical failure.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗

Colliding localized, lumpy holographic shocks with a granular nuclear structure

We apply a recent and simple technique which speeds up the calculation of localized collisions in holography to study more realistic models of the pre-hydrodynamic phase of heavy ion collisions using gauge/gravity duality. Our initial data reflects the lumpy nuclear structure of real heavy ions and our projectiles’ aspect ratio mimics the Lorentz contraction of nuclei during RHIC collisions. At the hydrodynamization time of the central region of the quark gluon plasma developed during the collision, we find that most of the system’s vorticity is located well outside the hydrodynamized part of the plasma. Only the relativistic corrections to the thermal vorticity within the hydrodynamized region are non-negligible. We compare the transverse flow shortly after the collision with previous results which did not use granular initial conditions and determine the proper energy density and fluid velocity in the hydrodynamized subregion of the plasma.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrokinetic energy conversion using flow induced oscillations of single-cylinder with large passive turbulence control

Various types of flow-induced oscillations (FIOs) have been implemented in development of marine hydrokinetic (MHK) energy converters. With passive turbulence control (PTC), energy harvesting starts at a flow speed of about 0.5 m/s. However, there is worldwide MHK energy available in even slower currents. In the present study, the effect of damping on FIO and power extraction is investigated for a converter with large turbulence stimulation (PTC) consisting of straight strips with a height of 15% of the cylinder diameter and placed symmetrically on the cylinder surface. The oscillating amplitude decreases, as the damping ratio increases, with unchanged sinusoidal pattern of the displacement time-history. The frequency ratio is also affected by damping especially in the VIV initial branch and transition region between VIV and galloping. An important flow characteristic of the large-PTC cylinder is that a recirculation region is formed behind the PTC, causing appreciable disturbance to the flow past the cylinder. Power can be harvested in the whole FIO range and the harnessed power maximum appears at the largest inflow velocity tested. However, the optimum of harnessing efficiency is located at the beginning of the VIV upper branch. The gap between VIV and galloping is bridged when large PTC is used, eliminating the drop in power and efficiency even at higher damping, which would be a weakness of regular-PTC cylinder for energy harvesting. Finally, the mechanism behind the variation of harnessing efficiency with inflow velocity and damping ratio is revealed, and the optimality criterion for the converter design is discussed.

16 TIDAL AND WAVE POWER↗

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗