BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Barbara Wohlmuth

Barbara Wohlmuth contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.APComputational Engineering, Finance, and Sciencemath.OCphysics.geo-phComputationcond-mat.stat-mechmath.PRmath.STMathematical Softwarephysics.app-phStatistics TheoryTissues and Organs

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
14topics
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

17 published item(s)

preprint2026arXiv

Constrained Consensus-Based Optimization and Numerical Heuristics for the Few Particle Regime

Consensus-based optimization (CBO) is a versatile multi-particle optimization method for performing nonconvex and nonsmooth global optimizations in high dimensions. Proofs of global convergence in probability have been achieved for a broad class of objective functions in unconstrained optimizations. In this work we adapt the algorithm for solving constrained optimizations on compact and unbounded domains with boundary by leveraging emerging reflective boundary conditions. In particular, we close a relevant gap in the literature by providing a global convergence proof for the many-particle regime comprehensive of convergence rates. On the one hand, for the sake of minimizing running cost, it is desirable to keep the number of particles small. On the other hand, reducing the number of particles implies a diminished capability of exploration of the algorithm. Hence numerical heuristics are needed to ensure convergence of CBO in the few-particle regime. In this work, we also significantly improve the convergence and complexity of CBO by utilizing an adaptive region control mechanism and by choosing geometry-specific random noise. In particular, by combining a hierarchical noise structure with a multigrid finite element method, we are able to compute global minimizers for a constrained $p$-Allen-Cahn problem with obstacles, a very challenging variational problem.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Directional flow in perivascular networks: Mixed finite elements for reduced-dimensional models on graphs

The flow of cerebrospinal fluid through the perivascular spaces of the brain is believed to play a crucial role in eliminating toxic waste proteins. While the driving forces of this flow have been enigmatic, experiments have shown that arterial wall motion is central. In this work, we present a network model for simulating pulsatile fluid flow in perivascular networks. We establish the well-posedness of this model in the primal and dual mixed variational settings, and show how it can be discretized using mixed finite elements. Further, we utilize this model to investigate fundamental questions concerning the physical mechanisms governing perivascular fluid flow. Notably, our findings reveal that arterial pulsations can induce directional flow in branching perivascular networks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solving time-fractional differential equation via rational approximation

Fractional differential equations (FDEs) describe subdiffusion behavior of dynamical systems. Its non-local structure requires taking into account the whole evolution history during the time integration, which then possibly causes additional memory use to store the history, growing in time. An alternative to a quadrature for the history integral is to approximate the fractional kernel with the sum of exponentials, which is equivalent to considering the FDE solution as a sum of solutions to a system of ODEs. One possibility to construct this system is to approximate the Laplace spectrum of the fractional kernel with a rational function. In this paper, we use the adaptive Antoulas--Anderson (AAA) algorithm for the rational approximation of the kernel spectrum which yields only a small number of real valued poles. We propose a numerical scheme based on this idea and study its stability and convergence properties. In addition, we apply the algorithm to a time-fractional Cahn-Hilliard problem.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Determining kernels in linear viscoelasticity

In this work, we investigate the inverse problem of determining the kernel functions that best describe the mechanical behavior of a complex medium modeled by a general nonlocal viscoelastic wave equation. To this end, we minimize a tracking-type data misfit function under this PDE constraint. We perform the well-posedness analysis of the state and adjoint problems and, using these results, rigorously derive the first-order sensitivities. Numerical experiments in a three-dimensional setting illustrate the method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elasto-acoustic modelling and simulation for the seismic response of structures: The case of the Tahtalı dam in the 2020 İzmir earthquake

As a mean to assess the risk dam structures are exposed to during earthquakes, we employ an abstract mathematical, three dimensional, elasto-acoustic coupled wave-propagation model taking into account (i) the dam structure itself, embedded into (ii) its surrounding topography, (iii) different material soil layers, (iv) the seismic source as well as (v) the reservoir lake filled with water treated as an acoustic medium. As a case study for extensive numerical simulations we consider the magnitude 7 seismic event of the 30$^{\rm th}$ of October 2020 taking place in the Icarian Sea (Greece) and the Tahtali dam around 30 km from there (Turkey). A challenging task is to resolve the multiple length scales that are present due to the huge differences in size between the dam building structure and the area of interest, considered for the propagation of the earthquake. Interfaces between structures and highly non-conforming meshes on different scales are resolved by means of a discontinuous Galerkin approach. The seismic source is modeled using inversion data about the real fault plane. Ultimately, we perform a real data driven, multi-scale, full source-to-site, physics based simulation based on the discontinuous Galerkin spectral element method, which allows to precisely validate the ground motion experienced along the Tahtali dam, comparing the synthetic seismograms against actually observed ones. A comparison with a more classical computational method, using a plane wave with data from a deconvolved seismogram reading as an input, is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Impacts of peak-flow events on hyporheic denitrification potential

Subsurface flows, particularly hyporheic exchange fluxes, driven by streambed topography, permeability, channel gradient and dynamic flow conditions provide prominent ecological services such as nitrate removal from streams and aquifers. Stream flow dynamics cause strongly nonlinear and often episodic contributions of nutrient concentrations in river-aquifer systems. Using a fully coupled transient flow and reactive transport model, we investigated the denitrification potential of hyporheic zones during peak-flow events. The effects of streambed permeability, channel gradient and bedform amplitude on the spatio-temporal distribution of nitrate and dissolved organic carbon in streambeds and the associated denitrification potential were explored. Distinct peak-flow events with different intensity, duration and hydrograph shape were selected to represent a wide range of peak-flow scenarios. Our results indicated that the specific hydrodynamic characteristics of individual flow events largely determine the average positive or negative nitrate removal capacity of hyporheic zones, however the magnitude of this capacity is controlled by geomorphological settings (i.e. channel slope, streambed permeability and bedform amplitude). Specifically, events with longer duration and higher intensity were shown to promote higher nitrate removal efficiency with higher magnitude of removal efficiency in the scenarios with higher slope and permeability values. These results are essential for better assessment of the subsurface nitrate removal capacity under the influence of flow dynamics and particularly peak-flow events in order to provide tailored solutions for effective restoration of interconnected river-aquifer systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A 3D-1D coupled blood flow and oxygen transport model to generate microvascular networks

In this work, we introduce an algorithmic approach to generate microvascular networks starting from larger vessels that can be reconstructed without noticeable segmentation errors. Contrary to larger vessels, the reconstruction of fine-scale components of microvascular networks shows significant segmentation errors, and an accurate mapping is time and cost intense. Thus there is a need for fast and reliable reconstruction algorithms yielding surrogate networks having similar stochastic properties as the original ones. The microvascular networks are constructed in a marching way by adding vessels to the outlets of the vascular tree from the previous step. To optimise the structure of the vascular trees, we use Murray's law to determine the radii of the vessels and bifurcation angles. In each step, we compute the local gradient of the partial pressure of oxygen and adapt the orientation of the new vessels to this gradient. At the same time, we use the partial pressure of oxygen to check whether the considered tissue block is supplied sufficiently with oxygen. Computing the partial pressure of oxygen, we use a 3D-1D coupled model for blood flow and oxygen transport. To decrease the complexity of a fully coupled 3D model, we reduce the blood vessel network to a 1D graph structure and use a bi-directional coupling with the tissue which is described by a 3D homogeneous porous medium. The resulting surrogate networks are analysed with respect to morphological and physiological aspects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Generalized bounds for active subspaces

In this article, we consider scenarios in which traditional estimates for the active subspace method based on probabilistic Poincaré inequalities are not valid due to unbounded Poincaré constants. Consequently, we propose a framework that allows to derive generalized estimates in the sense that it enables to control the trade-off between the size of the Poincaré constant and a weaker order of the final error bound. In particular, we investigate independently exponentially distributed random variables in dimension two or larger and give explicit expressions for corresponding Poincaré constants showing their dependence on the dimension of the problem. Finally, we suggest possibilities for future work that aim for extending the class of distributions applicable to the active subspace method as we regard this as an opportunity to enlarge its usability.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stencil scaling for vector-valued PDEs on hybrid grids with applications to generalized Newtonian fluids

Matrix-free finite element implementations for large applications provide an attractive alternative to standard sparse matrix data formats due to the significantly reduced memory consumption. Here, we show that they are also competitive with respect to the run time in the low order case if combined with suitable stencil scaling techniques. We focus on variable coefficient vector-valued partial differential equations as they arise in many physical applications. The presented method is based on scaling constant reference stencils originating from a linear finite element discretization instead of evaluating the bilinear forms on-the-fly. This method assumes the usage of hierarchical hybrid grids, and it may be applied to vector-valued second-order elliptic partial differential equations directly or as a part of more complicated problems. We provide theoretical and experimental performance estimates showing the advantages of this new approach compared to the traditional on-the-fly integration and stored matrix approaches. In our numerical experiments, we consider two specific mathematical models. Namely, linear elastostatics and incompressible Stokes flow. The final example considers a non-linear shear-thinning generalized Newtonian fluid. For this type of non-linearity, we present an efficient approach to compute a regularized strain rate which is then used to define the node-wise viscosity. Depending on the compute architecture, we could observe maximum speedups of 64% and 122% compared to the on-the-fly integration. The largest considered example involved solving a Stokes problem with 12288 compute cores on the state of the art supercomputer SuperMUC-NG.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The surrogate matrix methodology: Accelerating isogeometric analysis of waves

The surrogate matrix methodology delivers low-cost approximations of matrices (i.e., surrogate matrices) which are normally computed in Galerkin methods via element-scale quadrature formulas. In this paper, the methodology is applied to a number of model problems in wave mechanics treated in the Galerkin isogeometic setting. Herein, the resulting surrogate methods are shown to significantly reduce the assembly time in high frequency wave propagation problems. In particular, the assembly time is reduced with negligible loss in solution accuracy. This paper also extends the scope of previous articles in its series by considering multi-patch discretizations of time-harmonic, transient, and nonlinear PDEs as particular use cases of the methodology. Our a priori error analysis for the Helmholtz equation demonstrates that the additional consistency error introduced by the presence of surrogate matrices is independent of the wave number. In addition, our floating point analysis establishes that the computational complexity of the methodology compares favorably to other contemporary fast assembly techniques for isogeometric methods. Our numerical experiments demonstrate clear performance gains for time-harmonic problems, both with and without the presence of perfectly matched layers. Notable speed-ups are also presented for a transient problem with a compressible neo-Hookean material.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A high-order discontinuous Galerkin method for nonlinear sound waves

We propose a high-order discontinuous Galerkin scheme for nonlinear acoustic waves on polytopic meshes. To model sound propagation with and without losses, we use Westervelt's nonlinear wave equation with and without strong damping. Challenges in the numerical analysis lie in handling the nonlinearity in the model, which involves the derivatives in time of the acoustic velocity potential, and in preventing the equation from degenerating. We rely in our approach on the Banach fixed-point theorem combined with a stability and convergence analysis of a linear wave equation with a variable coefficient in front of the second time derivative. By doing so, we derive an a priori error estimate for Westervelt's equation in a suitable energy norm for the polynomial degree $p \geq 2$. Numerical experiments carried out in two-dimensional settings illustrate the theoretical convergence results. In addition, we demonstrate efficiency of the method in a three-dimensional domain with varying medium parameters, where we use the discontinuous Galerkin approach in a hybrid way.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A statistical framework for generating microstructures of two-phase random materials: application to fatigue analysis

Random microstructures of heterogeneous materials play a crucial role in the material macroscopic behavior and in predictions of its effective properties. A common approach to modeling random multiphase materials is to develop so-called surrogate models approximating statistical features of the material. However, the surrogate models used in fatigue analysis usually employ simple microstructure, consisting of ideal geometries such as ellipsoidal inclusions, which generally does not capture complex geometries. In this paper, we introduce a simple but flexible surrogate microstructure model for two-phase materials through a level-cut of a Gaussian random field with covariance of Matérn class. Such parametrization of the covariance function allows for the representation of a few key design parameters while representing the geometry of inclusions in a more general setting for a large class of random heterogeneous two-phase media. In addition to the traditional morphology descriptors such as porosity, size and aspect ratio, it provides control of the regularity of the inclusions interface and sphericity. These parameters are estimated from a small number of real material images using Bayesian inversion. An efficient process of evaluating the samples, based on the Fast Fourier Transform, makes possible the use of Monte-Carlo methods to estimate statistical properties for the quantities of interest in a given material class. We demonstrate the overall framework of the use of the surrogate material model in application to the uncertainty quantification in fatigue analysis, its feasibility and efficiency, and its role in the microstructure design.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Local and nonlocal phase-field models of tumor growth and invasion due to ECM degradation

We present and analyze new multi-species phase-field mathematical models of tumor growth and ECM invasion. The local and nonlocal mathematical models describe the evolution of volume fractions of tumor cells, viable cells (proliferative and hypoxic cells), necrotic cells, and the evolution of MDE and ECM, together with chemotaxis, haptotaxis, apoptosis, nutrient distribution, and cell-to-matrix adhesion. We provide a rigorous proof of the existence of solutions of the coupled system with gradient-based and adhesion-based haptotaxis effects. In addition, we discuss finite element discretizations of the model, and we present the results of numerical experiments designed to show the relative importance and roles of various effects, including cell mobility, proliferation, necrosis, hypoxia, and nutrient concentration on the generation of MDEs and the degradation of the ECM.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the unsteady Darcy-Forchheimer-Brinkman equation in local and nonlocal tumor growth models

A mathematical analysis of local and nonlocal phase-field models of tumor growth is presented that includes time-dependent Darcy-Forchheimer-Brinkman models of convective velocity fields and models of long-range cell interactions. A complete existence analysis is provided. In addition, a parameter-sensitivity analysis is described that quantifies the sensitivity of key quantities of interest to changes in parameter values. Two sensitivity analyses are examined; one employing statistical variances of model outputs and another employing the notion of active subspaces based on existing observational data. Remarkably, the two approaches yield very similar conclusions on sensitivity for certain quantities of interest. The work concludes with the presentation of numerical approximations of solutions of the governing equations and results of numerical experiments on tumor growth produced using finite element discretizations of the full tumor model for representative cases.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The surrogate matrix methodology: a priori error estimation

We give the first mathematically rigorous analysis of an emerging approach to finite element analysis (see, e.g., Bauer et al. [Appl. Numer. Math., 2017]), which we hereby refer to as the surrogate matrix methodology. This methodology is based on the piece-wise smooth approximation of the matrices involved in a standard finite element discretization. In particular, it relies on the projection of smooth so-called stencil functions onto high-order polynomial subspaces. The performance advantage of the surrogate matrix methodology is seen in constructions where each stencil function uniquely determines the values of a significant collection of matrix entries. Such constructions are shown to be widely achievable through the use of locally-structured meshes. Therefore, this methodology can be applied to a wide variety of physically meaningful problems, including nonlinear problems and problems with curvilinear geometries. Rigorous a priori error analysis certifies the convergence of a novel surrogate method for the variable coefficient Poisson equation. The flexibility of the methodology is also demonstrated through the construction of novel methods for linear elasticity and nonlinear diffusion problems. In numerous numerical experiments, we demonstrate the efficacy of these new methods in a matrix-free environment with geometric multigrid solvers. In our experiments, up to a twenty-fold decrease in computation time is witnessed over the classical method with an otherwise identical implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The surrogate matrix methodology: A reference implementation for low-cost assembly in isogeometric analysis

A reference implementation of a new method in isogeometric analysis (IGA) is presented. It delivers low-cost variable-scale approximations (surrogates) of the matrices which IGA conventionally requires to be computed by element-scale quadrature. To generate surrogate matrices, quadrature must only be performed on a fraction of the elements in the computational domain. In this way, quadrature determines only a subset of the entries in the final matrix. The remaining matrix entries are computed by a simple B-spline interpolation procedure. We present the modifications and extensions required for a reference implementation in the open-source IGA software library GeoPDEs. The exposition is fashioned to help facilitate similar modifications in other contemporary software libraries.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The surrogate matrix methodology: Low-cost assembly for isogeometric analysis

A new methodology in isogeometric analysis (IGA) is presented. This methodology delivers low-cost variable-scale approximations (surrogates) of the matrices which IGA conventionally requires to be computed from element-scale quadrature formulas. To generate surrogate matrices, quadrature must only be performed on certain elements in the computational domain. This, in turn, determines only a subset of the entries in the final matrix. The remaining matrix entries are computed by a simple B-spline interpolation procedure. Poisson's equation, membrane vibration, plate bending, and Stokes' flow problems are studied. In these problems, the use of surrogate matrices has a negligible impact on solution accuracy. Because only a small fraction of the original quadrature must be performed, we are able to report beyond a fifty-fold reduction in overall assembly time in the same software. The capacity for even further speed-ups is clearly demonstrated. The implementation used here was achieved by a small number of modifications to the open-source IGA software library GeoPDEs. Similar modifications could be made to other present-day software libraries.

0 citations0 reviewsNo code linkedOpen access