Source author record

Joachim Mathiesen

Joachim Mathiesen 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

25works
17topics
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

25 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.

preprint2016arXiv

Evaluation of the Finite Element Lattice Boltzmann Method for Binary Fluid Flows

In contrast to the commonly used lattice Boltzmann method, off-lattice Boltzmann methods decouple the velocity discretization from the underlying spatial grid, thus allowing for more efficient geometric representations of complex boundaries. The current work combines characteristic-based integration of the streaming step with the free-energy based multiphase model by Lee et. al. [Journal of Computational Physics, 206 (1), 2005 ]. This allows for simulation time steps more than an order of magnitude larger than the relaxation time. Unlike previous work by Wardle et. al. [Computers and Mathematics with Applications, 65 (2), 2013 ] that integrated intermolecular forcing terms in the advection term, the current scheme applies collision and forcing terms locally for a simpler finite element formulation. A series of thorough benchmark studies reveal that this does not compromise stability and that the scheme is able to accurately simulate flows at large density and viscosity contrasts.

preprint2016arXiv

Measure of Node Similarity in Multilayer Networks

The weight of links in a network is often related to the similarity of the nodes. Here, we introduce a simple tunable measure for analysing the similarity of nodes across different link weights. In particular, we use the measure to analyze homophily in a group of 659 freshman students at a large university. Our analysis is based on data obtained using smartphones equipped with custom data collection software, complemented by questionnaire-based data. The network of social contacts is represented as a weighted multilayer network constructed from different channels of telecommunication as well as data on face-to-face contacts. We find that even strongly connected individuals are not more similar with respect to basic personality traits than randomly chosen pairs of individuals. In contrast, several socio-demographics variables have a significant degree of similarity. We further observe that similarity might be present in one layer of the multilayer network and simultaneously be absent in the other layers. For a variable such as gender, our measure reveals a transition from similarity between nodes connected with links of relatively low weight to dis-similarity for the nodes connected by the strongest links. We finally analyze the overlap between layers in the network for different levels of acquaintanceships.

preprint2015arXiv

Emergent user behavior on Twitter modelled by a stochastic differential equation

Data from the social-media site, Twitter, is used to study the fluctuations in tweet rates of brand names. The tweet rates are the result of a strongly correlated user behavior, which leads to bursty collective dynamics with a characteristic 1/f noise. Here we use the aggregated "user interest" in a brand name to model collective human dynamics by a stochastic differential equation with multiplicative noise. The model is supported by a detailed analysis of the tweet rate fluctuations and it reproduces both the exact bursty dynamics found in the data and the 1/f noise.

preprint2015arXiv

Rotation-limited growth of three dimensional body-centered cubic crystals

According to classical grain growth laws, grain growth is driven by the minimization of surface energy and will continue until a single grain prevails. These laws do not take into account the lattice anisotropy and the details of the microscopic rearrangement of mass between grains. Here we consider coarsening of body-centered cubic polycrystalline materials in three dimensions using the phase field crystal model. We observe as function of the quenching depth, a cross over between a state where grain rotation halts and the growth stagnates and a state where grains coarsen rapidly by coalescence through rotation and alignment of the lattices of neighboring grains. We show that the grain rotation per volume change of a grain follows a power law with an exponent of $-1.25$. The scaling exponent is consistent with theoretical considerations based on the conservation of dislocations.

preprint2014arXiv

Detailed analysis of the lattice Boltzmann method on unstructured grids

The lattice Boltzmann method has become a standard for efficiently solving problems in fluid dynamics. While unstructured grids allow for a more efficient geometrical representation of complex boundaries, the lattice Boltzmann methods is often implemented using regular grids. Here we analyze two implementations of the lattice Boltzmann method on unstructured grids, the standard forward Euler method and the operator splitting method. We derive the evolution of the macroscopic variables by means of the Chapman-Enskog expansion, and we prove that it yields the Navier-Stokes equation and is first order accurate in terms of the temporal discretization and second order in terms of the spatial discretization. Relations between the kinetic viscosity and the integration time step are derived for both the Euler method and the operator splitting method. Finally we suggest an improved version of the bounce-back boundary condition. We test our implementations in both standard benchmark geometries and in the pore network of a real sample of a porous rock.

preprint2014arXiv

Intermittent dislocation density fluctuations in crystal plasticity from a phase-field crystal model

Plastic deformation mediated by collective dislocation dynamics is investigated in the two-dimensional phase-field crystal model of sheared single crystals. We find that intermittent fluctuations in the dislocation population number accompany bursts in the plastic strain-rate fluctuations. Dislocation number fluctuations exhibit a power-law spectral density $1/f^2$ at high frequencies $f$. The probability distribution of number fluctuations becomes bimodal at low driving rates corresponding to a scenario where low density of defects alternate at irregular times with high population of defects. We propose a simple stochastic model of dislocation reaction kinetics that is able to capture these statistical properties of the dislocation density fluctuations as a function of shear rate.

preprint2014arXiv

Statistics of co-occurring keywords on Twitter

Online social media such as the micro-blogging site Twitter has become a rich source of real-time data on online human behaviors. Here we analyze the occurrence and co-occurrence frequency of keywords in user posts on Twitter. From the occurrence rate of major international brand names, we provide examples on predictions of brand-user behaviors. From the co-occurrence rates, we further analyze the user-perceived relationships between international brand names and construct the corresponding relationship networks. In general the user activity on Twitter is highly intermittent and we show that the occurrence rate of brand names forms a highly correlated time signal.

preprint2013arXiv

Excitable human dynamics driven by extrinsic events in massive communities

Using empirical data from a social media site (Twitter) and on trading volumes of financial securities, we analyze the correlated human activity in massive social organizations. The activity, typically excited by real-world events and measured by the occurrence rate of international brand names and trading volumes, is characterized by intermittent fluctuations with bursts of high activity separated by quiescent periods. These fluctuations are broadly distributed with an inverse cubic tail and have long-range temporal correlations with a $1/f$ power spectrum. We describe the activity by a stochastic point process and derive the distribution of activity levels from the corresponding stochastic differential equation. The distribution and the corresponding power spectrum are fully consistent with the empirical observations.

preprint2013arXiv

Modular networks of word correlations on Twitter

Complex networks are important tools for analyzing the information flow in many aspects of nature and human society. Using data from the microblogging service Twitter, we study networks of correlations in the appearance of words from three different categories, international brands, nouns and US major cities. We create networks where the strength of links is determined by a similarity measure based on the rate of coappearance of words. In comparison with the null model, where words are assumed to be uncorrelated, the heavy-tailed distribution of pair correlations is shown to be a consequence of modules of words representing similar entities.

preprint2013arXiv

Rotation induced grain growth and stagnation in phase field crystal models

We consider the grain growth and stagnation in polycrystalline microstructures. From the phase field crystal modelling of the coarsening dynamics, we identify a transition from a grain-growth stagnation upon deep quenching below the melting temperature $T_m$ to a continuous coarsening at shallower quenching near $T_m$. The grain evolution is mediated by local grain rotations. In the deep quenching regime, the grain assembly typically reaches a metastable state where the kinetic barrier for recrystallization across boundaries is too large and grain rotation with subsequent coalescence is infeasible. For quenching near $T_m$, we find that the grain growth depends on the average rate of grain rotation, and follows a power law behavior with time with a scaling exponent that depends on the quenching depth.

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

Communication dynamics in finite capacity social networks

In communication networks structure and dynamics are tightly coupled. The structure controls the flow of information and is itself shaped by the dynamical process of information exchanged between nodes. In order to reconcile structure and dynamics, a generic model, based on the local interaction between nodes, is considered for the communication in large social networks. In agreement with data from a large human organization, we show that the flow is non-Markovian and controlled by the temporal limitations of individuals. We confirm the versatility of our model by predicting simultaneously the degree-dependent node activity, the balance between information input and output of nodes and the degree distribution. Finally, we quantify the limitations to network analysis when it is based on data sampled over a finite period of time.

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.

preprint2012arXiv

Labyrinthine clustering in a spatial rock-paper-scissors ecosystem

The spatial rock-paper-scissors ecosystem, where three species interact cyclically, is a model example of how spatial structure can maintain biodiversity. We here consider such a system for a broad range of interaction rates. When one species grows very slowly, this species and its prey dominate the system by self-organizing into a labyrinthine configuration in which the third species propagates. The cluster size distributions of the two dominating species have heavy tails and the configuration is stabilized through a complex, spatial feedback loop. We introduce a new statistical measure that quantifies the amount of clustering in the spatial system by comparison with its mean field approximation. Hereby, we are able to quantitatively explain how the labyrinthine configuration slows down the dynamics and stabilizes the system.

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.

preprint2009arXiv

Drying Patterns: Sensitivity to Residual Stresses

Volume alteration in solid materials is a common cause of material failure. Here we investigate the crack formation in thin elastic layers attached to a substrate. We show that small variations in the volume contraction and substrate restraint can produce widely different crack patterns ranging from spirals to complex hierarchical networks. The networks are formed when there is no prevailing gradient in material contraction whereas spirals are formed in the presence of a radial gradient in the contraction of a thin elastic layer.

preprint2004arXiv

Dynamics of Crystal Formation in the Greenland NorthGRIP Ice Core

The North Greenland Ice Core Project (NorthGRIP) provides paleoclimatic information back to at about 120 kyr before present (Dahl-Jensen and others, 2002). Each year, precipitation on the ice sheet covers it with a new layer of snow, which gradually transforms into ice crystals as the layer sinks into the ice sheet. The size distribution of ice crystals has been measured at selected depths in the upper 880 m of the NorthGRIP ice core (Svensson and others, 2003b), which covers a time span of 5300 years. The distributions change with time toward a universal curve, indicating a common underlying physical process in the formation of crystals. We identify this process as an interplay between fragmentation of the crystals and diffusion of their grain boundaries. The process is described by a two-parameter differential equation to which we obtain the exact solution. The solution is in excellent agreement with the experimentally observed distributions.

preprint2004arXiv

Selection by pairwise comparisons with limited resources

We analyze different methods of sorting and selecting a set of objects by their intrinsic value, via pairwise comparisons whose outcome is uncertain. After discussing the limits of repeated Round Robins, two new methods are presented: The {\it ran-fil} requires no previous knowledge on the set under consideration, yet displaying good performances even in the least favorable case. The {\it min-ent} method sets a benchmark for optimal dynamic tournaments design.