Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA model for reactive porous transport during re-wetting of hardened concrete
271
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A mathematical model is developed that captures the transport of liquid water in hardened concrete, as well as the chemical reactions that occur between the imbibed water and the residual calcium silicate compounds residing in the porous concrete matrix. The main hypothesis in this model is that the reaction product -- calcium silicate hydrate gel -- clogs the pores within the concrete thereby hindering water transport. Numerical simulations are employed to determine the sensitivity of the model solution to changes in various physical parameters, and compare to experimental results available in the literature.
rate research
Read More
Multi-phase reactive transport processes are ubiquitous in igneous systems. A challenging aspect of modelling igneous phenomena is that they range from solid-dominated porous to liquid-dominated suspension flows and therefore entail a wide spectrum of rheological conditions, flow speeds, and length scales. Most previous models have been restricted to the two-phase limits of porous melt transport in deforming, partially molten rock and crystal settling in convecting magma bodies. The goal of this paper is to develop a framework that can capture igneous system from source to surface at all phase proportions including not only rock and melt but also an exsolved volatile phase. Here, we derive an n-phase reactive transport model building on the concepts of Mixture Theory, along with principles of Rational Thermodynamics and procedures of Non-equilibrium Thermodynamics. Our model operates at the macroscopic system scale and requires constitutive relations for fluxes within and transfers between phases, which are the processes that together give rise to reactive transport phenomena. We introduce a phase- and process-wise symmetrical formulation for fluxes and transfers of entropy, mass, momentum, and volume, and propose phenomenological coefficient closures that determine how fluxes and transfers respond to mechanical and thermodynamic forces. Finally, we demonstrate that the known limits of two-phase porous and suspension flow emerge as special cases of our general model and discuss some ramifications for modelling pertinent two- and three-phase flow problems in igneous systems.
We study the asymptotic behaviour of sharp front solutions arising from the nonlinear diffusion equation theta_t = (D(theta)theta_x)_x, where the diffusivity is an exponential function D({theta}) = D_o exp(betatheta). This problem arises for example in the study of unsaturated flow in porous media where {theta} represents the liquid saturation. For physical parameters corresponding to actual porous media, the diffusivity at the residual saturation is D(0) = D_o << 1 so that the diffusion problem is nearly degenerate. Such problems are characterised by wetting fronts that sharply delineate regions of saturated and unsaturated flow, and that propagate with a well-defined speed. Using matched asymptotic expansions in the limit of large {beta}, we derive an analytical description of the solution that is uniformly valid throughout the wetting front. This is in contrast with most other related analyses that instead truncate the solution at some specific wetting front location, which is then calculated as part of the solution, and beyond that location the solution is undefined. Our asymptotic analysis demonstrates that the solution has a four-layer structure, and by matching through the adjacent layers we obtain an estimate of the wetting front location in terms of the material parameters describing the porous medium. Using numerical simulations of the original nonlinear diffusion equation, we demonstrate that the first few terms in our series solution provide approximations of physical quantities such as wetting front location and speed of propagation that are more accurate (over a wide range of admissible {beta} values) than other asymptotic approximations reported in the literature.
We analyze the dynamics of liquid filling in a thin, slightly inflated rectangular channel driven by capillary forces. We show that although the amount of liquid $m$ in the channel increases in time following the classical Lucas-Washburn law, $m propto t^{1/2}$, the prefactor is very sensitive to the deformation of the channel because the filling takes place by the growth of two parts, the bulk part (where the cross-section is completely filled by the liquid), and the finger part (where the cross-section is partially filled). We calculate the time dependence of $m$ accounting for the coupling between the two parts and show that the prefactor for the filling can be reduced significantly by a slight deformation of the rectangular channel, e.g., the prefactor is reduced 50% for a strain of 0.1%. This offers an explanation for the large deviation in the value of the prefactor reported previously.
A reactive fluid dissolving the surrounding rock matrix can trigger an instability in the dissolution front, leading to spontaneous formation of pronounced channels or wormholes. Theoretical investigations of this instability have typically focused on a steadily propagating dissolution front that separates regions of high and low porosity. In this paper we show that this is not the only possible dissolutional instability in porous rocks; there is another instability that operates instantaneously on any initial porosity field, including an entirely uniform one. The relative importance of the two mechanisms depends on the ratio of the porosity increase to the initial porosity. We show that the inlet instability is likely to be important in limestone formations where the initial porosity is small and there is the possibility of a large increase in permeability. In quartz-rich sandstones, where the proportion of easily soluble material (e.g. carbonate cements) is small, the instability in the steady-state equations is dominant.
Cable subsystems characterized by long, slender, and flexible structural elements are featured in numerous engineering systems. In each of them, interaction between an individual cable and the surrounding fluid is inevitable. Such a Fluid-Structure Interaction (FSI) has received little attention in the literature, possibly due to the inherent complexity associated with fluid and structural semi-discretizations of disparate spatial dimensions. This paper proposes an embedded boundary approach for filling this gap, where the dynamics of the cable are captured by a standard finite element representation $mathcal C$ of its centerline, while its geometry is represented by a discrete surface $Sigma_h$ that is embedded in the fluid mesh. The proposed approach is built on master-slave kinematics between $mathcal C$ and $Sigma_h$, a simple algorithm for computing the motion/deformation of $Sigma_h$ based on the dynamic state of $mathcal C$, and an energy-conserving method for transferring to $mathcal C$ the loads computed on $Sigma_h$. Its effectiveness is demonstrated for two highly nonlinear applications featuring large deformations and/or motions of a cable subsystem and turbulent flows: an aerial refueling model problem, and a challenging supersonic parachute inflation problem. The proposed approach is verified using numerical data, and validated using real flight data.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
