Drying stress formation by inhomogeneous moisture and temperature distribution

The field of moisture concentration and the field of temperature are the two factors that induce stresses in dried materials. The main aim of this paper is to estimate the influence of these two factors on the formation of drying-induced stresses. The considerations are based on the model elaborated by the authors and the analysis of the drying-induced stresses is carried out on a convectively dried prismatic bar.     

Numerical Modelling of Swelling Pressure in Unsaturated Expansive Elasto-Plastic Porous Media

The focus of this work is to provide a new concept for accessing the swelling stress in expansive porous media, especially in highly compacted bentonite. The key to the new approach is the simulation with a chemical swelling model of an infinitesimal volume change followed by a back compaction Process. Free extension is allowed in the first step, to calculate the interlayer porosity change (micro) and the induced volume change potential (macro). The object-oriented FEM simulator GeoSys/RockFlow allows the combination of different processes. The hydro-mechanic/chemical (H²M/C) model takes into consideration two phase flow and deformation, as well as chemical swelling effects. The negative displacements on each boundary after the free extension simulation are taken as Dirichlet boundary conditions of the back compaction problem. The deformation step is simulated in the context of elasto-plasticity using the modified Cam-Clay model. The stresses obtained by back compaction represent the swelling pressure. A 2D example of compacted bentonite is analyzed with the new H²M/C model.     

Pressure Distribution in a Curvilinear Thrust Bearing with One Porous Wall Lubricated by a Bingham Fluid

The influence of wall porosity on the pressure distribution of a Bingham fluid flowing in the clearance of a curvilinear thrust bearing is considered. The formulae expressing the pressure distribution are obtained for two cases, namely: an externally pressurized bearing and a squeeze film bearing. The example of a squeeze film between parallel disks is discussed in detail.     

A Film-Flow Model to Describe Free Water Transport during Drying of a Hygroscopic Capillary Porous Medium

In this work, a two-phase film-flow model in a hygroscopic capillary tube is developed and extended to describe the two-phase capillary viscous transport in a network of parallel capillary tubes in terms of relative permeabilities. This film-flow approach is further considered to predict the longitudinal moisture transport in oak wood during drying. Numerical results obtained from this prediction are compared with data of convective drying experiments performed on samples of this wood. The comparison seems to confirm the physical relevance of a film-flow model to correctly represent the moisture transfer until the hygroscopic regime is reached.     

Heat-Mass-Transport and Thermal Stresses in Porous Charring Materials

The aim of this paper is to develop a method of asymptotic averaging for processes occurring in porous charring materials under high temperatures. The advantage of the method is the ability to calculate not only averaged macrocharacteristics of the processes, namely internal gas generation, filtration and deforming processes, but also microcharacteristics, such as microstresses in phases of charring material, gas velocity in a pore, etc. To determine microcharacteristics, the method allows us to formulate special mathematical problems on a periodic cell. To calculate macrocharacteristics, such as pore pressure of filtrating gas, rate of charring and macrostresses, with the help of asymptotic averaging method, averaged global equations are formulated. Here effective characteristics of porous medium (gas permeability coefficient, rate of charring, elasticity modulus, thermal expansion coefficient) are determined not empirically, as in most works on porous materials, but on the basis of solving the local problems. Solution of these problems over the periodic cell allows us to derive analytically the law of the Darcy type for a gas phase flow in porous media, to obtain an expression for intensive mass transfer between solid and gas phases, to set the form of constitutive relations for charring porous media, and also to calculate microstresses in a vicinity of a growing pore. As an example of solving a global averaged problem, the problem on one-sided high-temperature heating of a plate made of epoxy binder has been solved numerically.     

Coupled Solvent and Heat Transport of a Mixture of Swelling Porous Particles and Fluids: Single Time-Scale Problem

A three-spatial scale, single time-scale model for both moisture and heat transport is developed for an unsaturated swelling porous media from first principles within a mixture theoretic framework. On the smallest (micro) scale, the system consists of macromolecules (clay particles, polymers, etc.) and a solvating liquid (vicinal fluid), each of which are viewed as individual phases or nonoverlapping continua occupying distinct regions of space and satisfying the classical field equations. These equations are homogenized forming overlaying continua on the intermediate (meso) scale via hybrid mixture theory (HMT). On the mesoscale the homogenized swelling particles consisting of the homogenized vicinal fluid and colloid are then mixed with two bulk phase fluids: the bulk solvent and its vapor. At this scale, there exists three nonoverlapping continua occupying distinct regions of space. On the largest (macro) scale the saturated homogenized particles, bulk liquid and vapor solvent, are again homogenized forming four overlaying continua: doubly homogenized vicinal fluid, doubly homogenized macromolecules, and singly homogenized bulk liquid and vapor phases. Two constitutive theories are developed, one at the mesoscale and the other at the macroscale. Both are developed via the Coleman and Noll method of exploiting the entropy inequality coupled with linearization about equilibrium. The macroscale constitutive theory does not rely upon the mesoscale theory as is common in other upscaling methods. The energy equation on either the mesoscale or macroscale generalizes de Vries classical theory of heat and moisture transport. The momentum balance allows for flow of fluid via volume fraction gradients, pressure gradients, external force fields, and temperature gradients.     

Moisture Removal from a Two-Layer Porous Media: A Conceptual Model and Experimental Results

Moisture removal from a two-layer porous media in which air is circulated through one layer and moisture is removed from the second has not been well studied due to the emphasis given to single-layer systems. This two-layer configuration is common in natural and engineered systems and can be used as a means to create a barrier to downward migrating fluids and to remove liquids and gases that may be present in the finer layer. However, there is little data on moisture removal from a two-layer porous media in which air is circulated through one layer parallel to the interface and moisture is removed from the finer second layer by evaporation. A conceptual model of the moisture removal from a two-layer porous media system was developed and compared to experimental moisture removal rates from laboratory scale dry barriers. The limited experimental data agrees well with the results predicted by the conceptual model, providing an initial validation.     

Mixed Convection in a Vertical Porous Channel

A numerical study of mixed convection in a vertical channel filled with a porous medium including the effect of inertial forces is studied by taking into account the effect of viscous and Darcy dissipations. The flow is modeled using the Brinkman–Forchheimer-extended Darcy equations. The two boundaries are considered as isothermal–isothermal, isoflux–isothermal and isothermal–isoflux for the left and right walls of the channel and kept either at equal or at different temperatures. The governing equations are solved numerically by finite difference method with Southwell–Over–Relaxation technique for extended Darcy model and analytically using perturbation series method for Darcian model. The velocity and temperature fields are obtained for various porous parameter, inertia effect, product of Brinkman number and Grashof number and the ratio of Grashof number and Reynolds number for equal and different wall temperatures. Nusselt number at the walls is also determined for three types of thermal boundary conditions. The viscous dissipation enhances the flow reversal in the case of downward flow while it counters the flow in the case of upward flow. The Darcy and inertial drag terms suppress the flow. It is found that analytical and numerical solutions agree very well for the Darcian model. An erratum to this article is available at .     

The Laminar Boundary Layer over a Permeable Wall

An analysis is given of the laminar boundary layer over a permeable/porous wall. The porous wall is passive in the sense that no suction or blowing velocity is imposed. To describe the flow inside and above the porous wall a continuum approach is employed based on the Volume-Averaging Method (S. Whitaker The method of volume averaging). With help of an order-of-magnitude analysis the boundary-layer equations are derived. The analysis is constrained by: (a) a low wall permeability; (b) a low Reynolds number for the flow inside the porous wall; (c) a sufficiently high Reynolds number for the freestream flow above the porous wall. Two boundary layers lying on top of each other can be distinguished: the Prandtl boundary layer above the porous wall, and the Brinkman boundary layer inside the porous wall. Based on the analytical solution for the Brinkman boundary layer in combination with the momentum transfer model of Ochoa-Tapia and Whitaker (Int. J. Heat Mass Transfer 38 (1995) 2635). for the interface region, a closed set of equations is derived for the Prandtl boundary layer. For the stream function a power series expansion in the perturbation parameter is adopted, where is proportional to ratio of the Brinkman to the Prandtl boundary-layer thickness. A generalization of the Falkner–Skan equation for boundary-layer flow past a wedge is derived, in which wall permeability is incorporated. Numerical solutions of the Falkner–Skan equation for various wedge angles are presented. Up to the first order in wall permeability causes a positive streamwise velocity at the interface and inside the porous wall, but a wall-normal interface velocity is a second-order effect. Furthermore, wall permeability causes a decrease in the wall shear stress when the freestream flow accelerates, but an increase in the wall shear stress when the freestream flow decelerates. From the latter it follows that separation, as indicated by zero wall shear stress, is delayed to a larger positive pressure gradient.     

Swelling Induced Finite Strain Flexure in a Rectangular Block of an Isotropic Elastic Material

The deformation of a rectangular block into an annular wedge is studied with respect to the state of swelling interior to the block. Nonuniform swelling fields are shown to generate these flexure deformations in the absence of resultant forces and bending moments. Analytical expressions for the deformation fields demonstrate these effects for both incompressible and compressible generalizations of conventional hyperelastic materials. Existing results in the absence of a swelling agent are recovered as special cases.     

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two-dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left-hand vertical wall has temperature T _h and the right-hand vertical wall is maintained at temperature T _c (T _h > T _c) and the horizontal walls are insulated. At time t > 0, the left-hand vertical wall temperature is suddenly raised to which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.     

Internal Stresses in a Solid with Dislocations

Formulas for calculating internal stresses in a material, generated by continuously distributed dislocations, are found on the basis of the gauge theory of defects. It is shown that internal stresses are selfbalanced and satisfy the equilibrium equations and boundary conditions in the absence of external loads.     

Free Convection in Spherical Annular Sectors Filled with a Porous Medium

Calculated free convection flows and heat transfer are presented for concentric spherical annular sectors, filled with a porous medium. Two isothermal walls and an adiabatic radial wall at the sector angle define the sectors. The governing equations (in the stream function and temperature formulation) are solved numerically using ADI (alternative direction implicit) finite-difference method. Over the range of geometric parameters examined, the obtained results for spherical annuli and low Rayleigh number Ra. As Ra increases, multicellular flows develop for small values of the aspect ratio parameter . In addition, analytical solutions of the governing equations were obtained for small values of Ra (1) and it was shown that these solutions agree well with those obtained numerically. Significant differences in the local heat transfer rates on the inner and outer walls of the spherical annuli were observed from these solutions.     

Determination of Interlaminar Stresses in Bending Problems for a Stack of Transversely Isotropic Plates

A mechanical and mathematical bending model for a stack of transversely isotropic plates is developed. The resolving equations for deflections and tangential displacements are supplemented with a system of differential equations for normal and tangential contact stresses. It is demonstrated that for stacks consisting of an arbitrary number of identical plates with no friction between them, the initial system of equations for contact stresses can be reduced to Helmholtz equations. This transition allows obtaining the complete eigenvalue spectrum for the Laplasian of the problem and, in special cases, eigenfunctions. They are Krylov functions when bending is cylindrical and Bessel functions when bending is axisymmetric     

Mixed Convection Boundary Layer Flow Past a Horizontal Circular Cylinder Embedded in a Bidisperse Porous Medium

Adopting a two-temperature and two-velocity model, appropriate to a bidisperse porous medium (BDPM) proposed by Nield and Kuznetsov (2008), the classical steady, mixed convection boundary layer flow about a horizontal, isothermal circular cylinder embedded in a porous medium has been theoretically studied in this article. It is shown that the boundary layer analysis leads to expressions for the flow and heat transfer characteristics in terms of an inter-phase momentum parameter, a thermal diffusivity ratio, a thermal conductivity ratio, a permeability ratio, a modified thermal capacity ratio, and a buoyancy or mixed convection parameter. The transformed partial differential equations governing the flow and heat transfer in the f-phase (the macro-pores) and the p-phase (the remainder of the structure) are solved numerically using a very efficient implicit finite-difference technique known as Keller-box method. A good agreement is observed between the present results and those known from the open literature in the special case of a traditional Darcy formulation (monodisperse system).     

Temperature and pressure field visualizations in a porous medium dried in superheated steam

This work focuses on heat and mass transfers with phase change in porous media. The experimental analysis was carried out in an aerodynamic return-flow wind tunnel, with very small cylinders of cellular concrete. For the local analysis, the samples were fitted with thermocouples and pressure sensors. A method of temperature and pressure field visualization is developed to highlight the dynamics of transport phenomena. We show two specific phenomena: (1) liquid outflow generated by the overpressure and (2) vaporization of the water inside a two-phase zone that progressively pervades the sample.     

Linear Stability and Convection in a Gravity Modulated Porous Layer Heated from Below: Transition from Synchronous to Subharmonic Solutions

The linear stability theory is used to investigate analytically the effects of gravity modulation on convection in a homogenous porous layer heated from below. The linear stability results are presented for both the synchronous and subharmonic solutions and the exact point for the transition from synchronous to subharmonic solutions is computed. It is also demonstrated that increasing the excitation frequency rapidly stabilizes the convection up to the transition point from synchronous to subharmonic convection. Beyond the transition point, the effect of increasing the frequency is to slowly destabilize the convection.     

Moisture transport in initially fully saturated concrete during drying

Moisture content changes during drying were investigated in the present work. Particular emphasis was placed on the initial stage of drying of saturated concrete, where moisture contents are high. For this stage of drying, experimental data are lacking, and no comprehensive theory exists to describe it. The present investigation was performed experimentally and numerically for drying of cylinders with one exposed end, made of normal weight and lightweight concrete with varying water to cement ratio (w/c). The gravimetric technique was employed to obtain the spatial distribution of moisture content. The experimental results obtained indicate that drying of concrete becomes diffusion controlled when the average moisture content decreases below 70 to 80% of the initial saturation. Typical drying rates are in the order of magnitude of 0.18 kg/day/m² and 0.02 kg/day/m² for the first and the second stage of drying, respectively. The lightweight concrete cylinders as compared to those made of normal weight concrete exhibited higher levels of moisture content throughout the process. At high w/c ratios, the moisture profiles for both types of cylinders, as expected, show steeper changes with time. Large, constant drying rates were observed both experimentally and numerically in the beginning of the drying. The numerical model developed is based on a generalized mathematical formulation for mass and heat transfer in porous media, and its predictions are in agreement with the experimental data within the uncertainty range of the input data.