Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located...Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located in the middle of the cavity with inlet and outlet ports with different configurations of inlet/outlet ports.Convective drying performance is numerically assessed for different values of Reynolds numbers(between 50 and 250),dry air inlet temperature(between 40 and 80℃)and different locations of the port.It is observed that changing the port locations has significant impacts on the flow recirculaitons inside the triangular chamber while convective drying performance is highly affected.The moisture content reduces with longer time and for higher Reynolds number(Re)values.Case P4 where inlet and outlet ports are in the middle of the walls provides the most effective configuration in terms of convective drying performance while the worst case is seen for P1 case where inlet and outlet are closer to the corners of the chamber.There is a 192% difference between the moisture reduction of these two cases at Re=250,T=80℃ and t=120 min.展开更多
Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid...Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.展开更多
Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc....Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc.Motivated by such applications,in this article,a numerical study of entropy generation impacts on the heat and momentum transfer in time-dependent laminar incompressible boundary layer flow of a Casson viscoplastic fluid over a uniformly heated vertical cylinder embedded in a porous medium is presented.Darcy’s law is used to simulate bulk drag effects at low Reynolds number for an isotropic,homogenous porous medium.Heat line visualization is also included.The mathematical model is derived and normalized using appropriate transformation variables.The resulting non-linear time-dependent coupled governing equations with associated boundary conditions are solved via an implicit finite difference method which is efficient and unconditionally stable.The outcomes show that entropy generation and Bejan number are both elevated with increasing values of Darcy number,Casson fluid parameter,group parameter and Grashof number.To analyze the heat transfer process in a two-dimensional domain,plotting heat lines provides an excellent approach in addition to streamlines and isotherms.It is remarked that as the Darcy number increases,the deviations of heat lines from the hot wall are reduced.展开更多
The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper ...The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.展开更多
A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering...A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering ocean waves of deploying process.The discretization of the RANS equations is achieved by a finite volume approach(FV).The volume of fluid method(VOF) is employed to track the complicated free surface.A numerical wave tank is built to generate the ocean waves which are suitable for deploying seafloor mining tool.A typical deploying condition is employed to reflect the process of flatted-bottom body impacting with waves,and the pressure distribution of bottom is also presented.Four different lowering velocities are applied to obtain the time histories of maximum pressure of bottom,and it can be concluded that the pressure coefficient decreases with water velocity increasing,which is similar with ordinary water entry case.The numerical results clearly demonstrate the characteristics of flatted-bottom body entering ocean waves.展开更多
Coalbed gas non-Darcy flow has been observed in high permeable fracture systems,and some mathematical and numerical models have been proposed to study the effects of non-Darcy flow using Forchheimer non-Darcy model.Ho...Coalbed gas non-Darcy flow has been observed in high permeable fracture systems,and some mathematical and numerical models have been proposed to study the effects of non-Darcy flow using Forchheimer non-Darcy model.However,experimental results show that the assumption of a constant Forchheimer factor may cause some limitations in using Forchheimer model to describe non-Darcy flow in porous media.In order to investigate the effects of non-Darcy flow on coalbed methane production,this work presents a more general coalbed gas non-Darcy flow model according to Barree-Conway equation,which could describe the entire range of relationships between flow velocity and pressure gradient from low to high flow velocity.An expanded mixed finite element method is introduced to solve the coalbed gas non-Darcy flow model,in which the gas pressure and velocity can be approximated simultaneously.Error estimate results indicate that pressure and velocity could achieve first-order convergence rate.Non-Darcy simulation results indicate that the non-Darcy effect is significant in the zone near the wellbore,and with the distance from the wellbore increasing,the non-Darcy effect becomes weak gradually.From simulation results,we have also found that the non-Darcy effect is more significant at a lower bottom-hole pressure,and the gas production from non-Darcy flow is lower than the production from Darcy flow under the same permeable condition.展开更多
文摘Heat and mass transfer of a circular-shaped porous moist object inside a two-dimensional triangle cavity is investigated by using finite element method.The porous object is considered to be a moist food sample,located in the middle of the cavity with inlet and outlet ports with different configurations of inlet/outlet ports.Convective drying performance is numerically assessed for different values of Reynolds numbers(between 50 and 250),dry air inlet temperature(between 40 and 80℃)and different locations of the port.It is observed that changing the port locations has significant impacts on the flow recirculaitons inside the triangular chamber while convective drying performance is highly affected.The moisture content reduces with longer time and for higher Reynolds number(Re)values.Case P4 where inlet and outlet ports are in the middle of the walls provides the most effective configuration in terms of convective drying performance while the worst case is seen for P1 case where inlet and outlet are closer to the corners of the chamber.There is a 192% difference between the moisture reduction of these two cases at Re=250,T=80℃ and t=120 min.
基金Projects(2006AA06Z105, 2007AA06Z134) supported by the National High-Tech Research and Development Program of ChinaProjects(2007, 2008) supported by China Scholarship Council (CSC)
文摘Based on the fact that 3-D model discretization by artificial could not always be successfully implemented especially for large-scaled problems when high accuracy and efficiency were required, a new adaptive multigrid finite element method was proposed. In this algorithm, a-posteriori error estimator was employed to generate adaptively refined mesh on a given initial mesh. On these iterative meshes, V-cycle based multigrid method was adopted to fast solve each linear equation with each initial iterative term interpolated from last mesh. With this error estimator, the unknowns were nearly optimally distributed on the final mesh which guaranteed the accuracy. The numerical results show that the multigrid solver is faster and more stable compared with ICCG solver. Meanwhile, the numerical results obtained from the final model discretization approximate the analytical solutions with maximal relative errors less than 1%, which remarkably validates this algorithm.
基金DST-INSPIRE (Code No. IF160028) for the grant of research fellowship
文摘Thermal transport in porous media has stimulated substantial interest in engineering sciences due to increasing applications in filtration systems,porous bearings,porous layer insulation,biomechanics,geomechanics etc.Motivated by such applications,in this article,a numerical study of entropy generation impacts on the heat and momentum transfer in time-dependent laminar incompressible boundary layer flow of a Casson viscoplastic fluid over a uniformly heated vertical cylinder embedded in a porous medium is presented.Darcy’s law is used to simulate bulk drag effects at low Reynolds number for an isotropic,homogenous porous medium.Heat line visualization is also included.The mathematical model is derived and normalized using appropriate transformation variables.The resulting non-linear time-dependent coupled governing equations with associated boundary conditions are solved via an implicit finite difference method which is efficient and unconditionally stable.The outcomes show that entropy generation and Bejan number are both elevated with increasing values of Darcy number,Casson fluid parameter,group parameter and Grashof number.To analyze the heat transfer process in a two-dimensional domain,plotting heat lines provides an excellent approach in addition to streamlines and isotherms.It is remarked that as the Darcy number increases,the deviations of heat lines from the hot wall are reduced.
基金Project(10872219) supported by the National Natural Science Foundation of China
文摘The finite element method was used to solve fluid dynamic interaction problems between the crust and mantle of the Earth. To consider different mechanical behaviours, the lithosphere consisting of the crust and upper mantle was simulated as fluid-saturated porous rocks, while the upper aesthenospheric part of the mantle was simulated as viscous fluids. Since the whole lithosphere was computationally simulated, the dynamic interaction between the crust and the upper mantle was appropriately considered. In particular, the mixing of mantle fluids and crustal fluids was simulated in the corresponding computational model. The related computational simulation results from an example problem demonstrate that the mantle fluids can flow into the crust and mix with the crustal fluids due to the resulting convective flows in the crust-mantle system. Likewise, the crustal fluids can also flow into the upper mantle and mix with the mantle fluids. This kind of fluids mixing and exchange is very important to the better understanding of the governing processes that control the ore body formation and mineralization in the upper crust of the Earth.
基金Project(51305463)supported by National Natural Science Foundation of ChinaProject(2012QNZT01601005125)supported by Free Exploration Plan of Central South University,ChinaProject supported by Postdoctoral Foundation of Central South university,China
文摘A numerical wave load model based on two-phase(water-air) Reynolds-Averaged Navier Stokes(RANS) type equations is used to evaluate hydrodynamic forces exerted on flatted-bottom seafloor mining tool during its entering ocean waves of deploying process.The discretization of the RANS equations is achieved by a finite volume approach(FV).The volume of fluid method(VOF) is employed to track the complicated free surface.A numerical wave tank is built to generate the ocean waves which are suitable for deploying seafloor mining tool.A typical deploying condition is employed to reflect the process of flatted-bottom body impacting with waves,and the pressure distribution of bottom is also presented.Four different lowering velocities are applied to obtain the time histories of maximum pressure of bottom,and it can be concluded that the pressure coefficient decreases with water velocity increasing,which is similar with ordinary water entry case.The numerical results clearly demonstrate the characteristics of flatted-bottom body entering ocean waves.
基金Projects(91330106,11171190)supported by the National Natural Science Foundation of ChinaProjects(15CX05065A,15CX05003A)supported by the Fundamental Research Funds for the Central Universities,China
文摘Coalbed gas non-Darcy flow has been observed in high permeable fracture systems,and some mathematical and numerical models have been proposed to study the effects of non-Darcy flow using Forchheimer non-Darcy model.However,experimental results show that the assumption of a constant Forchheimer factor may cause some limitations in using Forchheimer model to describe non-Darcy flow in porous media.In order to investigate the effects of non-Darcy flow on coalbed methane production,this work presents a more general coalbed gas non-Darcy flow model according to Barree-Conway equation,which could describe the entire range of relationships between flow velocity and pressure gradient from low to high flow velocity.An expanded mixed finite element method is introduced to solve the coalbed gas non-Darcy flow model,in which the gas pressure and velocity can be approximated simultaneously.Error estimate results indicate that pressure and velocity could achieve first-order convergence rate.Non-Darcy simulation results indicate that the non-Darcy effect is significant in the zone near the wellbore,and with the distance from the wellbore increasing,the non-Darcy effect becomes weak gradually.From simulation results,we have also found that the non-Darcy effect is more significant at a lower bottom-hole pressure,and the gas production from non-Darcy flow is lower than the production from Darcy flow under the same permeable condition.
