Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluids”

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 595 records · Page 33

Seismicity and fault aseismic deformation caused by fluid injection in decametric in-situ experiments

Seismicity induced by fluid perturbations became an important societal concern since felt earthquakes ( M w up to 6) occurred after anthropogenic activities. In order to mitigate the risks associated with undesired seismicity, as well as to be able to use the micro-seismicity as a probe for in-depth investigation of fluid-driven processes, it is of crucial importance to understand the links between seismicity, fluid pressure and flow. We have developed a series of in-situ, decameter-scale experiments of fault zone reactivation by controlled fluid injection, in order to improve the near-source geophysical and hydromechanical observations. The deployed geophysical monitoring close to the injection allows one to cover the full frequency range of the fault responses from the static deformation to the very high-frequency seismic emissions (up to 4 kHz). Here, we focus on the microseismicity ( M w ~ –4 to –3) recorded during two fluid injection experiments in low-permeable shale and highly-fractured limestone formations. In both experiments, the spatio-temporal distribution of the seismic events, the energy balance, and the seismic velocity changes of the fractured medium show that most of the deformation does not actually emit seismic signals. The induced deformation is mainly aseismic. Based on these high-resolution multiparametric observations in the near-field, we therefore proposed a new model for injection-induced seismicity: the seismicity is not directly induced by the increasing fluid pressure, but it is rather triggered by the stress perturbations transferred from the aseismic motion caused by the injection.

58 GEOSCIENCES↗

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗

Experimental procedures to mitigate electron beam induced artifacts during in situ fluid imaging of nanomaterials

Currently there are few standardized experimental practices in the field of fluid stage transmission electron microscopy. To obtain consistency in this emerging field, a crucial step is to establish the common artifacts and electron beam-sample interactions that can occur. Recently many unintended phenomena have been observed during in situ fluid stage scanning transmission electron microscopy (STEM) experiments, including growth of crystals on the fluid stage windows, repulsion of particles from the irradiated area, bubble formation, and the loss of atomic information during prolonged imaging of individual nanoparticles. Here we provide a comprehensive review of these fluid stage artifacts, and we present new experimental evidence that sheds light on their origins in terms of experimental apparatus issues and indirect electron beam sample interactions with the fluid layer. The results here will provide a methodology for minimizing fluid stage imaging artifacts and acquiring quantitative in situ observations of nanomaterial behavior in a liquid environment.

in situ fluid↗

The Effects of Earthquakes and Fluids on the Metamorphism of the Lower Continental Crust

Rock rheology and density have first-order effects on the lithosphere's response to plate tectonic forces at plate boundaries. Changes in these rock properties are controlled by metamorphic transformation processes that are critically dependent on the presence of fluids. At the onset of a continental collision, the lower crust is in most cases dry and strong. However, if exposed to internally produced or externally supplied fluids, the thickened crust will react and be converted into a mechanically weaker lithology by fluid-driven metamorphic reactions. Fluid introduction is often associated with deep crustal earthquakes. Microstructural evidence, suggest that in strong highly stressed rocks, seismic slip may be initiated by brittle deformation and that wall-rock damage caused by dynamic ruptures plays a very important role in allowing fluids to enter into contact with dry and highly reactive lower crustal rocks. The resulting metamorphism produces weaker rocks which subsequently deform by viscous creep. Volumes of weak rocks contained in a highly stressed environment of strong rocks may experience significant excursions toward higher pressure without any associated burial. Slow and highly localized creep processes in a velocity strengthening regime may produce mylonitic shear zones along faults initially characterized by earthquake-generated frictional melting and wall rock damage. However, stress pulses from earthquakes in the shallower brittle regime may kick start new episodes of seismic slip at velocity weakening conditions. These processes indicate that the evolution of the lower crust during continental collisions is controlled by the transient interplay between brittle deformation, fluid-rock interactions, and creep flow.

58 GEOSCIENCES↗

Analytical and computational study of the ideal full two-fluid plasma model and asymptotic approximations for Hall-magnetohydrodynamics

The 5-moment two-fluid plasma model uses Euler equations to describe the ion and electron fluids and Maxwell's equations to describe the electric and magnetic fields. Two-fluid physics becomes significant when the characteristic spatial scales are on the order of the ion skin depth and characteristic time scales are on the order of the ion cyclotron period. The full two-fluid plasma model has disparate characteristic speeds ranging from the ion and electron speeds of sound to the speed of light. Two asymptotic approximations are applied to the full two-fluid plasma to arrive at the Hall-MHD model, namely negligible electron inertia and infinite speed of light. The full two-fluid plasma model and the Hall-MHD model are studied for applications to an electromagnetic plasma shock, geospace environmental modeling (GEM challenge) magnetic reconnection, an axisymmetric Z-pinch, and an axisymmetric field reversed configuration (FRC).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determination of the thermodynamic correction factor of fluids confined in nano-metric slit pores from molecular simulation

The multi-component diffusive mass transport is generally quantified by means of the Maxwell-Stefan diffusion coefficients when using molecular simulations. These coefficients can be related to the Fick diffusion coefficients using the thermodynamic correction factor matrix, which requires to run several simulations to estimate all the elements of the matrix. In a recent work, Schnell et al. [“Thermodynamics of small systems embedded in a reservoir: A detailed analysis of finite size effects,” Mol. Phys. 110, 1069–1079 (2012)] developed an approach to determine the full matrix of thermodynamic factors from a single simulation in bulk. This approach relies on finite size effects of small systems on the density fluctuations. We present here an extension of their work for inhomogeneous Lennard Jones fluids confined in slit pores. We first verified this extension by cross validating the results obtained from this approach with the results obtained from the simulated adsorption isotherms, which allows to determine the thermodynamic factor in porous medium. We then studied the effects of the pore width (from 1 to 15 molecular sizes), of the solid-fluid interaction potential (Lennard Jones 9-3, hard wall potential) and of the reduced fluid density (from 0.1 to 0.7 at a reduced temperature T* = 2) on the thermodynamic factor. The deviation of the thermodynamic factor compared to its equivalent bulk value decreases when increasing the pore width and becomes insignificant for reduced pore width above 15. We also found that the thermodynamic factor is sensitive to the magnitude of the fluid-fluid and solid-fluid interactions, which softens or exacerbates the density fluctuations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cosmology with moving bimetric fluids

We study cosmological implications of bigravity and massive gravity solutions with non-simultaneously diagonal metrics by considering the generalized Gordon and Kerr-Schild ansatzes. The scenario that we obtain is equivalent to that of General Relativity with additional non-comoving perfect fluids. We show that the most general ghost-free bimetric theory generates three kinds of effective fluids whose equations of state are fixed by a function of the ansatz. Different choices of such function allow to reproduce the behaviour of different dark fluids. In particular, the Gordon ansatz is suitable for the description of various kinds of slowly-moving fluids, whereas the Kerr-Schild one is shown to describe a null dark energy component. The motion of those dark fluids with respect to the CMB is shown to generate, in turn, a relative motion of baryonic matter with respect to radition which contributes to the CMB anisotropies. CMB dipole observations are able to set stringent limits on the dark sector described by the effective bimetric fluid.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep Eutectic Salt Formulations Suitable as Advanced Heat Transfer Fluids

Concentrating solar power (CSP) facilities are comprised of many miles of fluid-filled pipes arranged in large grids with reflective mirrors used to capture radiation from the sun. Solar radiation heats the fluid which is used to produce steam necessary to power large electricity generation turbines. Currently, organic, oil-based fluid in the pipes has a maximum temperature threshold of 400 °C, allowing for the production of electricity at approximately 15 cents per kilowatt hour. The DOE hopes to foster the development of an advanced heat transfer fluid that can operate within higher temperature ranges. The new heat transfer fluid, when used with other advanced technologies, could significantly decrease solar electricity cost. Lower costs would make solar thermal electricity competitive with gas and coal and would offer a clean, renewable source of energy. Molten salts exhibit many desirable heat transfer qualities within the range of the project objectives. Halotechnics developed advanced heat transfer fluids (HTFs) for application in solar thermal power generation. This project focused on complex mixtures of inorganic salts that exhibited a high thermal stability, a low melting point, and other favorable characteristics. A high-throughput combinatorial research and development program was conducted in order to achieve the project objective. Over 19,000 candidate formulations were screened. The workflow developed to screen various chemical systems to discover salt formulations led to mixtures suitable for use as HTFs in both parabolic trough and heliostat CSP plants. Furthermore, salt mixtures which will not interfere with fertilizer based nitrates were discovered. In addition for use in CSP, the discovered salt mixtures can be applied to electricity storage, heat treatment of alloys and other industrial processes.

25 ENERGY STORAGE↗

Development of an Environmentally Acceptable Hydraulic Fluid

Hydroelectric plants operators like the United States Army Corps of Engineers seek their traditional petroleum-based fluids with environmentally acceptable hydraulic fluids (EAHF) to reduce their impact and liability related to spills and leakages to the water and environment. Artesion, Inc. and its sub-contractor Washington State University (WSU) sought to optimize WSU’s existing EAHF blend into an improved product that would demonstrate long-term stability while meeting the EAHF requirements of the EPA’s Vessel General Permit (VGP). Multiple rounds of Design of Experiments testing were performed from WSU’s base formulation, modifying mixing ratios to generate multiple test fluids. These fluids were then tested for hydraulic fluid properties. It was found that there was a significant tradeoff between oxidative stability and biodegradability; highly stable materials were not sufficiently biodegradable, while the highly biodegradable materials quickly failed in stability tests. After four rounds of Design of Experiments, AI and WSU had found a base fluid which met the VGP, but only had a Rotating Pressure Vessel Oxidative Test (RPVOT) of 284-minute RPVOT stability, far short of the 500-minute target that would be needed for hydroelectric plants. It was determined that WSU’s route had shown initial promise but needed more refinement to be able to produce an acceptable EAHF.

36 MATERIALS SCIENCE↗

Systematic study on the phase transition of confined fluid mixture up to the critical region

Systematic experimental study on the phase transition of confined fluids up to the critical region has been performed, in which many new insights have been obtained. For confined pure fluids, compared to the bulk phase, the critical point shifts to lower pressure and temperature, but the relative shift of critical pressure is much larger than that of critical temperature. The behavior of the shift, i.e., the shift as a function of pore size, is dependent on the pore type. For confined fluid mixtures, compared to the bulk phase, the critical point also shifts to lower pressure and temperature, and the relative shift of critical pressure is also much larger than that of critical temperature. Confined fluid mixture behaves similarly to confined pure fluid, i.e., it exhibits hysteresis at lower temperature range and does not exhibit phase coexistence or phase envelope. An equation of state for confined pure fluids, based on the Generalized van der Waals partition function, has also been developed and demonstrated to accurately predict the capillary condensation of simple molecules in MCM-41. The parameters of the equation of state are derived from experimental critical point.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Effect of fluid ingestion on orthostatic responses following acute exercise

Orthostatic tolerance is impaired following an acute bout of exercise. This study examined the effect of fluid ingestion following treadmill exercise in restoring the cardiovascular responses to an orthostatic stress. Five men (age, 29.6 +/- 3.4 yrs) were exposed to a graded lower body negative (LBNP) pressure protocol (0 to -50 mmHg) during euhydration without exercise (C), 20 minutes after exercise dehydration (D), 20 minutes after exercise and fluid ingestion (FI20), and 60 minutes after exercise and fluid ingestion (FI60). Fluid ingestion (mean +/- SE) consisted of water-ingestion equivalent to 50% of the body weight lost during exercise (520 +/- 15 ml). Exercise dehydration resulted in significantly higher heart rates (119 +/- 8 vs 82 +/- 7 bpm), lower systolic blood pressures (95 +/- 1.7 vs 108 +/- 2.3 mmHg), a smaller increase in leg circumference (3.7 +/- 4 vs 6.9 +/- 1.0 mm), and an attenuated increase in total peripheral resistance (2.58 +/- 1.2 vs 4.28 +/- 0.9 mmHg/L/min) at -50 mmHg LBNP compared to the C condition. Fluid ingestion (both 20 and 60), partially restored the heart rate, systolic blood pressure, and total peripheral resistance responses to LBNP, but did not influence the change in leg circumference during LBNP (4 +/- 0.3 for R20 and 2.8 +/- 0.4 mm for R60). These data illustrate the effectiveness of fluid ingestion on improving orthostatic responses following exercise, and suggest that dehydration is a contributing factor to orthostatic intolerance following exercise.

NASA Center JSC↗

NASA Fluid Lensing & MiDAR: Next-Generation Remote Sensing Technologies for Aquatic Remote Sensing

We present two recent instrument technology developments at NASA, Fluid Lensing and MiDAR, and their application to remote sensing of Earth's aquatic systems. Fluid Lensing is the first remote sensing technology capable of imaging through ocean waves in 3D at sub-cm resolutions. MiDAR is a next-generation active hyperspectral remote sensing and optical communications instrument capable of active fluid lensing. Fluid Lensing has been used to provide 3D multispectral imagery of shallow marine systems from unmanned aerial vehicles (UAVs, or drones), including coral reefs in American Samoa and stromatolite reefs in Hamelin Pool, Western Australia. MiDAR is being deployed on aircraft and underwater remotely operated vehicles (ROVs) to enable a new method for remote sensing of living and nonliving structures in extreme environments. MiDAR images targets with high-intensity narrowband structured optical radiation to measure an objectâ€"TM"s non-linear spectral reflectance, image through fluid interfaces such as ocean waves with active fluid lensing, and simultaneously transmit high-bandwidth data. As an active instrument, MiDAR is capable of remotely sensing reflectance at the centimeter (cm) spatial scale with a signal-to-noise ratio (SNR) multiple orders of magnitude higher than passive airborne and spaceborne remote sensing systems with significantly reduced integration time. This allows for rapid video-frame-rate hyperspectral sensing into the far ultraviolet and VNIR wavelengths. Previously, MiDAR was developed into a TRL 2 laboratory instrument capable of imaging in thirty-two narrowband channels across the VNIR spectrum (400-950nm). Recently, MiDAR UV was raised to TRL4 and expanded to include five ultraviolet bands from 280-400nm, permitting UV remote sensing capabilities in UV A, B, and C bands and enabling mineral identification and stimulated fluorescence measurements of organic proteins and compounds, such as green fluorescent proteins in terrestrial and aquatic organics.

fluid lensing↗

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗

THE POTENTIAL OF NANOPARTICLE ENHANCED IONIC LIQUIDS (NEILS) AS ADVANCED HEAT TRANSFER FLUIDS

Interest in capturing the energy of the sun is rising as demands for renewable energy sources increase. One area of developing research is the use of concentrating solar power (CSP), where the solar energy is concentrated by using mirrors to direct the sunlight towards a collector filled with a heat transfer fluid (HTF). The HTF transfers the collected energy into pressurized steam, which is used to generate energy. The greater the energy collected by the HTF, the more efficent the electrical energy production is, thus the overall efficiency is controlled by the thermal fluid. Commercial HTFs such as Therminol{reg_sign} (VP-1), which is a blend of biphenyl and diphenyl oxide, have a significant vapor pressure, especially at elevated temperatures. In order for these volatile compounds to be used in CSP systems, the system either has to be engineered to prevent the phase change (i.e., volatilization and condensation) through pressurization of the system, or operate across the phase change. Over thirty years ago, a class of low-melting organic compounds were developed with negligible vapor pressure. These compounds are referred to as ionic liquids (ILs), which are organic-based compounds with discrete charges that cause a significant decrease in their vapor pressure. As a class, ILs are molten salts with a melting point below 100 C and can have a liquidus range approaching 400 C, and in several cases freezing points being below 0 C. Due to the lack of an appreciable vapor pressure, volatilization of an IL is not possible at atmospheric pressure, which would lead to a simplification of the design if used as a thermal fluid and for energy storage materials. Though the lack of a vapor pressure does not make the use of ILs a better HTF, the lack of a vapor pressure is a compliment to their higher heat capacity, higher volummetric density, and thus higher volumetric heat capacity. These favorable physical properties give ILs a pontential advantage over the current commerically used thermal fluids. Also within the past decade nanofluids have gained attention for thermal conductivity enhancment of fluids, but little analysis has been completed on the heat capacity effects of the nanoparticle addition. The idea of ILs or nanofluids as a HTF is not new, as there are several references that have proposed the idea. However, the use of ionic liquid nanofluids containing nanomaterials other than carbon nanotubes has never before been studied. Here, for the first time, nano-particle enhanced ILs (NEILs) have been shown to increase the heat capacity of the IL with no adverse side effects to the ILs thermal stability and, only at high nanoparticle loading, are the IL physical properties affected. An increase of volumetric heat capacity translates into a better heat transfer fluid as more energy is stored per volumetric unit in the solar concentrating section, thus more efficency in increased steam pressure. Results show that the properties of the NEIL are highly dependant on the suspended nanomaterial and careful materials selection is required to fully optimize the nanofluid properties.

14 SOLAR ENERGY↗

Fluid-Structure Interaction Modeling of High-Aspect Ratio Nuclear Fuel Plates Using COMSOL

The High Flux Isotope Reactor at the Oak Ridge National Lab is in the research stage of converting its fuel from high-enriched uranium to low-enriched uranium. Due to different physical properties of the new fuel and changes to the internal fuel plate design, the current safety basis must be re-evaluated through rigorous computational analyses. One of the areas being explored is the fluid-structure interaction phenomenon due to the interaction of thin fuel plates (50 mils thickness) and the cooling fluid (water). Detailed computational fluid dynamics and fluid-structure interaction simulations have only recently become feasible due to improved numerical algorithms and advancements in computing technology. For many reasons including the already built-in fluid-structure interaction module, COMSOL has been chosen for this complex problem. COMSOL's ability to solve multiphysics problems using a fully-coupled and implicit solution algorithm is crucial in obtaining a stable and accurate solution. Our initial findings show that COMSOL can accurately model such problems due to its ability to closely couple the fluid dynamics and the structural dynamics problems.

Fluid-Structure Interaction↗

A Study of Physics-Informed Deep Learning for System Fluid Dynamics Closures

The sub-grid-scale (SGS) physics models, or so-called closure relations (CR), are essential in thermal-fluid modeling and simulation codes, ranging from large-eddy simulation to system-level thermo-fluid dynamic problems. While details vary, all thermo-fluid codes solve coupled mass-momentum-energy conservation partial differential equations (PDE). Traditionally, the analytical form of closure models has to be known, and then we calibrate the model using data from a set of applicable experiments. Characteristically, it takes long effort to gain insights and mechanistic understanding needed for developing closure models. Such a model is a compact form of data. The long time required for new model development limits the pace of model application while dealing with new geometries, new coolants or system conditions, particularly in newly designed systems. Alternative to mechanistic and semi-analytical models, it is possible to discover the pattern of closure relations by applying the modern machine learning methods to the body of relevant experimental data. Deep learning (DL) is a relatively new branch in the machine learning family, and a major breakthrough has been achieved since Hinton introduces a fast algorithm to efficiently train a feed-forward neural network (NN) with multilayer perceptrons. Generally, any NN with more than two layers is considered as deep learning, and it contains a large number of hyper-parameters to achieve pattern recognition or regression for complex datasets. Hornik shows that multilayer NN is compatible with universal approximation theorem. There is no constraint for a multilayer NN to capture the properties of any measurable information. This capability makes DL attractive for closure models in a system thermal-hydraulic code. As system geometry and/or operating conditions change, new physical experiments (increasingly sophisticated) are conducted to support evaluation and revision of closure models. Such a so-called data-driven modeling relaxes the structural form of system codes to directly assimilate data instead of hard-wiring the model to pre-calibrated correlations. The objective of this research is to enable cost-effective development of a DL-based closure that reflects the data. The previous work on applying machine learning techniques to fluid dynamics focused on statistical learning with basic NN from direct numerical simulations (DNS) dataset. In the present work, we show that DL-based closure models can cause the fluid calculation to diverge if multilayer NN is used. To address the issue, we investigate conditions by which the DL-based closure models work compatibly, stably, and effectively with PDE-constrained forward prediction problems. Optimally meeting the compatibility conditions constitutes a physics-informed deep learning (PIDL) strategy. We refer to this process as PDE-constrained machine learning or PDE-ML. To develop the technical basis for PDE-ML using DL-based fluid dynamics closure, numerical experiments are performed for manufactured problems of increasing complexity. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tabulated Fluid Properties Research Report

The Multiphysics Object-Oriented Simulation Environment (MOOSE) enables a wide range of advanced nuclear reactor simulations.[6] Under the guidance of MOOSE’s Thermal Hydraulics Team,I worked to expand the capabilities of Tabulated Fluid Properties (TFP) in the fluid properties module. The fluid properties module allows the user to determine a variety of fluid properties by interpolating points between tabulated data. I implemented the ability to use bilinear interpolation instead of bicubic interpolation for interpolating tabulated data. I also changed the method of variable set inversions to use a 2-dimensional Newton’s Method utility that I created. Variable set inversions are often done from (v,e) to (p,T), where v is specific volume, e is specific internal energy, p is pressure and T is temperature. New routines have also been added into TFP such that it can be used with more applications, such as the Navier Stokes and Thermal Hydraulics modules in MOOSE for Pronghorn[5] and RELAP-7[1] respectively. This work was spurred by interest from NASA in testing a Nuclear Thermal Propulsion (NTP) engine system. NTP engines have drastically different fluid properties throughout the engine and Tabulated Fluid Properties provides the flexibility needed to properly simulate and test these engines.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An adaptive lattice Boltzmann scheme for modeling two-fluid-phase flow in porous medium systems

In this paper, we formulate a multiple-relaxation-time (MRT) lattice-Boltzmann method (LBM) to simulate two-fluid-phase flow in porous medium systems. The MRT LBM is applied to simulate the displacement of a wetting fluid by a nonwetting fluid in a system corresponding to a microfluidic cell. Analysis of the simulation shows widely varying time scales for the dynamics of fluid pressures, fluid saturations, and interfacial curvatures that are typical characteristics of such systems. Displacement phenomena include Haines jumps, which are relatively short duration isolated events of rapid fluid displacement driven by capillary instability. An adaptive algorithm is advanced using a level-set method to locate interfaces and estimate their rate of advancement. Because the displacement dynamics are confined to the interfacial regions for a majority of the relaxation time, the computational effort is focused on these regions. The proposed algorithm is shown to reduce computational effort by an order of magnitude, while yielding essentially identical solutions to a conventional fully coupled approach. Finally, the challenges posed by Haines jumps are also resolved by the adaptive algorithm. Possible extensions to the advanced method are discussed.

42 ENGINEERING↗