Source author record

Alex Hansen

Alex Hansen 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

37works
15topics
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

37 published item(s)

preprint2025arXiv

A new kind of science

We discuss whether science is in the process of being transformed from a quest for causality to a quest for correlation in light of the recent development in artificial intelligence. We observe that while a blind trust in the most seductive promises of AI is surely to be avoided, a judicious combination of computer simulation based on physical insight and the machine learning ability to explore ultra-dimensional spaces, holds potential for transformative progress in the way science is going to be pursued in the years to come.

preprint2022arXiv

Steady-state two-phase flow of compressible and incompressible fluids in a capillary tube of varying radius

We study immiscible two-phase flow of a compressible and an incompressible fluid inside a capillary tube of varying radius under steady-state conditions. The incompressible fluid is Newtonian and the compressible fluid is an inviscid ideal gas. The surface tension associated with the interfaces between the two fluids introduces capillary forces that vary along the tube due to the variation in the tube radius. The interplay between effects due to the capillary forces and the compressibility results in a set of properties that are different from incompressible two-phase flow. As the fluids move towards the outlet, the bubbles of the compressible fluid grow in volume due to the decrease in pressure. The volumetric growth of the compressible bubbles makes the volumetric flow rate at the outlet higher than at the inlet. The growth is not only a function of the pressure drop across the tube, but also of the ambient pressure. Furthermore, the capillary forces create an effective threshold below which there is no flow. Above the threshold, the system shows a weak non-linearity between the flow rates and the effective pressure drop, where the non-linearity also depends on the absolute pressures across the tube.

preprint2020arXiv

Crack Localization and the Interplay between Stress Enhancement and Thermal Noise

We study the competition between thermal fluctuations and stress enhancement in the failure process of a disordered system by using a local load sharing fiber bundle model. The thermal noise is introduced by defining a failure probability that constitutes the temperature and elastic energy of the fibers. We observe that at a finite temperature and low disorder strength, the failure process, which nucleate in the absence of any thermal fluctuation, becomes spatially uncorrelated when the applied stress is sufficiently low. The dynamics of the model in this limit lies closely to the universality class of ordinary percolation. When applied stress is increased beyond a threshold value, localized fractures appear in the system that grow with time. We identify the boundary between the localized and random failure process in the space of temperature and applied stress, and find that the threshold of stress corresponding to the onset of localized crack growth increases with the increase of temperature.

preprint2020arXiv

Fatigue and failure of a polymer chain under tension

The rupture of a polymer chain maintained at temperature $T$ under fixed tension is prototypical to a wide array of systems failing under constant external strain and random perturbations. Past research focused on analytic and numerical studies of the mean rate of collapse of such a chain. Surprisingly, an analytic calculation of the probability distribution function (PDF) of collapse rates appears to be lacking. Since rare events of rapid collapse can be important and even catastrophic, we present here a theory of this distribution, with a stress on its tail of fast rates. We show that the tail of the PDF is a power law with a {\em universal} exponent that is theoretically determined. Extensive numerics validate the offered theory. Lessons pertaining to other problems of the same type are drawn.

preprint2020arXiv

Onsager-Symmetry Obeyed in Athermal Mesoscopic Systems: Two-Phase Flow in Porous Media

We compute the fluid flow time-correlation functions of incompressible, immiscible two-phase flow in porous media using a 2D network model. Given a properly chosen representative elementary volume, the flow rate distributions are Gaussian and the integrals of time correlation functions of the flows are found to converge to a finite value. The integrated cross-correlations become symmetric, obeying Onsager's reciprocal relations. These findings support the proposal of a non-equilibrium thermodynamic description for two-phase flow in porous media.

preprint2020arXiv

The Three Extreme Value Distributions: An Introductory Review

The statistical distribution of the largest value drawn from a sample of a given size has only three possible shapes: it is either a Weibull, a Fréchet or a Gumbel extreme value distributions. I describe in this short review how to relate the statistical distribution followed by the numbers in the sample to the associate extreme value distribution followed by the largest value within the sample. Nothing I present here is new. However, from experience, I have found that a simple and compact guide on this matter written for the physics community is missing.

preprint2016arXiv

A Monte Carlo Algorithm for Immiscible Two-Phase Flow in Porous Media

We present a Markov Chain Monte Carlo algorithm based on the Metropolis algorithm for simulation of the flow of two immiscible fluids in a porous medium under macroscopic steady-state conditions using a dynamical pore network model that tracks the motion of the fluid interfaces. The Monte Carlo algorithm is based on the configuration probability, where a configuration is defined by the positions of all fluid interfaces. We show that the configuration probability is proportional to the inverse of the flow rate. Using a two-dimensional network, advancing the interfaces using time integration the computational time scales as the linear system size to the fourth power, whereas the Monte Carlo computational time scales as the linear size to the second power. We discuss the strengths and the weaknesses of the algorithm.

preprint2016arXiv

A new set of equations describing immiscible two-phase flow in isotropic porous media

Based on non-equilibrium thermodynamics we derive a set of general equations relating the partial volumetric flow rates to each other and to the total volumetric flow rate in immiscible two-phase flow in porous media. These equations together with the conservation of saturation reduces the immiscible two-phase flow problem to a single-phase flow problem of a complex fluid. We discuss the new equation in terms of the relative permeability equations. We test the equations on model systems, both analytically and numerically.

preprint2016arXiv

Effective Rheology of Two-phase Flow in Three-Dimensional Porous Media: Experiment and Simulation

We present an experimental and numerical study of immiscible two-phase flow in 3-dimensional (3D) porous media to find the relationship between the volumetric flow rate ($Q$) and the total pressure difference ($ΔP$) in the steady state. We show that in the regime where capillary forces compete with the viscous forces, the distribution of capillary barriers at the interfaces effectively creates a yield threshold, making the fluids reminiscent of a Bingham viscoplastic fluid in the porous medium, introducing a threshold pressure $P_t$. In this regime, $Q$ depends quadratically on an excess pressure drop ($ΔP-P_t$). While increasing the flow-rate, there is a transition, beyond which the flow is Newtonian and the relationship is linear. In our experiments, we build a model porous medium using a column of glass beads transporting two fluids -- de-ionized water and air. For the numerical study, reconstructed 3D pore-networks from real core samples are considered and the transport of wetting and non-wetting fluids through the network are modeled by tracking the fluid interfaces with time. We find agreement between our numerical and experimental results. Our results match the mean-field results reported earlier.

preprint2015arXiv

Cracks in random brittle solids: From fiber bundles to continuum mechanics

Statistical models are essential to get a better understanding of the role of disorder in brittle disordered solids. Fiber bundle models play a special role as a paradigm, with a very good balance of simplicity and non-trivial effects. We introduce here a variant of the fiber bundle model where the load is transferred among the fibers through a very compliant membrane. This Soft Membrane fiber bundle mode reduces to the classical Local Load Sharing fiber bundle model in 1D. Highlighting the continuum limit of the model allows to compute an equivalent toughness for the fiber bundle and hence discuss nucleation of a critical defect. The computation of the toughness allows for drawing a simple connection with crack front propagation (depinning) models.

preprint2015arXiv

Dynamic wettability alteration in immiscible two-phase flow in porous media: Effect on transport properties and critical slowing down

The change in contact angles due to the injection of low salinity water or any other wettability altering agent in an oil-rich porous medium is modeled by a network model of disordered pores transporting two immiscible fluids. We introduce a dynamic wettability altering mechanism, where the time dependent wetting property of each pore is determined by the cumulative flow of water through it. Simulations are performed to reach steady-state for different possible alterations in the wetting angle ($θ$). We find that deviation from oil-wet conditions re-mobilizes the stuck clusters and increases the oil fractional flow. However, the rate of increase in the fractional flow depends strongly on $θ$ and as $θ\to 90^\circ$, a critical angle, the system shows critical slowing down which is characterized by two dynamic critical exponents.

preprint2015arXiv

The Local load sharing fiber bundle model in higher dimensions

We consider the local load sharing fiber bundle model in one to five dimensions. Depending on the breaking threshold distribution of the fibers, there is a transition where the fracture process becomes localized. In the localized phase, the model behaves as the invasion percolation model. The difference between the local load sharing fiber bundle model and the equal load sharing fiber bundle model decreases with increasing dimensionality as a power law.

preprint2013arXiv

History-independence of steady-state in simultaneous two-phase flow through porous media

It is well known that the transient behavior during drainage or imbibition in multiphase flow in porous media strongly depends on the history and initial condition of the system. However, when the steady-state regime is reached and both drainage and imbibition take place at the pore level, the influence of the evolution history and initial preparation is an open question. Here, we present an extensive experimental and numerical work investigating the history dependence of simultaneous steady-state two-phase flow through porous media. Our experimental system consists of a Hele-Shaw cell filled with glass beads which we model numerically by a network of disordered pores transporting two immiscible fluids. From the measurements of global pressure evolution, histogram of saturation and cluster-size distributions, we find that when both phases are flowing through the porous medium, the steady state does not depend on the initial preparation of the system or on the way it has been reached.

preprint2013arXiv

Local dynamics of a randomly pinned crack front: A numerical study

We investigate numerically the dynamics of crack propagation along a weak plane using a model consisting of fibers connecting a soft and a hard clamp. This bottom-up model has previously been shown to contain the competition of two crack propagation mechanisms: coalescence of damage with the front on small scales and pinned elastic line motion on large scales. We investigate the dynamical scaling properties of the model, both on small and large scale. The model results compare favorable with experimental results.

preprint2013arXiv

Universality Classes in Constrained Crack Growth

Based on an extension of the fiber bundle model we investigate numerically the motion of the crack front through a weak plane separating a soft and an infinitely stiff block. We find that there are two regimes. At large scales the motion is consistent with the pinned elastic line model and we find a roughness exponent equal to 0.39 +/- 0.04 characterizing it. At smaller scales, coalescence of holes dominates the motion rather than pinning of the crack front. We find a roughness exponent consistent with 2/3, which is the gradient percolation value. The length of the crack front is fractal on these smaller scales. We determine its fractal dimension to be 1.77 +/- 0.02, consistent with the hull of percolation clusters, 7/4. This suggests that the crack front is described by two universality classes: on large scales, the pinned elastic line one and on small scales, the percolation universality class.

preprint2012arXiv

A model for stable interfacial crack growth

We present a model for stable crack growth in a constrained geometry. The morphology of such cracks show scaling properties consistent with self affinity. Recent experiments show that there are two distinct self-affine regimes, one on small scales whereas the other at large scales. It is believed that two different physical mechanisms are responsible for this. The model we introduce aims to investigate the two mechanisms in a single system. We do find two distinct scaling regimes in the model.

preprint2012arXiv

Effective Rheology of Bubbles Moving in a Capillary Tube

We calculate the average volumetric flux versus pressure drop of bubbles moving in a single capillary tube with varying diameter, finding a square-root relation from mapping the flow equations onto that of a driven overdamped pendulum. The calculation is based on a derivation of the equation of motion of a bubble train from considering the capillary forces and the entropy production associated with the viscous flow. We also calculate the configurational probability of the positions of the bubbles.

preprint2012arXiv

Effective Rheology of Immiscible Two-Phase Flow in Porous Media

We demonstrate through numerical simulations and a mean field calculation that immiscible two-phase flow in a porous medium behaves effectively as a Bingham viscoplastic fluid. This leads to a generalized Darcy equation where the volumetric flow rate depends quadratically on an excess pressure difference in the range of flow rates where the capillary forces compete with the viscous forces. At higher rates, the flow is Newtonian.

preprint2012arXiv

Local Wettability Reversal during Steady-State Two-Phase Flow in Porous Media

We study the effect of local wettability reversal on remobilizing immobile fluid clusters in steady-state two-phase flow in porous media. We consider a two dimensional network model for porous medium and introduce a wettability alteration mechanism. A qualitative change in the steady state flow patterns, destabilizing the percolating and trapped clusters, is observed as the system wettability is varied. When capillary forces are strong a finite wettability alteration is necessary to move the system from single-phase to two-phase flow regime. For the case of both phases being mobile we find a linear relationship between fractional flow and wettability alteration.

preprint2012arXiv

Onset of Localization in Heterogeneous Interfacial Failure

We study numerically the failure of an interface joining two elastic materials under load using a fiber bundle model connected to an elastic half space. We find that the breakdown process follows the equal load sharing fiber bundle model without any detectable spatial correlations between the positions of the failing fibers until localization sets in. The onset of localization is an instability, not a phase transition. Depending on the elastic constant describing the elastic half space, localization sets in before or after the critical load causing the interface to fail completely, is reached. There is a crossover between failure due to localization or failure without spatial correlations when tuning the elastic constant, not a phase transition. Contrary to earlier claims based on models different from ours, we find that a finite fraction of fibers must fail before the critical load is attained, even in the extreme localization regime, i.e.\ for very small elastic constant. We furthermore find that the critical load remains finite for all values of the elastic constant in the limit of an infinitely large system.

preprint2012arXiv

Relation Between First Arrival Time and Permeability in Self-Affine Fractures with Areas in Contact

We demonstrate that the first arrival time in dispersive processes in self-affine fractures are governed by the same length scale characterizing the fractures as that which controls their permeability. In one-dimensional channel flow this length scale is the aperture of the bottle neck, i.e., the region having the smallest aperture. In two dimensions, the concept of a bottle neck is generalized to that of a minimal path normal to the flow. The length scale is then the average aperture along this path. There is a linear relationship between the first arrival time and this length scale, even when there is strong overlap between the fracture surfaces creating areas with zero permeability. We express the first arrival time directly in terms of the permeability.

preprint2012arXiv

Topology of Fracture Networks

We propose a mapping from fracture systems consisting of intersecting fracture sheets in three dimensions to an abstract network consisting of nodes and links. This makes it possible to analyze fracture systems with the methods developed within modern network theory. We test the mapping for two-dimensional geological fracture outcrops and find that the equivalent networks are small-world and dissasortative. By anlayzing the Discrete Fracture Network model, which is used to generate artifical fracture networks, we also find small world networks. However, the networks turn out to be assortative.

preprint2011arXiv

Anomalous Scaling and Solitary Waves in Systems with Non-Linear Diffusion

We study a non-linear convective-diffusive equation, local in space and time, which has its background in the dynamics of the thickness of a wetting film. The presence of a non-linear diffusion predicts the existence of fronts as well as shock fronts. Despite the absence of memory effects, solutions in the case of pure non-linear diffusion exhibit an anomalous sub-diffusive scaling. Due to a balance between non-linear diffusion and convection we, in particular, show that solitary waves appear. For large times they merge into a single solitary wave exhibiting a topological stability. Even though our results concern a specific equation, numerical simulations supports the view that anomalous diffusion and the solitary waves disclosed will be general features in such non-linear convective-diffusive dynamics.

preprint2011arXiv

Heat diffusion in a two-dimensional thermal fuse model

We present numerical studies of electrical breakdown in disordered materials using a two-dimensional thermal fuse model with heat diffusion. A conducting fuse is heated locally by a Joule heating term. Heat diffuses to neighbouring fuses by a diffusion term. When the temperature reaches a given threshold, the fuse breaks and turns into an insulator. The time dynamics is governed by the time scales related to the two terms, in the presence of quenched disorder in the conductances of the fuses. For the two limiting domains, when one time scale is much smaller than the other, we find that the global breakdown time $t_r$ follows $t_r\sim I^2$ and $t_r\sim L^2$, where $I$ is the applied current, and $L$ is the system size. However, such power law does not apply in the intermediate domain where the competition between the two terms produces a subtle behaviour.

preprint2010arXiv

Model for Density Waves in Gravity-Driven Granular Flow in Narrow Pipes

A gravity-driven flow of grains through a narrow pipe in vacuum is studied by means of a one-dimensional model with two coefficients of restitution. Numerical simulations show clearly how density waves form when a strikingly simple criterion is fulfilled: that dissipation due to collisions between the grains and the walls of the pipe is greater per collision than that which stems from collisions between particles. Counterintuitively, the highest flow rate is observed when the number of grains per density wave grows large. We find strong indication that the number of grains per density wave always approaches a constant as the particle number tends to infinity, and that collapse to a single wave, which was often observed also in previous simulations, occurs because the number of grains is insufficient for multiple wave formation.

preprint2010arXiv

Permeability Estimates of Self-Affine Fracture Faults Based on Generalization of the Bottle Neck Concept

We propose a method for calculating the effective permeability of two-dimensional self-affine permeability fields based on generalizing the one-dimensional concept of a bottleneck. We test the method on fracture faults where the local permeability field is given by the cube of the aperture field. The method remains accurate even when there is substantial mechanical overlap between the two fracture surfaces. The computational efficiency of the method is comparable to calculating a simple average and is more than two orders of magnitude faster than solving the Reynolds equations using a finite-difference scheme.

preprint2009arXiv

Failure Processes in Elastic Fiber Bundles

The fiber bundle model describes a collection of elastic fibers under load. the fibers fail successively and for each failure, the load distribution among the surviving fibers change. Even though very simple, the model captures the essentials of failure processes in a large number of materials and settings. We present here a review of fiber bundle model with different load redistribution mechanism from the point of view of statistics and statistical physics rather than materials science, with a focus on concepts such as criticality, universality and fluctuations. We discuss the fiber bundle model as a tool for understanding phenomena such as creep, and fatigue, how it is used to describe the behavior of fiber reinforced composites as well as modelling e.g. network failure, traffic jams and earthquake dynamics.

preprint2006arXiv

Self-Affinity in the Gradient Percolation Problem

We study the scaling properties of the solid-on-solid front of the infinite cluster in two-dimensional gradient percolation. We show that such an object is self affine with a Hurst exponent equal to 2/3 up to a cutoff-length proportional to the gradient to the power (-4/7). Beyond this length scale, the front position has the character of uncorrelated noise. Importantly, the self-affine behavior is robust even after removing local jumps of the front. The previously observed multi affinity, is due to the dominance of overhangs at small distances in the structure function. This is a crossover effect.

preprint2000arXiv

Normal Stress Distribution of Rough Surfaces in Contact

We study numerically the stress distribution on the interface between two thick elastic media bounded by interfaces that include spatially correlated asperities. The interface roughness is described using the self-affine topography that is observed over a very wide range of scales from fractures to faults. We analyse the correlation properties of the normal stress distribution when the rough surfaces have been brought into full contact. The self affinity of the rough surfaces is described by a Hurst exponent H. We find that the normal stress field is also self affine, but with a Hurst exponent H-1. Fluctations of the normal stress are shown to be important, especially at local scales with anti-persistent correlations.

preprint1999arXiv

A Fast Algorithm for Generating Long Self-Affine Profiles

We introduce a fast algorithm for generating long self-affine profiles. The algorithm, which is based on the fast wavelet transform, is faster than the conventional Fourier filtering algorithm. In addition to increased performance for large systems, the algorithm, named the wavelet filtering algorithm, a priori gives rise to profiles for which the long-range correlation extends throughout the entire system independently of the length scale.

preprint1999arXiv

Pathways in Two-State Protein Folding

The thermodynamics of proteins indicate that folding/unfolding takes place either through stable intermediates or through a two-state process without intermediates. The rather short folding times of the two-state process indicate that folding is guided. We reconcile these two seemingly contradictory observations quantitatively in a schematic model of protein folding. We propose a new dynamical transition temperature which is lower than the thermodynamic one, in qualitative agreement with in vivo measurement of protein stability using E.coli. Finally we demonstrate that our framework is easily generalized to encompass cold unfolding, and make predictions that relate the sharpness of the cold and hot unfolding transitions.

preprint1998arXiv

Determination of the Hurst Exponent by Use of Wavelet Transforms

We propose a new method for (global) Hurst exponent determination based on wavelets. Using this method, we analyze synthetic data with predefined Hurst exponents, fracture surfaces and data from economy. The results are compared with those obtained from Fourier spectral analysis. When many samples are available, the wavelet and Fourier methods are comparable in accuracy. However, when one or only a few samples are available, the wavelet method outperforms the Fourier method by a large margin.

preprint1996arXiv

Can a local repulsive potential trap an electron?

We study the classical dynamics of a charged particle in two dimensions, under the influence of a perpendicular magnetic and an in-plane electric field. We prove the surprising fact that there is a finite region in phase space that corresponds to the otherwise drifting particle being trapped by a local repulsive potential. Our result is a direct consequence of KAM-theory and, in particular, of Moser's theorem. We illustrate it by numerical phase portraits and by an analytic approximation to invariant curves.