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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

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

Community Fabric: Visualizing communities and structure in dynamic networks

We present Community Fabric, a novel visualization technique for simultaneously visualizing communities and structure within dynamic networks. In dynamic networks, the structure of the network is continuously evolving throughout time and these underlying topological shifts tend to lead to communal changes. Community Fabric helps the viewer more easily interpret and understand the interplay of structural change and community evolution in dynamic graphs. To achieve this, we take a new approach, hybridizing two popular network and community visualizations. Community Fabric combines the likes of the Biofabric static network visualization method with traditional community alluvial flow diagrams to visualize communities in a dynamic network while also displaying the underlying network structure. Our approach improves upon existing state-of-the-art techniques in several key areas. We describe the methodologies of Community Fabric, implement the visualization using modern web-based tools, and apply our approach to three example data sets.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

MSASGCN : Multi-Head Self-Attention Spatiotemporal Graph Convolutional Network for Traffic Flow Forecasting

Traffic flow forecasting is an essential task of an intelligent transportation system (ITS), closely related to intelligent transportation management and resource scheduling. Dynamic spatial-temporal dependencies in traffic data make traffic flow forecasting to be a challenging task. Most existing research cannot model dynamic spatial and temporal correlations to achieve well-forecasting performance. The multi-head self-attention mechanism is a valuable method to capture dynamic spatial-temporal correlations, and combining it with graph convolutional networks is a promising solution. Therefore, we propose a multi-head self-attention spatiotemporal graph convolutional network (MSASGCN) model. It can effectively capture local correlations and potential global correlations of spatial structures, can handle dynamic evolution of the road network, and, in the time dimension, can effectively capture dynamic temporal correlations. Experiments on two real datasets verify the stability of our proposed model, obtaining a better prediction performance than the baseline algorithms. The correlation metrics get significantly reduced compared with traditional time series prediction methods and deep learning methods without using graph neural networks, according to MAE and RMSE results. Compared with advanced traffic flow forecasting methods, our model also has a performance improvement and a more stable prediction performance. We also discuss some problems and challenges in traffic forecasting.

Cao, Yang↗
DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ 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. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate 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 shown 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.

Kim, Jinyong↗
DOE OSTI2025

Flow-Driven Stretch Fluctuations Govern the Nonlinear Viscoelasticity of Elongating Associative Polymer Networks

We use nonequilibrium molecular dynamics simulations to verify recent tube-model predictions that associative polymer networks exhibit broad stretch fluctuations during elongational flow. Simulations further show that these fluctuating dynamics give rise to the rate-dependent extensional viscosity 𝜂 𝐸 measured in filament stretching experiments on H-bonding networks. Simulations model bivalent associative networks with a reactive bead-spring model for varying association strength and extensional strain rate. We observe that stretch fluctuations are driven by a new form of chain tumbling, where chains continually collapse and elongate as their associations break and reform within the advecting network. This produces a broad, nearly uniform distribution of chain stretch over a wide range of strain rates, manifesting as a rate-independent plateau in the extensional stress. Our results show that the nonlinear viscoelasticity of associative networks is dominated by large fluctuations in molecular response, which cannot be captured by current mean-field models.

biomimetic & bio-inspired materials↗
DOE OSTI2026

A methodology for domain overlapping coupling of thermal-hydraulic systems

Multi-scale coupling has increasingly drawn attention as a promising approach for modeling thermal systems. Traditional system codes provide fast and robust predictions at the plant scale, while high-fidelity computational fluid dynamics (CFD)-based tools resolve localized flow and heat transfer phenomena with greater accuracy. By combining these complementary strengths, co-simulations enable multi-scale analysis that would otherwise be computationally prohibitive for a standalone CFD code. Here, this work introduces a robust and problem-agnostic domain overlapping (DO) coupling between the system thermal-hydraulic (STH) code System Analysis Module (SAM) and the coarse-mesh CFD code Pronghorn. Both applications belong to the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite, a code suite in active development at the Idaho National Laboratory (INL), tailored for multi-physics analysis of advanced reactors. Unlike previous approaches, BlueCRAB supports an agnostic interface between codes based on different fidelity, while its coupling formulation can address arbitrary flow geometries with multiple inlets and outlets in coupled components. The implemented method leads to consistent pressure drops, enthalpies, and scalar concentrations between coupled SAM and Pronghorn simulations. The methodology is demonstrated through two verification tests, which ensure the numerical consistency and conservation across the codes, and through one validation test against experimental data. The proposed problems explore different physical aspects inherent to thermal systems, with particular attention given to nuclear reactor analysis. These include buoyancy-driven flows, complex flow patterns, and setups with multiple inlets and outlets, representing challenges in advanced reactor applications.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transition Zone Morphology Dynamics of Dissolved Oxygen (DO) in a Salinity-Impacted Hyporheic Zone

Freshwater salinization is a common yet underappreciated environmental problem in rivers, yet how this process impacts the transport of dissolved oxygen (DO) in the hyporheic zone remains unknown. Here, using flume experiments and numerical simulations, we have demonstrated that the morphological dynamics of the DO transition zone are controlled by the dimensionless number Da·Ra –2 , a measure of the rate of aerobic respiration to free convection. At low Da·Ra –2 , the oxic front breaks up into fingers which subsequently grow, and leads to accelerating of DO transport and shrinking of the anaerobic zone in the underlying sediment. As Da·Ra –2 increases, the growth rate of the instability decreases, and the DO plume fingers are suppressed and delayed compared with the saltwater-freshwater interface, the mixing area between saltwater and freshwater. These results indicate that the freshwater salinization syndrome can have significant impacts on the functioning of aquifer ecosystems unless regulated and managed effectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quantifying particle-scale 3D granular dynamics during rapid compaction from time-resolved in situ 2D x-ray images

Understanding the particle-scale dynamics of granular materials during rapid compaction and flow is of fundamental importance for manufacturing, planetary science, geology, and defense applications. Time-resolved 2D radiography and static 3D x-ray tomography are powerful in situ tools for studying particle-scale dynamics but provide detail only in 2D or with significant time-scale limitations, respectively. Here, we introduce a new method that uses 2D in situ x-ray imaging for determining time-resolved 3D particle-scale dynamics in rapidly compressed granular materials. The method employs initial particle packing structures obtained from x-ray tomography, a 2D x-ray image generation algorithm, and an optimization algorithm. We first describe and validate the method using finite element simulations. Furthermore, we then apply the technique to x-ray phase-contrast images obtained during rapid compaction of granular materials with varying particle morphology and sample thickness. The depth-resolved particle-scale dynamics reveal complex velocity and porosity fields evolving heterogeneously along and perpendicular to the compaction direction. We characterize these features, their fluctuations near the compaction front, and the compaction front thickness. Our technique can be applied to understanding granular dynamics during rapid compaction events, and rearrangements during slower, but non-quasi-static, flows.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗
DOE OSTI2020

Flow-Driven Branching in a Frangible Porous Medium

Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium. We provide a simple theoretical framework that embodies this feedback mechanism in a multiphase model for flow through a frangible porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. Furthermore, this allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.

36 MATERIALS SCIENCE↗
DOE OSTI2024

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 non-monotonic ratchet effect magnitude as a function of obstacle size, so that there is an optimal obstacle size for ratcheting at fixed activity.

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

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗
DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Area Localized Coupled Model for Analytical Mean Flow Prediction in Arbitrary Wind Farm Geometries

This work introduces the area localized coupled (ALC) model, which extends the applicability of approaches that couple classical wake superposition models and atmospheric boundary layer models to wind farms with arbitrary layouts. Coupling wake and top–down boundary layer models is particularly challenging since the latter requires averaging over planform areas associated with turbine-specific regions of the flow that need to be specified. The ALC model uses Voronoi tessellation to define this local area around each turbine. A top–down description of a developing internal boundary layer is then applied over Voronoi cells upstream of each turbine to estimate the local mean velocity profile. Coupling between the velocity at hub-height based on this localized top–down model and a wake model is achieved by enforcing a minimum least-square-error in mean velocity in each cell. The wake model in the present implementation takes into account variations in wind farm inflow velocity and represents the wake profile behind each turbine as a super-Gaussian function that smoothly transitions between a top-hat shape in the region immediately following the turbine to a Gaussian profile downstream. Detailed comparisons to large-eddy simulation (LES) data from two different wind farms demonstrate the efficacy of the model in accurately predicting both wind farm power output and local turbine hub-height velocity for different wind farm geometries. These validations using data generated from two different LES codes demonstrate the model's versatility with respect to capturing results from different simulation setups and wind farm configurations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗