Researcher profile

Joachim Mathiesen

Joachim Mathiesen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2025arXiv

A pseudo-dynamic phase-field model for brittle fracture

The enforcement of global energy conservation in phase-field fracture simulations has been an open problem for the last 25 years. Specifically, the occurrence of unstable fracture is accompanied by a loss in total potential energy, which suggests a violation of the energy conservation law. This phenomenon can occur even with purely quasi-static, displacement-driven loading conditions, where finite crack growth arises from an infinitesimal increase in load. While such behavior is typically seen in crack nucleation, it may also occur in other situations. Initial efforts to enforce energy conservation involved backtracking schemes based on global minimization, however in recent years it has become clearer that unstable fracture, being an inherently dynamic phenomenon, cannot be adequately resolved within a purely quasi-static framework. Despite this, it remains uncertain whether transitioning to a fully dynamic framework would sufficiently address the issue. In this work, we propose a pseudo-dynamic framework designed to enforce energy balance without relying on global minimization. This approach incorporates dynamic effects heuristically into an otherwise quasi-static model, allowing us to bypass solving the full dynamic linear momentum equation. It offers the flexibility to simulate crack evolution along a spectrum, ranging from full energy conservation at one extreme to maximal energy loss at the other. Using data from recent experiments, we demonstrate that our framework can closely replicate experimental load-displacement curves, achieving results that are unattainable with classical phase-field models.

preprint2020arXiv

Gender-specific behavior change following terror attacks

Terrorists use violence in pursuit of political goals. While terror often has severe consequences for victims, it remains an open question how terror attacks affect the general population. We study the behavioral response of citizens of cities affected by $7$ different terror attacks. We compare real-time mobile communication patterns in the first $24$ hours following a terror attack to the corresponding patterns on days with no terror attack. On ordinary days, the group of female and male participants have different activity patterns. Following a terror attack, however, we observe a significant increase of the gender differences. Knowledge about citizens' behavior response patterns following terror attacks may have important implications for the public response during and after an attack.

preprint2020arXiv

Onset of turbulence in channel flows with scale-invariant roughness

Using 3D direct numerical simulations of the Navier--Stokes equations, we study the effect of wall roughness on the onset of turbulence in channel flow. The dependence of the friction factor on the Reynolds number, $\mathrm{Re}$, is found to follow a generalized Forchheimer law, which interpolates between the laminar and inertial asymptotes. The transition between these two asymptotes occurs at a first critical $\mathrm{Re}$, $\mathrm{Re}_{\mathrm{c}}$, that depends nontrivially on the roughness amplitude. We identify the transition from subcritical to supercritical onset by looking at the dependence of the velocity fluctuations on Re for different roughness amplitudes. We find that this second critical $\mathrm{Re}$ is comparable in magnitude to $\mathrm{Re}_{\mathrm{c}}$, implying that transitional flow is an integral part of flow in open fractures when $\mathrm{Re}$ and the roughness amplitude are sufficiently high.

preprint2019arXiv

Substrate curvature governs texture orientation in thin films of smectic block copolymers

Self-assembly of ordered nanometer-scale patterns is interesting in itself, but its practical value depends on the ability to predict and control pattern formation. In this paper we demonstrate theoretically and numerically that engineering of extrinsic as well as intrinsic substrate geometry may provide such a controllable ordering mechanism for block copolymers films. We develop an effective two-dimensional model of thin films of striped-phase diblock copolymers on general curved substrates. The model is obtained as an expansion in the film thickness and thus takes the third dimension into account, which crucially allows us to predict the preferred orientations even in the absence of intrinsic curvature. We determine the minimum-energy textures on several curved surfaces and arrive at a general principle for using substrate curvature as an ordering field, namely that the stripes will tend to align along directions of maximal curvature.

preprint2019arXiv

Transient electrohydrodynamic flow with concentration dependent fluid properties: modelling and energy-stable numerical schemes

Transport of electrolytic solutions under influence of electric fields occurs in phenomena ranging from biology to geophysics. Here, we present a continuum model for single-phase electrohydrodynamic flow, which can be derived from fundamental thermodynamic principles. This results in a generalized Navier-Stokes-Poisson-Nernst-Planck system, where fluid properties such as density and permittivity depend on the ion concentration fields. We propose strategies for constructing numerical schemes for this set of equations, where solving the electrochemical and the hydrodynamic subproblems are decoupled at each time step. We provide time discretizations of the model that suffice to satisfy the same energy dissipation law as the continuous model. In particular, we propose both linear and non-linear discretizations of the electrochemical subproblem, along with a projection scheme for the fluid flow. The efficiency of the approach is demonstrated by numerical simulations using several of the proposed schemes.

preprint2018arXiv

Bernaise: A flexible framework for simulating two-phase electrohydrodynamic flows in complex domains

Bernaise (Binary ElectRohydrodyNAmIc SolvEr) is a flexible high-level finite element solver of two-phase electrohydrodynamic flow in complex geometries. Two-phase flow with electrolytes is relevant across a broad range of systems and scales, from 'lab-on-a-chip' devices for medical diagnostics to enhanced oil recovery at the reservoir scale. For the strongly coupled multi-physics problem, we employ a recently developed thermodynamically consistent model which combines a generalized Nernst-Planck equation for ion transport, the Poisson equation for electrostatics, the Cahn-Hilliard equation for the phase field (describing the interface separating the phases), and the Navier-Stokes equations for fluid flow. As an efficient alternative to solving the coupled system of partial differential equations in a monolithic manner, we present a linear, decoupled numerical scheme which sequentially solves the three sets of equations. The scheme is validated by comparison to limiting cases where analytical solutions are available, benchmark cases, and by the method of manufactured solution. The solver operates on unstructured meshes and is therefore well suited to handle arbitrarily shaped domains and problem set-ups where, e.g., very different resolutions are required in different parts of the domain. Bernaise is implemented in Python via the FEniCS framework, which effectively utilizes MPI and domain decomposition, and should therefore be suitable for large-scale/high-performance computing. Further, new solvers and problem set-ups can be specified and added with ease to the Bernaise framework by experienced Python users.

preprint2012arXiv

Clonal selection prevents tragedy of the commons when neighbors compete in a rock-paper-scissors game

The rock-paper-scissors game is a model example of the on-going cyclic turnover typical of many ecosystems, ranging from the terrestrial and aquatic to the microbial. Here we explore the evolution of a rock-paper-scissors system where three species compete for space. The species are allowed to mutate and change the speed by which they invade one another. In the case when all species have similar mutation rates, we observe a perpetual arms race where no single species prevails. When only two species mutate, their aggressions increase indefinitely until the ecosystem collapses and only the non-mutating species survives. Finally we show that when only one species mutates, group selection removes individual predators with the fastest growth rates, causing the growth rate of the species to stabilize. We explain this group selection quantitatively.

preprint2012arXiv

Emergence of diversity in a model ecosystem

The biological requirements for an ecosystem to develop and maintain species diversity are in general unknown. Here we consider a model ecosystem of sessile and mutually excluding organisms competing for space [Mathiesen et al. Phys. Rev. Lett. 107, 188101 (2011)]. The competition is controlled by an interaction network with fixed links chosen by a Bernoulli process. New species are introduced in the system at a predefined rate. In the limit of small introduction rates, the system becomes bistable and can undergo a phase transition from a state of low diversity to high diversity. We suggest that patches of isolated meta-populations formed by the collapse of cyclic relations are essential for the transition to the state of high diversity.

preprint2011arXiv

Ecosystems with mutually exclusive interactions self-organize to a state of high diversity

Ecological systems comprise an astonishing diversity of species that cooperate or compete with each other forming complex mutual dependencies. The minimum requirements to maintain a large species diversity on long time scales are in general unknown. Using lichen communities as an example, we propose a model for the evolution of mutually excluding organisms that compete for space. We suggest that chain-like or cyclic invasions involving three or more species open for creation of spatially separated sub-populations that subsequently can lead to increased diversity. In contrast to its non-spatial counterpart, our model predicts robust co-existence of a large number of species, in accordance with observations on lichen growth. It is demonstrated that large species diversity can be obtained on evolutionary timescales, provided that interactions between species have spatial constraints. In particular, a phase transition to a sustainable state of high diversity is identified.

preprint2010arXiv

Organizational structure and communication networks in a university environment

The ``six degrees of separation" between any two individuals on Earth has become emblematic of the 'small world' theme, even though the information conveyed via a chain of human encounters decays very rapidly with increasing chain length, and diffusion of information via this process may be very inefficient in large human organizations. The information flow on a communication network in a large organization, the University of Oslo, has been studied by analyzing e-mail records. The records allow for quantification of communication intensity across organizational levels and between organizational units (referred to as ``modules"). We find that the number of e-mails messages within modules scales with module size to the power of $1.29\pm .06$, and the frequency of communication between individuals decays exponentially with the number of links required upwards in the organizational hierarchy before they are connected. Our data also indicates that the number of messages sent by administrative units is proportional to the number of individuals at lower levels in the administrative hierarchy, and the ``divergence of information" within modules is associated with this linear relationship. The observed scaling is consistent with a hierarchical system in which individuals far apart in the organization interact little with each other and receive a disproportionate number of messages from higher levels in the administrative hierarchy.