Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lattice Boltzmann”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidating effects of catalyst loadings and porous transport layer morphologies on operation of proton exchange membrane water electrolyzers

Producing green hydrogen efficiently via proton exchange membrane water electrolysis (PEMWE) is the key for achieving decarbonization targets. Iridium catalyst is expensive, and it is important to minimize its use and to optimize interface between Ir and ionomer or water for higher utilization of catalyst in oxygen evolution reaction. In this work, x-ray computed tomography along with electrochemical and modeling techniques are used to characterize the interface for two different porous transport layers (PTLs) and catalyst layers at various loadings. We show that low porosity sintered PTLs exhibit higher interfacial contact with the catalyst and the membrane that results in improved kinetics. Radiography and modeling results indicate that oxygen taking multiple transport pathways through the PTL results in slug flow through the channels that reduces mass transport overpotential. Based on the results, we suggest design guidelines for high efficiency and durable PEMWE and their components.

08 HYDROGEN↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2023

Effect of constitutive law on the erythrocyte membrane response to large strains

Three constitutive laws, that is the Skalak, neo-Hookean and Yeoh laws, commonly employed for describing the erythrocyte membrane mechanics are theoretically analyzed and numerically investigated to assess their accuracy for capturing erythrocyte deformation characteristics and morphology. In this work, particular emphasis is given to the nonlinear deformation regime, where it is known that the discrepancies between constitutive laws are most prominent. Hence, the experiments of optical tweezers and micropipette aspiration are considered here, for which relationships between the individual shear elastic moduli of the constitutive laws can also be established through analysis of the tension-deformation relationship. All constitutive laws were found to adequately predict the axial and transverse deformations of a red blood cell subjected to stretching with optical tweezers for a constant shear elastic modulus value. As opposed to Skalak law, the neo-Hookean and Yeoh laws replicated the erythrocyte membrane folding, that has been experimentally observed, with the trade-off of sustaining significant area variations. For the micropipette aspiration, the suction pressure-aspiration length relationship could be excellently predicted for a fixed shear elastic modulus value only when Yeoh law was considered. Importantly, the neo-Hookean and Yeoh laws reproduced the membrane wrinkling at suction pressures close to those experimentally measured. None of the constitutive laws suffered from membrane area compressibility in the micropipette aspiration case.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Opening twisted polymer chains for simultaneously high printability and battery fast-charge

Fast-charging is key to the widespread adoption of battery-based electric vehicles. However, improving fast-charging through architecture optimization is expensive. To reduce costs and expand to commercialization, we applied roll-to-roll screen printing technology to create channels and decrease the tortuosity of electrodes. For the first time, this work successfully opened twisted polymer chains within high-solid-content inks to improve their screen printability and battery performance of as-printed electrodes. With LiNi 0.6 Mn 0.2 Co 0.2 O 2 as active materials, the 60% solid content ink presents superior screen printability after opening the twisted binder chains. As-printed electrode exhibits 33% higher charge capacity at 6 C than printed electrode with chains twisted ink at mass loading of 6.5 mg/cm 2 . Furthermore, coarse-grained molecular dynamics simulations are performed to study the underlying mechanism systematically. In conclusion, the new ink preparation procedure provides a scalable, effective strategy for manufacturing screen-printable battery ink and promotes screen-printed electrode technology.

25 ENERGY STORAGE↗
DOE OSTI2020

Controlled 3D nanoparticle deposition by drying of colloidal suspension in designed thin micro-porous architectures

Nanoparticle deposition by drying of colloidal suspension in thin micro-porous architectures has attracted a lot of attention in scientific research as well as industrial applications. However, the underlying mechanisms of such three-dimensional (3D) deposition are not yet fully revealed due to the complexity of the co-occurring processes of two-phase fluid flows, phase change and mass transport. Consequently, the control of 3D nanoparticle deposition remains a challenge. We use a combined experimental and numerical approach to achieve controlled 3D nanoparticle deposition by drying of colloidal suspension in two pillar-based thin micro-porous architectures. By the design of pillar layout, rectangular-spiral and circular-spiral deposition configurations are obtained globally. By varying the surface wettability, vertically symmetric and sloped nanoparticle depositions can be achieved locally. While the numerical modeling reveals the mechanisms of liquid internal flow, as well as the impact of local drying rate on nanoparticle transport, accumulation and final deposition, the experimental results of deposition configurations validate the controlling strategies. This combined experimental and numerical work provides a framework to achieve desired 3D nanoparticle deposition in thin micro-porous architectures, with a thorough understanding of the underlying mechanisms of two-phase fluid flows, phase change and mass transport.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-GPU immersed boundary method hemodynamics simulations

Large-scale simulations of blood flow that resolve the 3D deformation of each comprising cell are increasingly popular owing to algorithmic developments in conjunction with advances in compute capability. Among different approaches for modeling cell-resolved hemodynamics, fluid structure interaction (FSI) algorithms based on the immersed boundary method are frequently employed for coupling separate solvers for the background fluid and the cells within one framework. GPUs can accelerate these simulations; however, both current pre-exascale and future exascale CPU-GPU heterogeneous systems face communication challenges critical to performance and scalability. In this paper, we describe, to our knowledge, the largest distributed GPU-accelerated FSI simulations of high hematocrit cell-resolved flows with over 17 million red blood cells. We compare scaling on a fat node system with six GPUs per node and on a system with a single GPU per node. Through comparison between the CPU- and GPU-based implementations, we identify the costs of data movement in multiscale multi-grid FSI simulations on heterogeneous systems and show it to be the greatest performance bottleneck on the GPU.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Dataset of 3D Structural and Simulated Transport Properties of Complex Porous Media

Physical processes that occur within porous materials have wide-ranging applications including - but not limited to - carbon sequestration, battery technology, membranes, oil and gas, geothermal energy, nuclear waste disposal, water resource management. The equations that describe these physical processes have been studied extensively; however, approximating them numerically requires immense computational resources due to the complex behavior that arises from the geometrically-intricate solid boundary conditions in porous materials. Here, we introduce a new dataset of unprecedented scale and breadth, DRP-372: a catalog of 3D geometries, simulation results, and structural properties of samples hosted on the Digital Rocks Portal. The dataset includes 1736 flow and electrical simulation results on 217 samples, which required more than 500 core years of computation. This data can be used for many purposes, such as constructing empirical models, validating new simulation codes, and developing machine learning algorithms that closely match the extensive purely-physical simulation. This article offers a detailed description of the contents of the dataset including the data collection, simulation schemes, and data validation.

3D images↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2024

Learning a general model of single phase flow in complex 3D porous media

Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Low-G Slosh Workshop Results From 2023: State of the Art, Gaps and Forward Work

The NASA Engineering and Safety Center (NESC) sponsored a workshop to dis-cuss low-g slosh. Over a dozen different NASA centers and industry partners participated in this workshop. The main topics were to identify the state of the art regarding test data, Computational Fluid Dynamics (CFD) analysis resources, and GN&C models that are available for this problem. The primary intent was to determine deficiencies in the collective knowledge base and to identify priorities that should be established for proactive work in the areas of model development, validation, and model integration. Recommendations from workshop participants include establishing an online slosh modeling resource, conducting ground/flight testing, developing surrogate models, and establishing co-simulation techniques. This paper summarizes the workshop findings and recommendations for the larger community.

Lattice-Boltzmann↗
DOE OSTI2025

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗