• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.112-117
• /
• 2000

A numerical study on natural convection in a vertical square cavity filled with a porous medium is carried out with Brinkman-Forchheimer-extended Darcy flow model, and the validity of local thermodynamic equilibrium assumption is studied. The local thermodynamic equilibrium refers to the state in which a single temperature can be used to describe a heat transfer process in a multiphase system. With this assumption, the analysis is greatly simplified because only one equation is needed to describe the heat transfer process. But prior to using this assumption, it is necessary to know in what conditions the assumption can be used. The numerical results of this study reveal that large temperature difference between fluid phase and solid phase exists near wall region, paticularily when the convection becomes dominant over conduction. And the influence of flow parameters such as fluid Rayleigh number, fluid Prandtl number, dimensionless particle diameter and conductivity ratio are investigated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.12
• /
• pp.1635-1643
• /
• 2000

The present work investigates a heat transfer phenomenon at the interface between a porous medium and an impermeable wall. In an effort to appropriately describe the heat transfer phenomenon at the interface, the heat transfer at the interface between the microchannel heat sink, which is an ideally organized porous medium, and the finite-thickness substrate is examined. From the examination, it is clarified that the he heat flux distribution at the interface is not uniform for the impermeable wall with finite thickness. On the other hand, the first approach, based on the energy balance for the representative elementary volume in the porous medium, is physically reason able. When the first approach is applied to the thermal boundary condition, and additional boundary condition based on the local thermal equilibrium assumption at the interface is used. This additional boundary condition is applicable except for the very th in impermeable wall. Hence, for practical situations, the first approach in combination with the local thermal equilibrium assumption at the interface is suggested as an appropriate thermal boundary condition. In order to confirm our suggestion, convective flows both in a microchannel heat sink and in a sintered porous channel subject to a constant heat flux condition are analyzed. The analytically obtained thermal resistance of the microchannel heat sink and the numerically obtained overall Nusselt number for the sintered porous channel are shown to be in close agreement with available experimental results when our suggestion for the thermal boundary conditions is applied.

• The KSFM Journal of Fluid Machinery
• /
• v.20 no.2
• /
• pp.11-16
• /
• 2017

The accurate prediction of leakage flow through the brush element of brush seal at the steam turbine is important to find optimum design parameters for increasing an efficiency. In this study, CFD analysis method using commercial software FLUENT is proposed to predict leakage through the brush element. Since the brush element has a complex three-dimensional shape with many bristle assemblies, it is difficult to analyze the flow field. Therefore, if the brush element is assumed to be porous medium region, the analysis time can be shortened. Two determination methods of resistance coefficients of the Darcian porous medium equation are suggested. By comparing the 2D and 3D CFD analysis results for the leakage of the brush element using the two resistance coefficient determination methods, the effectiveness of the analysis for the porous medium assumption is proved.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.8
• /
• pp.1172-1182
• /
• 2002

In the present study a general criterion for local thermal equilibrium is presented in terms of parameters of engineering importance which include the Darcy number, the effective Prandtl number of fluid, and the Reynolds number. For this, an order of magnitude analysis is performed for the case when the effect of convection heat transfer is dominant in a porous structure. The criterion proposed in this study is more general than the previous criterion suggested by Carbonell and Whitaker, because the latter is applicable only when conduction is the dominant heat transfer mode in a porous medium while the former can be applied even when convection heat transfer prevails. In order to check the validity of the proposed criterion for local thermal equilibrium, the forced convection phenomena in a porous medium with a microchanneled structure subject to an impinging jet are studied using a similarity transformation. The proposed criterion is also validated with the existing experimental and numerical results for convection heat transfer in various porous materials that include some of the parameters used in the criterion such as a microchannel heat sink with a parallel flow, a packed bed, a cellular ceramic, and a sintered metal. It is shown that the criterion presented in this work well-predicts the validity of the assumption of local thermal equilibrium in a porous medium.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.7
• /
• pp.983-990
• /
• 2002

A numerical analysis has been carried out on forced convection heat transfer in the developing region of a porous channel. The channel is filled with an isotropic porous medium. At the channel walls, a uniform heat flux is given. Comprehensive numerical solutions are acquired to the Brinkman-Forchheimer extended Darcy equation and the LTNE model which does not employ the assumption of local thermal equilibrium between solid and fluid phases. Details of thermal fields in the developing region are examined over wide ranges of the thermal parameters. The numerical solutions at the fully developed region are compared with the previous analytical solutions. The correlation for predicting local Nusselt number in a porous channel is proposed.

• Steel and Composite Structures
• /
• v.45 no.3
• /
• pp.409-423
• /
• 2022

Nowadays, there is a high demand for great structural implementation and multifunctionality with excellent mechanical properties. The porous structures reinforced by graphene platelets (GPLs) having valuable properties, such as heat resistance, lightweight, and excellent energy absorption, have been considerably used in different engineering implementations. However, stiffness of porous structures reduces significantly, due to the internal cavities, by adding GPLs into porous medium, effective mechanical properties of the porous structure considerably enhance. This paper is relating to vibration analysis of fluidconveying cantilever porous graphene platelet reinforced (GPLR) pipe with fractional viscoelastic model resting on foundations. A dynamical model of cantilever porous GPLR pipes conveying fluid and resting on a foundation is proposed, and the vibration, natural frequencies and primary resonant of such a system are explored. The pipe body is considered to be composed of GPLR viscoelastic polymeric pipe with porosity in which Halpin-Tsai scheme in conjunction with the fractional viscoelastic model is used to govern the construction relation of nanocomposite pipe. Three different porosity distributions through the pipe thickness are introduced. The harmonic concentrated force is also applied to the pipe and the excitation frequency is close to the first natural frequency. The governing equation for transverse motions of the pipe is derived by the Hamilton principle and then discretized by the Galerkin procedure. In order to obtain the frequency-response equation, the differential equation is solved with the assumption of small displacement, damping coefficient, and excitation amplitude by the multiple scale method. A parametric sensitivity analysis is carried out to reveal the influence of different parameters, such as nanocomposite pipe properties, fluid velocity and nonlinear viscoelastic foundation coefficients, on the primary resonance and linear natural frequency. Results indicate that the GPLs weight fraction porosity coefficient, fractional derivative order and the retardation time have substantial influences on the dynamic response of the system.

• Journal of the Korean Society of Groundwater Environment
• /
• v.2 no.1
• /
• pp.22-29
• /
• 1995

Mobile bacterial particles can act as carriers and enhance the transport of hydrophobic contaminants in ground water by reducing retardation effects. Because of their colloidal size and favorable surface conditions, bacteria can act as efficient contaminant carriers. When such carriers exist in a porous medium, the system can be thought of as three phases: an aqueous phase, a carrier phase, and a stationary solid matrix phase. Contaminant can be present in either or all of these phases. In this study, a mathematical model based on mass balances is developed to describe the transport and fate of biodegradable contaminant in a porous medium. Bacterial mass transfer mechanism between aqueous and solid matrix phases, and contaminant mass transfer between aqueous and bacterial phases are represented by kinetic models. Governing equations are non-dimensionalized and solved to analyze the bacteria facilitated contaminant transport. The numerical results of the facilitation effect match favorably with experimental data reported in the literature. Results show that the contaminant transport can be described by local equilibrium assumption when Damkohler numbers are larger than 10. Significant sensitivities to model parameters, particularly bacterial growth rate and influent bacterial concentration, were discovered.

• Journal of Mechanical Science and Technology
• /
• v.15 no.6
• /
• pp.796-803
• /
• 2001

Numerical solutions for the convection-dominated melting in a rectangular cavity are presented. The enthalpy-porosity model is employed as the mathematical model. This model is applied in conjunction with the EIT method to detect boundary movement in a phase changing environment. The absorption and evolution of latent heat during the phase change is dealt with by the enthalpy-based energy equation. This seems to be more efficient than resolving the temperature-based energy equation. The velocity switch-off, which is required when solid changes into liquid, is modeled by the porous medium assumption. For efficiency and simplicity of the solutions procedure, this paper proposes a simple algorithm, which iterates the temperature and the liquid fraction of the cells comprising the front layer. The numerical results agree reasonably well with the experimental data and other previous works using the transformed-grid system.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.3
• /
• pp.271-277
• /
• 2014

Internal structures, especially those located in the upstream of a reactor core, may have a significant influence on the core inlet flow rate distribution depending on both their shapes and the relative distance between the internal structures and the core inlet. In this study, to examine the effect of the reactor internal structure geometry treatment method on the prediction accuracy for the scale-down APR+ flow distribution, simulations with real geometry modeling were conducted using ANSYS CFX R.14, a commercial computational fluid dynamics software, and the predicted results were compared with those of the porous medium assumption. It was concluded that the core inlet flow distribution could be predicted more accurately by considering the real geometry of the internal structures located in the upstream of the core inlet. Therefore, if sufficient computational resources are available, an exact representation of these internal structures, for example, lower support structure bottom plate and ICI nozzle support plate, is needed for the accurate simulation of the reactor internal flow.

• Korean Chemical Engineering Research
• /
• v.59 no.1
• /
• pp.138-151
• /
• 2021

The effect of a reactant ratio on the growth of a buoyancy-driven instability in an irreversible A+B→C reaction system is analyzed theoretically and numerically. Taking a non-stoichiometric reactant ratio into account, new linear stability equations are derived without the quasi-steady state assumption (QSSA) and solved analytically. It is found that the main parameters to explain the present system are the Damköhler number, the dimensionless density difference of chemical species and the ratio of reactants. The present initial grow rate analysis without QSSA shows that the system is initially unconditionally stable regardless of the parameter values; however, the previous initial growth rate analysis based on the QSSA predicted the system is unstable if the system is physically unstable. For time evolving cases, the present growth rates obtained from the spectral analysis and pseudo-spectral method support each other, but quite differently from that obtained under the conventional QSSA. Adopting the result of the linear stability analysis as an initial condition, fully nonlinear direct numerical simulations are conducted. Both the linear analysis and the nonlinear simulation show that the reactant ratio plays an important role in the onset and the growth of the instability motion.