Numerical Analysis of Coupled Heat and Mass Transfer Phenomena in Concrete at Elevated Temperatures

In this work, we present a novel methodology for incorporating the effect of fibre surface morphology on liquid water transport in polymer electrolyte membrane fuel cell gas diffusion layers (GDLs). Roughness features presented on the surface of the fibre are analysed using atomic force microscopy and are found to significantly impact the capillary pressure of liquid water pathways propagating through the GDL. A threshold capillary pressure was defined as the largest capillary pressure exhibited by the liquid water phase during the invasion of the throat. The threshold capillary pressures observed in the presence of roughness features are significantly greater than those in the absence of roughness features. Two-dimensional circumferential roughness models in cylindrical and converging-diverging throats are established, and an interfacial meniscus advancing algorithm is presented to determine the resulting threshold capillary pressures required for liquid water penetration. Revised Young–Laplace equations, which are particularly useful for pore network modeling, are suggested for calculating threshold capillary pressures that account for the effect of the roughness of throats with intrinsic contact angles greater than \(90^{\circ }\).     

Driving Potentials of Heat and Mass Transport in Porous Building Materials: A Comparison Between General Linear,Thermodynamic and Micromechanical Derivation Schemes

In systems of coupled transport processes the question of the appropriate driving potentials is a point of discussion. In this article, three different approaches to derive models for transport currents are systematically compared. According to a general linear approach, an arbitrary full set of independent state variables and material properties is sufficient to describe any transport current. This approach is derived here from a symmetry principle. Thermodynamic and micromechanical approaches are more complex and even less general, but they allow additional statements about the transport coefficients and they reduce the number of transport processes. In the thermodynamic approach the additional information stems from the calculation of the entropy production rate; the micromechanical approach involves a microphysical model of the considered porous system. As a practical example, the three derivation schemes are applied to the often-encountered case of non-hysteretic heat and moisture transport in homogeneous building materials. It is shown, how the general state variables of a porous system are reduced to only two. Then from the general linear approach it can be seen, that all equations for the moisture transport current using a main driving potential (e.g. moisture content, vapour pressure, chemical potential) and an independent secondary driving potential (e.g. temperature, liquid pressure) are equivalent, without recurrence to the thermodynamic or micromechanical approach. However, the transport coefficients are arbitrary phenomenological functions depending on the two state variables. Based on a literature survey it is shown, which additional statements can be made in the thermodynamic and in the micromechanical approach. The latter yields the pressure-driven model (vapour and liquid pressure as the two driving potentials). Finally it is shown, what is to be expected, if in more complex systems the number of state variables increases.     

Application of a 3-D FEM/FDM Overlapping Scheme to Heat and Mass Transfer and Moving-Boundary Problems

A 3-D FEM/FDM overlapping scheme for viscous, incompressible flow problems is presented that combines the finite element method, which is best suited to analyze flow in any arbitrarily shaped flow geometry, with the finite difference method, which is advantageous in both computing time and computer storage. The combination of both methods enables large-scale viscous flow to be analyzed, which is crucial both for detailed analysis of 3-D flows and for solving flow problems around moving bodies, A modified ABM AC method is used as the basic algorithm, to which a sophisticated time integration scheme, proposed by the present authors, has been applied. In this paper, some numerical results including 3-D heat and mass transfer problem and moving-boundary problems are presented.     

Contribution to tip leakage loss modeling in radial turbines based on 3D flow analysis and 1D characterization

The characterization of tip leakage flow plays an important role for one-dimensional loss modeling and design in radial turbine research. Tip leakage losses can be expressed as function of fluid momentum and mass flow passing through the tip gap. Friction-driven flow and contrariwise oriented pressure gradient-driven flow are highly coupled. However, these numbers are mostly unknown and dependent on tip gap geometry and turbine running condition. Based on a commonly used definition of a non-dimensional tip leakage momentum ratio, a novel correlation has been derived. This allows a consistent characterization for variable tip gap sizes over a wide range of operating conditions. The correlation has been validated by means of CFD data with high variety in reduced speed tip gap geometry and expansion ratios. Results of the novel number show significant improvements of quantitative and qualitative results over a wide range of running conditions in comparison to existing correlations. Furthermore, correlations for tip leakage velocities, that can easily be used in one-dimensional models, have been derived. Finally, it has been demonstrated, that the influence of inlet flow momentum on the tip leakage flow can be analyzed with presented correlations.     

Solutions using the Recursive Operator Method for the Dynamic Equations of Heat and Mass Transfer in Hardening Concrete under Heat Treatment in a Heating Chamber

A recursive operator method was used to obtain a general solution of the system of linearized differential equations of heat and mass transfer in hardening concrete in a heating chamber. The solution contains arbitrary analytic functions determined from the boundary and initial conditions.     

On A. D. Kovalenko's Research Works on the Thermomechanics of Coupled Fields in Materials and Structural Members and Its Further Development

A brief review is given to the research works of Academician A. D. Kovalenko that laid the foundations of the thermomechanics of coupled fields in inelastic materials and structural members under monoharmonic loading. The developments of Kovalenko's ideas by the associates of the Thermoelasticity Department founded by him are pointed out __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 9, pp. 16–25, September 2005.     

3D Micromechanical Modeling of Failure and Damage Evolution in Dual Phase Steel Based on a Real 2D Microstructure

A plate of dual phase steel was produced from low carbon steel with intercritical annealing treatment. Its optically determined surface microstructure was utilized to construct three different microstructural models. To describe the ductile damage in the ferritic matrix,the Gurson-Tvergaard-Needleman model was used with the failure in the martensite phase being ignored. The numerical results obtained for the mechanism of void initiation and coalescence were compared with the experimental observations. The numerical results obtained from the randomly extruded 3D model showed a significantly better agreement with the experimental ones than those obtained from the 2D model or the uniformly extruded 3D model.     

Investigation of the Influence of Different Heterogeneous Recombination Mechanisms on the Heat Fluxes to a Catalytic Surface in Dissociated Carbon Dioxide

A kinetic model of heterogeneous recombination in dissociated carbon dioxide on high-temperature heat-shield coatings is developed; the model takes into account the nonequilibrium adsorption-desorption reactions of oxygen atoms and their recombination in the Eley-Rideal and Langmuir-Hinshelwood reactions. On the basis of a comparison of the calculated heat fluxes in dissociated carbon dioxide with those measured in the VGU-3 plasma generator of the Institute for Problems in Mechanics of the Russian Academy of Sciences (IPM RAS) and the available literature data, the parameters of the catalysis model are chosen for the glassy coating of the Buran orbiter tile heat shield based on the SiO₂–B₂O₃–SiB₄ system. The effects of heterogeneous recombination proceeding in accordance with the Langmuir-Hinshelwood mechanism, as well as the processes involving carbon atoms and those involving physically adsorbed oxygen atoms, on the heat fluxes to the glassy coating are analyzed on the surface temperature range from 300 to 2000 K.     

Upscaling and Simulation of Waterflooding in Heterogeneous Reservoirs Using Wavelet Transformations: Application to the SPE-10 Model

We have developed an accurate and highly efficient method for upscaling and simulation of immiscible displacements in three-dimensional (3D) heterogeneous reservoirs, which is an extension of the technique that we developed previously for 2D systems. The method utilizes wavelet transformations (WTs) to upscale the geological model of a reservoir, based on the spatial distribution of the single-phase permeabilities and the locations of the wells in the reservoir. It generates a non-uniform grid in which the resolved structure of the fine grid around the wells, as well as in the high-permeability sectors, are preserved, but the rest of the grid is upscaled. A robust uplayering procedure is used to reduce the number of the layers, and the WTs are used to upscale each layer areally. To demonstrate the method’s accuracy and efficiency, we have applied it to the geological model of a highly heterogeneous reservoir put forward in the tenth Society of Petroleum Engineers comparative solution project (the SPE-10 model), and carried out simulation of waterflooding in the upscaled model. Various upscaling scenarios were examined, and although some of them resulted in efficient simulations and accurate predictions, the results when non-uniform upscaling is used based on the WT technique are in excellent agreement with the solution of the same problem in the fine grid of the SPE-10 model. Most importantly, the speed-up factors that we obtain are several orders of magnitude. Hence, the method renders it unnecessary to use massively parallel computations for such problems.     

Drag of Bodies of Revolution in Supersonic Flow with Heat and Mass Addition to the Near Wake

It is shown experimentally that the base drag of bodies of revolution in Mach 1.15 to 3.0 flow can be completely eliminated using special techniques for injecting hydrogen and the products of incomplete combustion of pyrotechnic compounds into the near wake. The experimental data obtained are generalized.     

Simulation of Radon Transport through Building Materials: Influence of the Water Content on Radon Exhalation Rate

Radon transport in porous materials is strongly influenced by the presence of water. It is also necessary to be able to numerically control the effects of this parameter. The radon concentration and radon exhalation rate have been determined by simulation in various building materials containing an increasing water content. It has been proved that the presence of water does not involve the same variations in the concentration on the surface of the medium, according to its porosity. For porous media with low porosity like concrete or granite, ( < 0.2), the radon concentration and radon exhalation rate sharply increase with water until the volumetric water content becomes higher than 30%. At this point, radon emanation plays an important role, in relation to the molecular diffusion process. For materials with medium porosities (e.g., limestone, brick, cement: 0.3 < < 0.45), the concentration was observed to increase up to a volumetric water content of about 10% and then decreased from there. In this case, the molecular diffusion has a greater effect due to a greater quantity of pores in the material. For a small water content, this parameter tends to make the radon concentration decrease at the surface of the medium. These simulations have been compared with experimental analysis and are in strong accordance with the experimental results.     

Modelling Surface Runoff and Infiltration of Rain by an Elliptic-Parabolic Equation Coupled with a First‐Order Equation on the Boundary

. (Accepted April 9, 1998)     

Heat and Mass Transfer on the Mixed Convection of Non-Newtonian Fluids Over a Vertical Wedge with Soret/Dufour Effects and Internal Heat Generation: Variable Wall Temperature/Concentration

Transport in Porous Media - The presented work compares the mechanical behavior from standard unconfined compressive strength and indirect tensile strength tests of natural sandstone and artificial...     

Investigation on Heat and Mass Transfer in a Evaporator of a Capillary-Pumped Loop with the Lattice Boltzmann Method: Pore Scale Simulation

A two-dimensional transient numerical model based on the lattice Boltzmann method (LBM) for the global evaporator of a capillary-pumped loop (CPL) is proposed to describe heat and mass transfer with evaporation in the porous wick, heat conduction in the cover plate, and heat transfer in the vapor groove. To indicate the stochastic phase distribution characteristics of most porous wick, the quartet structure generation set (QSGS) is introduced for generating more realistic microstructures of porous media. By using the present lattice Boltzmann algorithm along with the porous structure, the heat and mass transfer of an evaporator on pore scale can be predicted without resorting to any empirical parameters determined case by case. The energy equations for entire evaporator are solved as a conjugate problem, which are solved by means of a spatially varying relaxation time in the lattice Boltzmann model and the liquid flow is driven via the interfacial mass flux. A convective boundary condition considering the latent heat during the evaporation on the interface is introduced into the lattice Boltzmann model based on the nonequilibrium extrapolation rule. Especially, the bounce-back rule and the equilibrium rule of the LBM are, respectively, introduced to deal with the momentum boundary conditions inside the porous wick and on the evaporation interface in order to ensure the stability and the efficiency of the LBM model. Numerical results corresponding to different working conditions and different working fluids are presented, which provide guidance for the evaporator design of a CPL system.     

LES Modeling of the Impact of Heat Losses and Differential Diffusion on Turbulent Stratified Flame Propagation: Application to the TU Darmstadt Stratified Flame

Common combustion chambers often exhibit turbulent flames propagating in partially-premixed mixtures. This propagation is generally governed by aerodynamics, unsteady mixing and chemical processes and may also be affected by conductive heat losses when the reactive zone develops close to the burner lips. The Filtered TAbulated Chemistry for Large Eddy Simulation (F-TACLES) model has been recently developed to include tabulated chemistry in Large Eddy Simulation (LES) of adiabatic stratified flames in flamelet regimes. The present article proposes a modeling approach to account for both differential diffusion and non-adiabatic effects on flame consumption speed following the F-TACLES formalism. The adiabatic F-TACLES model is first detailed using a generalized formalism for diffusive fluxes allowing either to account for differential diffusion or not. The F-TACLES model is then extended to non-adiabatic situations. A correction factor based on the non-adiabatic consumption rate is introduced to recover a realistic filtered flame consumption speed. The objective is here to tackle flame stabilization mechanisms when heat losses affect the reaction zone. The proposed approach is validated through the simulation of the unconfined stratified turbulent jet flame TSF-A for which stabilization process is affected by heat losses. Five simulations are performed for both adiabatic and non-adiabatic flow conditions comparing unity Lewis number and complex diffusion assumptions. The adiabatic F-TACLES model predicts a flame anchored at the burner lip disagreeing with experimental data. The non-adiabatic simulation exhibits local extinction due to heat losses near the burner exit. The flame is then lifted improving the comparison with experiments. Results also show a significant impact of molecular diffusion model on both mean flame consumption rate and angle.     

The Three-Stage Model Based on Strain Strength Distribution for the Tensile Failure Process of Rock and Concrete Materials

A three-stage model is introduced to describe the tensile failure process of rock and concrete materials.Failure of the material is defined to contain three stages in the model,which include elastic deformation stage,body damage stage and localization damage stage.The failure mode change from uniform body damage to localization damage is expressed.The heterogeneity of material is described with strain strength distribution.The fracture factor and intact factor,defined as the distribution function of strain strength,are used to express the fracture state in the failure process.And the distributive parameters can be determined through the experimental stress-strain curve.