SEARCH · Engineering Papers

Results for “Diffusion problems”

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,801 records · Page 100

All sources NASA NTRS DOE OSTI

NASA NTRS2004

The AstroBiology Explorer (ABE) Mission Concept

Infrared spectroscopy in the 2.5-16 micron range is a principle means by which organic compounds can be detected and identified in space via their vibrational transitions. Ground-based, airborne, and spaceborne IR spectral studies have already demonstrated that a significant fraction of the carbon in the interstellar medium (ISM) resides in the form of complex organic molecular species. Unfortunately, neither the distribution of these materials nor their genetic and evolutionary relationships with each other or their environments are well understood. The Astrobiology Explorer (ABE) is a MIDEX mission concept currently under study by a team of partners: NASA's Ames Research Center, Ball Aerospace and Technologies Corporation, and the Jet Propulsion Laboratory. ABE will conduct IR spectroscopic observations to address outstanding important problems in astrobiology, astrochemistry, and astrophysics. The core observational program would make fundamental scientific progress in understanding (1) The evolution of ices and organic matter in dense molecular clouds and young forming stellar systems, (2) The chemical evolution of organic molecules in the ISM as they transition from AGB outflows to planetary nebulae to the general diffuse ISM to HII regions and dense clouds, (3) The distribution of organics in the diffuse ISM, (4) The nature of organics in the Solar System (in comets, asteroids, satellites), and (5) The nature and distribution of organics in local galaxies. The technical considerations of achieving these science objectives in a MIDEX-sized mission will be presented.

Sandford, Scott A.↗

DOE OSTI2022

Approximate Bayesian Computation applied to the Diffuse Gamma-Ray Sky

ABSTRACT Many sources contribute to the diffuse gamma-ray background (DGRB), including star forming galaxies, active galactic nuclei, and cosmic ray interactions in the Milky Way. Exotic sources, such as dark matter annihilation, may also make some contribution. The photon counts-in-pixels distribution is a powerful tool for analysing the DGRB and determining the relative contributions of different sources. However, including photon energy information in a likelihood analysis of the counts-in-pixels distribution quickly becomes computationally intractable as the number of source types and energy bins increase. Here, we apply the likelihood-free method of approximate Bayesian computation (ABC) to the problem. We consider a mock analysis that includes contributions from dark matter annihilation in Galactic subhaloes as well as astrophysical backgrounds. We show that our results using ABC are consistent with the exact likelihood when energy information is discarded, and that significantly tighter parameter constraints can be obtained with ABC when energy information is included. ABC presents a powerful tool for analysing the DGRB and understanding its varied origins.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

DOE OSTI2018

A Weighted Least-Squares Transport Equation Compatible with Source Iteration and Voids

We present that second-order forms of the transport equation allow the use of continuous finite elements (CFEMs). This can be desired in multiphysics calculations where other physics require CFEM discretizations. Second-order transport operators are generally self-adjoint, yielding symmetric positive-definite (SPD) matrices, which allow the use of efficient linear algebra solvers with an enormous advantage in memory usage. Least-squares (LS) forms of the transport equation can circumvent the void problems of other second-order forms but are almost always nonconservative. Additionally, the standard LS form is not compatible with discrete ordinates method (SN) iterative solution techniques such as source iteration. A new form of the LS transport equation has recently been developed that is compatible with voids and standard SN iterative solution techniques. Performing nonlinear diffusion acceleration (NDA) using an independently differenced low-order equation enforces conservation for the whole system and makes this equation suitable for reactor physics calculations. In this context, “independent” means that both the transport and low-order solutions converge to the same scalar flux and current as the spatial mesh is refined, but for a given mesh, the solutions are not necessarily equal. In this paper we show that introducing a weight function into this LS equation improves issues with causality and can render our equation equal to the self-adjoint angular flux (SAAF) equation. Causality is a principle of the transport equation that states that information travels only downstream along characteristics. This principle can be violated numerically. We show how to limit the weight function in voids and demonstrate the effect of this limit on accuracy. Using the C5G7 benchmark, we compare our method to the SAAF formulation with a void treatment (SAAFτ) that is not self-adjoint and has a nonsymmetric coefficient matrix. Finallyl, we show that the weighted LS equation with NDA gives acceptable accuracy relative to the SAAFτ equation while maintaining a SPD system matrix.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

Senegalese land surface change analysis and biophysical parameter estimation using NOAA AVHRR spectral data

Surface biophysical estimates were derived from analysis of NOAA Advanced Very High Spectral Resolution (AVHRR) spectral data of the Senegalese area of west Africa. The parameters derived were of solar albedo, spectral visible and near-infrared band reflectance, spectral vegetative index, and ground temperature. Wet and dry linked AVHRR scenes from 1981 through 1985 in Senegal were analyzed for a semi-wet southerly site near Tambacounda and a predominantly dry northerly site near Podor. Related problems were studied to convert satellite derived radiance to biophysical estimates of the land surface. Problems studied were associated with sensor miscalibration, atmospheric and aerosol spatial variability, surface anisotropy of reflected radiation, narrow satellite band reflectance to broad solar band conversion, and ground emissivity correction. The middle-infrared reflectance was approximated with a visible AVHRR reflectance for improving solar albedo estimates. In addition, the spectral composition of solar irradiance (direct and diffuse radiation) between major spectral regions (i.e., ultraviolet, visible, near-infrared, and middle-infrared) was found to be insensitive to changes in the clear sky atmospheric optical depth in the narrow band to solar band conversion procedure. Solar albedo derived estimates for both sites were not found to change markedly with significant antecedent precipitation events or correspondingly from increases in green leaf vegetation density. The bright soil/substrate contributed to a high albedo for the dry related scenes, whereas the high internal leaf reflectance in green vegetation canopies in the near-infrared contributed to high solar albedo for the wet related scenes. The relationship between solar albedo and ground temperature was poor, indicating the solar albedo has little control of the ground temperature. The normalized difference vegetation index (NDVI) and the derived visible reflectance were more sensitive to antecedent rainfall amounts and green vegetation changes than were near-infrared changes. The information in the NDVI related to green leaf density changes primarily was from the visible reflectance. The contribution of the near-infrared reflectance to explaining green vegetation is largely reduced when there is a bright substrate.

Vukovich, Fred M.↗

NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

In this study, several variable turbulent Prandtl number formulations are examined for boundary layers, pipe flow, and axisymmetric jets. The model formulations include simple algebraic relations between the thermal diffusivity and turbulent viscosity as well as more complex models that solve transport equations for the thermal variance and its dissipation rate. Results are compared with available data for wall heat transfer and profile measurements of mean temperature, the root-mean-square (RMS) fluctuating temperature, turbulent heat flux and turbulent Prandtl number. For wall-bounded problems, the algebraic models are found to best predict the rise in turbulent Prandtl number near the wall as well as the log-layer temperature profile, while the thermal variance models provide a good representation of the RMS temperature fluctuations. In jet flows, the algebraic models provide no benefit over a constant turbulent Prandtl number approach. Application of the thermal variance models finds that some significantly overpredict the temperature variance in the plume and most underpredict the thermal growth rate of the jet. The models yield very similar fluctuating temperature intensities in jets from straight pipes and smooth contraction nozzles, in contrast to data that indicate the latter should have noticeably higher values. For the particular low subsonic heated jet cases examined, changes in the turbulent Prandtl number had no effect on the centerline velocity decay.

computational fluid dynamics↗

DOE OSTI2016

Group Structure Optimization Using the PyGroup Code

Group constants are used as coefficients in the energy domain discretized Linear Boltzmann Equation (LBE), in which the Continuous-Energy (CE) fluxes are also integrated over energy bins and transformed into group fluxes. This procedure, based on the Multi-Group (MG) theory, is widely applied in radiation transport and diffusion simulations Careful preparation of MG cross section data, or group constants, is thus a prerequisite for success in acquiring high-fidelity LBE solutions. Generation of MG cross section data itself is a complicated multi-step task, since group constants are indispensably coupled with the solution of the LBE, i.e. the flux spectrum in the problem to be solved. Theoretically, a MG cross section library is only valid for the problem that it is built upon. Nevertheless, in practice MG cross section libraries are prepared for similar type of problems, such as thermal reactors, fast reactors, etc., to facilitate the neutronic simulations. Generally, CE point-wise cross section data set, generated by both theoretical models and experimental data, is first converted to Fine-Group (FG) libraries using a pre-defined or precalculated flux spectra, then, Broad-Group (BG) libraries are collapsed from the fine-group libraries using problem-dependent flux spectra. Proper treatments of energy and spatial self-shielding and resonance cross section are required in the CE-to-MG converting process Previously, we have developed a computer code, YGROUP, for fine-to-broad group cross section collapsing. To maintain high accuracy of the resulting BG libraries, YGROUP provides users several options on collapsing strategies and weighting schemes for group binning and evaluation of BG cross section data in users' specific problems.. An approach based on Particle Swarm Optimization was applied to automatically determine the optimal group structure. We also developed a hybrid weighting scheme to further improve BG library performance for Special Nuclear Materials (SNM) detection applications. In this work, a new computer code, PyGroup, has been developed to implement the PSO approach with an improved simulation-driven fitness function design that includes both response and forward/adjoint flux errors. Together with YGROUP, and deterministic transport code TITAN/PENMSHXP, the PyGroup code package is tested on a Dual Range Coincident Counter (DRCC) model. In this paper, we demonstrated the capability of the PyGroup code for group structure optimization using the DRCC detector model. In the future, an iteration scheme can be implemented to determine the number of broad groups based on a user-defined tolerance. we will also apply PyGroup to more complicated reactor models, and investigate physics-based fitness function designs to further reduce computational cost. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

A mathematical and numerical framework for the analysis of compressible thermal convection in gases at very high temperatures

The relevance of non-equilibrium phenomena, nonlinear behavior, gravitational effects and fluid compressibility in a wide range of problems related to high-temperature gas-dynamics, especially in thermal, mechanical and nuclear engineering, calls for a concerted approach using the tools of the kinetic theory of gases, statistical physics, quantum mechanics, thermodynamics and mathematical modeling in synergy with advanced numerical strategies for the solution of the Navier–Stokes equations. The reason behind such a need is that in many instances of relevance in this field one witnesses a departure from canonical models and the resulting inadequacy of standard CFD approaches, especially those traditionally used to deal with thermal (buoyancy) convection problems. Starting from microscopic considerations and typical concepts of molecular dynamics, passing through the Boltzmann equation and its known solutions, we show how it is possible to remove past assumptions and elaborate an algorithm capable of targeting the broadest range of applications. Moving beyond the Boussinesq approximation, the Sutherland law and the principle of energy equipartition, the resulting method allows most of the fluid properties (density, viscosity, thermal conductivity, heat capacity and diffusivity, etc.) to be derived in a rational and natural way while keeping empirical contamination to the minimum. Special attention is deserved as well to the well-known pressure issue. With the application of the socalled multiple pressure variables concept and a projection-like numerical approach, difficulties with such a term in the momentum equation are circumvented by allowing the hydrodynamic pressure to decouple from its thermodynamic counterpart. The final result is a flexible and modular framework that on the one hand is able to account for all the molecule (translational, rotational and vibrational) degrees of freedom and their effective excitation, and on the other hand can guarantee adequate interplay between molecular and macroscopic-level entities and processes. Performances are demonstrated by computing some incompressible and compressible benchmark test cases for thermal (gravitational) convection, which are then extended to the high-temperature regime taking advantage of the newly developed features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

Finite difference methods for reducing numerical diffusion in TEACH-type calculations

A methodological evaluation for two-finite differencing schemes for computer-aided gas turbine design is presented. The two computational schemes include; a Bounded Skewed Finite Differencing Scheme (BSUDS); and a Quadratic Upwind Differencing Scheme (QSDS). In the evaluation, the derivations of the schemes were incorporated into two-dimensional and three-dimensional versions of the Teaching Axisymmetric Characteristics Heuristically (TEACH) computer code. Assessments were made according to performance criteria for the solution of problems of turbulent, laminar, and coannular turbulent flow. The specific performance criteria used in the evaluation were simplicity, accuracy, and computational economy. It is found that the BSUDS scheme performed better with respect to the criteria than the QUDS. Some of the reasons for the more successful performance BSUDS are discussed.

Syed, S. A.↗

NASA NTRS1992

[NASA/DOD Aerospace Knowledge Diffusion Research Project. Paper 18:] Scientific and Technical Information (STI) policy and the competitive position of the US aerospace industry

With its contribution to trade, its coupling with national security, and its symbolism of U.S. technological strength, the U.S. aerospace industry holds a unique position in the Nation's industrial structure. Federal science and technology policy and Federal scientific and technical information (STI) policy loom important as strategic contributions to the U.S. aerospace industry's leading competitive position. However, three fundamental policy problems exist. First, the United States lacks a coherent STI policy and a unified approach to the development of such a policy. Second, policymakers fail to understand the relationship of STI to science and technology policy. Third, STI is treated as a part of general information policy, without any recognition of its uniqueness. This paper provides an overview of the Federal information policy structure as it relates to STI and frames the policy issues that require resolution.

Hernon, Peter↗

NASA NTRS1992

NASA/DOD Aerospace Knowledge Diffusion Research Project. Paper 18: Scientific and Technical Information (STI) policy and the competitive position of the US aerospace industry

Hernon, Peter↗

DOE OSTI2016

Multi-scale diffuse interface modeling of multi-component two-phase flow with partial miscibility

In this paper, we introduce a diffuse interface model to simulate multi-component two-phase flow with partial miscibility based on a realistic equation of state (e.g. Peng–Robinson equation of state). Because of partial miscibility, thermodynamic relations are used to model not only interfacial properties but also bulk properties, including density, composition, pressure, and realistic viscosity. As far as we know, this effort is the first time to use diffuse interface modeling based on equation of state for modeling of multi-component two-phase flow with partial miscibility. In numerical simulation, the key issue is to resolve the high contrast of scales from the microscopic interface composition to macroscale bulk fluid motion since the interface has a nanoscale thickness only. To efficiently solve this challenging problem, we develop a multi-scale simulation method. At the microscopic scale, we deduce a reduced interfacial equation under reasonable assumptions, and then we propose a formulation of capillary pressure, which is consistent with macroscale flow equations. Moreover, we show that Young–Laplace equation is an approximation of this capillarity formulation, and this formulation is also consistent with the concept of Tolman length, which is a correction of Young–Laplace equation. At the macroscopical scale, the interfaces are treated as discontinuous surfaces separating two phases of fluids. Our approach differs from conventional sharp-interface two-phase flow model in that we use the capillary pressure directly instead of a combination of surface tension and Young–Laplace equation because capillarity can be calculated from our proposed capillarity formulation. A compatible condition is also derived for the pressure in flow equations. Furthermore, based on the proposed capillarity formulation, we design an efficient numerical method for directly computing the capillary pressure between two fluids composed of multiple components. Finally, numerical tests are carried out to verify the effectiveness of the proposed multi-scale method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗

DOE OSTI2016

Moisture ingress into electronics enclosures under isothermal conditions

The number of electronics used in outdoor environment is constantly growing. The humidity causes about 19 % of all electronics failures and, especially, moisture increases these problems due to the ongoing process of miniaturization and lower power consumption of electronic components. Moisture loads are still not understood well by design engineers, therefore this field has become one of the bottlenecks in the electronics system design. The objective of this paper is to model moisture ingress into an electronics enclosure under isothermal conditions. The moisture diffusion model is based on a 1D quasi-steady state (QSS) approximation for Fick’s second law. This QSS approach is also described with an electrical analogy which gives a fast tool in modelling of the moisture response. The same QSS method is applied to ambient water vapour variations. The obtained results are compared to an analytical solution and very good agreement is found.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Reacceleration of Galactic Cosmic Rays beyond the Knee at the Termination Shock of a Cosmic-Ray-driven Galactic Wind

The origin of cosmic rays (CRs) above the knee in the spectrum is an unsolved problem. We present a wind model in which interstellar gas flows along a nonrotating, expanding flux tube with a changing speed and cross-sectional area. CRs from Galactic sources, such as supernova remnants, which are coupled to the plasma via Alfvén waves, provide the main pressure source for driving this outflow. These CRs are then subject to diffusive shock reacceleration at the Galactic wind termination shock, which is located at a distance ~200 kpc. Some of the highest-energy reaccelerated particles propagate upstream against the wind and can contribute to the petaelectronvolt to exaelectronvolt range of the spectrum. We analyze the conditions under which efficient reacceleration can occur and find that rigidities ~10–40 PV can be obtained and that the termination shock may account for between 10% and 50% of the proton spectrum measured in IceCube/IceTop experiment. In our model, the termination shock is unable to fully explain the CR spectrum in the petaelectronvolt to exaelectronvolt range. The highest-energy particles that escape downstream from our termination shock, and similar shocks surrounding most galaxies, can be further accelerated by intergalactic shock fronts.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Direct Coupling Method for Time-Accurate Solution of Incompressible Navier-Stokes Equations

A noniterative finite difference numerical method is presented for the solution of the incompressible Navier-Stokes equations with second order accuracy in time and space. Explicit treatment of convection and diffusion terms and implicit treatment of the pressure gradient give a single pressure Poisson equation when the discretized momentum and continuity equations are combined. A pressure boundary condition is not needed on solid boundaries in the staggered mesh system. The solution of the pressure Poisson equation is obtained directly by Gaussian elimination. This method is tested on flow problems in a driven cavity and a curved duct.

Soh, Woo Y.↗

NASA NTRS1990

Evidence against a simple two-component model for the far-infrared emission from galaxies

Two of the first Infrared Astronomy Satellite (IRAS) results were that galaxies have a wide range of values for the ratio of 60 micron to 100 micron flux density (0.2 less than or equal to S sub 60/S sub 100 less than or equal to 1.0) and that this ratio is correlated with L sub fir, L sub b, L sub fir being the total far-infrared luminosity and L sub b being the luminosity at visible wavelengths (de Jong et al. 1984; Soifer et al. 1984). From these results arose the following simple model for the far-infrared emission from galaxies (de Jong et al. 1984), which has remained the standard model ever since. In this model, the far-infrared emission comes from two dust components: warm dust (T approx. equals 50 K) intermingled with, and heated by, young massive OB stars in molecular clouds and HII regions, and colder dust (T approx. equals 20 K) associated with the diffuse atomic hydrogen in the interstellar medium and heated by the general interstellar radiation field. As the number of young stars in a galaxy increases, S sub 60/S sub 100 increases, because there is a greater proportion of warm dust, and so does L sub fir/L sub b, because most of the radiation from the young stars is absorbed by the dust, leading to a swifter increase in far-infrared emission than in visible light. Although this model explains the basic IRAS results, it is inelegant. It uses two free parameters to fit two data (the 60 and 100 micron flux densities)-and there are now several observations that contradict it. Despite these major problems with the two-component model, it is not clear what should be put in its place. When considering possible models for the far-infrared emission from galaxies, the observational evidence for our own galaxy must be considered. Researchers suspect that the study by Boulanger and Perault (1988) of the far-infrared properties of the local interstellar medium may be particularly relevant. They showed that molecular clouds are leaky - that most of the light from OB stars in molecular clouds does not heat the dust in the clouds, but instead leaks out. The consequence of this is that that while most of the far-infrared emission from the solar neighborhood is from dust associated with diffuse HI, this dust is mostly heated by young stars.

Eales, Stephen A.↗

DOE OSTI2017

Pore-network model of evaporation-induced salt precipitation in porous media: The effect of correlations and heterogeneity

Salt transport and precipitation in porous media constitute a set of complex and fascinating phenomena that are of considerable interest to several important problems, ranging from storage of CO 2 in geological formations, to soil fertility, and protection of pavements and roads, as well as historical monuments. The phenomena occur at the pore scale and are greatly influenced by the heterogeneity of the pore space morphology. In this paper, we present a pore-network (PN) model to study the phenomena. Vapor diffusion, capillary effect at the brine-vapor interface, flow of brine, and transport of salt and its precipitation in the pores that plug the pores partially or completely are all accounted for. The drying process is modeled by the invasion percolation, while transport of salt in brine is accounted for by the convective-diffusion equation. We demonstrate that the drying patterns, the clustering and connectivity of the pore throats in which salt precipitation occurs, the saturation distribution, and the drying rate are all strongly dependent upon the pore-size distribution, the correlations among the pore sizes, and the anisotropy of the pore space caused by stratification that most natural porous media contain. In particular, if the strata are more or less parallel to the direction of injection of the gas that dries out the pore space (air, for example) and/or causes salt precipitation (CO 2 , for example), the drying rate increases significantly. Moreover, salt tends to precipitate in clusters of neighboring pores that are parallel to the open surface of the porous medium.

42 ENGINEERING↗

Refine your search to explore more results.