Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nonlinear dynamics, avalanches, and noise for driven Wigner crystals

Here we consider the driven dynamics of Wigner crystals interacting with random disorder. Using numerical simulations, we find a rich variety of transport phenomena as a function of charge density, drive, and pinning strength. For weak pinning, the system forms a defect-free crystal that depins elastically. When the pinning is stronger, a pinned glass phase appears that depins into a filamentary flow state, transitions at higher drives into a disordered flow phase, and finally forms a moving smectic. Within the filamentary flow phase, the conduction curves can show switching dynamics as well as negative differential conductivity in which switching events cause some flow channels to be blocked. The velocity noise in the filamentary flow regime exhibits narrow band characteristics due to the one-dimensional nature of the motion, while the moving smectic has narrow band velocity noise with a washboard frequency. In the disordered flow state, the noise power reaches a peak value and the noise has a 1/$\mathcal{f}$ character. Our transport results are consistent with recent experimental transport studies in systems where Wigner crystal states are believed to be occurring. Below the conduction threshold, we find that avalanches with a power-law size distribution appear when there are sudden local rearrangements of charges from one pinned configuration to another.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

On steady alternate bars forced by a localized asymmetric drag distribution in erodible channels

Studying the effect of different in-stream fluvial turbines siting on river morphodynamics allowed us to witness the onset of a time-averaged, large-scale, alternate distortion of bed elevations, which could not be exclusively related to the turbine rotor blockage. The longitudinal profiles of this two-dimensional bathymetric perturbation resemble those of steady fluvial bars. In this contribution we generalize the problem addressing a spatially impulsive, asymmetric distribution of drag force in the channel cross-section. Furthermore, this is experimentally investigated through the deployment of differently sized grids perpendicular to the flow, and analytically explored as a finite perturbation of an open channel flow over an erodible sediment layer, as described by a coupled flow–sediment shallow water equation. The steady solutions of this fluvial morphodynamic problem, physically represented by alternate bars scaling with the channel width, highlight the importance of the resonant conditions in defining the spatial extent of the bed deformation. The equations further suggest that in very shallow flows any asymmetric obstruction may lead to an upstream propagation of the steady bars, consistent with previous studies on the effects of channel curvature. In broad terms, this study provides the preliminary framework to control the onset of river meandering through imposed finite perturbations of the cross-section. In a more applied sense, it provides a tool to predict non-local scour–deposition patterns associated with the deployment of energy converters or other flow obstructions.

42 ENGINEERING↗
DOE OSTI2023

Hydrodynamic spin-orbit coupling in asynchronous optically driven micro-rotors

Abstract Vortical flows of rotating particles describe interactions ranging from molecular machines to atmospheric dynamics. Yet to date, direct observation of the hydrodynamic coupling between artificial micro-rotors has been restricted by the details of the chosen drive, either through synchronization (using external magnetic fields) or confinement (using optical tweezers). Here we present a new active system that illuminates the interplay of rotation and translation in free rotors. We develop a non-tweezing circularly polarized beam that simultaneously rotates hundreds of silica-coated birefringent colloids. The particles rotate asynchronously in the optical torque field while freely diffusing in the plane. We observe that neighboring particles orbit each other with an angular velocity that depends on their spins. We derive an analytical model in the Stokes limit for pairs of spheres that quantitatively explains the observed dynamics. We then find that the geometrical nature of the low Reynolds fluid flow results in a universal hydrodynamic spin-orbit coupling. Our findings are of significance for the understanding and development of far-from-equilibrium materials.

47 OTHER INSTRUMENTATION↗
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 OSTI2021

Temperature Homogeneity under Selective and Localized Microwave Heating in Structured Flow Reactors

Selective heating of different phases of multiphase systems via microwaves can result in energy savings and suppression of side reactions. However, materials properties and operating conditions that maximize temperature gradients are poorly understood. Here we utilize computational fluid dynamics (CFD) computations and temperature measurements in structured flow reactors (monoliths) in a monomodal microwave cavity to assess the temperature difference between the walls and the fluid and develop a simple lumped model to estimate when temperature gradients exist. We also explore the material’s thermal and electrical properties of structured reactors for isothermal catalyst conditions. We propose that CFD simulations can be used as a nonintrusive, predictive tool of temperature homogeneity. Importantly, we demonstrate that localized heating in the bed under several conditions rather than selective heating is responsible for the selectivity enhancement. Furthermore, our results indicate that structured beds made of high thermal conductivity materials avoid arcing and enable temperature homogeneity and low electrical conductivity materials allow microwaves to penetrate the domain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Mesh and model requirements for capturing deep-stall aerodynamics in low-Mach-number flows

Here, this work presents a comprehensive computational fluid dynamics investigation of the effects of grid resolution and turbulence-model choice for capturing the unsteady three-dimensional aerodynamic performance of NACA 0012 and 0021 airfoils, with specific focus on the deep-stall regime. At high angles of attack (α), wind turbine blades routinely experience vortex-induced vibrations, which can cause significant structural damages. Accurate predictions of post-stall aerodynamics can identify the frequencies at which such vibrations maybe triggered. In this context, the NACA 0012 airfoil simulations are conducted at a chord-based Reynolds number, Re c =2×10 6 , with the k-ω Shear-Stress Transport Reynolds-Averaged Navier-Stokes (RANS) and Improved Delayed Detached Eddy Simulation (IDDES) hybrid RANS-Large Eddy Simulation turbulence models. The effect of mesh resolution both in the wall-normal and spanwise directions is investigated. Only the IDDES model with a minimum spanwise resolution of 24 cells per chord length correctly predicts the aerodynamic forces. Spectral analysis shows the peak primary shedding frequency at α=30°, which signifies the end of the stall region. In the post-stall regime, both lift and drag frequencies drop asymptotically with increasing α. The Strouhal number, based on normalised chord length, remains nearly constant in this region. Based on this study, NACA 0021 airfoil runs are performed with IDDES for Re c =2.7×10 5 and 2.0×10 6 on the finest wall-normal mesh and three spanwise grids. Simulations conducted on the finer spanwise grids demonstrate grid independence and show good agreement with experiments. The effect of varying Rec on the airfoil frequency statistics is investigated. Additionally, comparison studies are presented to investigate the impact of airfoil thickness on the frequency content at Re c =2.0×10 6 . The results from the study provide guidance on the choice of mesh resolution with the IDDES model to accurately capture aerodynamic quantities for complex industrial applications.

17 WIND ENERGY↗
DOE OSTI2026

Assessment and Improvement Proposals for High-Pressure Rinsing Processes during Cleanroom Preparation of Superconducting Radio-Frequency Cavities for Particle Accelerators

The goal of this internship is to develop a fluid dynamics characterization of the high-pressure flow used for cleaning the internal surface of the SRF cavities. Initially, the High Pressure Rinsing (HPR) process is introduced, focusing on the problems of the SRF cavities and how this method can be used to increase their acceleration efficiency and performance. The characterization of the problem is then presented, including an introduction to the multiphase flow modeling and the analysis of the Fermilab facility for the HPR. To accomplish the goal, a preliminary estimation of the variables was performed using MATLAB. Then, multiple 2D simulations were developed with Ansys Fluent, and the influence of the main parameter was investigated. The simulation results are then presented and analysed, with a focus on the possible future developments of the HPR process fluid dynamic characterization.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Pseudo-viscous modeling of transport in dense granular flows for thermal energy storage applications

Dense, granular flows were examined to effectively capture and model bulk viscous properties in thin packed beds. A modified Couette cell with particle image velocimetry was used to experimentally determine pseudo-viscosity properties of four particulate media with varying morphologies: (1) iron oxide-coated SiO 2 particles, (2) CARBOBEAD CP30-60 particles, (3) CARBOBEAD CP40-100 particles, and (4) Al 2 O 3 beads. The pseudo-viscosity functions were fitted using a power law to correlate the measured shear stress as a function of measured shear rate. The pseudo-viscous functions were used as inputs to computation fluid dynamics models for a single-phase viscous fluid to predict granular flow profiles. Steady-state free surface velocity profiles at angular velocities <7 rad/s predicted by the model were in good agreement with the experimental particle image velocimetry measurements, resulting in Pearson correlation coefficients of 0.97 for iron-oxide coated SiO 2 particles and 0.95 for CP30-60 particles. As a result, this alternative approach to measuring pseudo-viscous properties under shearing and modeling bulk transport behavior of granular flow using computation fluid dynamics model offered significant reduction in computational load compared to discrete element methods.

14 SOLAR ENERGY↗
DOE OSTI2020

Molecular dynamics and continuum analyses of the electrokinetic zeta potential in nanostructured slit channels

This work presents a theoretical and numerical study of electrokinetic flow and the zeta potential for the case of slit channels with nanoscale surface features of dimensions comparable to the Debye length, by employing molecular dynamics simulations and continuum-level analyses. Additionally, a simple analytical model for considering the average effect of such surface features producing nanoscale roughness is proposed by employing matched asymptotic solutions for the charge density and fluid flow field, and matching conditions that satisfy electroneutrality and the Onsager reciprocal relation between the electroosmotic flow rate and streaming current. The proposed analytical model quantitatively accounts for results from molecular dynamics simulations that consider the presence of ion solvation shells and surface hydration layers. Our analysis indicates that a simultaneous knowledge of the electroosmotic and pressure-driven flow rate or streaming current can be instrumental to unambiguously determine the zeta potential in the presence of nanoscale surface features.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Hourly dynamic line ratings for existing transmission across the contiguous United States (preliminary results) [Slides]

The thermal rating of a transmission line depends on the rate at which heat generated by current flow through the line can be dissipated into the environment by convective and radiative cooling, keeping the conductor below a maximum operating temperature. Static or seasonal line ratings assume nearly-worst-case weather conditions (high ambient air temperature, high solar irradiance, and low wind speed), resulting in a low risk of "overrating" a line (i.e., allowing more current to flow than the line is designed for) but a high fraction of time when the line is "underrated" (i.e., able to flow more current while staying below the maximum temperature). Dynamic transmission line ratings (DLR) use real-time measurements of weather conditions, line sag, or line temperature to determine the thermal rating of transmission lines at hourly or higher time resolution, allowing for increased ratings during some (cold, windy, and/or dark) periods and correcting for overrating during other (hot, calm, and/or sunny) periods. Many studies have explored the benefits of DLR, but typically for small collections of transmission lines under a limited sample of weather conditions, leaving open the question of how widely applicable the benefits of DLR might be over diverse geographies and weather conditions. Here, we combine geospatial data on tens of thousands of transmission line routes across the contiguous United States with high-resolution hourly historical weather data for 2007-2013 to estimate the impact of DLR and other weather-adjusted rating systems. We compare static and dynamic line ratings to intermediate ratings adjusted by ambient temperature, day/night irradiance, clear-sky irradiance, and/or measured irradiance, illustrating the benefits of incorporating progressively more detailed weather data into the calculated line ratings. Moving from nighttime-adjusted ratings to clear-sky-adjusted ratings provides more benefit than moving from clear-sky-adjusted ratings to measured-irradiance-adjusted ratings. Depending on the assumed static windspeed for irradiance-adjusted-ratings, DLR either corrects for overrated hours or provides substantially increased ratings in the majority of hours on the majority of lines.

14 SOLAR ENERGY↗
DOE OSTI2023

Control methods for mitigating flow oscillations in a supercritical CO 2 recompression closed Brayton cycle

A dynamic model of a 10 MWe supercritical CO 2 (sCO 2 ) recompression closed Brayton power generation cycle is used to investigate control methods for mitigating oscillations in flow conditions and power demand load, especially as the load ramp rate is increased and the load setpoint is more closely tracked. The focus is on control of the main compressor inlet temperature (MCIT) since oscillations are prompted by the strong nonlinearities in sCO 2 fluid properties near the critical point and by changes in system inventory for efficient load management. In this study, the oscillatory behavior during load ramping is first shown and then two potential control solutions for substantially reducing oscillations are given: 1) the use of the main compressor inlet guide vanes (IGV) and 2) bypassing a portion of the sCO 2 around the cooler to maintain the MCIT. These control methods reduce the effect of cycle feedback from interactions between the MCIT and load. When implementing either control method, a load setpoint ramp of 7.5%/min was maintainable while achieving closer load setpoint tracking with limited oscillations. Finally, the impact of sensor noise on the control was examined since derivative action was used in the IGV and cooler sCO 2 bypass control. Finally, while some effect was evident, signal noise was not problematic to the control.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗
DOE OSTI2023

Precision Tests of the Nonlinear Mode Coupling of Anisotropic Flow via High-Energy Collisions of Isobars

Valuable information on dynamics of expanding fluids can be inferred from the response of such systems to perturbations in their initial geometry. We apply this technique in high-energy 96 Ru+ 96 Ru and 96 Zr+ 96 Zr collisions to scrutinize the expansion dynamics of the quark-gluon plasma, where the initial geometry perturbations are sourced by the differences in deformations and radial profiles between 96 Ru and 96 Zr, and the collective response is captured by the change in anisotropic flow V n between the two collision systems. Using a transport model, we analyze how the nonlinear coupling between lower-order flow harmonics V 2 and V 3 to the higher-order flow harmonics V 4 and V 5 , expected to scale as V 4 NL = χ 4 V 2 2 and V 5NL = χ 5 V 2 V 3 , gets modified as one moves from 96 Ru+ 96 Ru to 96 Zr+ 96 Zr systems. We find that these scaling relations are valid to high precision: variations of order 20% in V 4NL and V 5NL due to differences in quadrupole deformation, octupole deformation, and nuclear skin modify χ 4 and χ 5 by about 1–2%. Percent-level deviations are however larger than the expected experimental uncertainties and could be measured. Therefore, collisions of isobars with different nuclear structures are a unique tool to isolate subtle nonlinear effects in the expansion of the quark-gluon plasma that would be otherwise impossible to access in a single collision system.

Physics↗
DOE OSTI2022

Modeling atomically mixed graded density impactors

Graded density impactors (GDIs) are multi-material composite impactors used in gas gun experiments to tailor the drive conditions imparted to a sample test material. Previous graded density impactors generally rely on thin, but discrete, layers of different materials. The thinner and the greater number of layers will result in smoother compression. Taken to the limit of very thin layers would be pure material 1 at one surface, such as the front surface of an impactor, smoothly transitioning at the atomic scale to pure material 2 on the back surface. Such an impactor can initially shock, then smoothly compresses a material during a dynamic experiment. This type of experiment can serve to explore a larger region of thermodynamic space than a single or even multi-shock experiments. An overview of how graded density impactors are made is reviewed and sample results are given. Additionally, a strategy for modeling these kinds of impactors is presented. The length scales of constituent mixing are given from the experimental build through electrochemical-deposition. Equation of state models for pure constituents and their subsequent mixtures are presented. It is demonstrated that the time scales for pressure and temperature equilibration, for atomically mixed GDIs, are short enough to be a justifiable closure for the resulting multiphase flow. Furthermore, we present simulation results of dynamic shock followed by a ramp compression, utilizing a silver/gold graded density impactor, onto a tantalum sample.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
DOE OSTI2022

In situ feature analysis for large-scale multiphase flow simulations

The study of multiphase flow is essential for designing chemical reactors such as fluidized bed reactors (FBR), as a detailed understanding of hydrodynamics is critical for optimizing reactor performance and stability. An FBR allows scientists to conduct different types of chemical reactions involving multiphase materials, especially interaction between gas and solids. During such complex chemical processes, the formation of void regions in the reactor, generally termed as bubbles, is an important phenomenon. The study of these bubbles has a deep implication in predicting the reactor’s overall efficiency. But physical experiments needed to understand bubble dynamics are costly and non-trivial due to the technical difficulties involved and harsh working conditions of the reactors. Therefore, to study such chemical processes and bubble dynamics, a state-of-the-art computational simulation MFIX-Exa is being developed. Despite the proven accuracy of MFIX-Exa in modeling bubbling phenomena, the large-scale output data prohibits the use of traditional post hoc analysis capabilities in both storage and I/O time. Herein, to address these issues and allow the application scientists to explore the bubble dynamics in an efficient and timely manner, we have developed an end-to-end analytics pipeline that enables in situ detection of bubbles, followed by a flexible post hoc visual exploration methodology of bubble dynamics. The proposed method enables interactive analysis of bubbles, along with quantification of several bubble characteristics, enabling experts to understand the bubble interactions in detail. Positive feedback from the experts has indicated the efficacy of the proposed approach for exploring bubble dynamics in very-large-scale multiphase flow simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

In-situ L-TEM observations of dynamics of nanometric skyrmions and antiskyrmions

Nanometer-scale magnetic skyrmions and antiskyrmions exhibit unique dynamical behaviors in response to external stimuli, which are critical for their applications in low-power-consumption spintronic devices. This review discusses recent advancements in in-situ Lorentz transmission electron microscopy (L-TEM) observations of skyrmion and antiskyrmion dynamics, and demonstrates the manipulation and evolution of these textures in various magnetic materials under electric, magnetic, and thermal stimuli. Specifically, the motion tracking of single skyrmions and their clusters, and the deformation and transformation of skyrmions has been demonstrated in chiral helimagnets FeGe, Co 9 Zn 9 Mn 2 , and Co 10 Zn 10 with precise application of electric currents. Skyrmions can undergo dynamic transitions in current-driven skyrmion motions, from pinned states to linear flows, and even exhibit deformation into elliptical shapes, underscoring their topological robustness and dynamic flexibility. In addition, the manipulation of single antiskyrmions and antiskyrmion-lattice phases in (Fe 0.63 Ni 0.3 Pd 0.07 ) 3 P with S 4 symmetry is discussed, highlighting their high mobility and unique sliding capabilities along stripe domains at room temperature, facilitated by nanosecond pulsed electric currents. Finally, the temperature gradient-driven motion and topological transformation of elliptical skyrmions and antiskyrmions in this same material are investigated. In conclusion, the comprehensive insights gained from the L-TEM imaging technique are pivotal in advancing the design and functionality of next-generation skyrmion/antiskyrmion-based spintronic devices.

(Anti)skyrmion↗
DOE OSTI2023

Highly-scalable GPU-accelerated compressible reacting flow solver for modeling high-speed flows

Emerging supercomputing systems utilize a combination of central processing units (CPUs) and graphics processing units (GPUs) in an effort to reach exascale capabilities while minimizing the energy footprint of operating such systems. Such heterogeneous machines introduce new challenges for fluids solvers because the hardware architecture and operation of a GPU are fundamentally different from conventional CPUs. In this work, a general approach for efficient implementation of finite-volume based reacting flow solvers on such heterogeneous systems is presented. Three main challenges, namely, data access pattern, thread divergence, and thread safety, are addressed. Since compressible reacting flows require special methods to deal with chemical reactions, hyperbolic and nonlinear convection terms, and the presence of turbulence, specific algorithms that ensure GPU-based efficiency are developed. The approach is demonstrated on the widely available OpenFOAM open source software by modifying core algorithms for GPU accessibility. The scalability of the resulting solver, is demonstrated using practical test cases, including flow through a scramjet engine and the dynamics of a rotating detonation engine. Here, the solver provides near-ideal scaleup on a large number of GPUs (>3000), and extremely efficient use of the GPUs, with throughput nearly a constant even when processing a large number of control volumes.

42 ENGINEERING↗