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

Thomas Richter

Thomas Richter contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysiseess.IVphysics.comp-phphysics.flu-dynComputational Engineering, Finance, and Sciencecond-mat.soft

Trust snapshot

Quick read

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

21 published item(s)

preprint2023arXiv

Centroid adapted frequency selective extrapolation for reconstruction of lost image areas

Lost image areas with different size and arbitrary shape can occur in many scenarios such as error-prone communication, depth-based image rendering or motion compensated wavelet lifting. The goal of image reconstruction is to restore these lost image areas as close to the original as possible. Frequency selective extrapolation is a block-based method for efficiently reconstructing lost areas in images. So far, the actual shape of the lost area is not considered directly. We propose a centroid adaption to enhance the existing frequency selective extrapolation algorithm that takes the shape of lost areas into account. To enlarge the test set for evaluation we further propose a method to generate arbitrarily shaped lost areas. On our large test set, we obtain an average reconstruction gain of 1.29 dB.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Efficient lossless coding of highpass bands from block-based motion compensated wavelet lifting using JPEG 2000

Lossless image coding is a crucial task especially in the medical area, e.g., for volumes from Computed Tomography or Magnetic Resonance Tomography. Besides lossless coding, compensated wavelet lifting offers a scalable representation of such huge volumes. While compensation methods increase the details in the lowpass band, they also vary the characteristics of the wavelet coefficients, so an adaption of the coefficient coder should be considered. We propose a simple invertible extension for JPEG 2000 that can reduce the filesize for lossless coding of the highpass band by 0.8% on average with peak rate saving of 1.1%.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reconstruction of images taken by a pair of non-regular sampling sensors using correlation based matching

Multi-view image acquisition systems with two or more cameras can be rather costly due to the number of high resolution image sensors that are required. Recently, it has been shown that by covering a low resolution sensor with a non-regular sampling mask and by using an efficient algorithm for image reconstruction, a high resolution image can be obtained. In this paper, a stereo image reconstruction setup for multi-view scenarios is proposed. A scene is captured by a pair of non-regular sampling sensors and by incorporating information from the adjacent view, the reconstruction quality can be increased. Compared to a state-of-the-art single-view reconstruction algorithm, this leads to a visually noticeable average gain in PSNR of 0.74 dB.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reconstruction of Videos Taken by a Non-Regular Sampling Sensor

Recently, it has been shown that a high resolution image can be obtained without the usage of a high resolution sensor. The main idea has been that a low resolution sensor is covered with a non-regular sampling mask followed by a reconstruction of the incomplete high resolution image captured this way. In this paper, a multi-frame reconstruction approach is proposed where a video is taken by a non-regular sampling sensor and fully reconstructed afterwards. By utilizing the temporal correlation between neighboring frames, the reconstruction quality can be further enhanced. Compared to a state-of-the-art single-frame reconstruction approach, this leads to a visually noticeable gain in PSNR of up to 1.19 dB on average.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reducing Randomness of Non-Regular Sampling Masks for Image Reconstruction

Increasing spatial image resolution is an often required, yet challenging task in image acquisition. Recently, it has been shown that it is possible to obtain a high resolution image by covering a low resolution sensor with a non-regular sampling mask. Due to the masking, however, some pixel information in the resulting high resolution image is not available and has to be reconstructed by an efficient image reconstruction algorithm in order to get a fully reconstructed high resolution image. In this paper, the influence of different sampling masks with a reduced randomness of the non-regularity on the image reconstruction process is evaluated. Simulation results show that it is sufficient to use sampling masks that are non-regular only on a smaller scale. These sampling masks lead to a visually noticeable gain in PSNR compared to arbitrary chosen sampling masks which are non-regular over the whole image sensor size. At the same time, they simplify the manufacturing process and allow for efficient storage.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatio-temporal prediction in video coding by non-local means refined motion compensation

The prediction step is a very important part of hybrid video codecs. In this contribution, a novel spatio-temporal prediction algorithm is introduced. For this, the prediction is carried out in two steps. Firstly, a preliminary temporal prediction is conducted by motion compensation. Afterwards, spatial refinement is carried out for incorporating spatial redundancies from already decoded neighboring blocks. Thereby, the spatial refinement is achieved by applying Non-Local Means de-noising to the union of the motion compensated block and the already decoded blocks. Including the spatial refinement into H.264/AVC, a rate reduction of up to 14 % or respectively a gain of up to 0.7 dB PSNR compared to unrefined motion compensated prediction can be achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A neural network multigrid solver for the Navier-Stokes equations

We present the deep neural network multigrid solver (DNN-MG) that we develop for the instationary Navier-Stokes equations. DNN-MG improves computational efficiency using a judicious combination of a geometric multigrid solver and a recurrent neural network with memory. DNN-MG uses the multi-grid method to classically solve on coarse levels while the neural network corrects interpolated solutions on fine ones, thus avoiding the increasingly expensive computations that would have to be performed there. This results in a reduction in computation time through DNN-MG's highly compact neural network. The compactness results from its design for local patches and the available coarse multigrid solutions that provides a "guide" for the corrections. A compact neural network with a small number of parameters also reduces training time and data. Furthermore, the network's locality facilitates generalizability and allows one to use DNN-MG trained on one mesh domain also on different ones. We demonstrate the efficacy of DNN-MG for variations of the 2D laminar flow around an obstacle. For these, our method significantly improves the solutions as well as lift and drag functionals while requiring only about half the computation time of a full multigrid solution. We also show that DNN-MG trained for the configuration with one obstacle can be generalized to other time dependent problems that can be solved efficiently using a geometric multigrid method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A goal oriented error estimator and mesh adaptivity for sea ice simulations

For the first time we introduce an error estimator for the numerical approximation of the equations describing the dynamics of sea ice. The idea of the estimator is to identify different error contributions coming from spatial and temporal discretization as well as from the splitting in time of the ice momentum equations from further parts of the coupled system. The novelty of the error estimator lies in the consideration of the splitting error, which turns out to be dominant with increasing mesh resolution. Errors are measured in user specified functional outputs like the total sea ice extent. The error estimator is based on the dual weighted residual method that asks for the solution of an additional dual problem for obtaining sensitivity information. Estimated errors can be used to validate the accuracy of the solution and, more relevant, to reduce the discretization error by guiding an adaptive algorithm that optimally balances the mesh size and the time step size to increase the efficiency of the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Low-rank Approach for Nonlinear Parameter-dependent Fluid-structure Interaction Problems

Parameter-dependent discretizations of linear fluid-structure interaction problems can be approached with low-rank methods. When discretizing with respect to a set of parameters, the resulting equations can be translated to a matrix equation since all operators involved are linear. If nonlinear FSI problems are considered, a direct translation to a matrix equation is not possible. We present a method that splits the parameter set into disjoint subsets and, on each subset, computes an approximation of the problem related to the upper median parameter by means of the Newton iteration. This approximation is then used as initial guess for one Newton step on a subset of problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Adaptive time-step control for a monolithic multirate scheme coupling the heat and wave equation

We consider the dynamics of a parabolic and a hyperbolic equation coupled on a common interface and develop time-stepping schemes that can use different time-step sizes for each of the subproblems. The problem is formulated in a strongly coupled (monolithic) space-time framework. Coupling two different step sizes monolithically gives rise to large algebraic systems of equations where multiple states of the subproblems must be solved at once. For efficiently solving these algebraic systems, we inherit ideas from the partitioned regime and present two decoupling methods, namely a partitioned relaxation scheme and a shooting method. Furthermore, we develop an a posteriori error estimator serving as a mean for an adaptive time-stepping procedure. The goal is to optimally balance the time step sizes of the two subproblems. The error estimator is based on the dual weighted residual method and relies on the space-time Galerkin formulation of the coupled problem. As an example, we take a linear set-up with the heat equation coupled to the wave equation. We formulate the problem in a monolithic manner using the space-time framework. In numerical test cases, we demonstrate the efficiency of the solution process and we also validate the accuracy of the a posteriori error estimator and its use for controlling the time step sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An unfitted Eulerian finite element method for the time-dependent Stokes problem on moving domains

We analyse a Eulerian Finite Element method, combining a Eulerian time-stepping scheme applied to the time-dependent Stokes equations using the CutFEM approach with inf-sup stable Taylor-Hood elements for the spatial discretisation. This is based on the method introduced by Lehrenfeld \& Olshanskii [ESAIM: M2AN 53(2):585--614] in the context of a scalar convection-diffusion problems on moving domains, and extended to the non-stationary Stokes problem on moving domains by Burman, Frei \& Massing [arXiv:1910.03054 [math.NA]] using stabilised equal-order elements. The analysis includes the geometrical error made by integrating over approximated levelset domains in the discrete CutFEM setting. The method is implemented and the theoretical results are illustrated using numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient approximation of flow problems with multiple scales in time

In this article we address flow problems that carry a multiscale character in time. In particular we consider the Navier-Stokes flow in a channel on a fast scale that influences the movement of the boundary which undergoes a deformation on a slow scale in time. We derive an averaging scheme that is of first order with respect to the ratio of time-scales $ε$. In order to cope with the problem of unknown initial data for the fast scale problem, we assume near-periodicity in time. Moreover, we construct a second-order accurate time discretisation scheme and derive a complete error analysis for a corresponding simplified ODE system. The resulting multiscale scheme does not ask for the continuous simulation of the fast scale variable and shows powerful speed-ups up to 1:10000 compared to a resolved simulation. Finally, we present some numerical examples for the full Navier-Stokes system to illustrate the convergence and performance of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Falling balls in a viscous fluid with contact: Comparing numerical simulations with experimental data

We evaluate a number of different finite element approaches for fluid-structure (contact) interaction problems against data from physical experiments. For this we take the data from experiments by Hagemeier [Mendeley Data, doi: 10.17632/mf27c92nc3.1]. This consists of trajectories of single particles falling through a highly viscous fluid and rebounding off the bottom fluid tank wall. The resulting flow is in the transitional regime between creeping and turbulent flows. This type of configuration is particularly challenging for numerical methods due to the large change of the fluid domain and the contact between the wall and particle. In the numerical simulations we consider both rigid body and linear elasticity models for the falling particles. In the first case, we compare results obtained with the well established Arbitrary Lagrangian Eulerian (ALE) approach and a moving domain CutFEM method together with a simple and common approach for contact avoidance. For the full fluid-structure interaction (FSI) problem with contact, we use a fully Eulerian approach in combination with a unified FSI-contact treatment using Nitsche's method. For higher computational efficiency we use the geometrical symmetry of the experimental set up to reformulate the FSI system into two spatial dimensions. Finally, we show full three dimensional ALE computations to study the effects of small perturbations in the initial state of the particle to investigate deviations from a perfectly vertical fall observed in the experiment. The methods are implemented in open-source finite element libraries and the results are made freely available to aide reproducibility.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite Element Error Estimates on Geometrically Perturbed Domains

We develop error estimates for the finite element approximation of elliptic partial differential equations on perturbed domains, i.e. when the computational domain does not match the real geometry. The result shows that the error related to the domain can be a dominating factor in the finite element discretization error. The main result consists of $H^1-$ and $L_2-$ error estimates for the Laplace problem. Theoretical considerations are validated by a computational example.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low-rank Linear Fluid-structure Interaction Discretizations

Fluid-structure interaction models involve parameters that describe the solid and the fluid behavior. In simulations, there often is a need to vary these parameters to examine the behavior of a fluid-structure interaction model for different solids and different fluids. For instance, a shipping company wants to know how the material, a ship's hull is made of, interacts with fluids at different Reynolds and Strouhal numbers before the building process takes place. Also, the behavior of such models for solids with different properties is considered before the prototype phase. A parameter-dependent linear fluid-structure interaction discretization provides approximations for a bundle of different parameters at one step. Such a discretization with respect to different material parameters leads to a big block-diagonal system matrix that is equivalent to a matrix equation as discussed in [KressnerTobler 2011]. The unknown is then a matrix which can be approximated using a low-rank approach that represents the iterate by a tensor. This paper discusses a low-rank GMRES variant and a truncated variant of the Chebyshev iteration. Bounds for the error resulting from the truncation operations are derived. Numerical experiments show that such truncated methods applied to parameter-dependent discretizations provide approximations with relative residual norms smaller than $10^{-8}$ within a twentieth of the time used by individual standard approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the Impact of Fluid Structure Interaction in Blood Flow Simulations: Stenotic Coronary Artery Benchmark

We study the impact of using fluid-structure interactions (FSI) to simulate blood flow in a large stenosed artery. We compare typical flow configurations using Navier-Stokes in a rigid geometry setting to a fully coupled FSI model. The relevance of vascular elasticity is investigated with respect to several questions of clinical importance. Namely, we study the effect of using FSI on the wall shear stress distribution, on the Fractional Flow Reserve and on the damping effect of a stenosis on the pressure amplitude during the pulsatory cycle. The coupled problem is described in a monolithic variational formulation based on Arbitrary Lagrangian Eulerian (ALE) coordinates. For comparison, we perform pure Navier-Stokes simulations on a prestressed geometry to give a good matching of both configurations. A series of numerical simulations that cover important hemodynamical factors are presented and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Settling of spherical particles in the transitional regime

The settling process and wall impact of large spherical particles in a stagnant, highly viscous fluid has been observed by means of high-speed shadow imaging. The particles included in this study vary in size and material properties: steel, polytetrafluorethylen (PTFE), polyoxymethylen (POM), or rubber. The corresponding terminal Reynolds numbers range from 333 to 4012, covering in principle the transitional and Newton regime for drag forces. For collision studies, the Stokes number just before impact is also relevant, and lies in the range 50 < St < 2250. The settling curves obtained experimentally are compared with numerical and analytical solutions. The latter has been derived on the basis of nominal terminal velocity and relaxation time for the Stokes and Newton regimes. A deviation between experimental results and analytical solution was observed in all cases where the terminal Reynolds number is larger than 300 and smaller than 1100. In this flow regime, the settling process of the spheres is already affected at a long distance from the impact object, leading to an early but significant deceleration. Moreover, a reduced settling velocity was observed along the whole trajectory for the PTFE particles with the lowest terminal Reynolds number. All processed datasets are available via the Mendeley Data repository. Two additional effects have been observed during this study. First, a bright region was detected around all PTFE spheres. It finally was found to be due to total light reflection around the sphere, but, to the best of our knowledge, this peculiarity has never been reported before. Being purely optical, this does not effect the settling behavior for PTFE. On the other hand, partial absorption of liquid at the particle surface was observed for rubber, leading to a reduced sedimentation velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Structure Preservation for the Deep Neural Network Multigrid Solver

The simulation of partial differential equations is a central subject of numerical analysis and an indispensable tool in science, engineering and related fields. Existing approaches, such as finite elements, provide (highly) efficient tools but deep neural network-based techniques emerged in the last few years as an alternative with very promising results. We investigate the combination of both approaches for the approximation of the Navier-Stokes equations and to what extent structural properties such as divergence freedom can and should be respected. Our work is based on DNN-MG, a deep neural network multigrid technique, that we introduced recently and which uses a neural network to represent fine grid fluctuations not resolved by a geometric multigrid finite element solver. Although DNN-MG provides solutions with very good accuracy and is computationally highly efficient, we noticed that the neural network-based corrections substantially violate the divergence freedom of the velocity vector field. In this contribution, we discuss these findings and analyze three approaches to address the problem: a penalty term to encourage divergence freedom of the network output; a penalty term for the corrected velocity field; and a network that learns the stream function, i.e. the scalar potential of the divergence free velocity vector field and which hence yields by construction divergence free corrections. Our experimental results show that the third approach based on the stream function outperforms the other two and not only improves the divergence freedom but in particular also the overall fidelity of the simulation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The candy wrapper problem -- a temporal multiscale approach for pde/pde systems

We describe a temporal multiscale approach for the simulation of long-term processes with short-term influences involving partial differential equations. The specific problem under consideration is a growth process in blood vessels. The \emph{Candy Wrapper Process} describes a restenosis in a vessel that has previously be widened by inserting a stent. The development of a new stenosis takes place on a long time horizon (months) while the acting forces are mainly given by the pulsating blood flow. We describe a coupled pde model and a finite element simulation that is used as basis for our multiscale approach, which is based on averaging the long scale equation and approximating the fast scale impact by localized periodic-in-time problems. Numerical test cases in prototypical 3d configurations demonstrate the power of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Parallel time-stepping for fluid-structure interactions

We present a parallel time-stepping method for fluid-structure interactions. The interaction between the incompressible Navier-Stokes equations and a hyperelastic solid is formulated in a fully monolithic framework. Discretization in space is based on equal order finite element for all variables and a variant of the Crank-Nicolson scheme is used as second order time integrator. To accelerate the solution of the systems, we analyze a parallel-in time method. For different numerical test cases in 2d and in 3d we present the efficiency of the resulting solution approach. We also discuss some special challenges and limitations that are connected to the special structure of fluid-structure interaction problem.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Second order pressure estimates for the Crank-Nicolson discretization of the incompressible Navier-Stokes Equations

We provide optimal order pressure error estimates for the Crank-Nicolson semidiscretization of the incompressible Navier-Stokes equations. Second order estimates for the velocity error are long known, we prove that the pressure error is of the same order if considered at interval midpoints, confirming previous numerical evidence. For simplicity we first give a proof under high regularity assumptions that include nonlocal compatibility conditions for the initial data, then use smoothing techniques for a proof under reduced assumptions based on standard local conditions only.

0 citations0 reviewsNo code linkedOpen access