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 1,045 records · Page 58

Probabilistic learning on manifolds constrained by nonlinear partial differential equations for small datasets

A novel extension of the Probabilistic Learning on Manifolds (PLoM) is presented. It makes it possible to synthesize solutions to a wide range of nonlinear stochastic boundary value problems described by partial differential equations (PDEs) for which a stochastic computational model (SCM) is available and which depend on a vector-valued random control parameter. The cost of a single numerical evaluation of this SCM is assumed to be such that only a limited number of points can be computed for constructing the training dataset (small data). Each point of the training dataset is made up of realizations from a vector-valued stochastic process (the stochastic solution) and the associated random control parameter on which it depends. The presented PLoM constrained by PDE allows for generating a large number of learned realizations of the stochastic process and its corresponding random control parameter. These learned realizations are generated so as to minimize the vector-valued random residual of the PDE in the mean-square sense. Appropriate novel methods are developed to solve this challenging problem. Three applications are presented. The first one is a simple uncertain nonlinear dynamical system with a nonstationary stochastic excitation. The second one concerns the 2D nonlinear unsteady Navier–Stokes equations for incompressible flows in which the Reynolds number is the random control parameter. Here, the last one deals with the nonlinear dynamics of a 3D elastic structure with uncertainties. The results obtained make it possible to validate the PLoM constrained by stochastic PDE but also provide further validation of the PLoM without constraint.

Machine learning↗

River Dynamics Control Transit Time Distributions and Biogeochemical Reactions in a Dam-Regulated River Corridor

Residence Time Distributions (RTDs) exert an important control on biogeochemical translation in the hydrological systems. Previous tracer studies have revealed that RTDs often followed time-invariant exponential, lognormal, or heavy-tailed RTDs that have power-law behavior for long tails in headwater or low-order streams. However, there is a recognition that RTDs can be more complicated and time-variable in response to dynamic hydrological forcing. Here, we use particle tracking to estimate RTDs along the Hanford Reach of the Columbia River and to quantify the influences of river stage fluctuations. Particle tracking is conducted using the velocity field from high-resolution 3D groundwater flow simulations. The effects of hydrological forcing on the residence time distribution were evaluated by varying river flow boundary conditions and releasing particles in different time windows. Our results revealed that dynamic stage fluctuations created rapidly changing losing-gaining conditions in the river and led to highly transient RTDs, which contributes to multiple modes of RTDs. Dam-induced high-frequency (sub-daily) flow variation contributes to the short-time (sub-daily) component of the RTDs. Deviation of the reactant consumption under the single-mode assumption compared to the multimodal RTDs is relatively small (~5%) and appears when the Damköhler number was close to one.More specifically, high-frequency (daily) dam-induced stage variations have higher impacts on biogeochemical reactions with faster reaction rate, and can potentially increase the reactant consumption by 26.68%. Based on these findings, we suggest that current river basin models could be improved by including bank storage and more complex RTDs influenced by both short and long term river stage fluctuations.

54 ENVIRONMENTAL SCIENCES↗

A DATA EFFICIENT SPARSE MODELING FRAMEWORK FOR POWER ESTIMATION IN WATER TREATMENT SENSING OPERATIONS

With increasing freshwater scarcity, advanced process design mechanisms such as Closed-Circuit Reverse Osmosis (CCRO) and Digital/Physical Twin systems are gaining traction in water treatment and reuse operations. While digital and physical twin models enable improved system insight and control, their development is often expensive and computationally intensive, requiring large volumes of synthetic or experimental data to characterize underlying process dynamics. This work introduces a sparse surrogate modeling framework to estimate power consumption from measured flow and pressure variables, along with their nonlinear polynomial and interaction expansions. To ensure model reliability and reduce overfitting, a two-stage pipeline is proposed. First, a dynamic data filtering algorithm is employed to remove uninformative observations and transient operational states. Second, a sparse penalized regression technique is applied to select a minimal set of parsimonious features. The proposed model achieves high sparsity, retaining only 7 out of 34 candidate features (≈79.41% sparsity) while delivering a root mean square error (RMSE) of 0.072 on the test dataset.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

Temperature outweighs light and flow as the predominant driver of dissolved oxygen in US rivers

The concentration of dissolved oxygen (DO), an important measure of water quality and river metabolism, varies tremendously in time and space. Riverine DO is commonly perceived as regulated by interacting and competing drivers (light, temperature and flow) that define rivers’ climate. Its continental-scale drivers, however, have remained elusive, partly due to the scarcity and spatio-temporal inconsistency of water quality data. Here we show, via a deep learning model (long short-term memory) trained using data from 580 rivers, that temperature predominantly drives daily DO dynamics in the contiguous United States. Light comes a close second, whereas flow imparts minimal influence. This work showcases the promise of using deep learning models for data filling that enables large-scale systematic analysis of patterns and drivers. Results show fairly accurate prediction of DO by temperature alone, and declining DO in warming rivers, which has important implications for water security and ecosystem health in the future climate.

54 ENVIRONMENTAL SCIENCES↗

Energy Flows in Static and Programmable Catalysts

Programmable catalysts that change on the time scale of a catalytic cycle provide a new opportunity to control the flow of energy to reactants and products to promote faster and more selective chemistry. While traditional chemical manufacturing processes consume energy to achieve favorable reaction conditions, programmable catalysts aim to dynamically add or remove energy to catalytic cycles through perturbations of the catalytic surface via strain, charge, or light. These surface energy flows are quantified by the changes in adsorbate binding energy with time, and the overall efficiency relating energy inputs to catalytic performance is defined by the characteristics of the undulating catalytic surface. Here, understanding and quantifying energy flows in programmable catalysts provides baseline definitions and metrics for comparing dynamic conditions and identifying optimal catalytic performance for more efficient chemical manufacturing.

Binding energy↗

Scaling and Complexity of Stress Fluctuations Associated with Smooth and Jerky Flow in FeCoNiTiAl High-Entropy Alloy

Recent observations of jerky flow in high-entropy alloys (HEA) revealed a high role of self-organization of dislocations in their plasticity. The present work reports the first results of the investigation of stress fluctuations during plastic deformation of an FeCoNiTiAl alloy, examined in a wide temperature range covering both smooth and jerky flow. These fluctuations, which accompany the overall deformation behavior representing an essentially slower stress evolution controlled by the work hardening, were processed using complementary approaches comprising Fourier spectral analysis, refined composite multiscale entropy, and multifractal formalisms. The joint analysis at distinct scales testified that even a macroscopically smooth plastic flow is accompanied by nonrandom fluctuations, disclosing the self-organized dynamics of dislocations. Qualitative changes in such a fine-scale “noise” were found with varying temperature. The observed diversity is significant for understanding the relationships between different scales of plasticity of HEAs and crystal materials in general.

36 MATERIALS SCIENCE↗

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantifying Energy Flows in PV Circularity Processes

As sustainable deployment and end-of-life management become a hot topic to timely address in the PV community, a dynamic comparative evaluation of the benefits of circular pathways such as reuse, and remanufacturing, recycling has not been performed holistically beyond material flows or LCA analysis. Energy flows are critical for evaluating energy generation technologies. Previously they have been used to compare renewables to fossil generation and then between PV technologies. This paper quantifies energy flows to evaluate circular pathways for PV. The energy flows tracking manufacturing, generation, and losses complementary to the mass flows of silicon are quantified leveraging the PV ICE framework.

circular economy↗

A modified model parametrization algorithm for solving a special type of heat and mass transfer systems

A new method for solving nonlinear heat and mass transfer design tasks was considered. Systems using the Number of Transfer Units (NTU) method are a special type of mathematical model of heat and mass exchangers. It was observed, that the NTU models in a form of differential-algebraic equations (DAEs) cannot be directly solved with higher values of NTU. The requirements for consistent initial conditions, as well as numerical limitations of DAEs solvers, result, that the solution to the considered design problems that cannot be obtained by a classical direct shooting procedure. To overcome the presented difficulties, the αDAE model optimization algorithm was adjusted for solving NTU-based models. The new approach consists of 3 main steps: 1) task discretization by a multiple-shooting approach, 2) design an appropriate function $f_{NTU}$(α) to effectively influence the variability of the state variables described by dynamical relations, 3) the iterative numerical optimization algorithm for the new parametrized system. Moreover, computations can be performed by a chosen numerical optimization approach, which can be communicated with an available outer procedure for solving differential-algebraic equations. The presented algorithm was implemented and applied to solve the design task with the NTU model of a counter-flow exchanger. Here, the new approach was used to modify the system dynamics to influence the difficulty of the considered problem. Finally, the presented method enabled failure-free numerical computations for the higher values of the NTU parameter.

97 MATHEMATICS AND COMPUTING↗

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗

Air quality and comfort constrained energy efficient operation of multi-zone buildings

Maintaining indoor air quality (IAQ) through effective ventilation is essential for the well-being and productivity of building occupants. Control strategies aimed at improving the efficiency of heating, ventilation and air conditioning (HVAC) systems must jointly determine ventilation and heating and cooling processes. Here, in this paper, we study the problem of minimizing the energy consumption of the HVAC system in a multi-zone building, while meeting thermal comfort and IAQ requirements. We first perform a steady state analysis of the zonal carbon dioxide (CO 2 ) concentration and the temperature dynamics. The resulting expressions are convex in the zonal mass flow rates and zonal temperatures. Guided by the steady state solutions for meeting the thermal comfort constraints, we develop two control policies for improving the energy efficiency of building HVAC systems while jointly satisfying indoor temperature and IAQ constraints. We compare the performance of our proposed approaches with those of multiple baseline approaches which implement separate regimes for controlling zonal temperature and IAQ for a typical work-day in a multi-zone campus building. We have evaluated the performance of our proposed approaches under varying levels of flexibility in zonal temperatures. We have shown that zonal temperature flexibility can result in energy savings up to 32% (for the same control strategies) as compared to the case where no such flexibility is permitted. Our proposed approaches were seen to offer potential savings of nearly 29% compared to the baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Distinct Gas-Particle Partitioning and Viscosity Characteristics of Secondary Organic Aerosols Derived from α-Pinene versus Ocimene

Secondary organic aerosols (SOA) have complex, multicomponent composition that controls particle viscosity and gas-particle partitioning, key factors to their atmospheric evolution. This study investigates the chemical composition, volatility and viscosity of SOA formed by ozonolysis of cyclic α-pinene (PSOA) and acyclic ocimene (OSOA) monoterpenes. Using Temperature-Programmed Desorption combined with Direct Analysis in Real-Time ionization and High-Resolution Mass Spectrometry, we determined the molecular composition and saturation mass concentration of individual SOA constituents. These data enabled gas-particle partitioning and viscosity estimates under varied atmospheric conditions. PSOA, composed of higher molecular weight and less oxidized species, shows higher condensability and viscosity under high total organic mass (tOM) loadings. Here, in contrast, OSOA, consisting of more oxidized, lower molecular weight species, exhibits greater sensitivity to tOM, with viscosity increasing significantly upon dilution. Poke-flow experiments support this trend, indicating that OSOA undergoes more dynamic compositional and phase changes during atmospheric aging. These observations reveal distinct dynamic trends in the atmospheric transformations and reactivity of SOA from cyclic and acyclic monoterpenes, with the latter showing greater compositional changes during aging that alter viscosity and diffusion. This highlights the importance of incorporating such dynamic transformations into atmospheric models to improve predictions of SOA atmospheric loadings, lifetimes, and impacts.

cyclic and acyclic monoterpenes↗

Simulation of Microcapsule Transport in Fractured Media Using Coupled CFD‐DEM

Geothermal energy is sustainable and gaining momentum as a solution to energy crises and environmental issues. However, challenges like production temperature and thermal breakthrough can impact geothermal project efficiency. One innovative solution to alleviate the thermal breakthrough is to inject polymer-based materials that are encapsulated in microcapsules into fractures to modify fracture permeability and prevent preferential flow. In our study, we utilized a coupled computational fluid dynamics and discrete element method to simulate the transport of microcapsules under various scenarios controlled by microcapsule size, microcapsule concentration, and fracture roughness. For a smooth fracture, the results indicate that small microcapsules can travel through a smooth fracture regardless of their concentrations. Large microcapsules can transport through a smooth fracture when present in lower concentrations. However, medium and mixed-size microcapsules tend to cause the sealing of a smooth fracture, irrespective of their concentrations. For a rough fracture, the transport of microcapsules is complicated by their interactions with the rough fracture walls. The presence of two sealing positions in a rough fracture adds further complexity to this transport phenomenon. The size and concentration of microcapsules control one sealing location, while the rough fracture walls determine the other sealing location. The rough walls substantially affect microcapsule transport, rendering the role of microcapsule size and concentration less significant. The simulation results suggest that complex fracture surfaces significantly elevate the occurrence of sealing behavior. To mitigate sealing behavior within more complex fractures, it would be beneficial to use smaller and lower concentrations of microcapsules.

15 GEOTHERMAL ENERGY↗

In situ multimodal imaging for nanoscale visualization of tribofilm formation

A third of the energy from fuel combustion in passenger car gasoline engines is lost due to friction. Carefully designed engine lubricants can recover some of these losses by reducing friction and wear by forming a nanometer-scale chemico-physico tribofilm between surfaces. Accordingly, attention has focused on developing oil formulations that form low-friction tribofilms. However, analyses of resultant tribofilms are typically conducted after tribo-tests with conventional characterization tools and do not offer insights into tribofilm formation and evolution, precluding information critical to tuning tribofilm properties. In this work, we developed a unique multimodal methodology based on Atomic Force Microscopy (AFM) with local probe heating for in situ tribological studies that activates friction modifiers and simultaneously captures the evolution of friction and surface roughness, with nanometer resolution. As a platform to demonstrate the ability of this methodology to visualize dynamics of tribofilm formation in situ, we apply it to molybdenum-based friction modifiers to distinguish key factors in their functionality and correlate nanoscale AFM and Friction Force Microscopy data to bench tribo-tests used in the industry. To decode the formation mechanisms observed in situ and underlying chemistry of tribofilms, we performed ab initio Molecular Dynamics (AIMD) simulations at comparable conditions. AIMD simulations confirmed both nanoscale and bench tribo-test results and showed deviations in molecular organization in tribofilms that are formed from different molybdenum-based friction modifiers that corroborates with surface functionality. With this innovative methodology, we demonstrate proof-of-principle in situ formation of molybdenum-based tribofilms directly on steel surfaces that could be applied generally to studying tribofilm formation.

36 MATERIALS SCIENCE↗

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗