Source author record

Jan M. Nordbotten

Jan M. Nordbotten appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2023arXiv

PoroTwin: A digital twin for a FluidFlower rig

We present a framework for integrated experiments and simulations of tracer transport in heterogeneous porous media using digital twin technology. The physical asset in our setup is a meter-scale FluidFlower rig. The digital twin consists of a traditional physics-based forward simulation tool and a correction technique which compensates for mismatches between simulation results and observations. The latter augments the range of the physics-based simulation and allows us to bridge the gap between simulation and experiments in a quantitative sense. We describe the setup of the physical and digital twin, including data transfer protocols using cloud technology. The accuracy of the digital twin is demonstrated on a case with artificially high diffusion that must be compensated by the correction approach, as well as by simulations in geologically complex media. The digital twin is then applied to control tracer transport by manipulating fluid injection and production in the experimental rig, thereby enabling two-way coupling between the physical and digital twins.

preprint2022arXiv

An Adaptive Penalty Method for Inequality Constrained Minimization Problems

The primal-dual active set method is observed to be the limit of a sequence of penalty formulations. Using this perspective, we propose a penalty method that adaptively becomes the active set method as the residual of the iterate decreases. The adaptive penalty method (APM) therewith combines the main advantages of both methods, namely the ease of implementation of penalty methods and the exact imposition of inequality constraints inherent to the active set method. The scheme can be considered a quasi-Newton method in which the Jacobian is approximated using a penalty parameter. This spatially varying parameter is chosen at each iteration by solving an auxiliary problem.

preprint2021arXiv

Modelling and discretization of flow in porous media with thin, full-tensor permeability inclusions

When modelling fluid flow in fractured reservoirs, it is common to represent the fracturesas lower-dimensional inclusions embedded in the host medium. Existing discretizationsof flow in porous media with thin inclusions assume that the principal directions of theinclusion permeability tensor are aligned with the inclusion orientation. While this mod-elling assumption works well with tensile fractures, it may fail in the context of faults,where the damage zone surrounding the main slip surface may introduce anisotropy thatis not aligned with the main fault orientation. In this paper, we introduce a generalizeddimensional reduced model which preserves full-tensor permeability effects also in theout-of-plane direction of the inclusion. The governing equations of flow for the lower-dimensional objects are obtained through vertical averaging. We present a framework fordiscretization of the resulting mixed-dimensional problem, aimed at easy adaptation ofexisting simulation tools. We give numerical examples that show the failure of existingformulations when applied to anisotropic faulted porous media, and go on to show theconvergence of our method in both 2D and 3D

preprint2019arXiv

Functional Analysis and Exterior Calculus on Mixed-Dimensional Geometries

We are interested in differential forms on mixed-dimensional geometries, in the sense of a domain containing sets of $d$-dimensional manifolds, structured hierarchically so that each $d$-dimensional manifold is contained in the boundary of one or more $d + 1$ dimensional manifolds. On any given $d$-dimensional manifold, we then consider differential operators tangent to the manifold as well as discrete differential operators (jumps) normal to the manifold. The combined action of these operators leads to the notion of a semi-discrete differential operator coupling manifolds of different dimensions. We refer to the resulting systems of equations as mixed-dimensional, which have become a popular modeling technique for physical applications including fractured and composite materials. We establish analytical tools in the mixed-dimensional setting, including suitable inner products, differential and codifferential operators, Poincaré lemma, and Poincaré--Friedrichs inequality. The manuscript is concluded by defining the mixed-dimensional minimization problem corresponding to the Hodge-Laplacian, and we show that this minimization problem is well-posed.

preprint2018arXiv

A Non-standard Model for Microbial Enhanced Oil Recovery Including the Oil-water Interfacial Area

In this work we present a non-standard model for microbial enhanced oil recovery including the oil-water interfacial area. Including the interfacial area in the model, we eliminate the hysteresis in the capillary pressure relationship. One of the characteristics that a surfactant should have, it is biological production at the oil-water interface. Therefore, we consider the production rate of surfactants not only as a function of the nutrient concentration, but also the interfacial area. To solve the model equations, we use an efficient and robust linearization scheme that considers a linear approximation of the capillary pressure gradient. A comprehensive, 1D implementation based on two-point flux approximation of the model is achieved. We consider different parameterizations for the interfacial tension and residual oil saturation reduction. Illustrative numerical simulations are presented, where we study the spatial distribution and evolution in time of the average pressure, water saturation, interfacial area, capillary pressure, residual oil saturation and bacterial, nutrient and surfactant concentrations. Inclusion of the interfacial area in the model leads to different predictions of oil recovery. The model can also be used to design new experiments contributing to a better understanding and optimization of MEOR.

preprint2016arXiv

Effects of Spatial Heterogeneity in Rainfall and Vegetation Type on Soil Moisture and Evapotranspiration

Nonlinear plant-scale interactions controlling the soil-water balance are generally not valid at larger spatial scales due to spatial heterogeneity in rainfall and vegetation type. The relationships between spatially averaged variables are hysteretic even when unique relationships are imposed at the plant scale. The characteristics of these hysteretic relationships depend on the size of the averaging area and the spatial properties of the soil, vegetation, and rainfall. We upscale the plant-scale relationships to the scale of a regional land-surface model based on simulation data obtained through explicit representation of spatial heterogeneity in rainfall and vegetation type. The proposed upscaled function improves predictions of spatially averaged soil moisture and evapotranspiration relative to the effective-parameter approach for a water-limited Texas shrubland. The degree of improvement is a function of the scales of heterogeneity and the size of the averaging area. We also find that single-valued functions fail to predict spatially averaged leakage accurately. Furthermore, the spatial heterogeneity results in scale-dependent hysteretic relationships for the statistical-dynamic and Montaldo & Albertson approaches.

preprint2016arXiv

Modeling and simulation of microbial enhanced oil recovery including interfacial area

The focus of this paper is the derivation of a non-standard model for microbial enhanced oil recovery (MEOR) that includes the interfacial area (IFA) between the oil and water. We consider the continuity equations for water and oil, a balance equation for the oil-water interface and advective-dispersive transport equations for bacteria, nutrients and surfactants. Surfactants lower the interfacial tension (IFT), which improves the oil recovery. Therefore, we include in the model parameterizations of the IFT reduction and residual oil saturation as a function of the surfactant concentration. We consider for the first time in context of MEOR, the role of IFA in enhanced oil recovery (EOR). The motivation to include the IFA in the model is to reduce the hysteresis in the capillary pressure relationship, include the effects of observed bacteria migration towards the oil-water interface and biological production of surfactants at the oil-water interface. An efficient and robust linearization scheme was implemented, in which we use an implicit scheme that considers a linear approximation of the capillary pressure gradient, resulting in an efficient and stable scheme. A comprehensive, 2D implementation based on two-point flux approximation (TPFA) has been achieved. Illustrative numerical simulations are presented. We give an explanation of the differences in the oil recovery profiles obtained when we consider the IFA and MEOR effects. The model can also be used to design new experiments in order to gain a better understanding and optimization of MEOR.

preprint2012arXiv

Efficient Simulation of Geothermal Processes in Heterogeneous Porous Media based on the Exponential Rosenbrock-Euler and Rosenbrock-type Methods

Simulation of geothermal systems is challenging due to coupled physical processes in highly heterogeneous media. Combining the exponential Rosenbrock--Euler and Rosenbrock-type methods with control-volume (two-point flux approximation) space discretizations leads to efficient numerical techniques for simulating geothermal systems. In terms of efficiency and accuracy, the exponential Rosenbrock--Euler time integrator has advantages over standard time-dicretization schemes, which suffer from time-step restrictions or excessive numerical diffusion when advection processes are dominating. Based on linearization of the equation at each time step, we make use of matrix exponentials of the Jacobian from the spatial discretization, which provide the exact solution in time for the linearized equations. This is at the expense of computing the matrix exponentials of the stiff Jacobian matrix, together with propagating a linearized system. However, using a Krylov subspace or Leja points techniques make these computations efficient. The Rosenbrock-type methods use the appropriate rational functions of the Jacobian from the spatial discretization. The parameters in these schemes are found in consistency with the required order of convergence in time. As a result, these schemes are A-stable and only a few linear systems are solved at each time step. The efficiency of the methods compared to standard time-discretization techniques are demonstrated in numerical examples.