Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid characterization”

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 1,459 records · Page 81

Development of a CO 2 Chemical Sensor for Downhole CO 2 Monitoring in Carbon Sequestration

Geologic storage of carbon dioxide (CO 2 ) has been proposed as a viable means for reducing anthropogenic CO 2 emissions. The means for geological sequestration of CO 2 is injection of supercritical CO 2 underground, which requires the CO 2 to remain either supercritical, or in solution in the water/brine present in the underground formation. However, there are aspects of geologic sequestration that need further study, particularly in regards to safety. To date, none of the geologic sequestration locations have been tested for storage integrity under the changing stress conditions that apply to the sequestration of very large amounts of CO 2 . Establishing environmental safety and addressing public concerns require widespread monitoring of the process in the deep subsurface. In addition, studies of subsurface carbon sequestration such as flow simulations, models of underground reactions and transports require a comprehensive monitoring process to accurately characterize and understand the storage process. Real-time information about underground CO 2 movement and concentration change is highly helpful for: (1) better understanding the uncertainties present in CO 2 geologic storage; (2) improvement of simulation models; and (3) evaluation of the feasibility of geologic CO 2 storage. Current methods to monitor underground CO 2 storage include seismic, geoelectric, isotope and tracer methods, and fluid sampling analysis. However, these methods commonly resulted low resolution, high cost, and the inability to monitor continuously over the long time scales of the CO 2 storage process. A preferred way of monitoring in-situ underground CO 2 migration is to continuous measure CO 2 concentration change in brine during the carbon storage process. An approach to obtain the real time information on CO 2 concentration change in formation solution is highly demanded in carbon storage to understand the CO 2 migration subsurface and to answer the public safety problem. The objective of the study is to develop a downhole CO 2 sensor that can in-situ, continuously monitor CO 2 concentration change in deep saline. The sensor is a Severinghaus-type CO 2 sensor with small size, which renders it can be embedded in monitoring well casing or integrated with pressure/temperature transducers, enabling the development of “smart” wells. The studies included: (1) prepare and characterize metal-oxide electrodes. Test the electrodes response to pH change. Investigate different ions and brine concentration effects on the electrode’s performance. Study the stability of the electrode in brine solution; (2) fabricate a downhole CO 2 sensor with the metal-oxide electrodes prepared in the laboratory. Test the performance of the CO 2 sensor in brine solutions. Study high pressure effects on the performance of the sensor; (3) design and conduct CO 2 /brine coreflooding experiments with the CO2 sensor. Monitor CO 2 movement along the core and test the performance of the sensor in coreflooding tests. Develop a data acquisition system that can digitize the sensor’s output voltage. Our completed research has resulted in deep understanding of downhole CO 2 sensor development and CO 2 monitoring in CO 2 storage process. The developed downhole CO 2 sensor included a metal-oxide electrode, a gas-permeable membrane, a porous steel cup, and a bicarbonate-based internal electrolyte solution. Iridium oxide-based electrode was prepared and used for preparation the CO 2 sensor. The prepared iridium oxide-based electrode displayed a linearly response to pH change. Different factors such as different ions and ions concentration, temperature, and pressure effects on the electrode performance on pH response were investigated. The results indicated that the electrode exhibited a good performance even in high salt concentration of produced water. To improve the electrode performance under high pressure, IrO 2 nanoparticles with the particle size in the range of 1-2 nm were prepared and electrodeposited on stainless steel substrate by cyclic voltammetry. It was observed that the thin film of iridium oxide was formed on the substrate surface and such iridium oxide-based electrode displayed excellent performance under high pressure for longer term. A downhole CO 2 sensor with the iridium oxide-based electrode was prepared. The working principle of the CO 2 sensor is based on the measurement of the pH change of the internal electrolyte solution caused by the hydrolysis of CO 2 and then determination of the CO 2 concentration in water. The prepared downhole CO 2 sensor had the size of diameter of 0.7 in. and length of 1.5 in. The sensor was tested under the pressures of 500 psi, 2,000 psi, and 3,000 psi. A linear correlation was observed between the sensor potential change and dissolved CO 2 concentration in water. The response time of the CO 2 sensor was in the range of 60-100 minutes. Further tests indicated that the CO 2 sensor exhibited good reproducibility under high pressure. A CO 2 /brine coreflooding system was constructed to simulate the real-world CO 2 storage process. The prepared downhole CO 2 sensor was loaded in the system to monitor CO 2 movement during CO 2 /brine coreflooding test. The results indicated that the sensor could detect CO 2 movement in the tests. Further studies showed that the sensor could be recovered by brine flooding after CO 2 /brine flushed the core. The results of the coreflooding tests demonstrated that the sensor had potential application for CO 2 monitoring in carbon sequestration. A data acquisition system for the downhoe CO 2 sensor was developed and coded. The system converted the sensor output signal into digital data and transported the data from downhole to wellhead surface. The data acquisition system was tested and evaluated in the laboratory with the prepared sensor for data collection.

36 MATERIALS SCIENCE↗

Investigation of Shear-Thinning Behavior on Film Thickness and Friction Coefficient of Polyalphaolefin Base Fluids With Varying Olefin Copolymer Content

This study investigates the rheological properties, elastohydrodynamic (EHD) film-forming capability, and friction coefficients of low molecular mass poly-alpha-olefin (PAO) base stocks with varying contents of high molecular mass olefin copolymers (OCPs) to assess their shear stability and their potential for energy-efficient lubrication. Several PAO-OCP mixtures were blended in order to examine the relationship between their additive content and tribological performance. Gel permeation chromatography (GPC) and nuclear magnetic resonance (NMR) spectroscopy were used to characterize the molecular masses and structures, respectively. Density, viscosity, EHD film thickness, and friction were measured at 303 K, 348 K, and 398 K. Film thickness and friction were studied at entrainment speeds relevant to the boundary, mixed, and full-film lubrication regimes. The PAO-OCP mixtures underwent temporary shear-thinning resulting in decreases in film thickness and hydrodynamic friction. These results demonstrate that the shear characteristics of PAO-OCP mixtures can be tuned with the OCP content and provide insight into the effects of additives on EHD characteristics.

Elastohydrodynamic lubrication↗

DEEPEN Data Catalog for Magmatic Geothermal Systems in the United States

This data catalog contains information related to the Training Site Analysis for the Geothermica project "DE-risking Exploration of geothermal Plays in magmatic ENvironments (DEEPEN)." The DEEPEN project aims to reduce exploration risk for geothermal fluids in magmatic systems by developing improved an improved framework for interpretation of exploration data using the Play Fairway Analysis (PFA) methodology. The Training Site Analysis performed for DEEPEN leverages existing datasets to develop a customized PFA approach to exploration for multiple geothermal resource types in magmatic systems (conventional hydrothermal resources, supercritical fluid and superheated steam resources, and superhot EGS resources). This data catalog contains links to publicly available data files related to 8 training sites in the United States. US training sites are: the Cascades/Aleutians PFA project; the Hawaii PFA project, the Oregon Cascades PFA project, the Snake River Plain, Idaho PFA project, the Washington State PFA project, Newberry Volcano, Coso Geothermal Field, and the Geysers Geothermal field. This database contains an overview of these training sites, data sources, and links to publicly available exploration datasets. For the five PFA projects, details on exploration data related to PFA components (heat, fluid, permeability, sometimes seal) are provided, including a summary of data weighting methodologies.

15 GEOTHERMAL ENERGY↗

The Relative Importance of Saturated Silica Sand Interfacial and Pore Fluid Geochemistry on the Spectral Induced Polarization Response

Adsorption at the solid-pore fluid interface is a key mechanism controlling the mobility of nutrients and contaminants in subsurface soils and sediments. The spectral induced polarization (SIP) method has been shown to be sensitive to the quantity and type of adsorbed ions. Extending previous results, we investigated the relevance of pH, solution conductivity, and ion type on the SIP response of saturated silica sand. We also performed adsorption experiments to evaluate whether adsorption plays a relevant role on the effect of saturating solution conductivity and pH. Given their environmental relevance and different electrochemical characteristics, we focused on exploring the influence of Cu 2+ and Na + adsorption on the SIP signature. The adsorption results confirm the expected and modeled pH influence on the adsorption of both Cu 2+ and Na + . The measured quadrature conductivity spectra indicate that pH and solution conductivity control the electrical double layer electrochemical state and its capacitive behavior. On the contrary, no appreciable SIP signal changes are associated with ion substitution. The adsorption experiments highlight low values of site occupancy for Na and Cu, which suggests that the effects of pH and fluid conductivity are unrelated to their control on the ion adsorption. We interpret the solution conductivity as a proxy for ionic strength. The relative importance of pH and solution conductivity over ion type helps to further constrain the interpretation of SIP results in field geochemical and biogeochemical characterization and monitoring.

58 GEOSCIENCES↗

Wetland Soil Characterization and Methane Production Impacted by Nickel Addition, Argonne and Tims Branch Wetlands, September and October 2020

Abstract:Freshwater wetland soils are foci of biogeochemical cycling as they serve as key sources of methane to the atmosphere. An array of metalloenzymes is essential to anaerobic microbial carbon transformations. Nickel is notably recognized as playing key roles in the enzymatic pathways of methanogenesis. Low availability of trace metals limits microbial element cycling in laboratory studies, but the occurrence of such limitations in natural subsurface aquatic systems is poorly understood. Microcosm incubation studies were carried out using two distinct wetland soils, one from a marsh wetland and the second from a riparian wetland, to explore the effect of dissolved Ni concentrations on methane production. Data are provided for wetland soil characterization and soil incubation experiments using materials from marsh wetlands at Argonne National Laboratory and riparian wetlands in the Tims Branch watershed at Savannah River National Laboratory. The characterization data consists soil carbon, nitrogen, sulfur, and iron contents plus as well as the solid-phase concentrations of copper, nickel, cobalt, and zinc, bioessential trace metals that may limits microbial metabolic process if they have low availability. The data for the soil incubation experiments include fluid pH, fluid dissolved trace metal concentrations, and cumulative methane production. Three soil incubations are reported: marsh wetland soil with increasing nickel addition, marsh wetland soil in sulfate-free water with increasing nickel addition, and riparian wetland soil with increasing nickel addition. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units. Note that "u" is used in place of Greek lower case mu to indicate the micro prefix on units. A Table of Contents file (Yan_Soil_Incubations_2020_TOC.txt) provides an index for the data contained in the individual files.

54 ENVIRONMENTAL SCIENCES↗

Characterization of Reconsolidated Crushed Salt from the BAMBUS Site.

Observational petrofabrics, thermal, mechanical, and hydrological measurements were made on reconsolidated salt samples extracted from the field site in which a study called Backfilling and Sealing of Underground Repositories for Radioactive Waste in Salt was conducted. Similar characterization was completed more than a decade ago, so this work furthers previous measurements after sustained consolidation in situ . Porosity determined by traditional point-counting on polished sections and helium porosimeter methods ranged from 20-25% with consolidation governed by brittle processes, as evidence of fluid-aided, grain-boundary processes was rarely observed. Thermal conductivity in the range of 2.3 W /( m * K ) is consistent for granular halite in this porosity range. Gas flow measurements yielded permeability of the order of 5e -13 m 2 . Pressure-sensitive compressive strengths at 0.5, 1.0, and 2.0 MPa confining pressure were 8, 9, and 14 MPa, respectively, with apparent elastic moduli increase with deformation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Characterization of Reconsolidated Crushed Salt from the BAMBUS Site

Observational petrofabrics, thermal, mechanical, and hydrological measurements were made on reconsolidated salt samples extracted from the field site in which a study called Backfilling and Sealing of Underground Repositories for Radioactive Waste in Salt was conducted. Similar characterization was completed more than a decade ago, so this work furthers previous measurements after sustained consolidation in situ . Porosity determined by traditional point-counting on polished sections and helium porosimeter methods ranged from 20-25% with consolidation governed by brittle processes, as evidence of fluid-aided, grain-boundary processes was rarely observed. Thermal conductivity in the range of 2.3 W/(m∙K) is consistent for granular halite in this porosity range. Gas flow measurements yielded permeability of the order of 5e- 13 m 2 . Pressure-sensitive compressive strengths at 0.5, 1.0, and 2.0 MPa confining pressure were 8, 9, and 14 MPa, respectively, with apparent elastic moduli increase with deformation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Plasma core reactor applications

Analytical and experimental investigations were conducted to demonstrate the feasibility of fissioning uranium plasma core reactors and to characterize space and terrestrial applications for such reactors. Uranium hexafluoride fuel is injected into core cavities and confined away from the surface by argon buffer gas injected tangentially from the peripheral walls. Radiant heat transfer calculations were performed for a six-cavity reactor configuration. Axial working fluid channels are located along a fraction of each cavity peripheral wall. Results of calculations for outward-directed radiant energy fluxes corresponding to radiating temperatures of 2000 to 5000 K indicate total operating pressures from 80 to 650 atm, centerline temperatures from 6900 to 30,000 K, and total radiated powers from 25 to 2500 MW, respectively. Applications are described for this type of reactor such as (1) high-thrust, high specific impulse space propulsion, (2) highly efficient systems for generation of electricity, and (3) hydrogen or synthetic fuel production systems using the intense radiant energy fluxes.

Latham, T. S.↗

Cardiovascular responses to hypogravic environments

The cardiovascular deconditioning observed during and after space flight is characterized in a review of human space and simulation studies and animal simulations. The various simulation techniques (horizontal bed rest, head-down tilt, and water immersion in man, and immobilization of animals) are examined, and sample results are presented in graphs. Countermeasures such as exercise regimens, fluid replacement, drugs, venous pooling, G-suits, oscillating beds, electrostimulation of muscles, lower-body negative pressure, body-surface cooling, and hypoxia are reviewed and found to be generally ineffective or unreliable. The need for future space experimentation in both humans and animals is indicated.

Sandler, H.↗

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗

Leveraging Large Eddy Simulations to Characterize Aeroacoustic Environments on Launch Vehicles

NASA's Space Launch System (SLS) program has relied on multiple wind tunnel campaigns to predict expected aeroacoustic environments. Time, cost, Reynolds number matching, contamination from test section acoustics, and limitations in spatial resolution are all factors that impact the use experimental data for development of launch vehicle environments. Methodological and computational advances over the last decade are enabling scale resolving computational fluid dynamics (CFD) simulations to be completed that can address some the challenges with conventional datasets, enabling more complete/accurate modeling of dynamic loading events on ascent. These simulations are classified as so-called "Large-Eddy Simulations" (LES). As the name implies, these simulations seek to directly resolve the turbulence in the largest, most-energetic turbulent structures, or eddies, and the dissipation of the smallest eddies in the flow field, is modeled. Like anything else, these simulations still must undergo rigorous testing and validation for the results to ultimately be trusted in the analysis of a launch vehicle. EV33 has begun to conduct some of this early-stage testing and validation to better support NASA missions in the future.

Jonathan Boustani↗

Quasi‐dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

Abstract We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi‐linear viscoelastic stress–strain response of the parenchyma (alveoli) tissue. We adapted a continuum‐based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi‐dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

60 APPLIED LIFE SCIENCES↗

Combining selection valve and mixing chamber for nanoflow gradient generation: Toward developing a liquid chromatography cartridge coupled with mass spectrometer for protein and peptide analysis

In relation to developing a micro HPLC cartridge, we have recently built a high-pressure electroosmotic pump (EOP). However, we do not recommend people to use this pump to deliver an organic solvent directly, because it often makes the pump rate unstable. We have experimented several approaches to address this issue, but none of them are satisfactory. Here, we develop an innovative approach to address this issue. We first create an abruption (a dead-volume) within a fluid conduit. We then utilize an EOP to withdraw, via a selection valve, a train of eluent solutions having decreasing eluting power into the fluid conduit. When these solutions are further aspirated through the dead-volume, these solutions are partially mixed, smoothening concentration transitions between two adjacent eluent solutions. As these solutions are pushed back, through the dead-volume again, a smooth gradient profile is formed. In this work, we characterize this scheme for gradient formation, and we incorporate this approach with a high-pressure EOP, a nanoliter injection valve, and a capillary column, yielding a micro HPLC system. Lastly, we then couple this micro HPLC with an electrospray ionization – mass spectrometer for peptide and protein separations and identifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗

Solid–liquid separation of lignocellulosic sugars from biomass by rotating ceramic disc filtration

In many biomass conversion processes, the separation of cellulosic sugars from residual, lignin-rich solids is a critical step, and achieving high recovery yields of sugars by conventional tangential crossflow and vacuum filtration is challenged by the presence of biomass solids, which rapidly foul filters, resulting in decreased throughput. Considering the performance limitations of traditional filtration methods, dynamic filtration, which generates high shear at the membrane surface to decrease fouling, is emerging as a viable alternative for demanding solid–liquid separations. For high solids separations, there is little available information regarding the performance, limitations, and energy consumption of dynamic filtration. To that end, here we characterized the performance of a dynamic filtration module, specifically a rotating ceramic disc (RCD) filter, for the aseptic recovery of cellulosic sugars from biomass solids following pretreatment and enzymatic hydrolysis. We show how RCD rotational velocity and percent biomass solids impact the filter throughput. Additionally, we used computational fluid dynamics (CFD) simulations to estimate the shear rate at the membrane surface and to visualize hydrodynamic profiles within the module. With the combined CFD simulations and experimental results, we estimated the energy demand and operating expenses for a viable dynamic filtration system operating with a lignocellulosic feed slurry. Our results indicate that an RCD filter can achieve ≥95% recovery of sugars and produce a retentate slurry containing 12 wt% insoluble solids with low energy consumption (a 2.2-fold improvement over cross-flow filtration) and low operating costs ($\$$0.06 per kg sugars). These results show a viable path towards operationally reliable, energy efficient, and cost-effective separations of sterilized cellulosic sugars from biomass solids and highlight the potential of dynamic filtration systems for challenging solid–liquid separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Flow-through compression cell for small-angle and ultra-small-angle neutron scattering measurements

In situ measurements of geological materials under compression and with hydrostatic fluid pressure are important in understanding their behavior under field conditions, which in turn provides critical information for application-driven research. In particular, understanding the role of nano- to micro-scale porosity in the subsurface liquid and gas flow is critical for the high-fidelity characterization of the transport and more efficient extraction of the associated energy resources. In other applications, where parts are produced by the consolidation of powders by compression, the resulting porosity and crystallite orientation (texture) may affect its in-use characteristics. Small-angle neutron scattering (SANS) and ultra SANS are ideal probes for characterization of these porous structures over the nano to micro length scales. Here we show the design, realization, and performance of a novel neutron scattering sample environment, a specially designed compression cell, which provides compressive stress and hydrostatic pressures with effective stress up to 60 MPa, using the neutron beam to probe the effects of stress vectors parallel to the neutron beam. In conclusion, we demonstrate that the neutron optics is suitable for the experimental objectives and that the system is highly stable to the stress and pressure conditions of the measurements.

47 OTHER INSTRUMENTATION↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗