Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow rate ratio”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗
DOE OSTI2022

Uncertainty Quantification of Monolithic Tungsten Target Block Fatigue

The current STS spallation target design consists of 21 water cooled, tantalum clad, monolithic tungsten blocks. The target wheel assembly rotates such that each target block receives a direct proton pulse every 1.4 seconds (21 blocks / 15 Hz pulses). The quasi-steady heating induced by the repetitive beam pulses creates a significant mean stress distribution in the target. On top of this, the rapid energy deposition from a proton pulse induces dynamic stress waves that are of the same order of magnitude as the mean stresses. Over the design goal of 5,000 operational hours per year for a decade, each target must withstand over 130 million repetitions of these dynamic stress events. The fatigue life prediction of the target blocks is a complex sequence of simulations that depend upon a set of input parameters. For the tungsten block, the end goal is the fatigue factor-of-safety (FOS). This study documents a new structural analysis routine using software available to U.S. national laboratories (e.g. CUBIT, Sierra Multiphysics, Dakota) that enables more rapid target fatigue analysis. The increased throughput is amenable to parameter sensitivity and uncertainty quantification (UQ) studies. A set of input parameters (ultimate strength, fatigue strength, density, elastic modulus, Poisson’s ratio, thermal conductivity, coefficient of thermal expansion, specific heat, flow rate) with expected ranges are documented based on literature data, recent experimental testing, and engineering judgement. UQ bounds the resulting fatigue FOS at [0.5, 3.4] using epistemic interval analysis. These bounds are expected to be wider due to model effects (discretization) and beam effects (focus, offset) that were not included. Parameter sensitivity underscores the importance of ultimate strength and fatigue strength of irradiated tungsten on the fatigue calculations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Effect of Evaporation on Vapor Space Flammability in DWPF Process Vessels - 20278

The Savannah River Site (SRS) Defense Waste Processing Facility (DWPF) processes liquid high-level waste (HLW) and encapsulates it in glass within stainless-steel containers suitable for safe, long-term storage and disposal. Radiolytic and chemical decomposition of the liquid waste generates hydrogen and other flammable gases creating a hazard. The rate of flammable gas generation increases by an order of magnitude as temperature increases from sub-boiling to boiling. DWPF controls this hazard in process tanks by ensuring a minimum purge flow through each tank's head space. The minimum purge flow rate is established assuming dry air and neglecting evaporation. However, this study determined that steam from boiling in a representative process tank is sufficient to dilute flammable gases and, thereby, preclude flammability. For sub-boiling, water vapor is an insufficient diluent considering hydrogen's lower flammable limit (LFL) in air. However, the upper and lower flammable limits approach each other in a hydrogen-air-water-vapor environment as the water vapor concentration increases. This effect may preclude flammable conditions if water vapor concentration is about 60 volume percent or higher. A saturated, constant-pressure air-water system would achieve 60 volume percent water vapor at 86 deg. C. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Uniform Whole Wafer Anisotropic Etching of Structural Ta Thin Films

Tantalum (Ta) has promise as a structural material for micromechanical sensors and actuators. Anisotropic etching of alpha (α) phase Ta is required for micromachining applications. Uniform thickness and etch across a full wafer are desirable features. An experimental investigation is conducted to study plasma etching rate and anisotropy in etching α -Ta in relation to pressure and gas flow ratio. A comparatively low toxicity gas, carbon tetrafluoride (CF4), and argon (Ar) are used. Spectrometry by optical emission spectrometry (OES) and residual gas analysis (RGA) are employed to characterize the plasma to gain insight into the etch mechanisms. At low flow, the etch rate is slow due to an inadequate supply of the etching gas, CF4. The etch rate is also slowed at high flow due to a reduced CF4 residence time. Flow and pressure conditions to achieve a good etch rate and vertical sidewalls are identified and explained by means of a full factorial experiment coupled with emission spectra. Finally, with these conditions, uniform etching of 2.5 μm thick α-Ta across a 4-inch wafer is demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A novel 3-D pulsating heat pipe module for high heat-flux applications

This study proposed a novel pulsating heat pipes (PHP) for high-flux thermal module. According to the open literature, the maximum value of the heat flux is 23.7 W∙cm –2 handled by PHP, which is far behind industry requirements. Therefore, the present study proposes a novel 3D PHP suitable for high heat-flux applications. The PHP were experimentally tested for different working fluids, including methanol and water. Furthermore, the experiments were performed to analyze the influence of various operating parameters, such as filling ratio (25–60%), heat load (100–800 W), and air flow rate 100–250 cubic feet per minute (CFM). In this work, the results showed that this PHP can sustain a heat flux up to 70 W∙cm –2 with the maximum case temperature under 80 °C, which is approximately three times higher than the available literature. Furthermore, the performance of the PHP module was compared against the commercial module with wicked heat pipe. The PHP exhibited about 7% higher heat dissipation compared to the commercial module. Since the proposed module uses simple copper tubes in its construction over the wicked structure used in commercial module, the novel 3D PHP module provides additional economic and mass production capability.

3D pulsating heat pipe↗
DOE OSTI2022

An additively-manufactured molten salt-to-supercritical carbon di-oxide primary heat exchanger for solar thermal power generation – Design and techno-economic performance

The design and techno-economic performance of a compact additively manufactured (AM) molten salt (MS)-to-supercritical carbon di-oxide (sCO2) primary heat exchanger (PHE) for solar thermal application is described. The PHE design consists of sCO2 flow through an array of microscale pin fins while the MS flows through mm-scale rectangular channels. Constraints imposed by AM using laser powder bed fusion method are considered in the design. Structural and fluid flow simulations are performed to arrive at a viable design of the core and headers. A simplified one-dimensional steady state model for the PHE is developed including the impact of surface roughness from the AM process. A process-based cost model is used to determine the tradeoff between thermofluidic design and manufacturing cost. A parametric study is performed using the thermo-fluidic and cost models to determine the set of geometrical and flow variables that result in high power density and low cost, while restricting the pressure drop on the sCO2 side to less than 2% of line pressure. Flow rates of MS and sCO2 were varied over heat capacity rate ratios ranging from 0.2 to 1. Results indicate that it is possible to design a low-pressure drop AM PHE with an effectiveness of 90% and a power density in excess of 10 MW/m3 (including headers). Fabrication of representative nickel superalloy specimens are shown to demonstrate that low-porosity parts with the requisite dimensional tolerance of PHE core can be generated.

14 SOLAR ENERGY↗
DOE OSTI2024

Gaseous contaminant transfer in membrane-based air-to-air energy exchangers

Membrane-based air-to-air energy exchangers (M–AAEEs) transfer heat and moisture between building exhaust air and fresh ventilation air streams through a membrane, thereby reducing the energy required for conditioning the fresh ventilation air. Energy exchangers are typically used in buildings with relatively clean building exhaust air, such as office buildings and schools. Recently interest in using energy exchangers in a wider range of buildings has grown, to reduce the energy consumption associated with the heating, ventilating, and air-conditioning (HVAC) systems in these buildings. However, if the building exhaust air is not clean, as would be the case for laboratories or factories, new risks are encountered when using energy exchangers. It is possible that gaseous contaminants in the exhaust air may also transfer along with the moisture through the membranes, contaminating the incoming fresh ventilation air. Current test standards provide a test procedure to determine the contamination of the fresh ventilation air by measuring the transfer of an inert tracer gas in an energy exchanger. However, the tracer gas test may not represent the transfer of common indoor air contaminants due to differences in their transport properties. Therefore, in this study, an experimental facility is developed to determine the transfer of seven different contaminants through two membranes (porous and dense) at different flow rates. Contaminant transfer is quantified using a parameter called the exhaust contaminant transfer ratio (ECTR), which gives the fraction of the contaminants transferred from the exhaust air to the fresh ventilation air. A theoretical model based on the effectiveness-number of transfer units (ε-NTU) correlation and moisture transfer resistance of the membrane is presented to determine transfer through porous membranes and validated with experimental results. The major contribution of this paper is that it presents a simple method to predict the transfer of different contaminants through a porous membrane based on the moisture transfer resistance of the membrane at different operating conditions. It is found that as contaminant diffusivity decreases, ECTR generally also decreases, and as the flow rate increases, ECTR decreases, which is consistent with the predictions from the correlation.

42 ENGINEERING↗
DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗
DOE OSTI2022

Pattern formation and flocking for particles near the jamming transition on resource gradient substrates

Here, we numerically examine a bidisperse system of active and passive particles coupled to a resource substrate. The active particles deplete the resource at a fixed rate and move toward regions with higher resources, while all of the particles interact sterically with each other. We show that at high densities, this system exhibits a rich variety of pattern-forming phases along with directed motion or flocking as a function of the relative rates of resource absorption and consumption as well as the active to passive particle ratio. These include partial phase separation into rivers of active particles flowing through passive clusters, strongly phase separated states where the active particles induce crystallization of the passive particles, mixed jammed states, and fluctuating mixed fluid phases. For higher resource recovery rates, we demonstrate that the active particles can undergo motility-induced phase separation, while at high densities, there can be a coherent flock containing only active particles or a solid mixture of active and passive particles. The directed flocking motion typically shows a transient in which the flow switches among different directions before settling into one direction, and there is a critical density below which flocking does not occur. We map out the different phases as function of system density, resource absorption and recovery rates, and the ratio of active to passive particles.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2024

High-energy, low-jitter, narrowband ps probe laser for kHz-rate fs/ps coherent anti-Stokes Raman scattering

Hybrid fs/ps coherent anti-Stokes Raman scattering (CARS) thermometry often utilizes ps probe pulses derived from pulse shaping or spectrally filtering the primary laser source or by synchronization with a low repetition rate external laser. This results in limited energy, spectral resolution, and/or repetition rate of the ps probe. In this work, a master-oscillator power-amplifier (MOPA) laser was synchronized to the oscillator of a Ti:sapphire regenerative amplifier to achieve high-energy (600 µJ), narrowband (58 ps) probe pulses at kHz repetition rates. Furthermore, temporal filtering allows the pulse characteristics to be adjusted for each application. At 25 Torr, relevant to high-speed flows, the kHz-rate MOPA system generated signal-to-noise ratios 3× higher in nitrogen and had improved precision relative to a 10 ps probe derived from spectral filtering and the power-amplifier. The MOPA system also enabled single-shot ro-vibrational hybrid fs/ps CARS thermometry in 650 K heated air.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

On-line Measurement of Hydrogen Gas using Raman Spectroscopy for Process Gas Systems

Advanced reprocessing schemes such as the Zircex process are important for continued development of the nuclear fuel cycle. The hydrochlorination reaction of the Zircex process flow sheet for metallic fuels has a known off-gas stream of H 2 (g) and HCl. The flow sheet can be simplified through HCl recycling, which requires detection and quantification of any residual H 2 (g) for safety and purity. In this work, commercially available Raman spectroscopy was applied for low-level, on-line detection of H 2 (g) in process gas streams, and parameters were adjusted to optimize the H 2 (g) Raman signal. Increasing the number of scans and exposure time in the spectrometer increased the signal-to-noise ratio of the H 2 (g) Raman signal, while the gas flow rate was optimized at 2.3 L/min. The limit of detection for H 2 (g) was estimated to be 2,100 ppm H 2 (g) in an N 2 (g) background. This scoping study for the application of Raman spectroscopy for on-line measurements of H 2 (g) shows that H 2 (g) can be detected at low levels (approximately 5% of the lower explosive limit (LEL)) in flowing gas streams.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Criticality analysis of nuclear binding energy neural networks

Machine learning methods, in particular deep learning methods such as artificial neural networks (ANNs) with many layers, have become widespread and useful tools in nuclear physics. However, these ANNs are typically treated as ‘black boxes’, with their architecture (width, depth, and weight/bias initialization) and the training algorithm and parameters chosen empirically by optimizing learning based on limited exploration. We test a non-empirical approach to understanding and optimizing nuclear physics ANNs by adapting a criticality analysis based on renormalization group flows in terms of the hyperparameters for weight/bias initialization, training rates, and the ratio of depth to width. This treatment utilizes the statistical properties of neural network initialization to find a generating functional for network outputs at any layer, allowing for a path integral formulation of the ANN outputs as a Euclidean statistical field theory. We use a prototypical example to test the applicability of this approach: a simple ANN for nuclear binding energies. We find that with training using a stochastic gradient descent optimizer, the predicted criticality behavior is realized, and optimal performance is found with critical tuning. However, the use of an adaptive learning algorithm leads to somewhat superior results without concern for tuning and thus obscures the analysis. Nevertheless, the criticality analysis offers a way to look within the black box of ANNs, which is a first step towards potential improvements in network performance beyond using adaptive optimizers.

artificial neural network↗
DOE OSTI2020

A phenomenological rate of injection model for predicting fuel injection with application to mixture formation in light-duty diesel engines

Fuel injection rate laws are one of the most important pieces of information needed when modeling engine combustion with computational fluid dynamics. In this study, a simple phenomenological model of a common-rail injector was developed and calibrated for the Bosch CRI2.2 platform. The model requires three tunable parameter fits, making it relatively easy to calibrate and suitable for injector modeling when high-fidelity information about the internal injector’s geometry and electrical circuit details are not available. Each injection pulse is modeled as a sequence of up to four stages: an injection needle mechanical opening transient, a full-lift viscous flow inertial transient, a Bernoulli steady-state stage, and a needle descent transient. Parameters for each stage are obtained as polynomial fits from measured injection rate properties. The model enforces total injected mass, and the intermediate stages are only introduced if the injection pulse duration is long enough. Experimental rates of injection from two separate campaigns on the same injector were used to calibrate the model. The model was first validated against measured injection rate laws featuring pilot injections, short partially premixed combustion pulses, and conventional diesel combustion injection strategies. Then, it was employed as an input to engine computational fluid dynamics simulations, which were run to simulate experiments of mixture formation in an optically accessible light-duty diesel engine. It was found that, though simple, this model is capable of predicting both pilot and main injection pulse mass flow rates well: the simulations yielded accurate predictions of in-cylinder equivalence ratio distributions from injection strategies for both partially premixed combustion and pilot injections. Also, once calibrated, the model produced appropriate results for a wide range of injected mass and rail pressure values. Finally, it was observed that usage of such a relatively simple model can be a good choice when high-fidelity injection rate input and highly detailed information of the injector’s geometry and operation are not available, particularly as noticeable discrepancies can be present also among different experimental campaigns on similar hardware.

42 ENGINEERING↗
DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Comparison of Sub-Ppm Instrument Response Suggests Higher Detection Limits Could Be Used to Quantify Methane Emissions from Oil and Gas Infrastructure

Quantifying and controlling fugitive methane emissions from oil and gas facilities remains essential for addressing climate goals, but the costs associated with monitoring millions of production sites remain prohibitively expensive. Current thinking, supported by measurement and simple dispersion modelling, assumes single-digit parts-per-million instrumentation is required. To investigate instrument response, the inlets of three trace-methane (sub-ppm) analyzers were collocated on a facility designed to release gas of known composition at known flow rates between 0.4 and 5.2 kg CH 4 h –1 from simulated oil and gas infrastructure. Methane mixing ratios were measured by each instrument at 1 Hertz resolution over nine hours. While mixing ratios reported by a cavity ring-down spectrometer (CRDS)-based instrument were on average 10.0 ppm (range 1.8 to 83 ppm), a mid-infrared laser absorption spectroscopy (MIRA)-based instrument reported short-lived mixing ratios far larger than expected (range 1.8 to 779 ppm) with a similar nine-hour average to the CRDS (10.1 ppm). We suggest the peaks detected by the MIRA are likely caused by a micrometeorological phenomenon, where vortex shedding has resulted in heterogeneous methane plumes which only the MIRA can observe. Further analysis suggests an instrument like the MIRA (an optical-cavity-based instrument with cavity size ≤10 cm3 measuring at ≥2 Hz with air flow rates in the order of ≤0.3 slpm at distances of ≤20 m from the source) but with a higher detection limit (25 ppm) could detect enough of the high-concentration events to generate representative 20 min-average methane mixing ratios. Even though development of a lower-cost, high-precision, high-accuracy instrument with a 25 ppm detection threshold remains a significant problem, this has implications for the use of instrumentation with higher detection thresholds, resulting in the reduction in cost to measure methane emissions and providing a mechanism for the widespread deployment of effective leak detection and repair programs for all oil and gas infrastructure.

03 NATURAL GAS↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗