Numerical Investigation of Two-Phase Flow through a Fault

Imperial Users onl

A Force-Balanced Control Volume Finite Element Method for Multi-Phase Porous Media Flow Modelling

Dr D. Pavlidis would like to acknowledge the support from the following research grants: Innovate UK ‘Octopus’, EPSRC ‘Reactor Core-Structure Re-location Modelling for Severe Nuclear Accidents’) and Horizon 2020 ‘In-Vessel Melt Retention’. Funding for Dr P. Salinas from ExxonMobil is gratefully acknowledged. Dr Z. Xie is supported by EPSRC ‘Multi-Scale Exploration of Multi-phase Physics in Flows’. Part funding for Prof Jackson under the TOTAL Chairs programme at Imperial College is also acknowledged. The authors would also like to acknowledge Mr Y. Debbabi for supplying analytic solutions.Peer reviewedPublisher PD

A hierarchy of models for simulating experimental results from a 3D heterogeneous porous medium

In this work we examine the dispersion of conservative tracers (bromide and fluorescein) in an experimentally-constructed three-dimensional dual-porosity porous medium. The medium is highly heterogeneous ($\sigma_Y^2=5.7$), and consists of spherical, low-hydraulic-conductivity inclusions embedded in a high-hydraulic-conductivity matrix. The bi-modal medium was saturated with tracers, and then flushed with tracer-free fluid while the effluent breakthrough curves were measured. The focus for this work is to examine a hierarchy of four models (in the absence of adjustable parameters) with decreasing complexity to assess their ability to accurately represent the measured breakthrough curves. The most information-rich model was (1) a direct numerical simulation of the system in which the geometry, boundary and initial conditions, and medium properties were fully independently characterized experimentally with high fidelity. The reduced models included; (2) a simplified numerical model identical to the fully-resolved direct numerical simulation (DNS) model, but using a domain that was one-tenth the size; (3) an upscaled mobile-immobile model that allowed for a time-dependent mass-transfer coefficient; and, (4) an upscaled mobile-immobile model that assumed a space-time constant mass-transfer coefficient. The results illustrated that all four models provided accurate representations of the experimental breakthrough curves as measured by global RMS error. The primary component of error induced in the upscaled models appeared to arise from the neglect of convection within the inclusions. Interestingly, these results suggested that the conventional convection-dispersion equation, when applied in a way that resolves the heterogeneities, yields models with high fidelity without requiring the imposition of a more complex non-Fickian model.Comment: 27 pages, 9 Figure

Investigation of suitability of the method of volume averaging for the study of heat transfer in superconducting accelerator magnet cooled by superfluid helium.

In the field of applied superconductivity, there is a growing need to better understand heat transfers in superconducting accelerator magnets. Depending on the engineering point of view looked at, either 0-D, 1-D, 2D or 3D modeling may be needed. Because of the size of these magnets, alone or coupled together, it is yet, impossible to study this numerically for computational reasons alone without simplification in the description of the geometry and the physics. The main idea of this study is to consider the interior of a superconducting accelerator magnet as a porous medium and to apply methods used in the field of por-ous media physics to obtain the equations that model heat transfers of a superconducting accelerator magnet in different configurations (steady-state, beam losses, quench, etc.) with minimal compromises to the physics and geometry. Since the interior of a superconducting magnet is made of coils, collars and yoke filled with liquid helium, creating channels that interconnect the helium inside the magnet, an upscaling method provides models that describe heat transfer at the magnet scale and are suitable for numerical studies. This paper presents concisely the method and an example of application for super-conducting accelerator magnet cooled by superfluid helium in the steady-state regime in considering the thermal point of view