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 667 records · Page 37

Active nematic ratchet in asymmetric obstacle arrays

We numerically investigate the effect of an asymmetric periodic obstacle array in a two-dimensional active nematic. We find that activity in conjunction with the asymmetry leads to a ratchet effect or unidirectional flow of the fluid along the asymmetry direction. The directional flow is still present even in the active turbulent phase when the gap between obstacles is sufficiently small. We demonstrate that the dynamics of the topological defects transition from flow mirroring to smectic-like as the gap between obstacles is made smaller, and explain this transition in terms of the pinning of negative winding number defects between obstacles. This also leads to a nonmonotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

36 MATERIALS SCIENCE↗

A mountain-induced moist baroclinic wave test case for the dynamical cores of atmospheric general circulation models

Idealized test cases for the dynamical cores of atmospheric general circulation models are informative tools to assess the accuracy of the numerical designs and investigate the general characteristics of atmospheric motions. A new test case is introduced that is built upon a baroclinically unstable base state with an added orographic barrier. The topography is analytically prescribed and acts as a trigger of both baroclinic Rossby waves and inertia–gravity waves on a rotating, regular-sized planet. Both dry and idealized moist configurations are suggested. The latter utilizes the Kessler warm-rain precipitation scheme. The test case enhances the complexity of the existing test suite hierarchy and focuses on the impacts of two midlatitudinal mountain ridges on the circulation. Selected simulation examples from four dynamical cores are shown. These are the Spectral Element and Finite Volume dynamical cores, which are part of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM), versions 2.1.3 and 2.2, and the Cubed-Sphere Finite Volume dynamical cores, which is new to CESM version 2.2. In addition, the Model for Prediction Across Scales (MPAS) is tested. The overall flow patterns agree well in the four dynamical cores, but the details can vary greatly. The examples highlight the broad palette of use cases for the test case and reveal physics–dynamics coupling issues.

58 GEOSCIENCES↗

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗

Interpretable Models for Workflow Differentiation in High-Performance Scientific Networks

Scientific workflows in high-performance networks spawn hundreds of interdependent flows that must be managed collectively—yet existing network classifiers treat each flow in isolation, leading to fragmented QoS decisions and missed interflow patterns. We present a novel traffic classification solution that operates at the workflow level, distinguishing entire filetransfer operations from streaming analytics by capturing how concurrent flows interact and burst together. We introduce a workflow identification window (WIW) that ingests raw packet headers from parallel flows into unified tensors, preserving the spatial-temporal patterns that differentiate scientific workflows. This approach achieves 98.7% accuracy using CNN, LSTM, and hybrid architectures, while maintaining 84% accuracy on production traffic collected a week later—demonstrating robustness to temporal drift. By integrating SHAP and GradCAM explainability, we reveal that early-packet timing patterns and cross-flow correlations drive classification decisions, providing operators with interpretable insights. Our system enables coherent workflow-level QoS enforcement and dynamic bandwidth allocation in scientific networks, eliminating manual per-flow configuration while maintaining classification latency at millisecond level.

Giannakou, Anna [LBL, Berkeley]↗

Opportunistic Short‐Term Water Uptake Dynamics by Subalpine Trees Observed via In Situ Water Isotope Measurements

Abstract Variations in tree water sources are important to understand in semi‐arid ecosystems because climatic shifts towards lower snowpack and increased drought affect water availability in subalpine forests of the western US. Here, we use daily in situ measurements of stable isotopes ( 2 H & 18 O) in soil and tree stem water, soil matric potential and sap flow to study tree water uptake dynamics. We instrumented three soil profiles down to 90 cm, as well as three aspen and engelmann spruce trees near Gothic, Colorado, in the East River watershed. We observed the fate of natural isotopic variations in rainfall, soil, and plants from June to October 2022, and in August 2023 we conducted a 2 H labeled irrigation experiment. Our observations showed that all studied aspen trees compensated for water scarcity in the shallow soil by shifting the dominant water source at 60(±20) cm to ⅔ of uptake from 90 cm within a few days of a dry period. Both species relied on snowmelt stored in the subsoil to sustain transpiration. Intense rainfall caused the plant water uptake to shift partially to top soil layers within 2 days. Spruce transpiration was lower and relied more on snowmelt, because rainfall infiltration was low in the spruce stand due to high canopy interception. Our findings highlight the important role of snowmelt stored in the deep soil layers for subalpine forest drought response and the dominant fate of monsoonal rainfall to become transpiration rather than recharging groundwater and streams in the Upper Colorado River. Plain Language Summary There is a need to understand how trees in mountainous regions respond to dry conditions that lead to water scarcity, because climate projections suggest that such conditions will become more frequent in the future. Here we present a novel data set of measurements of daily stable isotopes of water across soil profiles and in tree stems of aspen and spruce. Our data show that when the upper soil dried out, aspen trees shifted to using water from deeper layers (beneath 60 cm) to keep transpiring. For spruce trees the uptake pattern is less clear, but both types of trees mainly used snowmelt stored in the deeper soil layers to survive the dry summer. After heavy rain, aspen and spruce trees switched to using water from the top 20 cm of soil. However, for spruce, only some rain reached the soil because the dense tree canopy intercepted it, so spruce trees stayed more dependent on snowmelt and used less water overall. This study shows how important deep snowmelt water is for helping forests survive dry periods and suggests that most summer rain is quickly used by trees rather than replenishing streams and groundwater in the headwaters of the Colorado River. Key Points Tree water resources changed within a few days from snow dominated to higher share of rainfall as soils wetted up after a dry period Compensatory plant water uptake by aspen from the deep layer (90 cm), while uptake from soil depths that became drier (60 cm) declined Strong differences between water sources and availability beneath aspen and spruce, respectively

Sprenger, Matthias↗

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗

Developing of Quaternary Pumped Storage Hydropower for Dynamic Studies

Quaternary pumped storage hydropower (Q-PSH) technology, as one of the new advanced-PSH technology, has been developed by taking advantage of Conventional-PSH (C-PSH) and Adjustable Speed-PSH (AS-PSH). By combining adjustable-speed pump unit and conventional hydropower turbine unit in the quaternary configuration, Q-PSH has the more competitive capability of providing fast power support in the future high renewable penetrated power system. Acting as energy storage (ES), Q-PSH provides promising power supply to deal with the uncertainty and variability from renewable energy generation. This paper focuses on the dynamic modeling of Q-PSH technology employing full-converter machine and the impact of Q-PSH on the frequency response in a system. The detailed model of Q-PSH is developed and implemented in the IEEE 14-bus system based on GE Positive Sequence Load Flow (PSLF) platform, which captures the dynamic of multiple operation modes, especially hydraulic short-circuit (HSC) operation mode. Several cases are set up to reveal the advantages of Q-PSH technology when power electronic based renewable energy generation is deployed in the system. Sensitivity studies of the controller in pump governor show the impact of parameters in pump response performance. The comparison case illustrates the impact of frequency response provided by the Q-PSH in the system.

dynamic modeling↗

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗

Verification and Validation of Spectral Element Code for Supercritical CO2 Flow in Vertical Heated Tubes

The investigation of heat transfer in supercritical CO2 (sCO2) has garnered considerable attention in recent decades, given sCO2's potential as a promising working fluid for advanced power conversion cycles. Despite previous research efforts, there are still gaps in our understanding of sCO2 heat transfer, particularly in conditions associated with heat transfer deterioration. To delve into sCO2 heat transfer more comprehensively, we propose employing the high-fidelity computational fluid dynamics code NekRS to simulate sCO2 flow using the large eddy simulation technique. Through graphics processing unit acceleration, NekRS achieves a higher computational speed than traditional CPU-based systems. However, before using NekRS in practical applications involving sCO2, it is imperative to perform verification and validation. Here, this paper presents our efforts to verify and validate the NekRS code's capability for simulating sCO2 using heated vertical tubes, where heat transfer deterioration usually happens. To accommodate the unique properties of sCO2, we have modified the NekRS code by integrating third-party property modules, such as REFPROP and PROPATH. Our simulations are compared with experimental and numerical data from the literature, instilling confidence in leveraging NekRS for future engineering applications. Our simulations also reveal that the accuracy of the property module significantly impacts the results, with REFPROP outperforming PROPATH for sCO2 properties. Additionally, we observed that, depending on the flow direction, buoyancy can either enhance or suppress turbulence in sCO2 flow. In upward flow, under certain conditions, the suppressed turbulence leads to heat transfer deterioration, resulting in elevated wall temperatures.

NekRS↗

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Exploring the effect of ELM and code-coupling frequencies on plasma and material modeling of dynamic recycling in divertors

Abstract Integrated modeling of plasma-surface interactions provides a comprehensive and self-consistent description of the system, moving the field closer to developing predictive and design capabilities for plasma facing components. One such workflow, including descriptions for the scrape-off-layer plasma, ion-surface interactions and the sub-surface evolution, was previously used to address steady-state scenarios and has recently been extended to incorporate time-dependence and two-way information flow. The new model can address dynamic recycling in transient scenarios, such as the application presented in this paper: the evolution of W samples pre-damaged by helium and exposed to ELMy H-mode plasmas in the DIII-D DiMES. A first set of simulations explored the effect of ELM frequency. This study was discussed in detail in this conference’s proceedings and is summarized here. The 2nd set of simulations, which is the focus of this paper, explores the effect of code-coupling frequency. These simulations include initial SOLPS solutions converged to the inter-ELM state, ion impact energy ( E in ) and angles ( A in ) calculated by hPIC2, and an improved heat transfer description in Xolotl. The model predicts increases in particle fluxes and decreases in heat fluxes by 10%–20% with the coupling time-step. Compared with the first set of simulations, the less shallow impact angle leads to smaller reflection rates and significant D implantation. The higher fraction of implanted flux (and deeper), in particular during ELMs, increases the accumulated D content in the W near-surface region. Future expansion of the workflow includes coupling to hPIC2 and GITR to ensure accurate descriptions of E in and A in , and W impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗

Discrete Element Method Analysis for Metal Powders Used in Additive Manufacturing, and DEM Simulation Tutorial Using LIGGGHTS-PUBLIC [PowerPoint and paper]

Discrete Element Method (DEM) is a method of analysis to evaluate the dynamic interactions between granular particles. This method has been used in the pharmaceutical industry to improve the powder compaction process for tablet manufacturing. There are also applications in agriculture, food industry, and manufacturing. Direct energy deposition is an additive manufacturing technique which uses metallic powders fed through a nozzle, melted using a directed laser, and transformed into a solid object layer by layer. One way of feeding metal particles into the system involves the use of a vibrating hopper. Given a specified amplitude and frequency input, the hopper will enable the powder to travel up a path, and inject through the system with assistance from a stream of gas. The mechanical properties of a printed object can vary, depending on the characteristics of the powder flow and the particles’ as-received properties. Improved understanding of dynamic interactions of flowing powders could enable additive manufacturing components with 2D or 3D variations in mechanical properties, e.g., density. This work uses DEM simulation software to investigate the effects of particle cohesion, friction, and density on the quality of the flow by performing an angle of repose simulation, which is often used as a metric to evaluate the flowability of powders.

36 MATERIALS SCIENCE↗

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗

The effect of dynamic near-wake modulation on utility-scale wind turbine wake development

High-resolution field-scale experiments using flow visualization with natural snowfall and high-fidelity large eddy simulations are combined to investigate the effect of dynamic turbine operation and atmospheric conditions on wind turbine wake mixing and recovery in the wake of a 2.5 MW wind turbine. Instantaneous near-wake expansion and deflection in response to changes in blade pitch and wind direction, termed dynamic wake modulation, is quantified using both techniques, demonstrating excellent agreement. The simulations are used to extend these results by calculating the energy flux into the wake 7 rotor diameters downstream, showing that dynamic turbine-atmospheric interactions enhance mixing in the far-wake. This finding is exhibited under both uniform and turbulent inflow conditions. Under turbulent flow, a synergistic relationship is also observed between dynamic wake modulation and wake meandering, as wake recovery can be further accelerated when the two phenomena occur together. The results of this study have implications for the development of more realistic far-wake models that include the significant impact of dynamic wake modulation on wake mixing and development. Additionally, the findings from the current study can be used to develop advanced control algorithms to speed up wake breakdown and recovery, further improving wind farm efficiency.

17 WIND ENERGY↗