A REDUCED MFE FORMULATION BASED ON POD FOR THE NON-STATIONARY CONDUCTION-CONVECTION PROBLEMS

In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.

Conduction-convection coupled heat transfer around a hollow cylinder under different buoyancy forces

In order to clarify the effect of a buoyancy force on conduction–convection coupled heat transfer in a hollow cylinder, the flow and thermal characteristics were analyzed using an RNG k-ε turbulence model. The Reynolds number was fixed at 1.014 × 10^(6), and the Rayleigh number varied from 1.122 × 10^(10)to 1.088 × 10^(11). Results have shown that, when considering the effect of an opposed buoyancy force, increasing the Rayleigh number has a positive impact on the rate of change and uniformity of the cylinder temperature. The temperature distributions along the axial and circumferential directions are similar for different Rayleigh numbers, but extreme values differ.Along the axial direction, the maximum temperature is obtained at the interface between the variable-diameter part and the constant-diameter part. The maximum dimensionless temperature value decreases to 0.12 when the Rayleigh number increases to 1.088 × 10^(11). Along the circumferential direction, the temperature distribution is affected by the buoyancy force, which results in the temperature of the upper part being higher than that of the lower part. After nondimensionalization of the temperature and time, a correlation was proposed to illustrate the transient heat transfer process quantitatively. The standard deviation of the maximum relative temperature, representing the temperature uniformity, was also calculated. It was found that the difference in the direction of the buoyancy force made a huge difference. Compared with the opposed buoyancy force, the maximum dimensionless temperature is almost two times higher with an assisted buoyancy force. Similarly, the heat transfer coefficient with an assisted buoyancy force is half of that with an opposed buoyancy force. Overall, an assisted buoyancy force plays a negative role in terms of thermal characteristics. The flow field around the hollow cylinder was also illustrated to reveal the mechanism of the buoyancy force on magnitude and direction aspects.

A Stabilized Finite Element Method for Non-Stationary Conduction-Convection Problems

This paper is concerned with a stabilized finite element method based on two local Gauss integrations for the two-dimensional non-stationary conductionconvection equations by using the lowest equal-order pairs of finite elements.This method only offsets the discrete pressure space by the residual of the simple and symmetry term at element level in order to circumvent the inf-sup condition.The stability of the discrete scheme is derived under some regularity assumptions.Optimal error estimates are obtained by applying the standard Galerkin techniques.Finally,the numerical illustrations agree completely with the theoretical expectations.

A defect-correction method for unsteady conduction convection problems Ⅰ:spatial discretization

In this paper, a semi-discrete defect-correction mixed finite element method (MFEM) for solving the non-stationary conduction-convection problems in two dimension is presented. In this method, we solve the nonlinear equations with an added artificial viscosity term on a finite element grid and correct this solutions on the same grid using a linearized defect-correction technique. The stability and the error analysis are derived. The theory analysis shows that our method is stable and has a good convergence property.

Comparative Study on Methods for Computing Soil Heat Storage and Energy Balance in Arid and Semi-Arid Areas

Observations collected in the Badan Jaran desert hinterland and edge during 19-23 August 2009 and in the Jinta Oasis during 12-16 June 2005 are used to assess three methods for calculating the heat storage of the5-20-cm soil layer.The methods evaluated include the harmonic method,the conduction-convection method,and the temperature integral method.Soil heat storage calculated using the harmonic method provides the closest match with measured values.The conduction-convection method underestimates nighttime soil heat storage.The temperature integral method best captures fluctuations in soil heat storage on sub-diurnal timescales,but overestimates the amplitude and peak values of the diurnal cycle.The relative performance of each method varies with the underlying land surface.The land surface energy balance is evaluated using observations of soil heat flux at 5-cm depth and estimates of ground heat flux adjusted to account for soil heat storage.The energy balance closure rate increases and energy balance is improved when the ground heat flux is adjusted to account for soil heat storage.The results achieved using the harmonic and temperature integral methods are superior to those achieved using the conduction-convection method.