Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortex 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.

1,500 records · Page 84

In Situ Chemical Monitoring and Imaging of Contents within Microfluidic Devices Having a Porous Membrane Wall Using Liquid Microjunction Surface Sampling Probe Mass Spectrometry

The ability to observe dynamic chemical processes (e.g., signaling, transport, etc.) in vivo or in situ using nondestructive chemical imaging opens a new door to understanding the complex dynamics of developing biological systems. With the advent of “biology-on-a-chip” devices has come the ability to monitor dynamic chemical processes in a controlled environment, using these engineered habitats to capture key features of natural systems while allowing visual observation of system development. Having the capability to spatially and temporally map the chemical signals within these devices may yield new insights into the forces that drive biosystem development. Here, a porous membrane sealed microfluidic device was designed to allow normal microfluidic operation while enabling continuous, location specific sampling and chemical characterization by liquid microjunction surface sampling probe mass spectrometry (LMJ-SSP MS). LMJ-SSP was used to extract fluids with nL-to-μL/min flow rates directly from selected areas of the microfluidic device without negatively impacting the device function. These extracts were subsequently characterized using MS. This technique was used to acquire MS images of the entirety of several multi-input microfluidic devices having different degrees of fluid mixing. LMJ-SSP MS imaging visualized the spatial distribution of chemical components within the microfluidic channels and could visualize chemical reactions occurring in the device. These microfluidic devices with a porous membrane wall are wholly compatible with the construction of biology-on-a-chip devices. This ultimately would enable correlation of biosystem physical structure with an evolving chemical environment

47 OTHER INSTRUMENTATION↗

A non-linear k-{epsilon} turbulent model sensitivity analysis for flow across a PWR spacer grid with mixing vanes

In this work we studied the sensitivity of different computational fluid dynamics (CFD) numerical solutions on turbulence model, domain discretization mesh type, and inlet boundary conditions. The case under examination is a turbulent one-phase flow (water) through a scaled nuclear pressurized water reactor (PWR) fuel assembly. The main feature of this geometrical configuration is the presence of a spacer grid with mixing vanes. The design of a spacer grid with mixing vanes plays a key role for the core heat transfer efficiency, increasing the coolant heat transfer coefficient. The main reason for this is the introduction of the top mixing vanes, which consist of an ordered repetitive succession of wings that generate secondary swirling flows. The increase of heat transfer efficiency also increases the critical heat flux (CHF), making it possible to increase the design and operational limits, far from boiling crises (departure of nucleate boiling in case of a PWR). The main focus of this work was to analyze the simulation capabilities of Reynolds-averaged Navier-Stokes (RANS) models (a cubic non-linear k-epsilon model formulation) for this complex flow geometry. The CFD code used was CD-Adapco STAR-CCM+.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High-order methods for wind turbine simulations and moving to next-generation platforms

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. This milestone represents an effort to increase the fidelity of Nalu-Wind at a fixed resolution through the implementation of a tensor-product based, matrix-free high order scheme. High order finite element methods have increased local work per datum communicated and have the potential to provide significantly more accurate solutions at a fixed number of degrees of freedom. Previous to this milestone, Nalu-Wind had an arbitrary order Control Volume Finite Element Method discretization as a solver option, but it required too much memory and was too slow to be of practical use. The work in this milestone addresses these issues by first implementing an implicit, high order solver that only partially assembles the global system. This reduces the memory footprint of the high-order scheme by orders of magnitude for higher polynomial orders. Second, a faster, tensor-product based method for evaluating the action of the left-hand side was implemented. This reduces the amount of computational work required by the scheme and dramatically enhanced the time-to-solution on example problems. Finally, this milestone is an evaluation of the value of high order methods in the wind application space. With the enhancements to memory and computational cost, accuracy vs. time-to-solution was evaluated for several resolutions on an under-resolved Taylor Green vortex test case. Results show that the high order scheme is cost-competitive with the production low-order schemes in Nalu-Wind, being moderately more expensive than the production edge-based vertex centered finite volume scheme. The evaluation of accuracy on the test case shows a potential benefit to high order at the highest resolution while not deteriorating accuracy on the lowest tested resolution. More work is needed to show value in the wind application, but positive strides have been made.

17 WIND ENERGY↗

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↗

Large Eddy Simulation of turbulent flow in rod bundle with spacer grids

Pressurized water nuclear reactor fuel assemblies are made of fuel rods held by spacer grids, which are the core components of the nuclear power plants. Hot fuel rods are cooled by the coolant water flow, and the spacer grids with mixing vanes distort the mainstream, producing the swirling flow which can enhance the heat transfer to improve CHF (Critical Heat Flux). Heat transfer between fuel rods and coolant is a crucial part of thermal-hydraulic analysis. For the purpose to improve the thermal hydraulic performance of the fuel assemblies, many experiments have been conducted to understand the flow structures and three-dimensional behaviors of the turbulence in the rod bundle sub-channels and then to optimize the design of the spacer grid. In the experimental work, which requires lots of efforts to get reliable data, LDA (Laser Doppler Anemometry) and PIV (Particle Image Velocimetry) are generally used to measure the velocity profile and turbulent intensity in the sub-channels. Computational Fluid Dynamics (CFD) is another useful way to study the characteristics of the detailed flow condition and heat transfer phenomena in the subchannels. The main issue of the single-phase CFD is to accurately predict the flow and turbulence characteristics. RANS method is generally used to simulate the induced lateral velocity and the heat transfer process. With the development of computer technology, some more precise but previously time-consuming turbulent models, such as LES (Large Eddy Simulation), are gradually used for the simulation of rod bundle flow with current advanced CPU resources. Ikeno and Kajishima explained the decay of swirling turbulent flow in the rod bundle using LES method. Benhamadouche acquired the turbulence induced forces on the central rod of the rod bundle through performing the LES of the flow downstream the spacer grid. Chang and Tavoularis conducted hybrid simulations with SAS SST and LES of a 5 x 5 rod bundle, and several vortices were identified in the vicinity of the vanes and the generation mechanisms, evolution progress and influence on the velocities of these vortices were clarified. In the paper we present CFD simulations of turbulent flow in rod bundle sub-channels with spacer grids of two different experiments. One is MATiS-H benchmark experiment with scaled spacer grids, which is conducted by KAERI, and is used as a CFD benchmark test, and the other is a 5 x 5 rod bundle with typical spacer grids using PIV technology, with the cooperation of Harbin Engineering University. The main objectives of the calculations are to validate LES methodology by the benchmark experiments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Kelvin–Helmholtz instability under stabilizing parallel magnetic field in nonhomogeneous compressible MHD flows

We study the Kelvin–Helmholtz instability (KHI) for the general case of a compressible, nonhomogeneous, magnetized plasma flow. The study is limited to a vortex sheet interface with an imposed parallel magnetic field. We introduce a new formalism based on a convective Mach number M c , a convective Alfvénic Mach number M Ac , and a total convective Mach number that combines the two. We derive an analytic expression of the KHI growth rate for a homogeneous flow (i.e., zero Atwood number, A=0) that converges toward both the expression for unmagnetized compressible flow and Chandrasekhar's expression for magnetized incompressible flow. Otherwise, the dispersion relation is solved numerically and allows deriving general stability diagrams of magnetized KHI for the triplet (A, M c , β −plasma) parameters. We show these parameters uniquely define all configurations for a parallel magnetic field. We also construct diagrams with respect to the convective Alfvénic Mach number, the β − plasma parameter, or the magnetic field showing which magnetic field strength is required for stabilizing a given shear flow. The theoretical growth rates are compared with 18 simulations made with the GAMERA code, currently used for 3D magnetospheric simulations. Finally, we apply our results to the analysis of a past KHI experiment performed at the OMEGA laser facility, showing linear theory succeeds to provide accurate estimates of the growth rate at early times. We further discuss how our results can inform future experiments in the high-Mach magnetized regime at the National Ignition Facility. Possible limitations of the study due to resistive, mixing, or turbulence effects are discussed.

compressible flows↗