共查询到20条相似文献，搜索用时 31 毫秒
1.
This study introduces the dispersive fluid flux of total fluid mass to the density-driven flow equation to improve thermohaline modeling of salt and heat transports in porous media. The dispersive fluid flux in the flow equation is derived to account for an additional fluid flux driven by the density gradient and mechanical dispersion. The coupled flow, salt transport and heat transport governing equations are numerically solved by a fully implicit finite difference method to investigate solution changes due to the dispersive fluid flux. The numerical solutions are verified by the Henry problem and the thermal Elder problem under a moderate density effect and by the brine Elder problem under a strong density effect. It is found that increment of the maximum ratio of the dispersive fluid flux to the advective fluid flux results in increasing dispersivity for the Henry problem and the brine Elder problem. The effects of the dispersive fluid flux on salt and heat transports under high density differences and high dispersivities are more noticeable than under low density differences and low dispersivities. Values of quantitative indicators such as the Nusselt number, mass flux, salt mass stored and maximum penetration depth in the brine Elder problem show noticeable changes by the dispersive fluid flux. In the thermohaline Elder problem, the dispersive fluid flux shows a considerable effect on the shape and the number of developed fingers and makes either an upwelling or a downwelling flow in the center of the domain. In conclusion, for the general case that involves strong density-driven flow and transport modeling in porous media, the dispersive fluid flux should be considered in the flow equation. 相似文献
2.
In this paper, a three-dimensional isopycnal approach is presented to simulate the dynamics of fluid mud covering the formation, development, transport, and disappearance of fluid mud. The basic assumption is the assignment of the fluid’s density as the indicating parameter for the rheological behavior. Considering stable stratification, as is usually the case for fluid mud, layers of constant density discretize the vertical domain. The non-Newtonian dynamics of fluid mud is simulated by solving the Cauchy equations for general continuum dynamics. Instead of using a turbulent viscosity approach, the viscosity is allowed to vary according to the rheological behavior of mud suspensions. This apparent viscosity can be determined for different rheological formulations in dependence of the volume solid fraction and the shear rate. An existing three-dimensional isopycnal hydrodynamic model was extended for vertical mass transport processes and was applied on a schematic system with hindered settling. For including the rheological behavior of fluid mud, the Worrall–Tuliani approach was parameterized and implemented. The resulting flow behavior is shown on a model application of fluid mud layers moving down an inclined plane. With these changes, it is demonstrated that the isopycnal model is capable of simulating fluid mud dynamics. 相似文献
3.
The results of 3D numerical simulation of coupled flow and transport processes of light-waste migration in heterogeneous aquifer
with brine are described. The study is based on the geological data and monitoring of deep-well injection of liquid waste
at Chepetsk Mechanical Plant (ChMZ, Glazov). The buoyancy of light waste is simulated for the post-injection period of 300
years using different conceptual models of heterogeneity. The procedure of 3D geostatistical simulation based on transition
probability analysis with TSIM code of the reservoir heterogeneity is presented. The results of the ascending waste migration
simulation using SEAWAT2000 code indicate that both the variable fluid density and the hydraulic heterogeneity must be considered
in the forecast models for determining the possible waste distribution during the injection and the post injection periods. 相似文献
4.
Fluid flow exerts a critical impact on the convection of thermal energy in geological media, whereas heat transport in turn affects fluid properties, including fluid dynamic viscosity and density. The interplay of flow and heat transport also affects solute transport. To unravel these complex coupled flow, heat, and solute transport processes, here, we present a theory for the idealized scale‐dependent Poiseuille flow model considering a constant temperature gradient (?T) along a single fracture, where fluid dynamic viscosity connects with temperature via an exponential function. The idealized scale‐dependent model is validated based on the solutions from direct numerical simulations. We find that the hydraulic conductivity (K) of the Poiseuille flow either increases or decreases with scales depending on ?T > 0°C/m or ?T < 0°C/m, respectively. Indeed, the degree of changes in K depends on the magnitude of ?T and fracture length. The scale‐dependent model provides an alternative explanation for the well‐known scale‐dependent transport problem, for example, the dispersion coefficient increases with travel distance when ?T > 0°C/m according to the Taylor dispersion theory, because K (or equivalently flux through fractures) scales with fracture length. The proposed theory unravels intertwined interactions between flow and transport processes, which might shed light on understanding many practical geophysical problems, for example, geothermal energy exploration. 相似文献
5.
Richard Hoffmann Jean-Christophe Maréchal Adrien Selles Alain Dassargues Pascal Goderniaux 《Ground water》2022,60(2):192-209
Heat as a tracer in fractured porous aquifers is more sensitive to fracture-matrix processes than a solute tracer. Temperature evolution as a function of time can be used to differentiate fracture and matrix characteristics. Experimental hot (50 °C) and cold (10 °C) water injections were performed in a weathered and fractured granite aquifer where the natural background temperature is 30 °C. The tailing of the hot and cold breakthrough curves, observed under different hydraulic conditions, was characterized in a log–log plot of time vs. normalized temperature difference, also converted to a residence time distribution (normalized). Dimensionless tail slopes close to 1.5 were observed for hot and cold breakthrough curves, compared to solute tracer tests showing slopes between 2 and 3. This stronger thermal diffusive behavior is explained by heat conduction. Using a process-based numerical model, the impact of heat conduction toward and from the porous rock matrix on groundwater heat transport was explored. Fracture aperture was adjusted depending on the actual hydraulic conditions. Water density and viscosity were considered temperature dependent. The model simulated the increase or reduction of the energy level in the fracture-matrix system and satisfactorily reproduced breakthrough curves tail slopes. This study shows the feasibility and utility of cold water tracer tests in hot fractured aquifers to boost and characterize the thermal matrix diffusion from the matrix toward the flowing groundwater in the fractures. This can be used as complementary information to solute tracer tests that are largely influenced by strong advection in the fractures. 相似文献
6.
This paper presents a study on suspended particle transport in porous medium with the aid of a sand layer transportation–deposition testing system to determine the kinetic characteristics of particles in porous medium under variable temperatures. Quartz sand and quartz powder were chosen as the porous medium and particle in the tests, respectively. Four size compositions and two operational modes, that is, temperature reduction mode (changing from 18°C to 5°C) and temperature increment mode (changing from 18°C to 35°C), were adopted. The turbidity and concentration of quartz powder were measured under various conditions. We observed a high temperature‐independent correlation between them. Breakthrough curves under different conditions were analysed using this testing system. The results showed that changes in temperature affected the particle transport process to some extent, and the degree of influence was closely related to the time moment of the temperature change onset. Moreover, we found a hysteresis phenomenon in the breakthrough curve under both temperature reduction and increment conditions. The results also indicated that the temperature effect was particularly significant for smaller particles. The typical curves to represent particle transport process under variable temperatures were put forward according to the results. To explain the test results, four factors, that is, water viscosity, adsorption effect, double layer force, and particle kinetic energy, were considered and categorized as promotion or constraining factors. 相似文献
7.
Ward E. Sanford Leonard F. Konikow Gary L. Rowe Jr. Susan L. Brantley 《Journal of Volcanology and Geothermal Research》1995,64(3-4)
Poa´s Volcano is an active stratovolcano in Costa Rica that has a lake in its active crater. The crater lake has high temperatures (50–90 °C), high acidity (pH ≈ 0.0), and a high dissolved-solids content (100 g/kg). The volcano has numerous freshwater springs on its flanks, but a few on the northwestern flank are highly acidic (pH = 1.6–2.5) and have high dissolved-solids concentrations (2–22 g/kg). This study analyzes the regional groundwater system at Poa´s and demonstrates the likelihood that the water discharging from the acidic springs in the Rio Agrio watershed originates at the acidic crater lake. Both heat and solute transport are analyzed on a regional scale through numerical simulations using the HST3D finite-difference model, which solves the coupled equations for fluid flow, heat transport, and solute transport. The code allows fluid viscosity and density to be functions of both temperature and solute concentration. The simulations use estimates for recharge to the mountain and a range of values and various distributions of permeability and porosity. Several sensitivity analyses are performed to test how the uncertainty in many of the model parameters affects the simulation results. These uncertainties yield an estimated range of travel times from the crater lake to the Rio Agrio springs of 1–30 years, which is in close agreement with the results of tritium analyses of the springs. Calculated groundwater fluxes into and out of the crater lake are both about several hundred kg/s. These fluxes must be accounted for in water budgets of the crater lake. 相似文献
8.
In two steady uniform flows at different physical scales in a small open channel, with variables characterizing flow, sediment, and fluid adjusted for dynamic similitude by means of four dimensionless modelling parameters (a Reynolds number, a Froude number, a density ratio, and a length ratio), measured frequency distributions of height, spacing, and migration rate of current ripples were almost identical when scaled, thus verifying that exact Reynolds-Froude modelling of loose-sediment transport is valid and workable. Modelling should be valid as well for a wide range of other transport conditions in the same kind of flow, because no additional kinds of forces or effects would be present in transport of loose grains in modes other than as ripples. In scaled-down modelling, a scale ratio of 2.5 is attainable without recourse to exotic fluids by use of water at 85°C to model natural flows at 10°C. 相似文献
9.
Numerical techniques for subsurface flow and transport modeling are often limited by computational limitations including fine mesh and small time steps to control artificial dispersion. Particle-tracking simulation offers a robust alternative for modeling solute transport in subsurface formations. However, the modeling scale usually differs substantially from the rock measurement scale, and the scale-up of measurements have to be made accounting for the pattern of spatial heterogeneity exhibited at different scales. Therefore, it is important to construct accurate coarse-scale simulations that are capable of capturing the uncertainties in reservoir and transport attributes due to scale-up. A statistical scale-up procedure developed in our previous work is extended by considering the effects of unresolved (residual) heterogeneity below the resolution of the finest modeling scale in 3D. First, a scale-up procedure based on the concept of volume variance is employed to construct realizations of permeability and porosity at the (coarse) transport modeling scale, at which flow or transport simulation is performed. Next, to compute various effective transport parameters, a series of realizations exhibiting detailed heterogeneities at the fine scale, whose domain size is the same as the transport modeling scale, are generated. These realizations are subjected to a hybrid particle-tracking simulation. Probabilistic transition time is considered, borrowing the idea from the continuous time random walk (CTRW) technique to account for any sub-scale heterogeneity at the fine scale level. The approach is validated against analytical solutions and general CTRW formulation. Finally, coarse-scale transport variables (i.e., dispersivities and parameterization of transition time distribution) are calibrated by minimizing the mismatch in effluent history with the equivalent averaged models. Construction of conditional probability distributions of effective parameters is facilitated by integrating the results over the entire suite of realizations. The proposed method is flexible, as it does not invoke any explicit assumption regarding the multivariate distribution of the heterogeneity. In contrast to other hierarchical CTRW formulation for modeling multi-scale heterogeneities, the proposed approach does not impose any length scale requirement regarding sub-grid heterogeneities. In fact, it aims to capture the uncertainty in effective reservoir and transport properties due to the presence of heterogeneity at the intermediate scale, which is larger than the finest resolution of heterogeneity but smaller than the representative elementary volume, but it is often comparable to the transport modeling scale. 相似文献
10.
The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree. 相似文献
11.
Joachim Bartsch 《Continental Shelf Research》1993,13(12)
A model system consisting of a three-dimensional circulation and transport model is used to simulate the dispersal of herring larvae in the North Sea. The driving forces of the circulation model are theM2-tide, time dependent wind stress and air pressure fields, as well as monthly climatological density fields interpolated on a daily basis. The transport model includes advection as well as diffusion and uses a tracer technique to pursue the larvae in the time and space domain. Furthermore, the transport model incorporates a simulation of active vertical movement of the larvae as vertical migration has a marked effect on the drift route of the larvae.Simulated North Sea currents are compared to observed currents gained during the Autumn Circulation Experiment (ACE) which took place from August 1987 to March 1988. The variability in observed currents at the location of the moorings is found to be larger than in the simulated currents. On a larger scale, the simulated circulation in the northern North Sea shows a close correspondence to the inferred circulation from a quasi-synoptic hydrographic survey. The quality of the predictions of larval transport and distribution patterns by the model system is tested using sequential larvae distributions observed during ACE. Common features and discrepancies of observed and simulated distributions are discussed and conclusions for further field investigations and modelling studies are drawn. 相似文献
12.
13.
Simulation of less‐mobile porosity dynamics in contrasting sediment water interface porous media 
下载免费PDF全文
下载免费PDF全文 Farzaneh MahmoodPoor Dehkordy Martin A. Briggs Frederick D. Day‐Lewis Amvrossios C. Bagtzoglou 《水文研究》2018,32(13):2030-2043
Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales. 相似文献
14.
Wolfgang Schmidt 《地球物理与天体物理流体动力学》2013,107(1-2):27-57
Abstract Models of a differentially rotating compressible convection zone are calculated, considering the inertial forces in the poloidal components of the equations of motion. Two driving mechanisms have been considered: latitude dependent heat transport and anisotropic viscosity. In the former case a meridional circulation is induced initially which in turn generates differential rotation, whereas in the latter case differential rotation is directly driven by the anisotropic viscosity, and the meridional circulation is a secondary effect. In the case of anisotropic viscosity the choice of boundary conditions has a big influence on the results: depending on whether or not the conditions of vanishing pressure perturbation are imposed at the bottom of the convection zone, one obtains differential rotation with a fast (≥ 10 ms?1) or a slow (～ 1 ms?1) circulation. In the latter case the rotation law is mainly a function of radius and the rotation rate increases inwards if the viscosity is larger in radial direction than in the horizontal directions. The models with latitude dependent heat transport exhibit a strong dependence on the Prandtl number. For values of the Prandtl number less than 0.2 the pole-equator temperature difference and the surface velocity of the meridional circulation are compatible with observations. For sufficiently small values of the Prandtl number the convection zone becomes globally unstable like a layer of fluid for which the critical Rayleigh number is exceeded. 相似文献
15.
The roles of surface heat flux and ocean heat transport convergence in determining Atlantic Ocean temperature variability 总被引：1，自引：0，他引：1
Jeremy P. Grist Simon A. Josey Robert Marsh Simon A. Good Andrew. C. Coward Beverly A. de Cuevas Steven G. Alderson Adrian L. New Gurvan Madec 《Ocean Dynamics》2010,60(4):771-790
The temperature variability of the Atlantic Ocean is investigated using an eddy-permitting (1/4°) global ocean model (ORCA-025) forced with historical surface meteorological fields from 1958 to 2001. The simulation of volume-averaged temperature and the vertical structure of the zonally averaged temperature trends are compared with those from observations. In regions with a high number of observations, in particular above a depth of 500 m and between 22° N and 65° N, the model simulation and the dataset are in good agreement. The relative contribution of variability in ocean heat transport (OHT) convergence and net surface heat flux to changes in ocean heat content is investigated with a focus on three regions: the subpolar and subtropical gyres and the tropics. The surface heat flux plays a relatively minor role in year-to-year changes in the subpolar and subtropical regions, but in the tropical North Atlantic, its role is of similar significance to the ocean heat transport convergence. The strongest signal during the study period is a cooling of the subpolar gyre between 1970 and 1990, which subsequently reversed as the mid-latitude OHT convergence transitioned from an anomalously weak to an anomalously strong state. We also explore whether model OHT anomalies can be linked to surface flux anomalies through a Hovmöller analysis of the Atlantic sector. At low latitudes, increased ocean heat gain coincides with anomalously strong northward transport, whereas at mid-high latitudes, reduced ocean heat loss is associated with anomalously weak heat transport. 相似文献
16.
二维水质模型在武汉东湖引水工程中的应用 总被引：2，自引：0，他引：2
以"大东湖生态水网"工程的实施为契机,应用MIKE21软件对武汉东湖引水工程进行数值模拟.利用东湖2008年6月和7月实测的水深、TN和TP数据率定了模型参数(床底摩擦力和涡粘系数),并利用2010年7月实测数据建立水动力学模型和对流扩散模型,模拟了四种引水工况下模型运行35 d之后东湖水体中TN、TP浓度分布情况,并比较四种工况的模拟结果.结果表明:在设置两个引水口以及考虑风速影响的情况下,东湖水体中TN、TP浓度分布最均匀.对"大东湖生态水网"工程的实施具有一定的指导意义. 相似文献
17.
D.L. Turcotte 《Earth and Planetary Science Letters》1980,48(1):53-58
The two principal contributions to the surface heat flow of the earth are the cooling of the earth and the heat production of radioactive isotopes. As the rate of heat production decreases with time the temperature of the interior of the earth also decreases. The rate of decrease is determined by the ability of solid-state mantle convection to transport the heat to the surface. The dominant effect is the exponential temperature dependence of the mantle viscosity. The non-dimensional mantle temperature can be parameterised in terms of the Rayleigh number for mantle convection. It is found that the mantle is currently cooling at a rate of 36°K/109 years and that three billion years ago the mean temperature was 150°K higher than it is today; 83% of the present surface heat flow is attributed to the decay of radioactive isotopes and 17% to the cooling of the earth. The corresponding mean concentration of uranium in the mantle is 32 ppb. 相似文献
18.
U.R. Christensen 《Physics of the Earth and Planetary Interiors》1984,35(4):264-282
The case is presented that the efficiency of variable viscosity convection in the Earth's mantle to remove heat may depend only very weakly on the internal viscosity or temperature. An extensive numerical study of the heat transport by 2-D steady state convection with free boundaries and temperature dependent viscosity was carried out. The range of Rayleigh numbers (Ra) is 104?107 and the viscosity contrast goes up to 250000. Although an absolute or relative maximum of the Nusselt number (Nu) is obtained at long wavelength in a certain parameter range, at sufficiently high Rayleigh number optimal heat transport is achieved by an aspect ratio close to or below one. The results for convection in a square box are presented in several ways. With the viscosity ratio fixed and the Rayleigh number defined with the viscosity at the mean of top and bottom temperature the increase of Nu with Ra is characterized by a logarithmic gradient β = ?ln(Nu)/? ln(Ra) in the range of 0.23–0.36, similar to constant viscosity convection. More appropriate for a cooling planetary body is a parameterization where the Rayleigh number is defined with the viscosity at the actual average temperature and the surface viscosity is fixed rather than the viscosity ratio. Now the logarithmic gradient β falls below 0.10 when the viscosity ratio exceeds 250, and the velocity of the surface layer becomes almost independent of Ra. In an end-member model for the Earth's thermal evolution it is assumed that the Nusselt number becomes virtually constant at high Rayleigh number. In the context of whole mantle convection this would imply that the present thermal state is still affected by the initial temperature, that only 25–50% of the present-day heat loss is balanced by radiogenic heat production, and the plate velocities were about the same during most of the Earth's history. 相似文献
19.
Richard W. Healy Sosina S. Haile David L. Parkhurst Scott R. Charlton 《Ground water》2018,56(5):810-815
Variably saturated groundwater flow, heat transport, and solute transport are important processes in environmental phenomena, such as the natural evolution of water chemistry of aquifers and streams, the storage of radioactive waste in a geologic repository, the contamination of water resources from acid‐rock drainage, and the geologic sequestration of carbon dioxide. Up to now, our ability to simulate these processes simultaneously with fully coupled reactive transport models has been limited to complex and often difficult‐to‐use models. To address the need for a simple and easy‐to‐use model, the VS2DRTI software package has been developed for simulating water flow, heat transport, and reactive solute transport through variably saturated porous media. The underlying numerical model, VS2DRT, was created by coupling the flow and transport capabilities of the VS2DT and VS2DH models with the equilibrium and kinetic reaction capabilities of PhreeqcRM. Flow capabilities include two‐dimensional, constant‐density, variably saturated flow; transport capabilities include both heat and multicomponent solute transport; and the reaction capabilities are a complete implementation of geochemical reactions of PHREEQC. The graphical user interface includes a preprocessor for building simulations and a postprocessor for visual display of simulation results. To demonstrate the simulation of multiple processes, the model is applied to a hypothetical example of injection of heated waste water to an aquifer with temperature‐dependent cation exchange. VS2DRTI is freely available public domain software. 相似文献
20.
The principle, that the rate of internal viscous dissipation is at a minimum, is incorrect when temperature and velocity fields are linked through temperature-dependent viscosity or density. This makes it inappropriate for the study of spreading ridge — transform fault systems or other plate-tectonic problems with large viscosity stratification resulting from large temperature gradients. Corrections to the principle are noted for cases without heat advection but with boundary tractions or non-linear materials. 相似文献