Search USGS⌕ Search

SEARCH · Search USGS

Results for “Numerical heat transfer”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 19 recordsLinked to original sources
USGS2016

The water content of recurring slope lineae on Mars

Observations of recurring slope lineae (RSL) from the High-Resolution Imaging Science Experiment have been interpreted as present-day, seasonally variable liquid water flows; however, orbital spectroscopy has not confirmed the presence of liquid H 2 O, only hydrated salts. Thermal Emission Imaging System (THEMIS) temperature data and a numerical heat transfer model definitively constrain the amount of water associated with RSL. Surface temperature differences between RSL-bearing and dry RSL-free terrains are consistent with no water associated with RSL and, based on measurement uncertainties, limit the water content of RSL to at most 0.5–3 wt %. In addition, distinct high thermal inertia regolith signatures expected with crust-forming evaporitic salt deposits from cyclical briny water flows are not observed, indicating low water salinity (if any) and/or low enough volumes to prevent their formation. Alternatively, observed salts may be preexisting in soils at low abundances (i.e., near or below detection limits) and largely immobile. These RSL-rich surfaces experience ~100 K diurnal temperature oscillations, possible freeze/thaw cycles and/or complete evaporation on time scales that challenge their habitability potential. The unique surface temperature measurements provided by THEMIS are consistent with a dry RSL hypothesis or at least significantly limit the water content of Martian RSL.

Geophysical Research Letters↗
USGS1985

Algorithm to reduce approximation error from the complex-variable boundary-element method applied to soil freezing.

An algorithm is presented for the numerical solution of the Laplace equation boundary-value problem, which is assumed to apply to soil freezing or thawing. The Laplace equation is numerically approximated by the complex-variable boundary-element method. The algorithm aids in reducing integrated relative error by providing a true measure of modeling error along the solution domain boundary. This measure of error can be used to select locations for adding, removing, or relocating nodal points on the boundary or to provide bounds for the integrated relative error of unknown nodal variable values along the boundary.

Numerical heat transfer↗
USGS2015

Temperature profile around a basaltic sill intruded into wet sediments

The transfer of heat into wet sediments from magmatic intrusions or lava flows is not well constrained from field data. Such field constraints on numerical models of heat transfer could significantly improve our understanding of water–lava interactions. We use experimentally calibrated pollen darkening to measure the temperature profile around a basaltic sill emplaced into wet lakebed sediments. It is well known that, upon heating, initially transparent palynomorphs darken progressively through golden, brown, and black shades before being destroyed; however, this approach to measuring temperature has not been applied to volcanological questions. We collected sediment samples from established Miocene fossil localities at Clarkia, Idaho. Fossils in the sediments include pollen from numerous tree and shrub species. We experimentally calibrated changes in the color of Clarkia sediment pollen and used this calibration to determine sediment temperatures around a Miocene basaltic sill emplaced in the sediments. Results indicated a flat temperature profile above and below the sill, with T > 325 °C within 1 cm of the basalt-sediment contact, near 300 °C at 1–2 cm from the contact, and ~ 250 °C at 1 m from the sill contact. This profile suggests that heat transport in the sediments was hydrothermally rather than conductively controlled. This information will be used to test numerical models of heat transfer in wet sediments on Earth and Mars.

Journal of Volcanology and Geothermal Research↗
USGS1985

Locating CVBEM collocation points for steady state heat transfer problems

The Complex Variable Boundary Element Method or CVBEM provides a highly accurate means of developing numerical solutions to steady state two-dimensional heat transfer problems. The numerical approach exactly solves the Laplace equation and satisfies the boundary conditions at specified points on the boundary by means of collocation. The accuracy of the approximation depends upon the nodal point distribution specified by the numerical analyst. In order to develop subsequent, refined approximation functions, four techniques for selecting additional collocation points are presented. The techniques are compared as to the governing theory, representation of the error of approximation on the problem boundary, the computational costs, and the ease of use by the numerical analyst.

Engineering Analysis↗
USGS1989

On numerical modeling of one-dimensional geothermal histories

Numerical models of one-dimensional geothermal histories are one way of understanding the relations between tectonics and transient thermal structure in the crust. Such models can be powerful tools for interpreting geochronologic and thermobarometric data. A flexible program to calculate these models on a microcomputer is available and examples of its use are presented. Potential problems with this approach include the simplifying assumptions that are made, limitations of the numerical techniques, and the neglect of convective heat transfer. ?? 1989.

Computers & Geosciences↗
USGS2021

A numerical model for the cooling of a lava sill with heat pipe effects

Understanding the cooling process of volcanic intrusions into wet sediments is a difficult but important problem, given the presence of extremely large temperature gradients and potentially complex water-magma interactions. This report presents a numerical model to study such interactions, including the effect of heat pipes on the cooling of volcanic intrusions. Udell (1985) has shown that heat pipes may develop in heated saturated granular media under laboratory conditions. In previous work, Baker and others (2015) calculated temperatures in the vicinity of a volcanic sill that intruded into wet sediment, showing an unexpected temperature profile in which peak temperatures remained near constant over a region extending a meter above and below the sill. This is challenging to explain with conduction or convection heating methods but is predicted if the heat transfer is performed primarily by a heat pipe. We have numerically modeled the cooling of a lava sill under similar circumstances, using the experimental findings of Udell (1985) to estimate the characteristics of the heat pipe. We have constructed a model using Microsoft C#.NET, complete with an intuitive graphical user interface. The model is available from the U.S. Geological Survey and is capable of being run on Microsoft Windows 7 and higher with modest hardware. We find that the resulting overall temperature profile has some key similarities to the profile inferred by Baker and others (2015). Future models including more detailed convective heat transfer physics will be necessary to fully reproduce the effects of boiling in sediments.

Techniques and Methods↗
USGS2015

Implications of climate change on the heat budget of lentic systems used for power station cooling: Case study Clinton Lake, Illinois

We use a numerical model to analyze the impact of climate change--in particular higher air temperatures--on a nuclear power station that recirculates the water from a reservoir for cooling. The model solves the hydrodynamics, the transfer of heat in the reservoir, and the energy balance at the surface. We use the numerical model to (i) quantify the heat budget in the reservoir and determine how this budget is affected by the combined effect of the power station and climate change and (ii) quantify the impact of climate change on both the downstream thermal pollution and the power station capacity. We consider four different scenarios of climate change. Results of simulations show that climate change will reduce the ability to dissipate heat to the atmosphere and therefore the cooling capacity of the reservoir. We observed an increase of 25% in the thermal load downstream of the reservoir, and a reduction in the capacity of the power station of 18% during the summer months for the worst-case climate change scenario tested. These results suggest that climate change is an important threat for both the downstream thermal pollution and the generation of electricity by power stations that use lentic systems for cooling.

Illinois↗
USGS1989

A numerical solution for the diffusion equation in hydrogeologic systems

The documentation of a computer code for the numerical solution of the linear diffusion equation in one or two dimensions in Cartesian or cylindrical coordinates is presented. Applications of the program include molecular diffusion, heat conduction, and fluid flow in confined systems. The flow media may be anisotropic and heterogeneous. The model is formulated by replacing the continuous linear diffusion equation by discrete finite-difference approximations at each node in a block-centered grid. The resulting matrix equation is solved by the method of preconditioned conjugate gradients. The conjugate gradient method does not require the estimation of iteration parameters and is guaranteed convergent in the absence of rounding error. The matrixes are preconditioned to decrease the steps to convergence. The model allows the specification of any number of boundary conditions for any number of stress periods, and the output of a summary table for selected nodes showing flux and the concentration of the flux quantity for each time step. The model is written in a modular format for ease of modification. The model was verified by comparison of numerical and analytical solutions for cases of molecular diffusion, two-dimensional heat transfer, and axisymmetric radial saturated fluid flow. Application of the model to a hypothetical two-dimensional field situation of gas diffusion in the unsaturated zone is demonstrated. The input and output files are included as a check on program installation. The definition of variables, input requirements, flow chart, and program listing are included in the attachments. (USGS)

Water-Resources Investigations Report↗
USGS1995

Duration of mineralization and fluid-flow history of the Upper Mississippi Valley zinc-lead district

Studies of fluid inclusions in sphalerite and biomarkers from the Upper Mississippi Valley zinc district show homogenization temperatures to be primarily between 90 and 150 ° C, yet show relatively low levels of thermal maturity. We use numerical calculations to simulate fluid and heat flow through fracture-controlled ore zones and heat transfer to the adjacent rocks. Combining a best-fit path through fluid-inclusion data with measured thermal alteration of biomarkers, we calculated the time interval during which mineralizing fluids circulated through the Upper Mississippi Valley district to be on the order of 200 ka. Cambrian and Ordovician aquifers underlying the district, principally the St. Peter and Mt. Simon Sandstones, were the source of the mineralizing fluid. The duration of mineralization thus reflects the fluid-flow history of these regional aquifers.

Geology↗
USGS2020

Timescales and processes of methane hydrate formation and breakdown, with application to geologic systems

Gas hydrate is an ice-like form of water and low molecular weight gas stable at temperatures of roughly -10ºC to 25ºC and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one-sixth of Earth’s methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean-atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean-atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth’s near-surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some fields studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

Journal of Geophysical Research-Solid Earth↗
USGS1987

Large-scale volcano-ground ice interactions on Mars

The process of volcano-ground ice interaction on Mars is investigated by thermodynamic calculations and observations of Viking Orbiter images. We develop a numerical model of volcano-ground ice interaction that includes heat transport by conduction, radiation from the surface, heat transfer to the atmosphere, and H 2 O phase changes in an ice-rich permafrost. We consider eruption of lava flows over permafrost, and intrusion of sills into permafrost. For eruption of lava over permafrost, most of the heat in the flow is lost by radiation and atmospheric effects. The amount of H 2 O liquid and vapor produced is small, and its removal would not be sufficient to cause collapse that would lower the surface of the lava flow below the surrounding terrain. For intrusion of a sill, most of the heat in the sill eventually goes into H 2 O phase changes, producing much larger amounts of water that could have profound geomorphic and geochemical effects. Approximate meltwater discharge rates are calculated for both extrusive and intrusive interactions. We examine two large regions of large-scale volcano-ground ice interactions. Near Aeolis Mensae, intrusion of a complex of dikes and sills into ice-rich ground has produced substantial melting, with mobilization and flow of material. This interaction probably also produced large quantities of palagonite tuff and breccia. Morphologic evidence for progressive fluidization implies that meltwater was stored beneath the surface for some time, and that most of the release of water and volcanic mudflow took place late in the interaction. Northeast of Hellas, several large channels emanate from the area near the volcano Hadriaca Patera. If genetically related to the volcanic activity, large collapse features at the sources of some channels must have originated due to heat from large buried magma bodies. A channel emerging directly from the base of Hadriaca Patera may have originated from release of heat from thick extruded material. Other small channels in the region results from heat released from surface lava flows. Inferred channel discharges may be compared to discharge rates calculated for lava-ground ice interactions. Such comparisons show that meltwater probably accumulated beneath the surface and then was released rapidly, with a discharge rate limited by soil permeability. Volcano-ground ice interaction has been a widespread and important geologic process on Mars, and may be the primary source of palagonites making up the ubiquitous Martian dust.

Icarus↗
USGS1997

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗
USGS1996

Fluid inclusions and biomarkers in the Upper Mississippi Valley zinc-lead district; implications for the fluid-flow and thermal history of the Illinois Basin

The Upper Mississippi Valley zinc-lead district is hosted by Ordovician carbonate rocks at the northern margin of the Illinois Basin. Fluid inclusion temperature measurements on Early Permian sphalerite ore from the district are predominantly between 90?C and I50?C. These temperatures are greater than can be explained by their reconstructed burial depth, which was a maximum of approximately 1 km at the time of mineralization. In contrast to the temperatures of mineral formation derived from fluid inclusions, biomarker maturities in the Upper Mississippi Valley district give an estimate of total thermal exposure integrated over time. Temperatures from fluid inclusions trapped during ore genesis with biomarker maturities were combined to construct an estimate of the district's overall thermal history and, by inference, the late Paleozoic thermal and hydrologic history of the Illinois Basin. Circulation of groundwater through regional aquifers, given sufficient flow rates, can redistribute heat from deep in a sedimentary basin to its shallower margins. Evidence for regional-scale circulation of fluids is provided by paleomagnetic studies, regionally correlated zoned dolomite, fluid inclusions, and thermal maturity of organic matter. Evidence for igneous acti vity contemporaneous with mineralization in the vicinity of the Upper Mississippi Valley district is absent. Regional fluid and heat circulation is the most likely explanation for the elevated fluid inclusion temperatures (relative to maximum estimated burial depth) in the Upper Mississippi Valley district. One plausible driving mechanism and flow path for the ore-forming fluids is groundwater recharge in the late Paleozoic Appalachian-Ouachita mountain belt and northward flow through the Reelfoot rift and the proto- Illinois Basin to the Upper Mississippi Valley district. Warm fluid flowing laterally through Cambrian and Ordovician aquifers would then move vertically upward through the fractures that control sphalerite mineralization in the Upper Mississippi Valley district. Biomarker reactant-product measurements on rock extracts from the Upper Mississippi Valley district define a relatively low level ofthermal maturity for the district, 0.353 for sterane and 0.577 for hopane. Recently published kinetic constants permit a time-temperature relationship to be determined from these biomarker maturities. Numerical calculations were made to simulate fluid heat flow through the fracture-controlled ore zones of the Thompson-Temperly mine and heat transfer to the adjacent rocks where biomarker samples were collected. Calculations that combine the fluid inclusion temperatures and the biomarker constraints on thermal maturity indicate that the time interval during which mineralizing fluids circulated through the Upper Mississippi Valley district is on the order of 200,000 years. Fluid inclusion measurements and thermal maturities from biomarkers in the district reflect the duration of peak temperatures resulting from regional fluid circulation. On the basis of thermal considerations, the timing of fluorite mineralization in southern Illinois, and the northward-decreasing pattern of fluorine enrichment in sediments, we hypothesize that the principal flow direction was northward through the Cambrian and Ordovician aquifers of the Illinois Basin. A basin-scale flow system would result in mass transport (hydrocarbon migration, transport of metals in solution) and energy (heat) transport, which would in turn drive chemical reactions (for example, maturation of organic matter, mineralization, diagenetic reactions) within the Illinois Basin and at its margins.

Bulletin↗
USGS1980

The influences of land use and land cover on climate; an analysis of the Washington-Baltimore area that couples remote sensing with numerical simulation

The Sun drives the atmospheric heat engine by warming the terrestrial surface which in turn warms the atmosphere above. Climate, therefore, is significantly controlled by complex interaction of energy flows near and at the terrestrial surface. When man alters this delicate energy balance by his use of the land, he may alter his climatic environment as well. Land use climatology has emerged as a discipline in which these energy interactions are studied; first, by viewing the spatial distributions of their surface manifestations, and second, by analyzing the energy exchange processes involved. Two new tools for accomplishing this study are presented: one that can interpret surface energy exchange processes from space, and another that can simulate the complex of energy transfers by a numerical simulation model. Use of a satellite-borne multispectral scanner as an imaging radiometer was made feasible by devising a gray-window model that corrects measurements made in space for the effects of the atmosphere in the optical path. The simulation model is a combination of mathematical models of energy transfer processes at or near the surface. Integration of these two analytical approaches was applied to the Washington-Baltimore area to coincide with the August 5, 1973, Skylab 3 overpass which provided data for constructing maps of the energy characteristics of the Earth's surface. The use of the two techniques provides insights into the relationship of climate to land use and land cover and in predicting alterations of climate that may result from alterations of the land surface.

Professional Paper↗
USGS2003

Constitutive relations associated with the Mott-Smith distribution function

It is shown that the distribution function assumed by Mott-Smith determines a unique relation between heat flux, stress, and fluid velocity given by q = (3/2) τ u, i.e., it provides a constitutive relation for heat flux, and it also determines a simple expression for this ratio of third-order central moments Q = (C 3 x) / C x C 2 . These expressions allow the equation of transfer for c x2 to be cast in a form that yields a nonlinear constitutive relation for stress. The results obtained from the Mott-Smith ansatz are compared with the theory of Baganoff and Nathenson and results from a numerical solution of the Boltzmann equation for shock-wave structure obtained by Hicks and Yen.

Physics of Fluids↗
USGS2010

Laboratory investigations of steam flow in a porous medium

Experiments were carried out in the laboratory to test a theory of transient flow of pure steam in a uniform porous medium. This theory is used in modeling pressure transient behavior in vapor dominated geothermal systems. Transient, superheated steam flow experiments were run by bringing a cylinder of porous material to a uniform initial pressure and then making a step increase in pressure at one end of the sample while monitoring the pressure transient breakthrough at the other end. It was found in experiments run at 100°, 125°, and 146°C that the time required for steam pressure transients to propagate through an unconsolidated material containing sand, silt, and clay was 10–25 times longer than predicted by conventional superheated steam flow theory. It is hypothesized that the delay in the steam pressure transient was caused by adsorption of steam in the porous sample. In order to account for steam adsorption, a sink term was included in the conservation of mass equation. In addition, energy transfer in the system has to be considered because latent heat is released when steam adsorption occurs, increasing the sample temperature by as much as 10°C. Finally, it was recognized that the steam pressure was a function of both the temperature and the amount of adsorption in the sample. This function was assumed to be an equilibrium adsorption isotherm, which was determined by experiment. By solving the modified mass and energy equations numerically, subject to the empirical adsorption isotherm relationship, excellent theoretical simulation of the experiments was achieved.

Water Resources Research↗
USGS1982

Laboratory investigations of the physics of steam flow in a porous medium

Experiments were carried out in the laboratory to test a theory of transient flow of pure steam in a uniform porous medium. This theory is used extensively in modeling pressure-transient behavior in vapor-dominated geothermal systems. Transient, superheated steam-flow experiments were run by bringing a cylinder of porous material to a uniform initial pressure, and then making a step increase in pressure at one end of the sample, while monitoring the pressure-transient breakthrough at the other end. It was found in experiments run at 100?, 125?, and 146?C that the time required for steam-pressure transients to propagate through an unconsolidated material containing sand, silt, and clay was 10 to 25 times longer than predicted by theory. It is hypothesized that the delay in the steam-pressure transient was caused by adsorption of steam in the porous sample. In order to account for steam adsorption, a sink term was included in the conservation of mass equation. In addition, energy transfer in the system has to be considered because latent heat is released when steam adsorption occurs, increasing the sample temperature by as much as 10?C. Finally, it was recognized that the steam pressure was a function of both the temperature and the amount of adsorption in the sample. For simplicity, this function was assumed to be in equilibrium adsorption isotherm, which was determined by experiment. By solving the modified mass and energy equations numerically, subject to the empirical adsorption isotherm relationship, excellent theoretical simulation of the experiments was achieved. The experiments support the hypothesis that adsorption of steam can strongly influence steam pressure-transient behavior in porous media; the results suggest that the modified steam-flow theory, which includes steam adsorption terms, should be used in modeling steam flow in vapor-dominated geothermal systems.

Open-File Report↗
USGS2012

Experimental determination of soil heat storage for the simulation of heat transport in a coastal wetland

Two physical experiments were developed to better define the thermal interaction of wetland water and the underlying soil layer. This information is important to numerical models of flow and heat transport that have been developed to support biological studies in the South Florida coastal wetland areas. The experimental apparatus consists of two 1.32 m diameter by 0.99 m tall, trailer-mounted, well-insulated tanks filled with soil and water. A peat–sand–soil mixture was used to represent the wetland soil, and artificial plants were used as a surrogate for emergent wetland vegetation based on size and density observed in the field. The tanks are instrumented with thermocouples to measure vertical and horizontal temperature variations and were placed in an outdoor environment subject to solar radiation, wind, and other factors affecting the heat transfer. Instruments also measure solar radiation, relative humidity, and wind speed. Tests indicate that heat transfer through the sides and bottoms of the tanks is negligible, so the experiments represent vertical heat transfer effects only. The temperature fluctuations measured in the vertical profile through the soil and water are used to calibrate a one-dimensional heat-transport model. The model was used to calculate the thermal conductivity of the soil. Additionally, the model was used to calculate the total heat stored in the soil. This information was then used in a lumped parameter model to calculate an effective depth of soil which provides the appropriate heat storage to be combined with the heat storage in the water column. An effective depth, in the model, of 5.1 cm of wetland soil represents the heat storage needed to match the data taken in the tank containing 55.9 cm of peat/sand/soil mix. The artificial low-density laboratory sawgrass reduced the solar energy absorbed by the 35.6 cm of water and 55.9 cm of soil at midday by less than 5%. The maximum heat transfer into the underlying peat–sand–soil mix lags behind maximum solar radiation by approximately 2 h. A slightly longer temperature lag was observed between the maximum solar radiation and maximum water temperature both with and without soil.

Journal of Hydrology↗