Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAsymptotic and numerical analysis of a porous medium model for transpiration-driven sap flow in trees
96
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We develop a 3D porous medium model for sap flow within a tree stem, which consists of a nonlinear parabolic partial differential equation with a suitable transpiration source term. Using an asymptotic analysis, we derive approximate series solutions for the liquid saturation and sap velocity for a general class of coefficient functions. Several important non-dimensional parameters are identified that can be used to characterize various flow regimes. We investigate the relative importance of stem aspect ratio versus anisotropy in the sapwood hydraulic conductivity, and how these two effects impact the radial and vertical components of sap velocity. The analytical results are validated by means of a second-order finite volume discretization of the governing equations, and comparisons are drawn to experimental results on Norway spruce trees.
rate research
Read More
Channel formation and branching is widely seen in physical systems where movement of fluid through a porous structure causes the spatiotemporal evolution of the medium in response to the flow, in turn causing flow pathways to evolve. We provide a simple theoretical framework that embodies this feedback mechanism in a multi-phase model for flow through a fragile porous medium with a dynamic permeability. Numerical simulations of the model show the emergence of branched networks whose topology is determined by the geometry of external flow forcing. This allows us to delineate the conditions under which splitting and/or coalescing branched network formation is favored, with potential implications for both understanding and controlling branching in soft frangible media.
We present and derive a novel double-continuum transport model based on pore-scale characteristics. Our approach relies on building a simplified unit cell made up of immobile and mobile continua. We employ a numerically resolved pore-scale velocity distribution to characterize the volume of each continuum and to define the velocity profile in the mobile continuum. Using the simplified unit cell, we derive a closed form model, which includes two effective parameters that need to be estimated: a characteristic length scale and a ratio of waiting times RD that lumps the effect of stagnant regions and escape process. To calibrate and validate our model, we rely on a set of pore-scale numerical simulation performed on a 2D disordered segregated periodic porous medium considering different initial solute distributions. Using a Global Sensitivity Analysis, we explore the impact of the two effective parameters on solute concentration profiles and thereby define a sensitivity analysis driven criterion for model calibration. The latter is compared to a classical calibration approach. Our results show that, depending on the initial condition, the mass exchange process between mobile and immobile continua impact on solute profile shape significantly. By introducing parameter RD we obtain a flexible transport model capable of interpreting both symmetric and highly skewed solute concentration profiles. We show that the effectiveness of the calibration of the two parameters closely depends on the content of information of calibration dataset and the selected objective function whose definition can be supported by of the implementation of model sensitivity analysis. By relying on a sensitivity analysis driven calibration, we are able to provide a good interpretation of the concentration profile evolution independent of the given initial condition relying on a unique set of effective parameter values.
Flows through porous media can carry suspended and dissolved materials. These sediments may deposit inside the pore-space and alter its geometry. In turn, the changing pore structure modifies the preferential flow paths, resulting in a strong coupling between structural modifications and transport characteristics. Here, we compare two different models that lead to channel obstruction as a result of subsequent deposition. The first model randomly obstructs pore-throats across the porous medium, while in the second model the pore-throat with the highest flow rate is always obstructed first. By subsequently closing pores, we find that the breakdown of the permeability follows a power-law scaling, whose exponent depends on the obstruction model. The pressure jumps that occur during the obstruction process also follow a power-law distribution with the same universal scaling exponent as the avalanche size distribution of invasion percolation, independent of the model. This result suggests that the clogging processes and invasion percolation may belong to the same universality class.
Gaseous thermal transpiration flows through a rectangular micro-channel are simulated by the direct simulation BGK (DSBGK) method. These flows are rarefied, within the slip and transitional flow regimes, which are beyond many traditional computational fluid dynamic simulation schemes, such as those based on the continuum flow assumption. The flows are very slow and thus many traditional particle simulation methods suffer large statistical noises. The adopted method is a combination of particle and gas kinetic methods and it can simulate micro-flows properly. The simulation results of mass flow rates have excellent agreement with experimental measurements. In another case of 2D channel, the DSBGK comparisons with the DSMC result and the solution of Shakhov equation are also in very good agreement. Another finding from this study is that numerical simulations by including two reservoirs at the channel ends lead to appreciable differences in simulation results of velocity and pressure distributions within the micro-channel. This is due to the inhaling and exhaling effects of reservoirs at the channel ends. Even though excluding those reservoirs may accelerate the simulations significantly by using a single channel in simulations, special attentions are needed because this treatment may over-simplify the problem, and some procedures and results may be questionable. One example is to determine the surface momentum accommodation coefficient by using analytical solution of the mass flow rate obtained in a single-channel problem without the confinement effect of reservoirs at the two ends.
Rapid and accurate simulation of cerebral aneurysm flow modifications by flow diverters (FDs) can help improving patient-specific intervention and predicting treatment outcome. However, with explicit FD devices being placed in patient-specific aneurysm model, the computational domain must be resolved around the thin stent wires, leading to high computational cost in computational fluid dynamics (CFD). Classic homogeneous porous medium (PM) methods cannot accurately predict the post-stenting aneurysmal flow field due to the inhomogeneous FD wire distributions on anatomic arteries. We propose a novel approach that models the FD flow modification as a thin inhomogeneous porous medium (iPM). It improves over classic PM approaches in that, first, FD is treated as a screen, which is more accurate than the classic Darcy-Forchheimer relation based on 3D PM. second, the pressure drop is calculated using local FD geometric parameters across an inhomogeneous PM, which is more realistic. To test its accuracy and speed, we applied the iPM technique to simulate the post stenting flow field in three patient-specific aneurysms and compared the results against CFD simulations with explicit FD devices. The iPM CFD ran 500% faster than the explicit CFD while achieving 94%-99% accuracy. Thus iPM is a promising clinical bedside modeling tool to assist endovascular interventions with FD and stents.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
