SEARCH · Engineering Papers

Results for “variable 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 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2023

FECM/NETL CO 2 Transport Cost Model (2023): Description and User’s Manual

The FECM/NETL CO 2 Transport Cost Model (CO 2 _T_COM) is an Excel spreadsheet model that calculates the cost of transporting CO 2 from the beginning to the end of a pipeline. This document provides two main functions. First, the document describes the equations and algorithms that are used by the model to calculate technical quantities (such as the minimum inner pipe diameter needed to transport a user-specified CO 2 mass flow rate a specified distance) and engineering-economic quantities (such as capital costs, operating costs and cash flows). Second, the document is a user’s manual for the model that describes the procedures the user must follow to run the model. The document also describes input variables and output variables (i.e., results) for the model. The model can be accessed at this URL: https://www.netl.doe.gov/energy-analysis/details?id=d3086f60-278d-4e97-a649-8e4d5ce5e93c

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2020

Seasonal to Decadal Variability in Focused Groundwater and Contaminant Discharge along a Channelized Stream

Abstract Focused groundwater discharge to streams is problematic at contaminated sites because high fluxes can limit natural attenuation in the hyporheic zone. However, information on location, spatial evolution, and temporal persistence of springs in unlithified sediments over multiyear time scales is limited. We examine discharge at point (~1‐m) to reach (~300‐m) scales along a stream that intercepts trichloroethene and technetium‐99 plumes from a Superfund site. During 2011 to 2012, we seasonally monitored stream and spring flow and contaminant concentrations, along with probing streambed temperatures on a grid in winter and summer, building on prior monitoring during 1999 to 2002. Baseflow measured by both gauging and dye dilution generally increased with distance downstream, and stream and spring discharge varied seasonally, from minima in October to January to maxima in February to June. Thermal anomalies identified by probing occupied approximately 3% to 6% of the reach and typically coincided with visible springs or seeps. Locations of anomalies were similar to those identified in summer 2002, although some orifices disappeared and others emerged. Vertical groundwater fluxes calculated from probing tended to be less than net fluxes calculated from stream discharge, perhaps in part because the assumption of one‐dimensional, steady‐state flow in calculating point fluxes was simplistic. Maximum contaminant concentrations and fluxes decreased between 1999 to 2001 and 2011 to 2012 as a result of partial capture by an upgradient pump‐and‐treat system. Our findings confirm that springs in unlithified sediments can remain stationary within a few meters over decadal time scales, and seasonal variability in discharge can be greater than decadal variability.

Tripathi, Ganesh N.↗

DOE OSTI2021

Holistic design principles for flow batteries: Cation dependent membrane resistance and active species solubility

Cation dependent resistance of a commercial cation exchange membrane, Nafion™ 212, as well as the solubility of select active materials are investigated, demonstrating practical consequences of cation choice on redox flow battery (RFB) performance. The conductivity of alkali and methyl- or -alcohol substituted ammonium cations through pretreated Nafion™ 212 reveal that cation size is the dominant factor affecting ionic transport. Among the alkali series the resistance increases with larger ionic radii and similarly, with increasing steric bulk through substitution on ammonia. The solubility of ferrocyanide, a commonly used redox flow battery posolyte (or catholyte), as well as iron, cobalt, and chromium metal-organic compounds are determined for the same series of cations. Full cells are assembled using various cation mixtures and membrane pre-treatments to demonstrate practical effects of each variable. Here, we highlight the need to carefully select the correct cations for efficient electrolyte design, since certain cations may lead to high conductivity but low solubility, or vice versa, and provide design principles for RFBs using Nafion™ 212.

25 ENERGY STORAGE↗

DOE OSTI2022

A Comparison of Linear Solvers for Resolving Flow in Three-Dimensional Discrete Fracture Networks

We compare various methods for resolving steady flow within three-dimensional discrete fracture networks, including direct methods, Krylov subspace methods with and without preconditioning, and multi-grid methods. We compared the performance of the methods based on compute times and scaling of the solution as a function of the number of grid nodes and log-variance of the hydraulic aperture. The methods are applied to three test cases: (a) variable density of networks with a truncated power-law distribution of fracture lengths, (b) a fixed network composed of monodisperse fracture sizes but varied permeability/aperture heterogeneity, (c) and a network based on field site in Nevada, US. We chose these cases to allow us to study the impact of the mesh size and flow properties, as well as to demonstrate our conclusions on a large-scale, realistic problem (more than 40 million mesh nodes). A direct solution using Cholesky factorization outperformed other methods for every example but was closely followed in performance by some algebraic multigrid (AMG) preconditioned Krylov subspace methods. Among the Krylov methods, conjugate gradients (CG) with an AMG preconditioner performs the best. Generally, Cholesky factorization is recommended, but CG with an AMG preconditioner may be suitable for very large problems beyond 40 million nodes where the entire linear system cannot reside in memory.

58 GEOSCIENCES↗

DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

DOE OSTI2024

Evaluation of a practical approach for field scale moisture flow modeling in heterogeneous media at a semiarid site

Abstract A practical approach for modeling field‐scale moisture flow in a highly heterogeneous unsaturated medium is described in this study. The validity of this approach is demonstrated through comparison of the numerical simulations with field observations at a semiarid site located in southcentral Washington State. The methodology is based on upscaling the core scale hydraulic properties and combining power‐law and tensorial connectivity‐tortuosity (PA‐TCT) approaches to derive macroscopic anisotropy parameters for each hydrostratigraphic unit (HSU) identified in the field. Each heterogeneous HSU is approximated by an equivalent homogeneous medium (EHM) model for which PA‐TCT parameters are used in the flow simulations. The available field data on moisture content and matric potential are compared with steady‐state flow simulations based on the mean form of Richards' equation. While the homogenization or averaging of heterogeneities, embedded in the EHM modeling approximation, cannot capture all of the field‐scale variability, the simulated steady‐state moisture and matric potential profiles capture well the central tendency of the field data. This approach is deemed practical for assessing the fate and transport of contaminants in highly heterogeneous unsaturated media at the transport scale of hundreds of meters.

Khaleel, Raziuddin↗

DOE OSTI2025

Semi-supervised permutation invariant particle-level anomaly detection

The development of analysis methods to distinguish potential beyond the Standard Model phenomena in a model-agnostic way can significantly enhance the discovery reach in collider experiments. However, the typical machine learning (ML) algorithms employed for this task require fixed length and ordered inputs that break the natural permutation invariance in collision events. To address this, a semi-supervised anomaly detection tool is presented that takes a variable number of particle-level inputs and leverages a signal model to encode this information into a permutation invariant, event-level representation via supervised training with a Particle Flow Network (PFN). Data events are then encoded into this representation and given as input to an autoencoder for unsupervised ANomaly deTEction on particLe flOw latent sPacE (ANTELOPE), classifying anomalous events based on a low-level and permutation invariant input modeling. Performance of the ANTELOPE architecture is evaluated on simulated samples of hadronic processes in a high energy collider experiment, showing good capability to distinguish disparate models of new physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Crystal plasticity model of residual stress in additive manufacturing using the element elimination and reactivation method

Selective laser melting is receiving increasing interest as an additive manufacturing technique. Residual stresses induced by the large temperature gradients and inhomogeneous cooling process can favour the generation of cracks. In this work, a crystal plasticity finite element model is developed to simulate the formation of residual stresses and to understand the correlation between plastic deformation, grain orientation and residual stresses in the additive manufacturing process. The temperature profile and grain structure from thermal-fluid flow and grain growth simulations are implemented into the crystal plasticity model. An element elimination and reactivation method is proposed to model the melting and solidification and to reinitialize state variables, such as the plastic deformation, in the reactivated elements. The accuracy of this method is judged against previous method based on the stiffness degradation of liquid regions by comparing the plastic deformation as a function of time induced by thermal stresses. The method is used to investigate residual stresses parallel and perpendicular to the laser scan direction, and the correlation with the maximum Schmid factor of the grains along those directions. The magnitude of the residual stress can be predicted as a function of the depth, grain orientation and position with respect to the molten pool. The simulation results are directly comparable to X-ray diffraction experiments and stress–strain curves.

316 stainless steel↗

DOE OSTI2023

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Moments of the Boltzmann collision operator for Coulomb interactions

Exact moments of the Boltzmann collision operator are calculated in the irreducible Hermitian moment expansion written in terms of the random-velocity variable of each species. The formulas are presented in closed, algebraic form and can be straightforwardly implemented in computer algebra systems. They are valid for two arbitrary masses, temperatures, and flow velocities, and hence include all other existing results derived for distribution functions expanded with respect to reference states of one temperature and flow velocity. In comparison, the Landau collisional moments are good approximations for large Coulomb logarithm and small relative flow velocity, but they fail to predict the correct behavior of most collisional moments for large relative flow even for weakly coupled plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Genotype–environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence

Abstract Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush ( Catharus guttatus ), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

An engineering perspective on evaluating mechanisms governing ductility in pure molybdenum

The United States lacks a stable domestic supply of 99 mTc, a critical medical imaging isotope generated from 99 Mo. Accelerator-based production using 100 Mo targets (aMo) introduces mechanical concerns due to the ductile-to-brittle transition temperature inherent to refractory metals like Mo. This study evaluates the tensile behavior of aMo targets from 25 to 1000 °C, compared to powder-metallurgy-processed natural Mo (PMo) and cast-and-rolled Mo (RMo), both as-received and after 5 ppm O 2 exposure in flowing He. RMo showed superior ductility (7.3% at 25 °C, 48% at 1000 °C) and strength, attributed to its fine, elongated grains and high geometrically necessary dislocation density. PMo exhibited variable ductility (up to 33%), while aMo remained brittle, with a maximum elongation of 10.4% at 600 °C. EBSD analysis revealed weak texture in PMo and aMo, but high defect density in aMo limited dislocation mobility. This work links processing, microstructure, and deformation mechanisms to guide fabrication of ductile refractory targets.

Hyer, Holden C. [Oak Ridge National Laboratory (OR↗

DOE OSTI2024

Energetic Particles - microturbulence interaction thrust (EP SIWG white paper)

This white paper outlines a strategic approach to Energetic Particle (EP) research, aimed at applying the results to the future burning devices such as ITER. Our goal is to deepen the understanding of the interplay between plasma background microturbulence and diverse species of energetic ions in burning plasma devices. This strategy is poised to advance our knowledge in this crucial area of fusion science, paving the way for significant breakthroughs in plasma physics and fusion energy research. We consider the confinement of energetic ions such as auxiliary heating beams, minority ions from ICRH, and fusion product alpha particles. Recent studies have underscored the potential of EP/microturbulence interaction to enhance fusion plasma performance by stabilizing microturbulence. However, it also poses risks, such as exacerbating Alfvénic eigenmode (AE) instabilities, which could compromise the sustainability of plasma discharge by causing EP losses. The interplay between EPs and background microturbulence is evident in the phenomenon of effective pitch angle scattering, a crucial aspect of the quasilinear (QL) theory. This theory is integral to developing numerically-efficient yet comprehensive and self-consistent approaches, recently employed to investigate the relaxation of energetic particle populations in the holistic modeling of fusion-grade plasmas. The primary objectives of this white paper thrust encompass a dual focus. First is that we need to delve into the mechanisms through which microturbulence engenders effective pitch angle scattering. This investigation will entail the examination of the presence of known or self-consistently predicted spectra of modes accountable for microturbulence. In these explorations, electrostatic microturbulence serves as a logical initial stride towards achieving the trust's objectives. This endeavor is anticipated to yield formulations expressing the parametric dependencies of the effective pitch angle scattering frequency on variables such as thermal plasma electron and/or ion temperatures, as well as their respective thermal conductivities. In our second objective, we aim to delve into the intricate formation of zonal flow (ZF) structures amidst the complex interplay of microturbulence and Alfvenic eigenmodes (AEs). This endeavor poses greater challenges as we seek to unravel the macroscopic manifestations influenced by microturbulence, termed as zonal structures (ZS), stemming from microturbulence-induced ZF. Drawing from this understanding, we anticipate employing a QL approach to yield comprehensive insights into the distribution function of energetic particles (EP) within phase space. This method entails resolving the dominant multidimensional phase space diffusion processes while effectively averaging over the rapid ballistic responses. Nevertheless, substantial strides remain imperative to realize a comprehensive whole-device modeling framework. This entails meticulous verification and validation exercises against experimental observations, as well as rigorous benchmarking against theoretical frameworks and numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Selective Dynamical Imaging of Interferometric Data

Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT’s (u, v)-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the EHT can contain regions of time over the course of a single observation that facilitate dynamical imaging. These optimal time regions typically have projected baseline distributions that are approximately angularly isotropic and radially homogeneous. We derive a metric of coverage quality based on baseline isotropy and density that is capable of ranking array configurations by their ability to produce accurate dynamical reconstructions. We compare this metric to existing metrics in the literature and investigate their utility by performing dynamical reconstructions on synthetic data from simulated EHT observations of sources with simple orbital variability. We then use these results to make recommendations for imaging the 2017 EHT Sgr A* data set.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Computational Risk Analysis of Propane Releases in Maintenance Facilities

Liquefied petroleum gas (LPG) is a viable, cleaner alternative to traditional diesel fuel used in busses and other heavy-duty vehicles and could play a role in helping the US meet its lower emission goals. While the LPG industry has focused efforts on developing vehicles and fueling infrastructure, we must also establish safe parameters for maintenance facilities which are servicing LPG fueled vehicles. Current safety standards aid in the design of maintenance facilities, but additional quantitative analysis is needed to prove safeguards are adequate and suggest improvements where needed. In this report we aim to quantify the amount of flammable mass associated with propane releases from vehicle mounted fuel vessels within enclosed garages. Furthermore, we seek to qualify harm mitigation with variable ventilations and facility layout. To accomplish this we leverage validated computational resources at Sandia National Laboratories to simulate various release scenarios representative of real world vehicles and maintenance facilities. Flow solvers are used to predict the dynamics of fuel systems as well as the evolution of propane during release events. From our simulated results we observe that both inflow and outflow ventilation locations play a critical role in reducing flammable cloud size and potential overpressure values during a possible combustion event.

03 NATURAL GAS↗

DOE OSTI2023

Developing Renewable Energy Applications for Water Treatment Technologies and Transmission

The need for smaller treatment facilities for more isolated or rural municipal entities and industrial water treatment facilities continues to grow. From environmental justice and energy needs in remote communities that have never been serviced by larger scale utilities to smaller onsite treatment requirements to meet discharge regulations, energy efficient treatment systems that can leverage renewable energy sources can be an important asset to owners. Climate change and population dynamics have changed the needs for distributed water treatment systems and the manner in which they function. New treatment technologies that can operate under dynamic conditions and leverage variable energy sources (e.g., renewable energy, RE) can provide substantial efficiencies and cost savings. The National Renewable Energy Laboratory (NREL) is embarking on research projects that will perform multi-scale studies to investigate treatment technologies on varying source water quality that can operate under unsteady state flow conditions (timing) in response to energy availability, supply, and demand. A pilot-scale system will be built to assess modularity and scaling of treatment system operating conditions and performance. Modeling will also be used to determine feasibility and integration into grid emulation platforms to understand the use and value of operation with RE. Optimization will be completed for the treatment and distribution systems to document strengths and weaknesses of using RE power sources. Results of this effort will advance dynamic water treatment understanding under variable conditions and integration of RE sources to power these systems. This will lead to further decarbonization of municipal and private treatment facilities while supporting a paradigm shift in water treatment process operation and system design to achieve net-zero goals and move toward energy-sensitive resilient water treatment systems.

ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A Computationally Efficient Linear Semi-Lagrangian Scheme for the Advection of Microphysical Variables in Cloud-Resolving Models

An efficient semi-Lagrangian (SL) numerical scheme for the advection of positive-definite microphysical variables in cloud-resolving models is proposed. The scheme is a linear combination of SL schemes of the first and second order. The proposed scheme is compared with two high-order, monotonic, nonlinear Eulerian schemes, in two 2D tests. In the first test, advection of a positive scalar field with large gradients is performed within complicated idealized steering flows. In the second test, a hail storm is simulated using the 2D Hebrew University Cloud Model, describing the interactions between 15 size distributions of different microphysical quantities (each distribution is described by 43 mass bins). The proposed scheme is computationally efficient, has a low numerical diffusion and a high level of mass conservation accuracy, and preserves the sum of multiple advected variables as well as the linear relationships between the advected scalar variables. Although the proposed SL scheme is not exactly positive definite, the negative values and their number, which may appear only in the case of very strong gradients, are negligible. The scheme produces results similar to those of nonlinear Eulerian advection schemes while reducing computation time by approximately a factor of 10.

54 ENVIRONMENTAL SCIENCES↗