Assumption Of Thermal Equilibrium Research Articles

Existence condition for detonations in condensed explosives with pressure–temperature equilibrium models

Most engineering tools dealing with detonation waves in condensed explosives are based on the reactive Euler equations, which involve a thermodynamic closure based on temperature and pressure equilibrium conditions. Indeed, the reactant and the detonation products are assumed to evolve in mechanical and thermal equilibrium. Although the assumption of thermal equilibrium is physically questionable, the reactive Euler equations are used for convenience. This choice is supported by several reasons. When considering thermal equilibrium, the assessment of thermal exchanges between the reactant and the detonation products is simplified, as no exchange coefficient needs to be determined. Furthermore, the reactive Euler equations are in conservative form and the shock relations are well-defined. This model appears to be the simplest option for addressing detonations in condensed explosives. However, this simplicity has limitations and conceals a subtle difficulty that can lead to pathological detonations. The present contribution investigates this difficulty and presents a global exothermic condition to preserve exothermic detonations in the frame of temperature–pressure equilibrium flow models. This condition depends on the equations of state only. It is independent of the kinetics used to model the decomposition of the explosive. By meticulously selecting the thermodynamic parameters of the equations of state for the reactant and detonation products, pathological detonations are prevented, thereby maintaining a globally exothermic reaction consistent with the Zeldovich–von Neumann–Döring theory.

Calculation of thermodynamic properties and transport coefficients of Ar/N2–H2–Si plasma

Abstract The composition, thermodynamic properties and transport coefficients of Ar – H 2 – Si and N 2 – H 2 – Si plasma within a temperature range of 300–30 000 K and pressure range of 0.1–10 atm are calculated under the assumptions of local thermal equilibrium (LTE) and local chemical equilibrium (LCE). Taking Debye–Hückel corrections into account, the chemical equilibrium composition and thermodynamic properties of these two plasma systems are derived using the mass action law and classical statistical thermodynamics respectively. The transport coefficients, including viscosity, conductivity, and thermal conductivity, are calculated using the Chapman–Enskog (C–E) method extended to a third-order approximation (second-order for viscosity and heavy particle translational thermal conductivity). Some of the results have been compared with those of other researchers, and there is a good level of agreement. The slight difference arises from the selection of interaction potential. The final calculation reveals that the introduction of silicon vapor significantly alters the thermodynamic properties and transport coefficients of Ar / N 2 – H 2 – Si plasma, even at a small concentration of silicon vapor (1%), the effect on the electrical conductivity cannot be ignored. Furthermore, our calculation results provide the fundamental data for numerical simulations of magnetohydrodynamics (MHD) for the synthesis of silicon nanoparticles and silicon composites.

Non-equilibrium thermal models of lithium batteries

Temperature fluctuations impact battery performance, safety, and health. Industry-standard cell-level models of these phenomena ignore thermal gradients within the electrodes’ active material, i.e., assume the latter to be in “thermal equilibrium”. We present a “non-equilibrium” thermal model that explicitly accounts for spatial variability of temperature with the active material (and the carbon-binder domain). We investigate the conditions, expressed in terms of the heat-generation rate and the thermal properties of a cell’s liquid (electrolyte) and solid (active material and CBD) phases, under which the thermal equilibrium assumption breaks down and our model should be used instead. The differences between these two thermal models are investigated further by coupling them with an industry-standard electrochemical model. The resulting thermal–electrochemical model demonstrates the importance of thermal gradients within the active material at high C-rates (discharge current densities) and for large grain sizes. Under these conditions, the equilibrium assumption underestimates internal temperature by as much as 50%. These two thermal models are then applied to a commercial NMC battery with multiple units. Our non-equilibrium model predicts the battery surface temperature that is in good agreement with measurements, while the equilibrium model underestimates the observed temperature.

Modelling and simulation on compressible multi-component gas-liquid flows with chemical reaction and phase transition effects

In this paper, mathematical models and numerical algorithm are conducted for simulating compressible two-phase multi-component reactive flows, coupling with the finite-rate-based chemistry and the instantaneous-equilibrium phase transition model. Distinctive from previous multi-component two-phase governing equations, all volume fraction transport equations, with non-conservative characteristics, of gaseous components are substituted with a thermal equilibrium assumption inside the gaseous mixture for closure. Moreover, based on the temperature polynomial (NASA-7) fitted gas properties and the stiffened gas equation of state (SG-EOS) described liquid properties, the instantaneous-equilibrium phase transition model is re-derived, and the solving strategy of the transient phase transition model is detailly provided using an algebraical approach rather than the iterative process. Having included these elements, the proposed model is validated on three canonical tests, first for pure reactive flows, second for pure two-phase flows, and third for complex compressible multi-component two-phase reactive flows with phase transition effects: hydrogen-oxygen detonation interacting with a water droplet. The cases have well demonstrated the proposed model capability and the expected fluid physics in both the multiphase and reactive flows are reproduced.

Linear and weakly nonlinear analyses of double‐diffusive convection in porous media with chemical reaction using LTNE model

AbstractThe onset of convection in a horizontal porous layer with chemical reaction and local thermal nonequilibrium is investigated. The nondimensional governing equations have been solved using the normal mode technique, which results in an eigenvalue problem. The analytical expressions for both stationary and oscillatory Rayleigh numbers are obtained. The effect of different parameters has been investigated and presented. The amplitude equation is derived using weakly nonlinear theory. Nusselt number is calculated using an amplitude equation to investigate heat transport. When modeling a fluid‐saturated porous medium, previous research on double‐diffusive convection has uniformly operated under the assumption of local thermal equilibrium (LTE) between the fluid and solid phases at all points within the medium. This standard practice assumes a minimal temperature gradient between the phases at any given location. However, in practical scenarios involving high‐speed flows or significant temperature differentials between the fluid and solid phases, the LTE assumption proves insufficient.

Natural convection of water-based nanofluids in three-dimensional enclosures

The single-phase approach, in which nanofluids are treated as pure fluids assuming that the solid and liquid phases are in local thermal and hydrodynamic equilibrium, has been used with good results to simulate forced convection flows. Conversely, its extension to natural convection has revealed to be more problematic, as the obtained numerical results are substantially different from the experimental data, mainly due to the effects of the slip motion occurring between the suspended nanoparticles and the base liquid, whose consequent nonuniform distribution of the solid phase concentration can significantly affect both heat and momentum transfer. In the present article, a two-phase model based on a double-diffusive approach is proposed to study the natural convection flows of water-based metal oxide nanofluids inside cylindrical and rectangular enclosures heated and cooled at their opposite sides with the scope to reproduce a number of experimental data-sets available in the literature. Given the assumption of local thermal equilibrium between nanoparticles and host liquid, and considering the Brownian and thermophoretic diffusion as responsible for causing significant relative velocity between the solid and liquid phases, the governing equations are solved by a control-volume numerical method implemented using the open-source platform OpenFOAM (Open Field Operation and Manipulation). It has been found that the nanoparticle diffusion gives the major contribution to the decrease of the heat transfer rate red usually observed experimentally in natural convection flows of nanofluids inside differentially-heated enclosures. This means that a two-phase model, in conjunction with proper correlations for the evaluation of the nanofluid effective physical properties, must be necessarily enforced to obtain reliable results for many engineering applications.

Numerical modeling and analysis of nanofluids subcooled boiling in a horizontal circular tube

Subcooled flow boiling, characterized by large heat transfer coefficient (HTC), is widely encountered and of great significance in energy industry. To enhance its thermal efficiency for compact applications, nanofluids are employed by improving the thermophysical properties. This article numerically investigates nanofluids subcooled boiling in a horizontal tube under wide-range, multi-factor impacts, including Reynolds number of 40,000–50,000, nanoparticle type, that is, Al2O3, TiO2, Cu, and volume fraction of 1%–5%. The thermal equilibrium assumption and Eulerian multiphase model are utilized after numerical validations. Results indicate that the vapor volume fraction augments along the tube while nanoparticles reduce the vapor production during boiling. As for thermal and flow performance, nanofluids can enhance the HTC at the expense of a rising pressure drop. In the case of Al2O3/Water nanofluid, the maximum of HTC increases by 10.56%. In contrast, the pressure drop increases by 2.1 times when the Al2O3 concentration changes from 1% to 5%. The performance evaluation criterion (PEC) in this article always exceeds 1 for nanofluids, demonstrating that the heat transfer enhancement outweighs the pressure drop penalty, and the Al2O3/Water nanofluid show superior performance with PEC achieving 1.1. This work can provide support for a profound understanding and utilization of nanofluids subcooled boiling.

Bubble-wall velocity in local thermal equilibrium: hydrodynamical simulations vs analytical treatment

We perform real-time hydrodynamical simulations of the growth of bubbles formed during cosmological first-order phase transitions under the assumption of local thermal equilibrium. We confirm that pure hydrodynamic backreaction can lead to steady-state expansion and that bubble-wall velocity in such case agrees very well with the analytical estimates. However, this is not the generic outcome. Instead, it is much more common to observe runaways, as the early-stage dynamics right after the nucleation allow the bubble walls to achieve supersonic velocities before the heated fluid shell in front of the bubble is formed. This effect is not captured by other methods of calculation of the bubble-wall velocity which assume stationary solutions to exist at all times and would have a crucial impact on the possible generation of both baryon asymmetry and gravitational wave signals.

Thermal lattice Boltzmann approach-based numerical study of nanofluid magnetohydrodynamic forced convection in a back-facing stepped porous channel

Unsteady laminar magnetohydrodynamic forced convection and entropy generation inside a backward-facing step (BFS) porous channel with Cu-H2O nanofluid under a tilted magnetic field is numerically handled using a thermal lattice Boltzmann method (TLBM) with two distribution functions. The Darcy-Brinkman-Forchheimer (DBF) equations filled up with the energy equation (with local thermal equilibrium assumption) have been drawn up as the governing equations model. Their numerical solutions unveiled streamlines, Nusselt number, pumping power, performance index and entropy generation, and impacts of influential parameters (Hartmann number (Ha = 0; 10; 20; 100, 200), Darcy number (Da = 10−3; 10−2; 10−1), magnetic tilt angle (γ = 0, π∕4, π∕2), nanoparticles’ volume fraction (φ = 1 − 4%), medium porosity (ε = 0.6, 0.7, 0.8). The irreversibility analysis is also addressed to highlight the flow behavior. Regarding the validation of the modeling approach deemed, a consensus has been achieved with results available in the literature. Computations revealed that the volume fraction provides a fair performance rating with moderate pumping power. Likewise, the average Nusselt number, the average entropy generation, the pumping power and the performance index have a direct relationship with the Hartmann number, the nanoparticles volume fraction and the magnetic field tilt. Moreover, it can be stated the Cu-H2O nanofluid that seems worthwhile in terms of heat transfer efficiency compared to other examined nanofluids. Entropy generation due to the magnetic field was identified as the main source of irreversibility processes in all cases. Based on the outcomes, it is thought that applying a magnetic field can help minimize entropy generation in practical applications.

Melting performance analysis of finned metal foam thermal energy storage tube under steady rotation

The latent heat thermal energy storage (LHTES) technology based on solid-liquid phase change material (PCM) is of great significance for the efficient utilization of thermal energy. To address the issues of slow thermal response and non-uniform melting of the LHTES technology, a hybrid heat transfer enhancement method combined with finned metal foam and steady rotation is proposed in this work. An enthalpy-porosity model considering non-Darcy porous effects and mechanical rotation is established based on the local thermal equilibrium assumption and the fixed grid system. Four different structures (uniform metal foam, graded metal foam, finned metal foam with uniform porosity, and finned metal foam with graded porosity) are investigated firstly to identify the optimal structural arrangement under the same volume of PCM. Numerical results demonstrate that the LHTES units strengthened by the finned metal foam with graded porosity achieve the shortest melting time and largest thermal energy storage rate (TESR). The graded porous structure reduces thermal resistance, rotation enhances flow and heat transfer inside the container, and the fins expand the heat source area. Moreover, the effect of different fin types, graded porosities, and rotational speeds are further considered. The findings suggest that the increase of these three parameters does not correspond to better thermal performance; instead, there is an optimal value. Compared with the base case, the optimal configuration (fin length of 20 mm, fin width of 1.5 mm, porosity gradient of 2%, and rotational speed of 0.5 rpm) can shorten the melting time by 46.68%, increase the TESR and Nu by 74.06% and 69.02%, respectively. This paper validates the feasibility of the hybrid heat transfer enhancement method with finned metal foam and steady rotation, which can offer new insights into the engineering practice of the LHTES technology.

A novel numerical method to evaluate the heat transfer characteristics of complicated CICC structures

In this article, we propose a numerical model for calculating the external heat transfer coefficient between the bundle region and its wall for cables in conduit conductor (CICC). With the assumption of local thermal equilibrium, one macro equation describing the steady heat transfer on the cross section of the CICC can be obtained. A key parameter, the effective transverse thermal conductivity keff, which takes the contribution of both strands and flowing fluid into consideration, is introduced. To calculate the effective transverse thermal conductivity keff, we first obtain the true distribution of different phases (fluid, copper and superconducting material) on a cross section of CICC with the help of the image recognition technique. Based on this, the value of the effective transverse thermal conductivity keff can be calculated numerically. Due to the large difference among the component thermal conductivities (at 4.5 K, typical values of the thermal conductivity of liquid helium, superconducting material and copper are kHe=0.024Wm−1K−1, kNb3Sn=0.04Wm−1K−1 and kCu=708Wm−1K−1; and the maximal ratio of thermal conductivity can be as high as about 30,000), the high-performance finite analytical method (FAM) is recommended to calculate keff. After obtaining the effective transverse thermal conductivity keff of the bundle region, the heat transfer coefficient can be calculated directly by solving a simple Poisson equation under proper boundary conditions. Two case studies are performed, including the JT-60 Super Advanced (JT-60SA) CICC and the dual channel International Thermonuclear Experimental Reactor, Toroidal Field Performance Sample (ITER-TFPS). The calculated values of heat transfer coefficient are consistent with the experimental results, which verifies our proposed model.

Computational fluid dynamic simulation of packed bed drying process: impact of particle properties, drying conditions, and lateral edge heating modes on drying kinetics

In this work, a two-dimensional computational fluid dynamic (CFD) model is developed to describe the drying process of a packed bed made of spherical particles. The volumetrical evaporation rate inside the bed is computed from the pressure difference between the particle surface and the airflow. By using the thermal equilibrium assumption, the heat conservation equation is derived. The CFD model is solved in the COMSOL Multiphysics environment. The obtained results indicate remarkable maldistributions of temperature and moisture content. These maldistributions can be explained by the impact of lateral edges on thermo-hydraulic behavior. Additionally, the impact of particle diameter, air velocity, and bed width on the spatial-temporal moisture content and temperature distribution is investigated. It shows that the CFD model can be simplified to the receding front drying model for a bed made of small particles. Furthermore, by changing the thermal boundary conditions at the lateral edges, the influence of the heating mode at the lateral edges on the drying behavior is explored. The results indicate that contact heating at the bed wall can help to accelerate the drying process significantly.

Local thermal equilibrium constraints for energy transport equations for thermal oil recovery processes

The mathematical framework of reservoir thermal simulators typically relies on the assumption of local thermal equilibrium (LTE) when modeling energy transport through porous media for petroleum engineering applications, including Thermal Enhanced Oil Recovery (TEOR) processes. The LTE presupposes instantaneous exchange of energy between the injected hot fluids and the reservoir's rock and fluids. Nevertheless, it is essential to acknowledge that LTE's validity becomes questionable in specific scenarios, potentially resulting in an overestimation of oil recovery, e.g., when high fluid velocities arise. Under such conditions, the local thermal non-equilibrium (LTNE) formulation may offer a more accurate representation of the energy transport physics. In this study, we present the physical constraints that must be satisfied for the application of two distinct mathematical formulations (1) the one-equation model for the fluids and rock as a whole (LTE is assumed), and (2) the two-equation model, one equation assigned to the fluids and another one to the rock (pseudo-LTNE is assumed). Such constraints were derived by examining conductive and convective heat transport properties via the theory of the method of volume averaging, specifically applied to an oil-water-gas-rock system at the pore scale. It was found that the one-equation model requires that 18 physical constraints to be satisfied, while the pseudo-LTNE model needs to fulfill 9 physical constraints. A sensitivity analysis was conducted based on factor upon which accumulation, conduction and convection constraint equations depend, aiming to identify the conditions under which LTNE occurs within the porous medium. These influencing factors encompass the characteristic time, the Peclet number, the presence of external heat sources, the distribution of phases at the pore-scale, and the disparity between rock and fluid properties. The combinations of these factors can induce the existence of local thermal non-equilibrium. One significant property is the heat transfer coefficient in the reservoir. We recommend that the pseudo-LTNE approach is reliable as long as the rock-fluid heat transfer coefficients are at least one order of magnitude larger than the fluid-fluid heat transfer coefficients (1000−10000W/(m3K)). Contrarily, the OEM should be employed only when the heat transfer coefficients surpass a threshold value of 10000W/(m3K). When injector and producer wells are hydraulically interconnected, resulting in heightened fluid velocities, it is encouraged the use of the two-equations model. The predictability and reliability of the energy models were verified by comparison of predicted oil recovery against data from laboratory experiments. The equations presented in this study consistently demonstrate their utility in assessing the reliability of energy transport models for TEOR scenarios.

Prediction and control of elevated temperatures within landfills under aeration and recirculation based on the thermal non-equilibrium model

Aeration is an effective approach to sustainable landfilling but may lead to elevated temperatures within landfills, resulting in landfill fires or explosions. Therefore, aeration is usually combined with leachate recirculation to control the elevated temperatures within landfills. To predict landfill temperatures during aeration and recirculation, a local thermal non-equilibrium model is developed considering the heat generation of biodegradation, the heat removal due to evaporation and leachate-gas flow, and the effects of the capillary. The solver is implemented in OpenFOAM based on the finite volume method and validated against a waste-column experiment and an in-situ aeration test. The simulation results demonstrate that the assumption of local thermal equilibrium will distinctly overestimate the temperature, maximally by 15 °C in the studied case. The model is then used to simulate a typical aerobic landfill unit to investigate the formation of explosive gas mixtures and elevated temperatures under different operating conditions. The simulation results of gas composition suggest that aeration will not result in explosive gas within landfills. A reasonable recirculation method for temperature control with corresponding operating parameters under a group of values of aeration pressure (2000–4000 Pa) and recirculation rate (0.0001–0.0008 m/s) are proposed, which can provide some guides for the design of an aeration and recirculation combined system. For a given total volume of added leachate, a higher recirculation rate does not always mean better cooling, and the cooling effect of continuous recirculation is better than that of intermittent recirculation.

Experimental investigation on the validity of the local thermal equilibrium assumption in ablative-material response models

Thermal Protection Systems (TPS) material response models rely on the assumption of local thermal equilibrium (LTE) between the solid phase and the gas phase. This assumption was challenged and investigated by several authors but a sufficiently precise knowledge of heat transfer coefficients in TPS materials was lacking to reach final conclusions. The objective of this work is to contribute to filling this gap by providing a literature review of available data in other communities (thermal energy storage, heat exchangers) and by performing an experimental characterization of Calcarb, a commercial carbon preform used for manufacturing thermal protection systems. Heat transfer within Calcarb was studied experimentally in the Through-Thickness and in the In-Plane directions for Reynolds numbers of 1 to 4 - representative of the TPS application - using the transient single-blow technique. Numerical parameter estimation was performed using the Porous material Analysis Toolbox based on OpenFoam (PATO) and the Design Analysis Kit for Optimization and Terascale Applications (DAKOTA). The heat transfer coefficient hv is found to be greater than or equal to 108 W/(m3K) and the LTE assumption is shown to be valid in the conditions of the experiment. To assess the validity of the LTE assumption for other conditions, the above bound of hv may now be used in combination with a local thermal non-equilibrium model.

A Spectroscopic Thermometer: Individual Vibrational Band Spectroscopy with the Example of OH in the Atmosphere of WASP-33b

Individual vibrational band spectroscopy presents an opportunity to examine exoplanet atmospheres in detail, by distinguishing where the vibrational state populations of molecules differ from the current assumption of a Boltzmann distribution. Here, retrieving vibrational bands of OH in exoplanet atmospheres is explored using the hot Jupiter WASP-33b as an example. We simulate low-resolution spectroscopic data for observations with the JWST's NIRSpec instrument and use high-resolution observational data obtained from the Subaru InfraRed Doppler instrument (IRD). Vibrational band–specific OH cross-section sets are constructed and used in retrievals on the (simulated) low- and (real) high-resolution data. Low-resolution observations are simulated for two WASP-33b emission scenarios: under the assumption of local thermal equilibrium (LTE) and with a toy non-LTE model for vibrational excitation of selected bands. We show that mixing ratios for individual bands can be retrieved with sufficient precision to allow the vibrational population distributions of the forward models to be reconstructed. A fit for the Boltzmann distribution in the LTE case shows that the vibrational temperature is recoverable in this manner. For high-resolution, cross-correlation applications, we apply the individual vibrational band analysis to an IRD spectrum of WASP-33b, applying an “unpeeling” technique. Individual detection significances for the two strongest bands are shown to be in line with Boltzmann-distributed vibrational state populations, consistent with the effective temperature of the WASP-33b atmosphere reported previously. We show the viability of this approach for analyzing the individual vibrational state populations behind observed and simulated spectra, including reconstructing state population distributions.

Entropy Generation Analysis on a Metal Foam in an Automotive Exhaust Line with Thermoelectric Generator

Abstract In the present paper, an entropy generation analysis on a 2-D steady state problem in convective regime of an aluminum foam partially and totally filled channel with an external TEG element is solved in numerical way. The numerical analyses are accomplished with the assumption of the local thermal equilibrium (LTE) hypothesis in order to model the metal foam and heat transfer inside the channel. The working fluid is exhaust gas characterized by the same properties of the air in correspondence to the TEG upper surface temperature. The TEG is considered as a solid component characterized by an internal energy generation. The thermophysical properties are assumed temperature independent. Ansys-Fluent code is employed in order to resolve the governing equations for exhaust gas, metal foam and TEG. Different exhaust gas mass flow rates on the inlet section are assumed. Several thicknesses of aluminum foam values are employed. The porous media are characterized by a porosity from 0.90 to 0.978 and number of pores per inch (PPI) equal to 5, 10, 20, 40. Results are given in terms of global entropy generations related to the thermal and viscous effects. The total global entropy generation increases with increasing of exhaust gas flow rate for all pore density and porosity values. Bejan number decreases with the increment of mass flow rate and thickness. It increases when the porosity value increases whereas at high mass flow rate and for assigned porosity the values present small difference for different pore density values.

Local thermal equilibrium analysis of complete phase change process inside porous diffuser using NanoFluids

This works implements an investigation on two-phase flow inside porous diffuser using Cu/water nanofluid instead of pure water. Interesting has been focused on how the using Cu/water affect the predictions of complete boiling process. Two-phase mixture model (TPMM), based on the assumptions of Local Thermal Equilibrium (LTE), has been converted to implement the actions of combining nano-fluid rather only water. The governing equations are solved by finite volume method. For the purpose of comparison, the solutions of temperature and liquid saturation, obtained employing Cu/water and only water as working fluid, have been displayed for a wide ranges of geometric parameters, operating parameters and properties of porous media to investigate impacts of using nanofluid. The findings illustrated that the adding Cu particle to water has a high influences on the heat diffusivity inside the diffuser as compared to the solutions of only water due to enhancing the axial diffusion. By employing Cu/water, the boiling and condensation fronts are considerably moved to downstream and upstream sections, respectively. The exit temperature is significantly reduced when utilizing Cu/water as compared to the only water solution. The findings of Cu/water show that the volume fraction of nano-fluid, operating conditions, geometric conditions and porous structure properties have large impacts on the locations of liquid and steam fronts. Accordingly, the nano-fluid is supposed as a safety key in the boiling applications since it improves the technique of heat transfer while simulating the complete boiling process.

Validity of Local Thermal Equilibrium Assumption for Heat Transfer in a Saturated Soil Layer

The assumption of local thermal equilibrium (LTE) between solid and fluid phases commonly is used for studies of heat transfer in saturated soil and is valid for a wide range of conditions. However, for certain conditions, the heat transfer process may give rise to local thermal nonequilibrium (LTNE) in which adjacent solid and fluid phases have different temperatures. This note presents the results of a numerical study of the validity of the LTE assumption for one-dimensional heat transfer in a saturated soil layer with fluid flow. For the conditions investigated, the LTE assumption holds for soil with particle sizes smaller than the gravel range. For soil with particles sizes in the gravel range and larger, the LTE assumption may be valid or invalid, depending primarily on fluid discharge velocity, with higher fluid velocity values more likely to produce LTNE conditions.

CO2 storage in porous media unsteady thermosolutal natural convection -Application in deep saline aquifer reservoirs

In this paper, the storage of CO2 is investigated numerically under supercritical conditions in deep saline aquifer reservoirs. The simulations were performed on a non-deformable, saturated porous material inside a vertical enclosure that was assumed to be impermeable and insulated on three sides. The porous medium is homogeneous and isotropic with constant thermo-physical properties except for the fluid density, which varies according to Boussinesq approximations. A dynamic model assumes that flow is two-dimensional and obeys Darcy's motion law. The thermal equilibrium assumption (LTE) is also considered. The set of conservation equations and the appropriate initial and boundary conditions have been resolved numerically by the classical finite volume method. Spatio-temporal variations of different state variables such as pressure, velocity, temperature, and concentration were numerically simulated and plotted versus controlling parameters, particularly the thermal and solutal Rayleigh numbers, RaT and RaS; the Lewis number, Le; and the reservoir's aspect ratio, A. According to regulating factors, and physical and geological qualities, we could predict the behavior of CO2 that would be stored in a geologic reservoir that is thought to be porous media. We could also estimate the time and direction of CO2 mass transfer based on the abovementioned factors. Great attention was paid to examining the sensitivity of fluid flow, heat, and mass transfer rates according to the reservoir form and the operating conditions.