BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Hannes Uecker

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

ResearcherUnclaimed source record
math.APnlin.PScond-mat.mes-hallmath-phmath.MPmath.NAmath.OCphysics.opticsastro-ph.EPastro-ph.SRcond-mat.dis-nncond-mat.softcond-mat.supr-conmath.DSnlin.SINumerical Analysisphysics.chem-phphysics.comp-phphysics.flu-dyn

Catalog footprint

What is connected

20works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

20 published item(s)

preprint2020arXiv

Bending and pinching of three-phase stripes: From secondary instabilities to morphological deformations in organic photovoltaics

Optimizing the properties of the mosaic morphology of bulk heterojunction (BHJ) organic photovoltaics (OPV) is not only challenging technologically but also intriguing from the mechanistic point of view. Among the recent breakthroughs is the identification and utilization of a three-phase (donor/mixed/acceptor) BHJ, where the (intermediate) mixed-phase can inhibit morphological changes, such as phase separation. Using a mean-field approach, we reveal and distinguish, between generic mechanisms that alter through transverse instabilities the evolution of stripes: the bending (zigzag mode) and the pinching (cross-roll mode) of the donor/acceptor domains. The results are summarized in a parameter plane spanned by the mixing energy and illumination, and show that donor-acceptor mixtures with higher mixing energy are more likely to develop pinching under charge-flux boundary conditions. The latter is notorious as it leads to the formation of disconnected domains and hence to loss of charge flux. We believe that these results provide a qualitative road-map for BHJ optimization, using mixed-phase composition and therefore, an essential step toward long-lasting OPV. More broadly, the results are also of relevance to study the coexistence of multiple-phase domains in material science, such as in ion-intercalated rechargeable batteries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pattern formation with pde2path -- a tutorial

We explain some pde2path setups for pattern formation in 1D, 2D and 3D. A focus is on new pde2path functions for branch switching at steady bifurcation points of higher multiplicity, typically due to discrete symmetries, but we also review general concepts of pattern formation and their handling in pde2path, including localized patterns and homoclinic snaking, again in 1D, 2D and 3D, based on the demo sh (Swift-Hohenberg equation). Next, the demos schnakpat (a Schnakenberg reaction-diffusion system) and chemtax (a quasilinear RD system with cross-diffusion from chemotaxis) simplify and unify previous results in a simple and concise way, CH (Cahn-Hilliard) deals with mass constraints, hexex deals with (multiple) branch points of higher degeneracy in a scalar problem on a hexagonal domain, and shgc illustrates some global coupling. The demos acS, actor, schnakS and schnaktor (the Allen-Cahn and Schnakenberg models on spheres and tori) consider pattern formation on curved surfaces, cpol considers a problem of cell polarization described by bulk-surface coupling, and bruosc (Brusselator) explains how to augment autonomous systems by a time periodic forcing. Along the way we also comment on the choice of meshes and mesh adaptation, on time integration, and we give some examples of branch point continuation and Hopf point continuation to approximate stability boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Snaking branches of planar BCC fronts in the 3D Brusselator

We present results of the application of the numerical continuation and bifurcation package pde2path to the 3D Brusselator model, focusing on snaking branches of planar fronts between body centered cubes (BCCs) and the spatial homogeneous solution, and on planar fronts between BCCs and tubes (also called square prisms). These solutions also yield approximations of localized BCCs, and of BCCs embedded in a background of tubes (or vice versa). Additionally, we compute some moving fronts between lamellas and tubes. To give some theoretical background, and to aid the numerics for the full system, we use the Maxwell points for the cubic amplitude system over the BCC lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stripes on finite domains: Why the zigzag instability is only a partial story

Stationary periodic patterns are widespread in natural sciences, ranging from nano-scale electrochemical and amphiphilic systems to mesoscale fluid, chemical and biological media and to macro-scale vegetation and cloud patterns. Their formation is usually due to a primary symmetry breaking of a uniform state to stripes, often followed by secondary instabilities to form zigzag and labyrinthine patterns. These secondary instabilities are well studied under idealized conditions of an infinite domain, however, on finite domains, the situation is more subtle since the unstable modes depend also on boundary conditions. Using two prototypical models, the Swift-Hohenberg equation and the forced complex Ginzburg-Landau equation, we consider bounded domains with no flux boundary conditions transversal to the stripes, and reveal a distinct mixed-mode instability that lies in between the classical zigzag and the Eckhaus lines. This explains the stability of stripes in the mildly zigzag unstable regime, and, after crossing the mixed-mode line, the evolution of zigzag stripes in the bulk of the domain and the formation of defects near the boundaries. The results are of particular importance for problems with large time scale separation, such as bulk-heterojunction deformations in organic photovoltaic and vegetation in semi-arid regions, where early temporal transients may play an important role.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

User guide on Hopf bifurcation and time periodic orbits with pde2path

We explain the setup for using the pde2path libraries for Hopf bifurcation and continuation of branches of periodic orbits and give implementation details of the associated demo directories. See [Uecker, Comm. in Comp. Phys., 2019] for a description of the basic algorithms and the mathematical background of the examples. Additionally we explain the treatment of Hopf bifurcations in systems with continuous symmetries, including the continuation of traveling waves and rotating waves in O(2) equivariant systems as relative equilibria, the continuation of Hopf bifurcation points via extended systems, and some simple setups for the bifurcation from periodic orbits associated to critical Floquet multipliers going through +-1.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimal harvesting and spatial patterns in a semi arid vegetation system

We consider an infinite time horizon spatially distributed optimal harvesting problem for a vegetation and soil water reaction diffusion system, with rainfall as the main external parameter. By Pontryagin's maximum principle we derive the associated four component canonical system, and numerically analyze this and hence the optimal control problem in two steps. First we numerically compute a rather rich bifurcation structure of flat (spatially homogeneous) and patterned canonical steady states (FCSS and PCSS, respectively), in 1D and 2D. Then we compute time dependent solutions of the canonical system that connect to some FCSS or PCSS. The method is efficient in dealing with non-unique canonical steady states, and thus also with multiple local maxima of the objective function. It turns out that over wide parameter regimes the FCSS, i.e., spatially uniform harvesting, are not optimal. Instead, controlling the system to a PCSS yields a higher profit. Moreover, compared to (a simple model of) private optimization, the social control gives a higher yield, and vegetation survives for much lower rainfall. In addition, the computation of the optimal (social) control gives an optimal tax to incorporate into the private optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sine-Gordon solitons in networks: Scattering and transmission at vertices

We consider the sine-Gordon equation on metric graphs with simple topologies and derive vertex boundary conditions from the fundamental conservation laws, such as energy and current conservation. Traveling wave solutions for star and tree graphs are obtained analytically in the form of kink, antikink and breather solitons for a special case. It is shown that these solutions provide reflectionless soliton transmission at the graph vertex. We find the sum rule for bond-dependent coefficients making the sine-Gordon equation embedded on the graph completely integrable. For the general case the problem is solved numerically and the vertex scattering is quantified. Applications of the obtained results to Josephson junction networks, DNA double helix and elastic fibre networks are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bifurcation of Nonlinear Bloch Waves from the Spectrum in the Gross-Pitaevskii Equation

We rigorously analyze the bifurcation of stationary so called nonlinear Bloch waves (NLBs) from the spectrum in the Gross-Pitaevskii (GP) equation with a periodic potential, in arbitrary space dimensions. These are solutions which can be expressed as finite sums of quasi-periodic functions, and which in a formal asymptotic expansion are obtained from solutions of the so called algebraic coupled mode equations. Here we justify this expansion by proving the existence of NLBs and estimating the error of the formal asymptotics. The analysis is illustrated by numerical bifurcation diagrams, mostly in 2D. In addition, we illustrate some relations of NLBs to other classes of solutions of the GP equation, in particular to so called out--of--gap solitons and truncated NLBs, and present some numerical experiments concerning the stability of these solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal management and spatial patterns in a distributed shallow lake model

We present a numerical framework to treat infinite time horizon spatially distributed optimal control problems via the associated canonical system derived by Pontryagin's Maximum Principle. The basic idea is to consider the canonical system in two steps. First we perform a bifurcation analysis of canonical steady states using the continuation and bifurcation package pde2path, yielding a number of so called flat and patterned canonical steady states. In a second step we link pde2path to the two point boundary value problem solver TOM to study time dependent canonical paths to steady states having the so called saddle point property. As an example we consider a shallow lake model with diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Soliton transport in tubular networks: transmission at vertices in the shrinking limit

Soliton transport in tube-like networks is studied by solving the nonlinear Schroedinger equation (NLSE) on finite thickness ("fat") graphs. The dependence of the solution and of the reflection at vertices on the graph thickness and on the angle between its bonds is studied and related to a special case considered in our previous work, in the limit when the thickness of the graph goes to zero. It is found that both the wave function and reflection coefficient reproduce the regime of reflectionless vertex transmission studied in our previous work.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The pde2path add-on toolbox p2pOC for solving infinite time-horizon spatially distributed optimal control problems - Quickstart Guide -

p2pOC is an add-on toolbox to the Matlab package pde2path. It is aimed at the numerical solution of optimal control (OC) problems with an infinite time horizon for parabolic systems of PDE over 1D or 2D spatial domains. The basic idea is to treat the OC problem via the associated canonical system in two steps. First we use pde2path to find branches of stationary solutions of the canonical system, also called canonical steady states (CSS). In a second step we use the results and the spatial discretization of the first step to calculate the objective values of time-dependent canonical paths ending at a CSS with the so called saddle point property. This is a (typically very high dimensional) boundary value problem (BVP) in time, which we solve by combining a modification of the BVP solver TOM with a continuation algorithm in the initial states. We explain the design and usage of the package via two example problems, namely the optimal management of a distributed shallow lake model, and of a semi-arid grazing system. Both show interesting bifurcations of so called patterned CSS, and in particular the latter also a variety of patterned optimal steady states. The package (library and demos) can be downloaded at www.staff.uni-oldenburg.de/hannes.uecker/pde2path.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

pde2path - version 2.0: faster FEM, multi-parameter continuation, nonlinear boundary conditions, and periodic domains - a short manual

pdepath 2.0 is an upgrade of the continuation/bifurcation package pde2path for elliptic systems of PDEs over bounded 2D domains, based on the Matlab pdetoolbox. The new features include a more efficient use of FEM, easier switching between different single parameter continuations, genuine multi-parameter continuation (e.g., fold continuation), more efficient implementation of nonlinear boundary conditions, cylinder and torus geometries (i.e., periodic boundary conditions), and a general interface for adding auxiliary equations like mass conservation or phase equations for continuation of traveling waves. The package (library, demos, manuals) can be downloaded at www.staff.uni-oldenburg.de/hannes.uecker/pde2path

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical results for snaking of patterns over patterns in some 2D Selkov-Schnakenberg Reaction-Diffusion systems

For a Selkov--Schnakenberg model as a prototype reaction-diffusion system on two dimensional domains we use the continuation and bifurcation software pde2path to numerically calculate branches of patterns embedded in patterns, for instance hexagons embedded in stripes and vice versa, with a planar interface between the two patterns. We use the Ginzburg-Landau reduction to approximate the locations of these branches by Maxwell points for the associated Ginzburg-Landau system. For our basic model, some but not all of these branches show a snaking behaviour in parameter space, over the given computational domains. The (numerical) non-snaking behaviour appears to be related to too narrow bistable ranges with rather small Ginzburg-Landau energy differences. This claim is illustrated by a suitable generalized model. Besides the localized patterns with planar interfaces we also give a number of examples of fully localized atterns over patterns, for instance hexagon patches embedded in radial stripes, and fully localized hexagon patches over straight stripes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

pde2path - A Matlab package for continuation and bifurcation in 2D elliptic systems

pde2path is a free and easy to use Matlab continuation/bifurcation package for elliptic systems of PDEs with arbitrary many components, on general two dimensional domains, and with rather general boundary conditions. The package is based on the FEM of the Matlab pdetoolbox, and is explained by a number of examples, including Bratu's problem, the Schnakenberg model, Rayleigh-Benard convection, and von Karman plate equations. These serve as templates to study new problems, for which the user has to provide, via Matlab function files, a description of the geometry, the boundary conditions, the coefficients of the PDE, and a rough initial guess of a solution. The basic algorithm is a one parameter arclength continuation with optional bifurcation detection and branch-switching. Stability calculations, error control and mesh-handling, and some elementary time-integration for the associated parabolic problem are also supported. The continuation, branch-switching, plotting etc are performed via Matlab command-line function calls guided by the AUTO style. The software can be downloaded from www.staff.uni-oldenburg.de/hannes.uecker/pde2path, where also an online documentation of the software is provided such that in this paper we focus more on the mathematics and the example systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Well-posedness of some initial-boundary-value problems for dynamo-generated poloidal magnetic fields

Given a bounded domain $G \subset \R^d$, $d\geq 3$, we study smooth solutions of a linear parabolic equation with non-constant coefficients in $G$, which at the boundary have to $C^1$-match with some harmonic function in $\R^d \setminus \ov{G}$ vanishing at spatial infinity. This problem arises in the framework of magnetohydrodynamics if certain dynamo-generated magnetic fields are considered: For example, in the case of axisymmetry or for non-radial flow fields, the poloidal scalar of the magnetic field solves the above problem. We first investigate the Poisson problem in $G$ with the above described boundary condition as well as the associated eigenvalue problem and prove the existence of smooth solutions. As a by-product we obtain the completeness of the well-known poloidal "free decay modes" in $\R^3$ if $G$ is a ball. Smooth solutions of the evolution problem are then obtained by Galerkin approximation based on these eigenfunctions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusive mixing of periodic wave trains in reaction-diffusion systems

We consider reaction-diffusion systems on the infinite line that exhibit a family of spectrally stable spatially periodic wave trains $u_0(kx-\om t;k)$ that are parameterized by the wave number $k$. We prove stable diffusive mixing of the asymptotic states $u_0(k x+ϕ_{\pm};k)$ as $x\ra \pm\infty$ with different phases $ϕ_-\neqϕ_+$ at infinity for solutions that initially converge to these states as $x\ra \pm\infty$. The proof is based on Bloch wave analysis, renormalization theory, and a rigorous decomposition of the perturbations of these wave solutions into a phase mode, which shows diffusive behavior, and an exponentially damped remainder. Depending on the dispersion relation, the asymptotic states mix linearly with a Gaussian profile at lowest order or with a nonsymmetric non-Gaussian profile given by Burgers equation, which is the amplitude equation of the diffusive modes in the case of a nontrivial dispersion relation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Perturbation theory for plasmonic eigenvalues

We develop a perturbative approach for calculating, within the quasistatic approximation, the shift of surface resonances in response to a deformation of a dielectric volume. Our strategy is based on the conversion of the homogeneous system for the potential which determines the plasmonic eigenvalues into an inhomogeneous system for the potential's derivative with respect to the deformation strength, and on the exploitation of the corresponding compatibility condition. The resulting general expression for the first-order shift is verified for two explicitly solvable cases, and for a realistic example of a deformed nanosphere. It can be used for scanning the huge parameter space of possible shape fluctuations with only quite small computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coupled Mode Equations and Gap Solitons for the 2D Gross-Pitaevskii equation with a non-separable periodic potential

Gap solitons near a band edge of a spatially periodic nonlinear PDE can be formally approximated by solutions of Coupled Mode Equations (CMEs). Here we study this approximation for the case of the 2D Periodic Nonlinear Schrödinger / Gross-Pitaevskii Equation with a non-separable potential of finite contrast. We show that unlike in the case of separable potentials [T. Dohnal, D. Pelinovsky, and G. Schneider, J. Nonlin. Sci. {\bf 19}, 95--131 (2009)] the CME derivation has to be carried out in Bloch rather than physical coordinates. Using the Lyapunov-Schmidt reduction we then give a rigorous justification of the CMEs as an asymptotic model for reversible non-degenerate gap solitons and even potentials and provide $H^s$ estimates for this approximation. The results are confirmed by numerical examples including some new families of CMEs and gap solitons absent for separable potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Statistics for surface modes of nanoparticles with shape fluctuations

We develop a numerical method for approximating the surface modes of sphere-like nanoparticles in the quasi-static limit, based on an expansion of (the angular part of) the potentials into spherical harmonics. Comparisons of the results obtained in this manner with exact solutions and with a perturbation ansatz prove that the scheme is accurate if the shape deviations from a sphere are not too large. The method allows fast calculations for large numbers of particles, and thus to obtain statistics for nanoparticles with random shape fluctuations. As an application we present some statistics for the distribution of resonances, polariziabilities, and dipole axes for particles with random perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

An integral boundary layer equation for film flow over inclined wavy bottoms

We study the flow of an incompressible liquid film down a wavy incline. Applying a Galerkin method with only one ansatz function to the Navier-Stokes equations we derive a second order weighted residual integral boundary layer equation, which in particular may be used to describe eddies in the troughs of the wavy bottom. We present numerical results which show that our model is qualitatively and quantitatively accurate in wide ranges of parameters, and we use the model to study some new phenomena, for instance the occurrence of a short wave instability for laminar flows which does not exist over flat bottom.

0 citations0 reviewsNo code linkedOpen access