Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porous 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.

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

At least 91 records · Page 5

Flow through very porous inclined screens

The steady, inviscid flow through and around a screen inclined at a uniform angle to the incoming flow was investigated. For a screen placed in an infinite flow field, an asymptotic analysis for small resistance coefficients was performed, and the effects of inclination were determined. The velocity at first order in the asymptotic expansion was nonuniform along the screen. This nonuniformity caused the wake behind the screen to contain distributed vorticity at second order. These effects therefore occurred at one order lower than for normal screens.

Muramoto, K. K.↗

Flow through very porous inclined screens

The steady, inviscid flow through and around a screen inclined at a uniform angle to the incoming flow was investigated. For a screen placed in an infinite flow field, an asymptotic analysis for small resistance coefficients was performed, and the effects of inclination were determined. The velocity at first order in the asymptotic expansion was nonuniform along the screen. This nonuniformity caused the wake behind the screen to contain distributed vorticity at second order. These effects therefore occurred at one order lower than for normal screens.

Muramoto, K. K.↗

Flow friction of the turbulent coolant flow in cryogenic porous cables

Considered are cryogenic power transmission cables with porous cores. Calculations of the turbulent coolant flow with injection or suction through the porous wall are presented within the framework of a two-layer model. Universal velocity profiles were obtained for the viscous sublayer and flow core. Integrating the velocity profile, the law of flow friction in the pipe with injection has been derived for the case when there is a tangential injection velocity component. The effect of tangential velocity on the relative law of flow friction is analyzed. The applicability of the Prandtl model to the problem under study is discussed. It is shown that the error due to the acceptance of the model increases with the injection parameter and at lower Reynolds numbers; under these circumstances, the influence of convective terms in the turbulent energy equation on the mechanism of turbulent transport should be taken into account.

Hendricks, R. C.↗

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun porous flow equation was applied to the brush seal data. The Ergun relation, which required heuristic information to characterize the coefficients, fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the (C(sub f)A x Constant) flow and area coefficients. Reasonable agreement between prediction and data was found, as expected, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

Preliminary experimental results for a cryogenic brush seal configuration

Preliminary fluid nitrogen flow data are reported for a five-brush, ceramic-coated-rub-runner brush seal system, where the brushes and the rub runner were placed at each end of a centrally pressurized multifunction tester ('back-to-back' set of brushes) and tested at rotor speeds of 0, 10, 18, and 22.5 krpm. After testing, both the brushes and the ceramic-coated rub runner appeared pristine. The coating withstood both the thermomechanical and dynamic loadings with minor wear track scarring. The bristle tips showed some indication of material shearing (smearing) wear. The Ergun (1952) porous flow equation was applied to the brush seal data. The Ergun relation fit the gaseous data but was in poor agreement with the fluid results. The brush seal exit conditions were two phase. Two-phase, choked-flow design charts were applied but required one data point at each rotor speed to define the flow and area coefficients. Reasonable agreement between prediction and data was found, but such methods are not to be construed as two-phase-flow brush seal analyses.

Carlile, J. A.↗

Multi-fidelity electrochemical modeling of thermally activated battery cells

Thermally activated batteries undergo a series of coupled physical changes during activation that influence battery performance. These processes include energetic material burning, heat transfer, electrolyte phase change, capillary-driven two-phase porous flow, ion transport, electrochemical reactions, and electrical transport. Several of these processes are strongly coupled and have a significant effect on battery performance, but others have minimal impact or may be suitably represented by reduced-order models. Additionally, assessing the relative importance of these phenomena must be based on comparisons to a high-fidelity model including all known processes. In this work, we first present and demonstrate a high-fidelity, multi-physics model of electrochemical performance. This novel multi-physics model enables predictions of how competing physical processes affect battery performance and provides unique insights into the difficult-to-measure processes that happen during battery activation. We introduce four categories of model fidelity that include different physical simplifications, assumptions, and reduced-order models to decouple or remove costly elements of the simulation. Using this approach, we show an order-of-magnitude reduction in computational cost while preserving all design-relevant quantities of interest within 5 percent. The validity of this approach and these model reductions is demonstrated by comparison between results from the full fidelity model and the different reduced models.

25 ENERGY STORAGE↗

Materials Engineering for High Performance and Durability Proton Exchange Membrane Water Electrolyzers

Proton exchange membrane water electrolyzers (PEMWEs) are expected to play a crucial role in the global green energy transition during the 21st century. They provide a versatile and sustainable solution for generating hydrogen with very high purity in combination with renewable energies, such as solar and wind. Despite their promise, PEMWEs face several critical problems, including high costs, performance limitations, and durability challenges, particularly at low iridium (Ir) loading on the anode. Advancing next-generation PEMWEs requires extensive work on materials engineering of all cell components, including the catalyst layer (CL), membrane, porous transport layer (PTL), bipolar plate (BPP), and gasket. This task must be performed with the complementary contribution of different modeling and characterization techniques. This review presents a critical perspective from academia, research centers, and industry, mapping main developments, remaining gaps, and strategic pathways to advance PEMWE technology. A focus is devoted to key aspects, such as operation at low Ir loading, membrane durability, multiscale transport layers, porous and non-porous flow fields, multiphysics modeling, and multipurpose characterization techniques, which are thoroughly discussed. By unifying these topics, this review provides readers with the essential knowledge to grasp current developments and tackle tomorrow's challenges in PEMWE engineering.

36 MATERIALS SCIENCE↗

Simplifying the creation of thermal decomposition mechanisms for designing mixed fibre composites

Among the most challenging aspects of simulating thermal decomposition of fibre reinforced polymers is the determination of appropriate reaction parameters. Thermogravimetric analysis is typically used to generate decomposition data, to which the reaction parameters are then fit. When designing mixed fibre materials (e.g. combinations of glass and carbon fibres), the number of TGA experiments needed to explore the entire design space may be intractable. Here, we demonstrate the creation of two candidate proxy mixed fibre mechanisms and compare them to a mechanism created by fitting parameters from TGA on the mixed fibre composite. These mechanisms are then demonstrated in a 2D axisymmetric numerical decomposition, heat transfer, and porous flow model. We find a maximum 11% uncertainty in mass and 4% in temperature difference when using a proxy mechanism.

Scott, Sarah N. [Sandia National Lab. (SNL-CA), Li↗

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗

HI-STORM Overpack and MPC-32 Thermal-Hydraulic Model with MOOSE Framework

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a MOOSE thermal model of the MPC-32 canister and HI-STORM overpack

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. The safe management of the spent nuclear fuel (SNF) is a key aspect of the back-end of the nuclear fuel cycle. Spent fuel dry storage systems are becoming a popular and effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution to maintain spent fuel for 60 years before final disposal. Dry cask storage has many characteristics that make it attractive. It fulfills the safety requirements of the Nuclear Regulatory Commission (NRC) while providing modularity and flexibility to contractors. The HI-STORM overpack and MPC-32 canister are the main parts of the HI-STORM 100 dry cask storage system. These components remove heat from the system using natural circulation, requiring no human intervention. This is the characteristic that provides passive heat removal and low maintenance features in dry cask storage systems. To develop a thermal model for a dry cask storage system, the physics behind the system should be defined clearly. There are two natural circulation loops in the system; circulation of helium cools down the nuclear assemblies in the MPC, while circulation of air cools down the walls of the MPC. This work aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. In this study, we will investigate and demonstrate the thermal-hydraulics modeling capabilities of the MOOSE framework, including natural circulation, heat transfer, and porous flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a MOOSE thermal model of the MPC-32 canister and HI-STORM overpack

Nuclear power is a significant source of electricity in the United States, but the average age of nuclear power plants is around 40 years old. Safe management of spent nuclear fuel (SNF) is a key aspect of the back end of the nuclear fuel cycle, and SNF dry storage systems are becoming a popular, effective solution in this area, given the absence of a final disposal system. The spent fuel cask system (dry cask method) provides a feasible solution for maintaining SNF (~60 years) prior to final disposal. This project aims to develop a thermal model of the MPC-32 canister and HI-STORM overpack, using the Multiphysics Object-Oriented Simulation Environment (MOOSE). MOOSE is an open-source framework developed by Idaho National Laboratory (INL) for multiscale, multiphysics simulations. This study will investigate and demonstrate the thermal-hydraulics capabilities of the MOOSE framework, including natural circulation, heat transfer, porous flows, etc. The ultimate goal of the project is to verify whether MOOSE tools (including Pronghorn) can be used to study the thermal performance of the SNF dry cask storage system. This study provides reliable and inclusive solving strategy for dry cask problems. The detailed information about the solving scheme and the governing equations related to the physics of the system is provided in the report. The results for thermal-hydraulic analysis of the HI-STORM system is produced with using open source modules of the MOOSE framework. This results highlights the flexibility and modularity of the MOOSE which makes it a unique candidate for the frameworks and code packages. Therefore, integration of the MOOSE to UNF ST&DARDS will improve the thermal-hydraulic capability of the system while providing distinctive features to users.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗

Surface characteristics of Venus derived from Pioneer Venus altimetry, roughness, and reflectivity measurements

Pioneer Venus data sets for global topography, surface roughness, and reflectivity are analyzed, and the geological significance of some of the distinctive surface types is locally calibrated by examining high-resolution images from the Venera lander spacecraft and the Arecibo Observatory. Regional rock and bedrock surfaces are found to cover the majority of the planet. The relatively smooth rock and bedrock surfaces probably represent lava flows. Porous and unconsolidated fine materials cover less than about 27 percent of the surface. The distribution of soil surfaces indicate that they are the result of local weathering and small amounts of lateral transport. Several hypotheses for the origin of the high-dielectric materials which cover a small portion of the surface are discussed.

Head, J. W., III↗

Melt migration in a silicate liquid-olivine system - An experimental test of compaction theory

An experimentally derived melt migration profile is compared with those obtained from compaction theory. A couple is prepared in which a source-for-melt (a disk composed of a silicate glass and olivine) was placed in contact with a sink-for-melt (a disk of polycrystalline olivine) in order to induce melt migration under laboratory conditions. It is noted that melt infiltrates into the sink disk along triple junctions driven by capillary forces. In order to analyze the melt migration profile in terms of compaction theory, the equations developed by McKenzie (1984) to describe porous flow of a liquid in a deformable matrix were modified. The bouyancy force term is replaced with a capillary force term and governing equations are then solved numerically with the initial boundary conditions specified by the experimental design. Also, it is determined that compaction theory provides a good description of the experimental results provided that the permeability of these solid-liquid materials increases linearly with increasing liquid fraction.

Riley, G. N., Jr.↗

Trace element and isotopic effects arising from magma migration beneath mid-ocean ridges

The trace element concentrations and isotopic ratios in the magma erupted on mid-ocean ridges may differ from those in the source material due to physical effects such as porous flow dispersion, exchange of trace elements between the fluid and solid phases during magma migration, and convective mixing in magma chambers. These differences are in addition to those produced by better known processes such as fractional crystallization and partial melting. The effects of the three former processes are described. It is predicted that magma typically reaches the subridge magma chambers with a spatial heterogeneity only slightly reduced from that of the source material, but with a subdued variation in time. Convective mixing then further reduces the spatial heterogeneity. Application of the results for convective mixing to a recent Fourier analysis of Sr-87/Sr-86 variations along the Mid-Atlantic Ridge suggests that the falloff in amplitude of variation observed with decreasing wavelength in the Mid-Atlantic Ridge data cannot be explained by convective mixing in magma chambers. Instead, it is postulated that this falloff is due to the mechanics of the production and/or the solid-state convective mixing of chemical and isotopic heterogeneities in the solid mantle.

Kenyon, Patricia M.↗