SEARCH · Engineering Papers

Results for “VISCOSITY”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,549 records · Page 86

All sources NASA NTRS DOE OSTI

NASA NTRS2009

A Hybrid Model for Leveed Lava Flows: Implications for Eruption Styles on Mars

Many cehannelized lava flows on the plains of Mars have substantial embanking margins and levees inferred to have been stationary while the central channel was active. Levee formation can be attributed to two end-member processes during emplacement; construction during passage of the flow front and growth along the entire length of the flow while it is active. It is shown here that the amount of lava that can be deposited by the flow front alone is limited. Estimates of the levee volume for many Mars plains flows exceed this limit and must have formed by processes that continued after the passage of the front. Experimental studies of analogous laboratory flows also indicate a combination of both modes of emplacement. A model that combines both modes of levee formation. is presented, including a method for estimating volumetric flow rate, eruption duration, and viscosity. Six lava flows on the plains of the Tharsis volcanic province are used as illustrative examples. Crustal thicknesses for the six flows examined range from 9 to 23 m. Estimated emplacement times required to cool crusts of these thicknesses range from I year to 10 years. Correspondini viscosities are on the order of 10 5-106 Pa s. Effusion rates range from 25 to 840 m 3 s - and are all within the range of terrestrial observations. Therefore, the large leveed plains flows on Mars are not dramatically different in eruption rate or lava viscosity from large terrestrial analogs.

Glaze, Lori S.↗

NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗

NASA NTRS2011

Venusian Applications of 3D Convection Modeling

This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.

Venus↗

NASA NTRS2014

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗

NASA NTRS2014

The Tethered Moon

A reasonable initial condition on Earth after the Moonforming impact is that it begins as a hot global magma ocean1,2. We therefore begin our study with the mantle as a liquid ocean with a surface temperature on the order of 3000- 4000 K at a time some 100-1000 years after the impact, by which point we can hope that early transients have settled down. A 2nd initial condition is a substantial atmosphere, 100-1000 bars of H2O and CO2, supplemented by smaller amounts of CO, H2, N2, various sulfur-containing gases, and a suite of geochemical volatiles evaporated from the magma. Third, we start the Moon with its current mass at the relevant Roche limit. The 4th initial condition is the angular momentum of the Earth-Moon system. Canonical models hold this constant, whilst some recent models begin with considerably more angular momentum than is present today. Here we present a ruthlessly simplified model of Earth's cooling magmasphere based on a full-featured atmosphere and including tidal heating by the newborn Moon. Thermal blanketing by H2O-CO2 atmospheres slows cooling of a magma ocean. Geochemical volatiles - chiefly S, Na, and Cl - raise the opacity of the magma ocean's atmosphere and slow cooling still more. We assume a uniform mantle with a single internal (potential) temperature and a global viscosity. The important "freezing point" is the sharp rheological transition between a fluid carrying suspended crystals and a solid matrix through which fluids percolate. Most tidal heating takes place at this "freezing point" in a gel that is both pliable and viscous. Parameterized convection links the cooling rate to the temperature and heat generation inside the Earth. Tidal heating is a major effect. Tidal dissipation in the magma ocean is described by viscosity. The Moon is entwined with Earth by the negative feedback between thermal blanketing and tidal heating that comes from the temperature-dependent viscosity of the magma ocean. Because of this feedback, the rate that the Moon's orbit evolves is limited by the modest radiative cooling rate of Earth's atmosphere, which in effect tethers the Moon to the Earth. Consequently the Moon's orbit evolves orders of magnitude more slowly than in conventional models. Slow orbital evolution promotes capture by orbital resonances that may have been important in the Earth-Moon system

Earth-Moon system↗

NASA NTRS2014

Thermophysical Property Measurements of Silicon-Transition Metal Alloys

Metals and metallic alloys often have high melting temperatures and highly reactive liquids. Processing reactive liquids in containers can result in significant contamination and limited undercooling. This is particularly true for molten silicon and it alloys. Silicon is commonly termed "the universal solvent". The viscosity, surface tension, and density of several silicon‐transition metal alloys were determined using the Electrostatic Levitator system at the Marshall Space Flight Center. The temperature dependence of the viscosity followed an Arrhenius dependence, and the surface tension followed a linear temperature dependence. The density of the melts, including the undercooled region, showed a linear behavior as well. Viscosity and surface tension values were obtain for several of the alloys in the undercooled region.

Banish, R. Michael↗

NASA NTRS2015

Hybrid Reynolds-Averaged/Large Eddy Simulation of a Cavity Flameholder; Assessment of Modeling Sensitivities

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit Reynolds stress model. Fortunately, the numerical error assessment at most of the axial stations used to compare with measurements clearly indicated that the scale-resolving simulations were improving (i.e. approaching the measured values) as the grid was refined. Hence, unlike a Reynolds-averaged simulation, the hybrid approach provides a mechanism to the end-user for reducing model-form errors.

Baurle, R. A.↗

NASA NTRS2016

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism↗

NASA NTRS2018

Thermophysical Properties of Nickel-Based Superalloys

The NASA Marshall Space Flight Center (MSFC) electrostatic levitation (ESL) laboratory has a long history of providing materials research and thermophysical property data. The lab can measure thermophysical properties such as density, surface tension, and viscosity of liquid materials, including elements, alloys, glasses, ceramics, and oxides. Nickel-based superalloys (e.g. Inconel, Hastelloy, and Waspaloy) have many high performance applications, including turbine engines for aerospace. Superalloy parts are typically manufactured by casting and forging. These processes generate both polycrystalline and monocrystalline parts. A relatively new method of fabrication of turbine disk materials is additive manufacturing, which is typically done for aerospace parts by powder-bed methods such as selective laser melting (SLM). Accurate modeling of casting and forging, as well as additively manufacturing processes, require thermophysical properties (density, surface tension, and viscosity). The thermophysical properties of liquid nickel-based superalloys, both conventional and additively manufactured, were measured at several temperatures in both undercooled and superheated condition. Surface tension and viscosity was measured using the oscillating drop method and density and thermal expansion was measured using edge detection methods.

SanSoucie, Michael P.↗

NASA NTRS2018

Numerical Study of Turbulent Separation Bubbles with Varying Pressure Gradient and Reynolds Number

A family of cases each containing a small separation bubble is treated by direct numerical simulation (DNS), varying two parameters: the severity of the pressure gradients, generated by suction and blowing across the opposite boundary, and the Reynolds number. Each flow contains a well-developed entry region with essentially zero pressure gradient, and all are adjusted to have the same value for the momentum thickness, extrapolated from the entry region to the centre of the separation bubble. Combined with fully defined boundary conditions this will make comparisons with other simulations and turbulence models rigorous; we present results for a set of eight Reynolds-averaged Navier–Stokes turbulence models. Even though the largest Reynolds number is approximately 5.5 times higher than in a similar DNS study we presented in 1997, the models have difficulties matching the DNS skin friction very closely even in the zero pressure gradient, which complicates their assessment. In the rest of the domain, the separation location per se is not particularly difficult to predict, and the most definite disagreement between DNS and models is near reattachment. Curiously, the better models tend to cluster together in their predictions of pressure and skin friction even when they deviate from the DNS, although their eddy-viscosity levels are widely different in the outer region near the bubble (or they do not rely on an eddy viscosity). Stratford’s square-root law is satisfied by the velocity profiles, both at separation and reattachment. The Reynolds-number range covers a factor of two, with the Reynolds number based on the extrapolated momentum thickness equal to approximately 1500 and 3000. This allows tentative estimates of the improvements that even higher values will bring to the model comparisons. The solutions are used to assess models through pressure, skin friction and other measures; the flow fields are also used to produce effective eddy-viscosity targets for the models, thus guiding turbulence-modelling work in each region of the flow.

Coleman, G. N.↗

NASA NTRSDate not supplied

Investigating Fifth Oxide Effect on CMXAS Glass Properties

Coatings for hot section aero turbine engine structural materials are subject to thermochemical degradation after interacting with calcium-magnesium-aluminosilicates (CMAS). Molten CMAS viscosity is the primary glass property of importance, correlating with coating infiltration and reaction depth. CMAS viscosity was measured using a spindle-dipped viscometer at temperatures between 1300-1550°C in 25°C increments. Experimental data were compared to Thermo-Calc and FactSage computational models, which showed good agreement with experimental trends. CMXAS compositions, where X denotes a fifth oxide, of Ca30.75-Mg9-X5-Al13-Si42.25 (single cation oxide mol%) for X = Fe2+, Fe3+, Gd3+, Hf4+, Na1+, Ti4+, Y3+, Yb3+, Zr4+ additions were used. Oxide precipitation was observed in X = Hf4+ , Zr4+ CMXAS specimens, indicating coating dissolution is an increased threat for T ≥ 1390°C. Molten CMAS and CMXAS viscosity measurements were shown to correlate with net cation field strength, offering an alternative in property prediction where thermodynamic data are unavailable.

Clark Avery Luckhardt↗

DOE OSTI2016

Process and Material Design for Micro-Encapsulated Ionic Liquids in Post-Combustion CO 2 Capture

Aprotic Heterocyclic Anion (AHA) Ionic Liquids (ILs) have been identified as promising new solvents for post-combustion carbon capture due to their high CO 2 uptake and the high tenability 1,2 of their binding energy with CO 2 . Some of these compounds change phase (solid to liquid) on absorption of CO 2 ; these Phase Change ILs (PCILs)3 offer the additional advantage that part of the heat needed to desorb the CO2 from the absorbent is provided by the heat of fusion as the PCIL solidifies upon release of CO 2 . However, the relatively high viscosity of AHA ILs and the occurrence of a phase change in PCILs present challenges for conventional absorption equipment. To overcome these challenges we are pursuing the use of new technology to micro-encapsulate the AHA ILs and PCILs. Our partners at Lawrence Livermore National Laboratory have successfully demonstrated this technology in the application of post-combustion carbon capture with sodium and potassium carbonate solutions,4 and have recently shown the feasibility of micro-encapsulation of an AHA IL for carbon capture.5 The large effective surface area and high CO 2 permeability of the micro-capsules is expected to offset the drawback of the high IL viscosity and to provide for a more efficient and cost-effective mass transfer operation involving AHA ILs and PCILs. These opportunities, however, present us with both process and materials design questions. For example, what is the target CO 2 absorption strength (enthalpy of chemical absorption) for the tunable AHA IL? What is the target for micro-capsule diameter in order to obtain a high mass transfer rate and good fluidization performance? What are the appropriate temperatures and pressures for the absorber and stripper? In order to address these and other questions, we have developed a rate-based model of a post-combustion CO 2 capture process using micro-encapsulated ILs. As a performance baseline, we have also developed a rate-based model of a standard packed bed absorber using an un-encapsulated AHA IL absorbent. Using such models we can determine optimal CO 2 capture performance and investigate the sensitivity of the optimum with respect to the key thermo-physical and transport properties of the IL (e.g., CO 2 binding energy, viscosity, etc.) and the micro-capsules (e.g. diameter, CO 2 permeability, etc.). Results of these process and material design studies will be presented, and the performance of this novel micro-encapsulation technology will be assessed.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2017

Internal Nozzle Flow Simulations of Gasoline-Like Fuels under Diesel Operating Conditions

Spray formation in internal combustion engines with direct injection is strictly correlated with internal nozzle flow characteristics, which are in turn influenced by fuel physical properties and injector needle motion. This paper pre-sents a series of 3D simulations that model the in-nozzle flow in a 5-hole mini-sac diesel injector. Two gasoline-like naphtha fuels, namely full-range and light naphtha, were tested under operating conditions typical of diesel applica-tions and were compared with n-dodecane, selected from a palette used as diesel surrogates. Validated methodolo-gies from our previous work were employed to account for realistic needle motion. The multi-phase nature of the problem was described by the mixture model assumption with the Volume of Fluid method. Cavitation effects were included by means of the Homogeneous Relaxation Model and turbulence closure was achieved with the Standard k-ε model in an Unsteady Reynolds-Averaged Navier-Stokes formulation. The results revealed that injector perfor-mance and propensity to cavitation are influenced by the fuel properties. Analyses of several physical quantities were carried out to highlight the fuel-to-fuel differences in terms of mass flow rate, discharge coefficients, and fuel vapor volume fraction inside the orifices. A series of parametric investigations was also performed to assess the fuel response to varied fuel injection temperature, injection pressure, and cross-sectional orifice area. For all cases, the strict correlation between cavitation magnitude and saturation pressure was confirmed. Owing to their higher volatil-ity, the two gasoline-like fuels were characterized by higher cavitation across all the simulated conditions. Occur-rence of cavitation was mostly found at the needle seat and at the orifice inlets during the injection event’s transient, when very small gaps exist between the needle and its seat. This behavior tended to disappear at maximum needle lift, where cavitation was absent for all fuels. Differences in mass flow rate between the naphtha fuels and n-dodecane were measured and ascribed to the different densities of the three fuels. Nevertheless, they were found to be smaller than expected, owing to the lower viscosity of the gasoline-like fuels. This beneficial influence of the lower viscosity was shown to be less effective at higher temperature, where the relative viscosity differences de-creased.

Torelli, R.↗

DOE OSTI2017

Sorting process of nanoparticles and applications of same

In one aspect of the present invention, a method for sorting nanoparticles includes preparing a high-viscosity density gradient medium filled in a container, dispersing nanoparticles into an aqueous solution to form a suspension of the nanoparticles, each nanoparticle having one or more cores and a shell encapsulating the one or more cores, layering the suspension of the nanoparticles on the top of the high-viscosity density gradient medium in the container, and centrifugating the layered suspension of the nanoparticles on the top of the high-viscosity density gradient medium in the container at a predetermined speed for a predetermined period of time to form a gradient of fractions of the nanoparticles along the container, where each fraction comprises nanoparticles in a respective one of aggregation states of the nanoparticles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2018

Compositions for, solutions for, and methods of use of siloxane based aromatic trisureas as viscosifiers

A compound represented by the following formula is provided: ##STR00001## Also provided is a solution including a compound disclosed herein, a volume of dense carbon dioxide (CO.sub.2), and a co-solvent, where the solution has an increased viscosity greater than the viscosity of dense CO.sub.2. Methods of increasing the viscosity of dense CO.sub.2 and natural gas liquids (NGLs) by, for example, dissolving a compound disclosed herein to form a solution, are also provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Foot-scale Evaluation of CO2-Responsive Polymer and CO2 Binary Fluid as an Alternative Fracturing Fluid for Enhanced Geothermal Systems

A foot-scale high-temperature true-triaxial fracturing apparatus was used to conduct hydraulic fracture tests on rock samples. The fracturing tests used different fluids including water, CO2, CO2 with water, and CO2 with a 1% polyallylamine (PAA)/water solution. The tests were conducted with both constant pressure injection and constant flow rate injection modes. CO2-based fracturing fluids were found to produce higher breakdown pressures, high transient flow rates, and produce higher-conductivity fractures as compared to water-based fracturing fluids. Additionally, faster pressurization rates with CO2-based fracturing fluids (obtained when fracturing at constant flow rate mode) are found to be associated with higher fracture conductivities. When fracturing with CO2-based fluids in the presence of a PAA aqueous solution, for example if the rock is saturated with 1wt% PAA in water, the volume expansion caused by CO2-induced crosslinking of PAA leads to a faster pressure increase due to the associated volume expansion and increase in viscosity. It was also found that CO2 as a fracturing fluid injected in hot dry rock (HDR) attain the highest fracture conductivity only when injected at very high flow rates, followed very closely by the CO2/PAA fracturing fluid system that generates fractures with, on average, similarly high conductivity values though independently of injection flow rate and using 1/6 of the mass of CO2 as compared to CO2 in HDR. Breakdown pressures were also similar for CO2 stimulation in HDR and CO2/PAA fluid system under identical injection flow rates. Finally, the well-known low viscosity of CO2 phase prevents the efficient transport of proppants while the reacted CO2/PAA fluid system is known to form a high viscosity binary fracturing fluid system with the potential to carry proppant in addition of being a fine fracturing fluid.

Enhanced Geothermal Systems, polyallylamine, Co2-r↗

DOE OSTI2022

Study on Emulsification/Demulsification Behavior and Mechanism of Produced Liquid from Polymer Flooding on Alaska North Slope

Heavy oil reservoirs on Alaska North Slope (ANS) are unconsolidated and contain abundant clay minerals, where the first-ever field pilot is currently implemented to validate the use of polymer floods for heavy oil enhanced oil recovery (EOR). The polymer molecules and/or fine clay particles carried with the produced liquid could potentially affect the oil/water separation, which is one of the major concerns for field operators. This dissertation aims to investigate the emulsification behavior of produced liquid, understand the emulsifying mechanism, and seek an adaptive and cost-effective method to treat the produced liquid from polymer flooding. Emulsions were prepared by mechanically mixing the actual heavy oil and the produced water from the pilot site, of which the stability was investigated by bottle test method or multiple light scattering method. Drop size distribution and interfacial properties were measured via microscope and pendant drop technique to probe the stability mechanism further. Results showed that oil-continuous or water-continuous emulsion could be generated depending on the water cut, clay types, clay concentration, and polymer concentration. In the crude oil/water system, the increasing water cut triggered the phase inversion of oil-continuous emulsion to water-continuous emulsion, resulting in faster separation and lower emulsion stability. Whereas, clay particles, no matter added to the oil or water phase, resulted in an unfavorable phase inversion from the loose water-continuous emulsion to the tight oil-continuous emulsion as clay concentration increased. For all four types of clay except Ca-montmorillonite, clay particles added to water led to an earlier phase inversion and higher emulsion stability than that added to the oil. The dual function of polymer on emulsion stability was observed. On the one hand, both sheared and unsheared polymer tended to convert the oil-continuous emulsion formed in either crude oil/water system or complex crude oil/water/clay particle system to the water-continuous emulsion, acting as a weak demulsifier to accelerate the oil/water separation. On the other hand, the addition of polymer to the water-continuous emulsion could result in enhanced emulsion stability, which is primarily attributed to the increased viscosity of the continuous phase and the decreased drop size of the dispersed oil phase. Particularly, the sheared polymer had a weaker ability to stabilize the o/w emulsion than the unsheared polymer due to the lower viscosity of the sheared polymer solution resulting from the breakdown of the polymer macromolecules. As for the chemical demulsification tests, the performance of demulsifiers showed a complex dependency upon the water cut, the shearing intensity, demulsifier type and dosage, and the polymer concentration. A compound emulsion breaker, E12+E18, exhibited the most satisfactory demulsification performance despite the varied test conditions. For severe water-continuous emulsions that might require a multi-fold dosage of demulsifier, a less expensive electrolyte, KCl, was proposed to be used in combination with demulsifier E12+E18 to improve the demulsification performance. In the proposed demulsifier formula, the effectiveness of the commercial demulsifier relied on its destructive effect on the interfacial film, while the efficacy of KCl was mainly dependent on its viscosity reduction effect on the continuous phase. This dissertation illustrates that intermediate layer elimination and water clarification are the major challenges for produced liquid treatment from polymer flooding. It also provides practical and theoretical guidance in advance for the demulsification strategy of the produced liquid from the ongoing first-ever polymer flooding pilot on ANS.

Chang, Hongli↗

DOE OSTI2024

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗