Search

In this article, a finite element model is implemented to analyze hydro-thermal convective flow in a porous medium. The mathematical model encompasses Darcy's law for incompressible fluid behavior, which is coupled with a convection-diffusion-type energy equation to characterize the temperature in the porous medium. The current investigation presents an efficient, stable, and accurate finite element discretization for the hydro-thermal convective flow model. The well-posedness of the proposed discrete Galerkin finite element formulation is guaranteed due to the decoupling property and the linearity of the numerical method. Computational experiments confirm the optimal convergence rates for a manufactured solution. Several numerical results are obtained for the variations of the hydraulic resistivity and thermal diffusivity. In the present study, the bottom wall is maintained at a constant higher hot temperature while side vertical walls are thermally insulated and the top wall is maintained at a constant cold temperature. Heat transfer rates at the heated bottom wall are presented in terms of local Nusselt number. A linear variation in hydraulic resistivity and a quadratic variation in thermal diffusivity show an increase in the heat transfer rate.

In this paper, we study the preferential stiffness and the crack-tip fields for an elastic porous solid of which material properties are dependent upon the density. Such a description is necessary to describe the failure that can be caused by damaged pores in many porous bodies such as ceramics, concrete and human bones. To that end, we revisit a new class of implicit constitutive relations under the assumption of small deformation. Although the constitutive relationship \textit{appears linear} in both the Cauchy stress and linearized strain, the governing equation bestowed from the balance of linear momentum results in a quasi-linear partial differential equation (PDE) system. For the linearization and obtaining a sequence of elliptic PDEs, we propose the solution algorithm comprise a \textit{Newton's method} coupled with a bilinear continuous Galerkin-type finite elements for the discretization. Our algorithm exhibits an optimal rate of convergence for a manufactured solution. In the numerical experiments, we set the boundary value problems (BVPs) with edge crack {under different modes of loading (i.e., the pure mode-I, II, and the mixed-mode). From the numerical results, we find that the density-dependent moduli model describes diverse phenomena that are not captured within the framework of classical linearized elasticity. In particular,numerical solutions clearly indicate that the nonlinear \textit{modeling} parameter depending on its sign and magnitude can control preferential mechanical stiffness along with the change of volumetric strain; larger the parameter is in the positive value}, the responses are such that the strength of porous solid gets weaker against the tensile loading while stiffer against the in-plane shear (or compressive) loading, which is vice versa for the negative value of it.