2015 ◽  
Vol 137 (7) ◽  
Author(s):  
James White

In order to better manage computational requirements in the study of thermal conduction with short-scale heterogeneous materials, one is motivated to arrange the thermal energy equation into an accurate and efficient form with averaged properties. This should then allow an averaged temperature solution to be determined with a moderate computational effort. That is the topic of this paper as it describes the development using multiple-scale analysis of an averaged thermal energy equation based on Fourier heat conduction for a heterogeneous material with isotropic properties. The averaged energy equation to be reported is appropriate for a stationary or moving solid and three-dimensional heat flow. Restrictions are that the solid must display its heterogeneous properties over short spatial and time scales that allow averages of its properties to be determined. One distinction of the approach taken is that all short-scale effects, both moving and stationary, are combined into a single function during the analytical development. The result is a self-contained form of the averaged energy equation. By eliminating the need for coupling the averaged energy equation with external local problem solutions, numerical solutions are simplified and made more efficient. Also, as a result of the approach taken, nine effective averaged thermal conductivity terms are identified for three-dimensional conduction (and four effective terms for two-dimensional conduction). These conductivity terms are defined with two types of averaging for the component material conductivities over the short-scales and in terms of the relative proportions of the short-scales. Numerical results are included and discussed.


1980 ◽  
Vol 20 (06) ◽  
pp. 501-507 ◽  
Author(s):  
M.H. Lee ◽  
L.D. Roberts

Abstract In a fracture acidizing treatment the acid reacts with the fracture faces. This acid/rock reaction generates heat that causes the acid temperature itself to increase. To predict accurately the temperature profile and acid spending rate of acid traveling down a hydraulically created fracture, this heat must be considered.Since the heat generated by reaction depends on the reaction rate, the thermal energy equation must be coupled with the acid spending equation. A model has been developed that, for the first time, examines the effect of the heat of reaction on fluid temperature and acid penetration in a fracture. Some sample calculations have also been made to illustrate the effects of the most important parameters on acid penetration in a fracture. Introduction Acid hydraulic fracturing is a common method of stimulating a reservoir. Acid selectively reacts with, and dissolves, portions of the fracture wall so that a finite fluid conductivity remains when the well is returned to production. An important aim in designing such fracturing treatments is determining the distance that live acid will penetrate down the hydraulically induced fracture. This distance is usually called the acid penetration distance and is essential to estimate the production improvement from a given treatment.Because of its importance in predicting stimulation ratio, acid penetration in fractures has been studied by numerous investigators. They assumed the temperature in the fracture was uniform. In real fractures, however, the temperature will vary from the wellbore to the tip of the fracture. Therefore, the assumption of constant temperature seems to be an oversimplification.Whitsitt and Dysart were among the first to study the temperature distribution in a fracture. They constructed a model but it could be applied only to a nonreacting fluid flowing in a fracture because the heat generated by an acid/rock reaction was not considered. In a fracture acidizing treatment, the acid is reacting with the rock walls. This acid/rock reaction generates heat, which causes the acid temperature itself to increase. To predict accurately the temperature profile along the fracture, this heat also must be considered. A model has been developed that, for the first time, examines the effect of the heat of reaction on fluid temperature and acid penetration distance. Mathematical Development The mathematical model is a modification of that introduced by Whitsitt and Dysart to allow for the heat of reaction in the energy-balance equation. Since the heat generated by the acid reaction also depends on the reaction rate, the thermal-energy equation is coupled with the mass-balance equation. These two equations must be solved simultaneously .The model for acid spending in a fractures is illustrated in Fig. 1. The fluid leakoff velocity Vw is assumed constant over the fracture length. Assuming steady-state flow in a vertical fracture and constant fluid properties, the mass-balance equation for acid flowing in a fracture is ................(1) SPEJ P. 501^


2009 ◽  
Vol 131 (11) ◽  
Author(s):  
James Sucec

An equilibrium thermal wake strength parameter is developed for a two-dimensional turbulent boundary layer flow and is then used in the combined thermal law of the wall and the wake to give an approximate temperature profile to insert into the integral form of the thermal energy equation. After the solution of the integral x momentum equation, the integral thermal energy equation is solved for the local Stanton number as a function of position x for accelerating turbulent boundary layers. A simple temperature distribution in the thermal “superlayer” is part of the present modeling. The analysis includes a dependence of the hydrodynamic and thermal wake strengths on the momentum thickness and enthalpy thickness Reynolds numbers, respectively. An approximate dependence of the turbulent Prandtl number, in the “log” region, on the strength of the favorable pressure gradient is proposed and incorporated into the solution. The resultant solution for the Stanton number distribution in accelerated turbulent flows is compared with experimental data in the literature. A comparison of the present predictions is also made to a finite difference solution, which uses the turbulent kinetic energy—turbulent dissipation model of turbulence, for a few cases of accelerating flows.


2016 ◽  
Vol 138 (6) ◽  
Author(s):  
James White

Numerical solution of heat conduction in a heterogeneous material with small spatial and time scales can lead to excessive compute times due to the dense computational grids required. This problem is avoided by averaging the energy equation over the small-scales, which removes the appearance of the short spatial and time scales while retaining their effect on the average temperature. Averaging does, however, increase the complexity of the resulting thermal energy equation by introducing mixed spatial derivatives and six different averaged conductivity terms for three-dimensional analysis. There is a need for a numerical method that efficiently and accurately handles these complexities as well as the other details of the averaged thermal energy equation. That is the topic of this paper as it describes a numerical solution for the averaged thermal energy equation based on Fourier conduction reported recently in the literature. The solution, based on finite difference techniques that are second-order time-accurate and noniterative, is appropriate for three-dimensional time-dependent and steady-state analysis. Speed of solution is obtained by spatially factoring the scheme into an alternating direction sequence at each time level. Numerical stability is enhanced by implicit algorithms that make use of the properties of tightly banded matrices. While accurately accounting for the nonlinearity introduced into the energy equation by temperature-dependent properties, the numerical solution algorithm requires only the consideration of linear systems of algebraic equations in advancing the solution from one time level to the next. Computed examples are included and compared with those for a homogeneous material.


2006 ◽  
Vol 14 (1) ◽  
pp. 27-40 ◽  
Author(s):  
Damian W.I. Rouson ◽  
Xiaofeng Xu ◽  
Karla Morris

The memory management rules for abstract data type calculus presented by Rouson, Morris & Xu [15] are recast as formal statements in the Object Constraint Language (OCL) and applied to the design of a thermal energy equation solver. One set of constraints eliminates memory leaks observed in composite overloaded expressions with three current Fortran 95/2003 compilers. A second set of constraints ensures economical memory recycling. The constraints are preconditions, postconditions and invariants on overloaded operators and the objects they receive and return. It is demonstrated that systematic run-time assertion checking inspired by the formal constraints facilitated the pinpointing of an exceptionally hard-to-reproduce compiler bug. It is further demonstrated that the interplay between OCL's modeling capabilities and Fortran's programming capabilities led to a conceptual breakthrough that greatly improved the readability of our code by facilitating operator overloading. The advantages and disadvantages of our memory management rules are discussed in light of other published solutions [11,19]. Finally, it is demonstrated that the run-time assertion checking has a negligible impact on performance.


2000 ◽  
Vol 64 (1) ◽  
pp. 1-11 ◽  
Author(s):  
L. FERRARIO

The hydrodynamic expansion of plasmas produced by laser beams focused on thin-foil targets is studied using a self-similar approach. The model is useful for studying long-scale-length plasmas, produced by uniform laser irradiation, which become fully underdense during the laser pulse. An essentially qualitative analysis of the hydrodynamics system in Cartesian geometry is given. In particular, a useful expression is obtained by integration of the thermal energy equation.


2005 ◽  
Vol 127 (2) ◽  
pp. 159-164 ◽  
Author(s):  
James Sucec

The combined thermal law of the wall and wake is used as the approximating sequence for the boundary layer temperature profile to solve an integral thermal energy equation for the local Stanton number distribution. The velocity profile in the turbulent boundary layer was taken to be the combined law of the wall and wake of Coles. This allows the solution of an integral form of the x-momentum equation to give the skin friction coefficient distribution. This, along with the velocity profile, is used to solve the thermal energy equation using inner coordinates. The strength of the thermal wake was found by analysis of earlier research results, in the literature, for equilibrium, constant property, turbulent boundary layers. Solutions for the Stanton number distribution with position are found for some adverse pressure gradient boundary layers as well as for those having zero pressure gradient. The zero pressure gradient results cover both fully heated plates and those with unheated starting lengths, including both isothermal surfaces and constant flux surfaces. Comparison of predictions of the present work is made with experimental data in the literature.


2001 ◽  
Vol 44 (12) ◽  
pp. 2379-2382 ◽  
Author(s):  
A.M.W. Wojcik ◽  
J.C. Jansen ◽  
Th. Maschmeyer

Sign in / Sign up

Export Citation Format

Share Document