SEARCH · Engineering Papers

Results for “oscillating flow”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗

DOE OSTI2023

FY23 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY23. As a continuation of the modeling effort from FY22, this year’s work focuses on improving the RELAP5-3D model developed previously for two-phase flow simulations. The RELAP5-3D model is updated to more accurately capture the heat loss experienced by the facility. The updated model is compared against experimental data for benchmarking purposes of the RELAP5-3D input model. By correctly accounting for heat loss, the updated RELAP5-3D model can now predict the two-phase baseline case more accurately. The onset and the duration of instability are captured well by the model. Furthermore, analyses are performed to better understand the instability mechanism experienced by the flow where the expansion of the boiling boundary in the chimney is studied in details and the fundamental frequencies of the oscillations are obtained. The updated RELAP5-3D model is further compared against four fault conditions, namely the reduction of riser header inlet flow area, depletion of system inventory, blocked riser channels, and static boiling scenario. For each fault condition, minor modifications and tuning are performed to improve the predictions of the model. The purpose of the analyses is to investigate the capability of RELAP5-3D in predicting complex two-phase flows in possible accident scenarios in actual Reactor Cavity Cooling System (RCCS). Overall, the model is able to capture the behaviors and trends of these fault conditions relatively well. Some discrepancies remain between the experimental data and the predictions, many of which are likely due to the differences in the predicted and experimental vapor generation rate. Future work will focus on the continued development of the current RELAP5-3D input model of the NSTF to both improve the accuracy of the model’s predictive capability and continue supporting the experimental program needs. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Neural-network learning of SPOD latent dynamics

Here, we aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower -dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Tool Chain for Simulations of Bi-Filar Coil Winding for Fast Quench Protection

The advancement of high-field magnets utilizing high-temperature superconductors (HTS) brings about complex challenges, especially in quench detection and protection. Traditional methods often fall short due to the inherently slow quench propagation in HTS materials. One promising approach to overcome this involves using a bifilar winding configuration, where two conductors are placed side by side. Under normal operation, they function in series, but during a quench event, they switch to an anti-parallel mode. This shift reduces the differential inductance of the coil to near zero, enabling rapid current oscillations through a capacitor discharge. The resulting high-frequency current flow leads to swift, uniform heating, triggering a full-coil quench within microseconds. Moreover, the strong mutual coupling between the two windings significantly reduces electrical noise in voltage measurements. In this work, we explore the viability of this concept by designing, constructing, and testing a REBCO bifilar racetrack coil in liquid nitrogen. We also present a validated simulation model that closely mirrors the coil's dynamic behavior under these conditions, aligning well with experimental observations.

Jayathilaka, Rehan [Northern Illinois U.]↗

DOE OSTI2024

Strain shift measured from stress-controlled oscillatory shear: Evidence for a continuous yielding transition and new techniques to determine recovery rheology measures

Understanding the yielding of complex fluids is an important rheological challenge that affects our ability to engineer and process materials for a wide variety of applications. Common theoretical understandings of yield stress fluids follow the Oldroyd–Prager formalism in which the material behavior below the yield stress is treated as solidlike, and above the yield stress as liquidlike, with an instantaneous transition between the two states. This formalism was built on a quasi-static approach to the yield stress, while most applications, ranging from material processing to end user applications, involve a transient approach to yielding over a finite timescale. Using stress-controlled oscillatory shear experiments, we show that yield stress fluids flow below their yield stresses. This is quantified through measuring the strain shift, which is the value about which the strain oscillates during a stress-controlled test and is a function of only the unrecoverable strain. Measurements of the strain shift are, therefore, measurements of flow having taken place. These experimental results are compared to the Herschel–Bulkley form of the Saramito model, which utilizes the Oldroyd–Prager formalism, and the recently published Kamani–Donley–Rogers (KDR) model, in which one constitutive equation represents the entire range of material responses. Scaling relationships are derived, which allow us to show why yield stress fluids will flow across all stresses, above and below their yield stress. Finally, derivations are presented that show strain shift can be used to determine average metrics previously attainable only through recovery rheology, and these are experimentally verified.

Mechanics↗

DOE OSTI2021

Experimental Investigation of Fuel Staging Effect on Modal Dynamics of Thermoacoustic Azimuthal Instabilities in a Multi-Nozzle Can Combustor

Abstract This paper analyzes the dynamics of unstable azimuthal thermoacoustic modes in a lean premixed combustor. Azimuthal modes can be decomposed into two counter rotating waves where they can either compete and potentially suppress one of them (spinning) or coexist (standing), depending on the operating conditions. This paper describes experimental results of the dynamical behaviors of these two waves. The experimental data were taken at different mass flow rates as well as different azimuthal fuel staging in a multi-nozzle can combustor. It is shown that at a low flow rate with uniform fuel distribution, the two waves have similar amplitudes, giving rise to a standing wave. However, the two amplitudes are slowly oscillating out of phase to each other, and the phase difference between the two waves also shows oscillatory behavior. For an intermediate flow rate, the dynamics show intermittency between standing and spinning waves, indicating that the system is bistable. In addition, the phase difference dramatically shifts when the mode switches between standing and spinning waves. For a high flow rate, the system stabilizes at a spinning wave most of the time. These experimental observations demonstrate that not only the amplitudes of two waves but also the phase difference plays an important role in the dynamics of azimuthal mode. For non-uniform azimuthal fuel staging, the modal dynamics exhibit only an oscillatory standing wave behavior regardless of the mass flow rate. Compared to the uniform fuel staging, however, the pressure magnitude is considerably reduced, which provides a potential strategy to mitigate and/or suppress the instabilities.

Kim, J.↗

DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Surface prediction and measurement for modulated tool path (MTP) turning

Here, this paper describes a time-domain simulation for predicting surface finish in modulated tool path (MTP) turning, which uses sinusoidal axis motions in the feed direction to produce discontinuous chips for ductile workpiece materials. The simulation includes: the low frequency and low amplitude tool oscillation in the feed direction; the time-varying chip thickness, cutting force, and tool displacement; and the plastic side flow effect used to calibrate the effective tool nose radius. Comparisons between predicted and measured surface profiles are presented for turning a 6061-T6 aluminum cylinder as a function of the MTP oscillation frequency and amplitude with discontinuous chip formation.

42 ENGINEERING↗

DOE OSTI2021

Coherent microwave scattering from resonance enhanced multi-photon ionization (Radar REMPI): A review

Coherent microwave scattering from laser-induced plasmas, including weakly ionized plasma, laser sparks, Multiphoton Ionization (MPI), and Resonance Enhanced Multi-Photon Ionizations (Radar REMPI) has achieved much successes in plasma, reactive and nonreactive flow diagnostics. Under illumination of microwaves (Radar), electrons inside the laser-induced plasma oscillate with the electric field of the microwave and re-radiate from the electrons forming coherent scattering. In the far-field approximation, the microwave scattering from the small volume plasma reflects the generation and evolution of unbounded electrons inside the plasma, when the microwave wavelength is much greater than the size of the plasma and the skin layer depth at the microwave frequency is larger than the size of the plasma. Laser excitation schemes, microwave detection methods, calibration of microwave scattering, and the novel applications of the technique have been significantly expanded and improved. Furthermore, this review paper summarizes physical principles, various REMPI excitation schemes for atomic and molecular species, and temperature measurements in plasma and reactive flows. Discussions on new research directions and applications are given at the end.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Effect of the divertor configuration on the JET edge radial electric field

The role of the divertor configuration and divertor plasma physics on the L–H transition is poorly understood, leading to large uncertainties in predicting the L–H power threshold in future devices. This contribution reports on edge perpendicular plasma flow measurements by Doppler backscattering in JET L–H transition experiments with the outer divertor strike-point at different positions: horizontal target (HT), vertical target (VT), and in the corner configuration (between the HT and VTs). The edge perpendicular flow was found to be significantly affected by changes in the divertor configuration in the region inside the separatrix. Our results do not show evidence for the existence of a critical edge flow shear needed to achieve H-mode for different divertor configurations, with a larger shear observed for the VT configuration. No significant change in the shear flow and in the density fluctuation level is measured preceding the L–H transition in the region just inside the separatrix. The dynamics of the L–H transitions are also influenced by the divertor configuration with divertor oscillations (DOs) observed only in the HT configuration. Interestingly, DOs are associated with marked changes in the edge perpendicular flow around the separatrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Convective Momentum Transport and Its Impact on the Madden‐Julian Oscillation in E3SM‐MMF

Abstract Convective momentum transport (CMT) is the process of vertical redistribution of horizontal momentum by small‐scale turbulent flows from moist convection. Traditional general circulation models (GCMs) and their multiscale modeling framework (MMF) versions poorly represent CMT due to insufficient information of subgrid‐scale flows at each GCM grid. Here the explicit scalar momentum transport (ESMT) scheme for representing CMT is implemented in the Energy Exascale Earth System Model‐Multiscale Modeling Framework (E3SM‐MMF) with embedded 2‐D cloud‐resolving models (CRMs), and verified against E3SM‐MMF simulations with 3‐D CRMs and observations. The results show that representing CMT by ESMT helps reduce climatological mean precipitation model bias over the western Pacific and the ITCZ regions, which is attributed to the weakened mean easterlies over the Pacific. Also, CMT from simulations with 2‐D and 3‐D CRMs impose a similar impact on Kelvin waves by reducing their variability and slowing down their phase speed, but opposite impacts on the Madden‐Julian Oscillation (MJO) variability. The ESMT scheme readily captures the climatological mean spatial patterns of the zonal and meridional components of CMT and their variability across multiple time scales, but shows some differences in estimating its magnitude. CMT mainly affects the MJO by decelerating its winds in the free troposphere, but accelerates its near‐surface winds. This study serves as a prototype for implementing CMT scheme in the MMF simulations, highlighting its crucial role in reducing model bias in mean state and spatiotemporal variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗

DOE OSTI2024

Mechanisms of Low-Level Jet Formation in the U.S. Mid-Atlantic Offshore

Abstract Low-level jets (LLJs), in which the wind speed attains a local maximum at low altitudes, have been found to occur in the U.S. mid-Atlantic offshore, a region of active wind energy deployment as of 2023. In contrast to widely studied regions such as the U.S. southern Great Plains and the California coastline, the mechanisms that underlie LLJs in the U.S. mid-Atlantic are poorly understood. This work analyzes floating lidar data from buoys deployed in the New York Bight to understand the characteristics and causes of coastal LLJs in the region. Application of the Hilbert–Huang transform, a frequency analysis technique, to LLJ case studies reveals that mid-Atlantic jets frequently occur during times of adjustment in synoptic-scale motions, such as large-scale temperature and pressure gradients or frontal passages, and that they do not coincide with motions at the native inertial oscillation frequency. Subsequent analysis with theoretical models of inertial oscillation and thermal winds further reveals that these jets can form in the stationary geostrophic wind profile from horizontal temperature gradients alone—in contrast to canonical LLJs, which arise from low-level inertial motions. Here, inertial oscillation can further modulate the intensity and altitude of the wind speed maximum. Statistical evidence indicates that these oscillations arise from stable stratification and the associated frictional decoupling due to warmer air flowing over a cold sea surface during the springtime land–sea breeze. These results improve our conceptual understanding of mid-Atlantic jets and may be used to better predict low-level wind speed maxima. Significance Statement The purpose of this work is to identify and characterize the atmospheric mechanisms that result in an occasional low-level maximum in the wind speed off the U.S. mid-Atlantic coastline. Our findings show that these low-level jets form due to horizontal temperature gradients arising from fronts and synoptic systems, as well as from the land–sea breeze that forces warmer air over the cold ocean surface. This work aids predictability of such jets, improves our understanding of this coastal environment, and has implications for future deployment of offshore wind energy in this region.

17 WIND ENERGY↗

DOE OSTI2023

Turbulence and E × B flow correlations across the L-H transition in DIII-D deuterium and hydrogen plasmas

The isotope dependence of the low- to high confinement-mode (L- to H-mode) transition power threshold P LH presents significant challenges for the initial (non-nuclear) hydrogen operations phase of the international thermonuclear experimental reactor (ITER). Here, we examine the isotope dependence of turbulence and E × B flow correlation properties in the L-mode edge plasma, leading up to the L-H transition. At marginal auxiliary power (near P LH ), turbulence is initially suppressed periodically during limit cycle oscillations (LCO) that precede the transition to sustained H-mode confinement. Here, we present evidence that the long-range (toroidal) correlation of the E × B edge plasma flow across the LCO phase is much weaker in hydrogen than in deuterium in the DIII-D tokamak, congruent with the higher threshold power P LH required to access LCO and H-mode in hydrogen. Concomitantly, the time required to initially quench edge turbulence via localized edge E × B flow shear is significantly longer in hydrogen (1–1.5 ms) than in deuterium (approx. 100 µs). No toroidal long-range correlation of the turbulence amplitude is observed, in agreement with expectations based on the relatively short poloidal turbulence correlation length. Radial edge turbulence and flow correlation lengths are longer in deuterium than in hydrogen plasmas as one would expect from ‘naïve’ gyro-Bohm isotope transport scaling, despite the substantially higher thermal flux across the last closed flux surface in hydrogen before the L-H transition. This article is part of a discussion meeting issue 'H-mode transition and pedestal studies in fusion plasmas'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Code Demonstration: Approximation of nearly-periodic symplectic maps via structure-preserving neural networks

In this notebook, we use the symplectic gyroceptron architecture from [Duruisseaux et al., 2022] to learn a surrogate map for the nearly-periodic symplectic flow map associated to a nearly-periodic Hamiltonian system composed of two nonlinearly coupled oscillators, where one of them oscillates significantly faster than the other: $\Big\lbrace$ $ {\dot{q}_1 = p_1 \hspace{5mm} \dot{p}_1 = -q_1 - \epsilon \partial_{q1}U(q_1 , q_2 )} \atop {\dot{q}_2 = \epsilon p_2 \hspace{4mm} \dot{p}_2 = -\epsilon q_2 - \epsilon \partial_{q2}U(q_1 , q_2 )}$. These equations of motion are the Hamilton's equations associated to the Hamiltonian $H_{\epsilon}(q_1,q_2,p_1,p_2)=\frac{1}{2}(q^2_1+p^2_1)+\frac{1}{2}\epsilon(q^2_2+p^2_2)+\epsilon$$U(q_1,q_2)$.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗