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.

SourceResearch connectionFrom yearThrough year

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

At least 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Computationally Aided Design of a Small-Scale Waste-Glass Melter - 20452

To provide mission support to the WTP for the vitrification of legacy nuclear tank waste at the Hanford site, a reduced-scale vitrification pilot system is being designed to process simulated and actual radioactive tank wastes. The waste stored in underground tanks will be separated into HLW and LAW fractions, where the majority is LAW by mass (approximately 90%) and HLW by activity (greater than 95%). The first tank waste to be processed at the WTP will be LAW. Because Hanford tank waste and the resultant melter feed compositions are known to vary widely, pilot-scale operations are essential to identify potential problems, provide needed data, confirm assumptions, and determine impacts to full-scale melters and off-gas systems from these different feeds. In addition, the reduced-scale melter system must be capable of quickly providing results to operations without the typically high cost of full-scale radioactive operations and to provide an engineering platform in close proximity to local operations staff. The objectives of this computational modeling effort are to produce process conditions like those encountered in the WTP LAW melters and to minimize the volume of radioactive waste necessary for the evaluations. To this end, melter design features are being explored to increase production without significantly increasing surface area or melter-glass volumes or compromising data quality. In support of these objectives, a set of computational fluid dynamic models was developed to provide insight into the flow patterns within these small melters and evaluate design features to increase production without significantly increasing glass inventory. Computed velocities were evaluated at the interface between the cold cap and the molten-glass pool and within the bulk-glass pool. Since melt rate correlates with the velocity underneath the cold cap, this information can be used to select design features that produce the highest gains. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Mass Transport in Membrane Systems: Flow Regime Identification by Fourier Analysis

The numerical calculation of local mass distributions in membrane systems by computational fluid dynamics (CFD) offers indispensable benefits. However, the concept to calculate such distributions in response to separate variations of operation conditions (OCs) makes it difficult to address overall, flow-physics-related questions, which require the consideration of the collective interaction of OCs. It is shown that such understanding-related relationships can be obtained by the analytical solution of the advection–diffusion equation considered. A Fourier series model (FSM) is presented, which provides exact solutions of an advection–diffusion equation for a wide range of OCs. On this basis, a new zeroth-order model is developed, which is very simple and as accurate as the complete FSM for all conditions of practical relevance. Advection-dominated blocked and diffusion-dominated unblocked flow regimes are identified (depending on a Péclet number which compares the flow geometry with a length scale imposed by the flow), which implies relevant requirements for the use of lab results for pilot- and full-scale applications. Analyses reveal the equivalence of variations of OCs, which offers a variety of options to accomplish desired flow regime changes.

Heinz, Stefan (ORCID:0000000248712416)↗

DOE OSTI2022

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Rectified Rotational Dynamics of Mobile Inclusions in Two-Dimensional Active Nematics

We investigate the dynamics of mobile inclusions embedded in 2D active nematics. The interplay between the inclusion shape, boundary-induced nematic order, and autonomous flows powers the inclusion motion. Disks and achiral gears exhibit unbiased rotational motion, but with distinct dynamics. In comparison, chiral gear-shaped inclusions exhibit long-term rectified rotation, which is correlated with dynamics and polarization of nearby +1/2 topological defects. Furthermore, the chirality of defect polarities and the active nematic texture around the inclusion correlate with the inclusion’s instantaneous rotation rate. Inclusions provide a promising tool for probing the rheological properties of active nematics and extracting ordered motion from their inherently chaotic motion.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗

DOE OSTI2020

Multi-Scale Simulation of Wind Farm Performance during a Frontal Passage

Predicting the response of wind farms to changing flow conditions is necessary for optimal design and operation. Here in this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine–flow interaction, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with 10 m grid spacing. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed background velocity, turbine wake structure, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD model framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗

DOE OSTI2023

Coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) for engineering flow measurements

Increasingly, fully 3D experimental measurements of flow in complex engineering geometries are required to validate computational fluid dynamics models that support and inform reactor design and licensing. One barrier to such measurements is the complexity of typical reactor components and subsequent lack of optical access in these systems. To overcome this, the deployment of coregistered positron emission particle tracking (PEPT) and X-ray computed tomography (CT) is explored for flow measurement in reactor thermal hydraulic components and model (scaled) systems. Through this methodology, fully 3D flow information (via PEPT) and detailed internal geometry (via CT) are captured in opaque systems such as pipes, rod bundles, packed beds, etc. The reconstructed flow field and geometry can then be overlain to reveal detailed flow features around internal structures within a given test section. This is enabled through the use of a combined preclinical PET/CT scanner with overlapping PET and CT fields of view. Such measurements are useful for characterizing flow inside such intricate nuclear thermal hydraulic components as core geometries and heat exchangers, among others, and providing valuable 3D validation data for CFD models. In this work, basic tests of this 3D flow/geometry mapping are presented, and the implications of such measurements are discussed. Further, preliminary measurements are made with both point sources and flow in a simple pipe flow geometry to evaluate the capabilities of this technique. PEPT and CT features are coregistered with up to 0.1 mm precision, and pipe flow mean velocity and Reynolds stresses are reconstructed with similar accuracy to previous PEPT demonstrations. The utility of PEPT/CT is shown herein, and suggestions for future measurements are made.

3D flow measurement↗

DOE OSTI2024

Code for the manuscript "Mori-Zwanzig Modal Decomposition"

We would like to create an open source repository in LANL's github on code written in Julia, in which we implement and extend the data-driven Mori-Zwanzig method for extracting large-scale spatio-temporal structures from data, which we call MZMD. This method is an extension of Dynamic Mode Decomposition (DMD) in which Mori-Zwanzig memory kernels are included into the associated companion matrix. In the code we would like to release, we apply MZMD to a flow over a cylinder with Reynolds number 100 rather than the much larger data set used in the associated manuscript. DMD is used extensively in the fluid dynamics community mainly for extracting large scale spatio-temporal structures (patters) from flow data. This is useful for understanding the key mechanisms that generate certain complex dynamical process relevant in engineering design. In MZMD, we improve upon DMD by adding the Mori-Zwanzig memory kernels, and show this improvement is especially important in strongly nonlinear regions of the flow.

Woodward, Michael↗

DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Covariant bit threads

We derive several new reformulations of the Hubeny-Rangamani-Takayanagi covariant holographic entanglement entropy formula. These include: (1) a minimax formula, which involves finding a maximal-area achronal surface on a timelike hypersurface homologous to D(A) (the boundary causal domain of the region A whose entropy we are calculating) and minimizing over the hypersurface; (2) a max V-flow formula, in which we maximize the flux through D(A) of a divergenceless bulk 1-form V subject to an upper bound on its norm that is non-local in time; and (3) a min U-flow formula, in which we minimize the flux over a bulk Cauchy slice of a divergenceless timelike 1-form U subject to a lower bound on its norm that is non-local in space. The two flow formulas define convex programs and are related to each other by Lagrange duality. For each program, the optimal configurations dynamically find the HRT surface and the entanglement wedges of A and its complement. The V-flow formula is the covariant version of the Freedman-Headrick bit thread reformulation of the Ryu-Takayanagi formula. We also introduce a measure-theoretic concept of a “thread distribution”, and explain how Riemannian flows, V-flows, and U-flows can be expressed in terms of thread distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Surfactant-Aided Low-Salinity Waterflooding for Low-Temperature Carbonate Reservoirs

Carbonate reservoirs tend to be oil-wet/mixed-wet and heterogeneous because of mineralogy and diagenesis. The objective of this study is to improve oil recovery in low-temperature dolomite reservoirs using low-salinity and surfactant-aided spontaneous imbibition. The low-salinity brine composition was optimized using ζ-potential measurements, contact-angle (CA) experiments, and a novel wettability-alteration measure. Significant wettability alteration was observed on dolomite rocks at a salinity of 2,500 ppm. We evaluated 37 surfactants by performing CA, interfacial-tension (IFT), and spontaneous-imbibition experiments. Here, three (quaternary ammonium) cationic and one (sulfonate) anionic surfactants showed significant wettability alteration and produced 43–63% of original oil in place (OOIP) by spontaneous imbibition. At a low temperature (35°C), oil recovery by low-salinity effect is small compared with that by wettability-altering surfactants. Coreflood tests were performed with a selected low-salinity cationic surfactant solution. A novel coreflood was proposed that modeled heterogeneity and dynamic imbibition into low-permeability regions. The results of the “heterogeneous” coreflood were consistent with that of spontaneous-imbibition tests. These experiments demonstrated that a combination of low-salinity brine and surfactants can make originally oil-wet dolomite rocks more water-wet and improve oil recovery from regions bypassed by waterflood at a low temperature of 35°C.

02 PETROLEUM↗

DOE OSTI2023

Effect of heat treatment paths on the aging and rejuvenation of metallic glasses

Here, this study explores the influence of heat treatment paths on the structural relaxation of metallic glasses (MGs) in the regime of fast dynamics. We create MG samples using various quenching rates, from 10 9 to 10 11 K/s, and expose these samples to near-T g heat treatments with assorted combinations of heating/cooling rates, from 5×10 10 to 5×10 12 K/s, and annealing durations, from 10 ps to 10 ns. Results show that the effect of the heating rate is intricately tied to the initial structure of the MG, while a decrease in the cooling rate invariably intensifies the aging process. Extending the annealing duration may induce either aging or rejuvenation, subject to the specific thermal history. To interpret these findings, we hypothesize that the memory effect, governed by the activation and annihilation of reversible and irreversible flow units, underpins the rejuvenation-aging competition in MGs. Moreover, we suggest viewing the cooling stage as an annealing phase controlled by gradient temperature and fine-tuning its rate to achieve the targeted active flow unit distributions. This study illuminates the role of fast dynamics during the glass relaxation process and offers practical strategies for tailoring heat treatments to optimize MG structures and mechanical performance.

36 MATERIALS SCIENCE↗

DOE OSTI2022

On the Viability of Quantum Annealers to Solve Fluid Flows

This paper explores the suitability of upcoming novel computing technologies, particularly adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. For our experiments, we start with a well-studied one-dimensional simple flow problem, and provide a framework to convert such problems in continuum to a form amenable for deployment on such quantum annealers. Since the DWave annealer returns multiple states sampling the energy landscape of the problem, we explore multiple solution selection strategies to approximate the solution of the problem. We analyze the continuum solutions obtained both qualitatively and quantitatively as well as their sensitivities to the particular solution selection scheme.

42 ENGINEERING↗

DOE OSTI2021

Distribution of blackouts in the power grid and the Motter and Lai model

Carreras, Dobson, and colleagues have studied empirical data on the sizes of the blackouts in real grids and modeled them with computer simulations using the direct current approximation. They have found that the resulting blackout sizes are distributed as a power law and suggested that this is because the grids are driven to the self-organized critical state. In contrast, more recent studies found that the distribution of cascades is bimodal resulting in either a very small blackout or a very large blackout, engulfing a finite fraction of the system. Here we reconcile the two approaches and investigate how the distribution of the blackouts changes with model parameters, including the tolerance criteria and the dynamic rules of failure of the overloaded lines during the cascade. Finally, we study the same problem for the Motter and Lai model and find similar results, suggesting that the physical laws of flow on the network are not as important as network topology, overload conditions, and dynamic rules of failure.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Novel Tube Design for Superheater Heat Exchanger Enabled Via Additive Manufacturing

Superheater tubes are critical boiler components that operate at relatively higher temperatures and pressure. Amongst the primary concerns for these tubes is the deposition of ash particles on the tube surface, leading to the reduced thickness of the tube due to material corrosion, consequently causing early creep failure of the component. In this research, a novel tube design has been proposed which resembles a teardrop or ogive shape to reduce the drag and concurrently improve the creep life of the superheater tubes. To administer the practicality of novel tubes, metal additive manufacturing (AM), for instance, laser-powder bed fusion (L-PBF), has been proposed. These unconventional designs were assessed and compared with the baseline circular tube design for mechanical design requirements (hoop stress and creep life) and the particle and flue gas flow characteristics around the differently shaped tubes. A thermomechanical finite element (FE) analysis was performed for hoop stress calculations. This study also emphasizes on effect of circumferential thermal variation on hoop stress distribution in tubes. Therefore, a detailed two-dimensional (2D) thermal simulation has been performed to report the circumferential thermal variation on the tube. A computational fluid dynamics (CFD) analysis coupled with particle tracing was performed for gas flow visualization and particle tracing around the proposed shapes and baseline circular-shaped tube design. The Schlieren optic setup was built and leveraged for qualitative validation of the proposed design. The complete design methodology established in the paper shows teardrop-shaped tubes better in terms of drag and creep life in contrast to the circular-shaped tube.

42 ENGINEERING↗

DOE OSTI2020

Computationally-Aided Design of a Small-Scale Radioactive Waste Glass Melter

To provide mission support to the Hanford Waste Immobilization and Treatment Plant for the vitrification of legacy nuclear tank waste, a reduced-scale vitrification pilot system is being designed to process simulated and actual radioactive tank wastes. The tank waste will be separated into high-level waste (HLW) and low-activity waste (LAW) fractions, where the majority by mass (~90%) is LAW and by activity (>95%) is HLW. The first tank waste to be processed at the WTP will be LAW. Since Hanford tank waste and the resultant melter feed compositions are known to vary widely, pilot-scale operations are essential to identify potential problems, provide needed data, confirm assumptions, and determine impacts to full-scale melters and off-gas systems from these different feeds. In addition, the reduced-scale melter system must be capable of quickly providing results to operations without the typically high costs of radioactive operations and provide an engineering platform in close proximity to local operations staff. The objective is to produce similar process conditions to those encountered in the full-scale LAW melters and minimize the volume of radioactive waste necessary for the evaluations. To this end, melter design features are being explored to increase production without significantly increasing surface area, melter glass volumes or compromising data quality. In support of these objectives, a set of computational fluid dynamic models have been developed to provide insight into the flow patterns within these small melters and evaluate design features to increase production without significantly increasing glass inventory. Computed velocities were evaluated at the interface between the cold cap and the molten glass pool and within the bulk glass pool. A method was developed to correlate results to an estimated melt rate and down-select design features that produce the highest gains.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Development and Validation of a Two-Phase Thermal-Hydraulic CFD Code NEK-2P

A project is underway to develop, verify and validate an advanced two-phase flow modeling capability for the highly-scalable, high-performance Computational Fluid Dynamics (CFD) code NEK5000. The goal of this work is to verify and validate the two-phase version of the NEK5000 code, named NEK-2P, to simulate the two-phase flow and heat transfer phenomena that occur in a Boiling Water Reactor (BWR) fuel bundle under various operating conditions. The NEK-2P two-phase flow models follow the approach used for the Extended Boiling Framework (EBF) previously developed at Argonne but include more fundamental physical models of boiling phenomena and advanced numerical algorithms for improved computational accuracy, robustness, and computational speed. The development of the NEK-2P two-phase solver and the implementation of the Extended Boiling Framework two-phase models were initially supported by Argonne National Laboratory (Argonne) through a Laboratory Directed Research and Development (LDRD) project during FY2014-2016. The development and validation of the two-phase models through analyses of selected two-phase boiling flow experiments was supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY2017-2020. This report focuses on verification and validation of the water-steam boiling model NEK-2P Two-Phase, CFD code. The NEK-2P was validated with Nuclear Power Engineering Corporation (NUPEC) Pressurized Water Reactor (PWR) Sub-channel and Bundle Test (PSBT) void distribution benchmark. Three different simulations were performed and analyzed for various operating conditions such as wall-heat flux and sub-cooled inlet temperatures. Reasonably good agreement with measured data was obtained in predicting the measured void distributions. Simulations were performed for Virginia Tech. (VT) 3x3 rod bundle geometry with and without spacers. The preliminary results were presented for Simplified Spacer Grid (SSG). In addition, the implementation of interface reconstruction model was tested with one of the Becker benchmark Critical Heat Flux (CHF) experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Computationally-Aided Design of a Small-Scale Radioactive Waste Glass Melter

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗