Publications

Monographs

  • E. Franck, J. Fuhrmann, V. Michel-Dansac, L. Navoret, eds., Finite Volumes for Complex Applications X -- Volume 1, Elliptic and Parabolic Problems: FVCA10, Strasbourg, France, October 30, 2023 -- November 03, 2023, Invited Contributions, 432 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, 396 pages, (Collection Published), DOI 10.1007/978-3-031-40864-9 .

  • E. Franck, J. Fuhrmann, V. Michel-Dansac, L. Navoret, eds., Finite Volumes for Complex Applications X -- Volume 2, Hyperbolic and Related Problems: FVCA10, Strasbourg, France, October 30, 2023 -- November 03, 2023, 433 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, 308 pages, (Collection Published), DOI 10.1007/978-3-031-40860-1 .

Articles in Refereed Journals

  • R. Araya, C. Cárcamo, A.H. Poza, E. Vino, An adaptive stabilized finite element method for the Stokes--Darcy coupled problem, Journal of Computational and Applied Mathematics, 443 (2024), pp. 115753/1--115753/24, DOI 10.1016/j.cam.2024.115753 .
    Abstract
    For the Stokes--Darcy coupled problem, which models a fluid that flows from a free medium into a porous medium, we introduce and analyze an adaptive stabilized finite element method using Lagrange equal order element to approximate the velocity and pressure of the fluid. The interface conditions between the free medium and the porous medium are given by mass conservation, the balance of normal forces, and the Beavers--Joseph--Saffman conditions. We prove the well-posedness of the discrete problem and present a convergence analysis with optimal error estimates in natural norms. Next, we introduce and analyze a residual-based a posteriori error estimator for the stabilized scheme. Finally, we present numerical examples to demonstrate the performance and effectiveness of our scheme.

  • G.R. Barrenechea, V. John, P. Knobloch, Finite element methods respecting the discrete maximum principle for convection-diffusion equations, SIAM Review, 66 (2024), pp. 1--86, DOI 10.1137/22M1488934 .
    Abstract
    Convection-diffusion-reaction equations model the conservation of scalar quantities. From the analytic point of view, solutions of these equations satisfy, under certain conditions, maximum principles, which represent physical bounds of the solution. That the same bounds are respected by numerical approximations of the solution is often of utmost importance in practice. The mathematical formulation of this property, which contributes to the physical consistency of a method, is called the discrete maximum principle (DMP). In many applications, convection dominates diffusion by several orders of magnitude. It is well known that standard discretizations typically do not satisfy the DMP in this convection-dominated regime. In fact, in this case it turns out to be a challenging problem to construct discretizations that, on the one hand, respect the DMP and, on the other hand, compute accurate solutions. This paper presents a survey on finite element methods, with the main focus on the convection-dominated regime, that satisfy a local or a global DMP. The concepts of the underlying numerical analysis are discussed. The survey reveals that for the steady-state problem there are only a few discretizations, all of them nonlinear, that at the same time both satisfy the DMP and compute reasonably accurate solutions, e.g., algebraically stabilized schemes. Moreover, most of these discretizations have been developed in recent years, showing the enormous progress that has been achieved lately. Similarly, methods based on algebraic stabilization, both nonlinear and linear, are currently the only finite element methods that combine the satisfaction of the global DMP and accurate numerical results for the evolutionary equations in the convection-dominated scenario.

  • J.P. Thiele, Th. Wick, Numerical modeling and open-source implementation of variational partition-of-unity localizations of space-time dual-weighted residual estimators for parabolic problems, Journal of Scientific Computing, 99 (2024), pp. 1--40, DOI 10.1007/s10915-024-02485-6 .
    Abstract
    In this work, we consider space-time goal-oriented a posteriori error estimation for parabolic problems. Temporal and spatial discretizations are based on Galerkin finite elements of continuous and discontinuous type. The main objectives are the development and analysis of space-time estimators, in which the localization is based on a weak form employing a partition-of-unity. The resulting error indicators are used for temporal and spatial adaptivity. Our developments are substantiated with several numerical examples.

  • D. Abdel, N.E. Courtier, P. Farrell, Volume exclusion effects in perovskite charge transport modeling, Optical and Quantum Electronics, 55 (2023), pp. 884/1--884/14, DOI 10.1007/s11082-023-05125-9 .
    Abstract
    Due to their flexible material properties, perovskite materials are a promising candidate for many semiconductor devices such as lasers, memristors, LEDs and solar cells. For example, perovskite-based solar cells have recently become one of the fastest growing photovoltaic technologies. Unfortunately, perovskite devices are far from commercialization due to challenges such as fast degradation. Mathematical models can be used as tools to explain the behavior of such devices, for example drift-diffusion equations portray the ionic and electric motion in perovskites. In this work, we take volume exclusion effects on ion migration within a perovskite crystal lattice into account. This results in the formulation of two different ionic current densities for such a drift-diffusion model -- treating either the mobility or the diffusivity as density-dependent while the other quantity remains constant. The influence of incorporating each current density description into a model for a typical perovskite solar cell configuration is investigated numerically, through simulations performed using two different open source tools.

  • D. Abdel, C. Chainais-Hillairet, P. Farrell, M. Herda, Numerical analysis of a finite volume scheme for charge transport in perovskite solar cells, IMA Journal of Numerical Analysis, published online on 10.6.2023, DOI 10.1093/imanum/drad034 .
    Abstract
    In this paper, we consider a drift-diffusion charge transport model for perovskite solar cells, where electrons and holes may diffuse linearly (Boltzmann approximation) or nonlinearly (e.g. due to Fermi-Dirac statistics). To incorporate volume exclusion effects, we rely on the Fermi-Dirac integral of order −1 when modeling moving anionic vacancies within the perovskite layer which is sandwiched between electron and hole transport layers. After non-dimensionalization, we first prove a continuous entropy-dissipation inequality for the model. Then, we formulate a corresponding two-point flux finite volume scheme on Voronoi meshes and show an analogous discrete entropy-dissipation inequality. This inequality helps us to show the existence of a discrete solution of the nonlinear discrete system with the help of a corollary of Brouwer's fixed point theorem and the minimization of a convex functional. Finally, we verify our theoretically proven properties numerically, simulate a realistic device setup and show exponential decay in time with respect to the L2 error as well as a physically and analytically meaningful relative entropy.

  • S. Katz, A. Caiazzo, V. John, Impact of viscosity modeling on the simulation of aortic blood flow, Journal of Computational and Applied Mathematics, 425 (2023), pp. 115036/1--115036/18, DOI 10.1016/j.cam.2022.115036 .
    Abstract
    Modeling issues for the simulation of blood flow in an aortic coarctation are studied in this paper. From the physical point of view, several viscosity models for non-Newtonian fluids as well as a Newtonian fluid model will be considered. From the numerical point of view, two different turbulence models are utilized in the simulations. The impact of both, the physical and the numerical modeling, on clinically relevant biomarkers is investigated and compared.

  • S. Katz, A. Caiazzo, B. Moreau, U. Wilbrandt, J. Brüning, L. Goubergrits, V. John, Impact of turbulence modeling on the simulation of blood flow in aortic coarctation, International Journal of Numerical Methods in Biomedical Engineering, 39 (2023), pp. e3695/1--e3695/36, DOI 10.1002/cnm.3695 .
    Abstract
    Numerical simulations of pulsatile blood flow in an aortic coarctation require the use of turbulence modeling. This paper considers three models from the class of large eddy simulation (LES) models (Smagorinsky, Vreman, -model) and one model from the class of variational multiscale models (residual-based) within a finite element framework. The influence of these models on the estimation of clinically relevant biomarkers used to assess the degree of severity of the pathological condition (pressure difference, secondary flow degree, normalized flow displacement, wall shear stress) is investigated in detail. The simulations show that most methods are consistent in terms of severity indicators such as pressure difference and stenotic velocity. Moreover, using second-order velocity finite elements, different turbulence models might lead to considerably different results concerning other clinically relevant quantities such as wall shear stresses. These differences may be attributed to differences in numerical dissipation introduced by the turbulence models.

  • F. Galarce Marín, K. Tabelow, J. Polzehl, Ch.P. Papanikas, V. Vavourakis, L. Lilaj, I. Sack, A. Caiazzo, Displacement and pressure reconstruction from magnetic resonance elastography images: Application to an in silico brain model, SIAM Journal on Imaging Sciences, 16 (2023), pp. 996--1027, DOI 10.1137/22M149363X .
    Abstract
    This paper investigates a data assimilation approach for non-invasive quantification of intracranial pressure from partial displacement data, acquired through magnetic resonance elastography. Data assimilation is based on a parametrized-background data weak methodology, in which the state of the physical system tissue displacements and pressure fields is reconstructed from partially available data assuming an underlying poroelastic biomechanics model. For this purpose, a physics-informed manifold is built by sampling the space of parameters describing the tissue model close to their physiological ranges, to simulate the corresponding poroelastic problem, and compute a reduced basis. Displacements and pressure reconstruction is sought in a reduced space after solving a minimization problem that encompasses both the structure of the reduced-order model and the available measurements. The proposed pipeline is validated using synthetic data obtained after simulating the poroelastic mechanics on a physiological brain. The numerical experiments demonstrate that the framework can exhibit accurate joint reconstructions of both displacement and pressure fields. The methodology can be formulated for an arbitrary resolution of available displacement data from pertinent images. It can also inherently handle uncertainty on the physical parameters of the mechanical model by enlarging the physics-informed manifold accordingly. Moreover, the framework can be used to characterize, in silico, biomarkers for pathological conditions, by appropriately training the reduced-order model. A first application for the estimation of ventricular pressure as an indicator of abnormal intracranial pressure is shown in this contribution.

  • C. Chainais-Hillairet, R. Eymard, J. Fuhrmann, A monotone numerical flux for quasilinear convection diffusion equations, Mathematics of Computation, 93 (2024), pp. 203--231 (published online in June 2023), DOI 10.1090/mcom/3870 .

  • R. Araya, C. Cárcamo, A.H. Poza, A stabilized finite element method for the Stokes--Temperature coupled problem, Applied Numerical Mathematics. An IMACS Journal, 187 (2023), pp. 24--49, DOI 10.1016/j.apnum.2023.02.002 .
    Abstract
    In this work, we introduce and analyze a new stabilized finite element scheme for the Stokes--Temperature coupled problem. This new scheme allows equal order of interpolation to approximate the quantities of interest, i.e. velocity, pressure, temperature, and stress. We analyze an equivalent variational formulation of the coupled problem inspired by the ideas proposed in [3]. The existence of the discrete solution is proved, decoupling the proposed stabilized scheme and using the help of continuous dependence results and Brouwer's theorem under the standard assumption of sufficiently small data. Optimal convergence is proved under classic regularity assumptions of the solution. Finally, we present some numerical examples to show the quality of our scheme, in particular, we compare our results with those coming from a standard reference in geosciences described in [38].

  • D. Budáč, V. Miloš, M. Carda, M. Paidar, J. Fuhrmann, K. Bouzek, Prediction of electrical conductivity of porous composites using a simplified Monte Carlo 3D equivalent electronic circuit network model: LSM--YSZ case study, Electrochimica Acta, 457 (2023), pp. 142512/1--142512/12, DOI 10.1016/j.electacta.2023.142512 .
    Abstract
    Multiphase electric charge conductors composed of materials with various properties are widely utilized in both research and industrial applications. The composite materials include porous electrodes and other components mainly applied in fuel cell and battery technologies. In this study, a simplified Monte Carlo equivalent electronic circuit (EEC) network model is presented. In comparison to similar models, the present EEC network model allows an accurate prediction of the electrical properties of such materials, thus saving time-consuming experimental determination. The distinct feature of this EEC network model is that it requires only experimentally easily obtainable data as the input parameters: phase composition, porosity and bulk electrical conductivity of the individual constituents. During its run, the model generates a large number of artificial cubically shaped specimens based on random distribution of individual phases according to the input composition. Each of the specimens generated was modelled by a corresponding EEC network. The EEC networks were solved using Kirchhoff's laws, resulting in impedance response simulation for the prediction of composite conductivity values. The EEC network model was validated using lanthanum strontium manganite mixed with yttria-stabilized zirconia. Excellent agreement was obtained between the experimentally determined and the calculated electrical conductivity for sample porosities of 0 to 60 %. Due to its variability, the EEC network model can be suitable for a wide range of practical applications. The presented approach has high potential to save an enormous amount of experimental effort, while maintaining sufficient accuracy, when designing corresponding multiphase electrode structures.

  • R. Finn, M. O'Donovan, P. Farrell, J. Moatti, T. Streckenbach, Th. Koprucki, S. Schulz, Theoretical study of the impact of alloy disorder on carrier transport and recombination processes in deep UV (Al,Ga)N light emitters, Applied Physics Letters, 122 (2023), pp. 241104/1--241104/7, DOI 10.1063/5.0148168 .
    Abstract
    Aluminum gallium nitride [(Al,Ga)N] has gained significant attention in recent years due to its potential for highly efficient light emitters operating in the deep ultra-violet (UV) range (<280 nm). However, given that current devices exhibit extremely low efficiencies, understand- ing the fundamental properties of (Al,Ga)N-based systems is of key importance. Here, using a multi-scale simulation framework, we study the impact of alloy disorder on carrier transport, radiative and non-radiative recombination processes in a c-plane Al 0.7 Ga0.3 N/Al0.8 Ga0.2 N quantum well embedded in a p-n junction. Our calculations reveal that alloy fluctuations can open "percolative" pathways that promote transport for the electrons and holes into the quantum well region. Such an effect is neglected in conventional and widely used transport sim- ulations. Moreover, we find that the resulting increased carrier density and alloy induced carrier localization effects significantly increase non-radiative Auger--Meitner recombination in comparison to the radiative process. Thus, to suppress such non-radiative process and poten- tially related material degradation, a careful design (wider well, multi-quantum wells) of the active region is required to improve the effi- ciency of deep UV light emitters.

  • B. García-Archilla, V. John, J. Novo, Second order error bounds for POD-ROM methods based on first order divided differences, Applied Mathematics Letters, 146 (2023), pp. 108836/1--108836/7, DOI 10.1016/j.aml.2023.108836 .
    Abstract
    This note proves for the heat equation that using BDF2 as time stepping scheme in POD-ROM methods with snapshots based on difference quotients gives both the optimal second order error bound in time and pointwise estimates.

  • B. García-Archilla, V. John, J. Novo, POD-ROMs for incompressible flows including snapshots of the temporal derivative of the full order solution, SIAM Journal on Numerical Analysis, 61 (2023), pp. 1340--1368, DOI 10.1137/22M1503853 .
    Abstract
    In this paper we study the influence of including snapshots that approach the velocity time derivative in the numerical approximation of the incompressible Navier--Stokes equations by means of proper orthogonal decomposition (POD) methods. Our set of snapshots includes the velocity approximation at the initial time from a full order mixed finite element method (FOM) together with approximations to the time derivative at different times. The approximation at the initial velocity can be replaced by the mean value of the velocities at the different times so that when implementing the method to the fluctuations, as done mostly in practice, only approximations to the time derivatives are included in the set of snapshots. For the POD method we study the differences between projecting onto L2 and H1. In both cases pointwise in time error bounds can be proved. Including grad-div stabilization in both the FOM and the POD methods, error bounds with constants independent of inverse powers of the viscosity can be obtained.

  • B. García-Archilla , V. John, S. Katz, J. Novo, POD-ROMs for incompressible flows including snapshots of the temporal derivative of the full order solution: Error bounds for the pressure, Journal of Numerical Mathematics, published online on 26.08.2023, DOI 10.1515/jnma-2023-0039 .
    Abstract
    Reduced order methods (ROMs) for the incompressible Navier?Stokes equations, based on proper orthogonal decomposition (POD), are studied that include snapshots which approach the temporal derivative of the velocity from a full order mixed finite element method (FOM). In addition, the set of snapshots contains the mean velocity of the FOM. Both the FOM and the POD-ROM are equipped with a grad-div stabilization. A velocity error analysis for this method can be found already in the literature. The present paper studies two different procedures to compute approximations to the pressure and proves error bounds for the pressure that are independent of inverse powers of the viscosity. Numerical studies support the analytic results and compare both methods.

  • A. Jha, V. John, P. Knobloch, Adaptive grids in the context of algebraic stabilizations for convection-diffusion-reaction equations, SIAM Journal on Scientific Computing, 45 (2023), pp. B564--B589, DOI 10.1137/21M1466360 .
    Abstract
    Three algebraically stabilized finite element schemes for discretizing convection-diffusion-reaction equations are studied on adaptively refined grids. These schemes are the algebraicflux correction (AFC) scheme with the Kuzmin limiter, the AFC scheme with the Barrenechea-John-Knobloch limiter, and the recently proposed monotone upwind--type algebraically stabilizedmethod. Both conforming closure of the refined grids and grids with hanging vertices are considered.A nonstandard algorithmic step becomes necessary before these schemes can be applied on gridswith hanging vertices. The assessment of the schemes is performed with respect to the satisfactionof the global discrete maximum principle, the accuracy, e.g., smearing of layers, and the efficiency insolving the corresponding nonlinear problems.

  • A. Jha, O. Pártl, N. Ahmed, D. Kuzmin, An assessment of solvers for algebraically stabilized discretizations of convection-diffusion-reaction equations, Journal of Numerical Mathematics, 31 (2023), pp. 79--103, DOI 10.1515/jnma-2021-0123 .
    Abstract
    We consider flux-corrected finite element discretizations of 3D convection-dominated transport problems and assess the computational efficiency of algorithms based on such approximations. The methods under investigation include flux-corrected transport schemes and monolithic limiters. We discretize in space using a continuous Galerkin method and P1 or Q1 finite elements. Time integration is performed using the Crank-Nicolson method or an explicit strong stability preserving Runge-Kutta method. Nonlinear systems are solved using a fixed-point iteration method, which requires solution of large linear systems at each iteration or time step. The great variety of options in the choice of discretization methods and solver components calls for a dedicated comparative study of existing approaches. To perform such a study, we define new 3D test problems for time dependent and stationary convection-diffusion-reaction equations. The results of our numerical experiments illustrate how the limiting technique, time discretization and solver impact on the overall performance.

  • P. Ral, A.K. Giri, V. John, Instantaneous gelation and non-existence of weak solutions for the Oort--Hulst--Safronov coagulation model, Proceedings of The Royal Society of London. Series A. Mathematical, Physical and Engineering Sciences, 479 (2023), pp. 20220385/1--20220385/13, DOI 10.1098/rspa.2022.0385 .
    Abstract
    The possible occurrence of instantaneous gelation to the Oort--Hulst--Safronov (OHS) coagulation equation is investigated for a certain class of unbounded coagulation kernels. The existence of instantaneous gelation is confirmed by showing the non-existence of mass-conserving weak solutions. Finally, it is shown that for such kernels, there is no weak solution to the OHS coagulation equation at any time interval.

  • B. Spetzler, D. Abdel, F. Schwierz, M. Ziegler, P. Farrell, The role of vacancy dynamics in two-dimensional memristive devices, Advanced Electronic Materials, published online on 08.11.2023, DOI 10.1002/aelm.202300635 .
    Abstract
    Two-dimensional (2D) layered transition metal dichalcogenides (TMDCs) are promising memristive materials for neuromorphic computing systems as they could solve the problem of the excessively high energy consumption of conventional von Neumann computer architectures. Despite extensive experimental work, the underlying switching mechanisms are still not understood, impeding progress in material and device functionality. This study reveals the dominant role of mobile defects in the switching dynamics of 2D TMDC materials. The switching process is governed by the formation and annihilation dynamics of a local vacancy depletion zone. Moreover, minor changes in the interface potential barriers cause fundamentally different device behavior previously thought to originate from multiple mechanisms. The key mechanisms are identified with a charge transport model for electrons, holes, and ionic point defects, including image-charge-induced Schottky barrier lowering (SBL). The model is validated by comparing simulations to measurements for various 2D MoS2-based devices, strongly corroborating the relevance of vacancies in TMDC devices and offering a new perspective on the switching mechanisms. The insights gained from this study can be used to extend the functional behavior of 2D TMDC memristive devices in future neuromorphic computing applications.

  • D. Frerichs-Mihov, L. Henning, V. John, Using deep neural networks for detecting spurious oscillations in discontinuous Galerkin solutions of convection-dominated convection-diffusion equations, Journal of Scientific Computing, 97 (2023), pp. 36/1--36/27, DOI 10.1007/s10915-023-02335-x .
    Abstract
    Standard discontinuous Galerkin (DG) finite element solutions to convection-dominated convection-diffusion equations usually possess sharp layers but also exhibit large spurious oscillations. Slope limiters are known as a post-processing technique to reduce these unphysical values. This paper studies the application of deep neural networks for detecting mesh cells on which slope limiters should be applied. The networks are trained with data obtained from simulations of a standard benchmark problem with linear finite elements. It is investigated how they perform when applied to discrete solutions obtained with higher order finite elements and to solutions for a different benchmark problem.

  • V. John, P. Knobloch, U. Wilbrandt, A posteriori optimization of parameters in stabilized methods for convection-diffusion problems -- Part II, Journal of Computational and Applied Mathematics, 428 (2023), pp. 115167/1--115167/17, DOI 10.1016/j.cam.2023.115167 .
    Abstract
    Extensions of algorithms for computing optimal stabilization parameters in finite element methods for convection-diffusion equations are presented. These extensions reduce the dimension of the control space, in comparison to available methods, and thus address the long computing times of these methods. One method is proposed that considers only relevant mesh cells, another method that uses groups of mesh cells, and the combination of both methods is also studied. The incorporation of these methods within a gradient-based optimization procedure, via solving an adjoint problem, is explained. Numerical studies provide impressions on the gain of efficiency as well as on the loss of accuracy if control spaces with reduced dimensions are utilized.

  • CH. Merdon, W. Wollner, Pressure-robustness in the context of optimal control, SIAM Journal on Control and Optimization, 61 (2023), pp. 342--360, DOI 10.1137/22M1482603 .
    Abstract
    This paper studies the benefits of pressure-robust discretizations in the scope of optimal control of incompressible flows. Gradient forces that may appear in the data can have a negative impact on the accuracy of state and control and can only be correctly balanced if their L2-orthogonality onto discretely divergence-free test functions is restored. Perfectly orthogonal divergence-free discretizations or divergence-free reconstructions of these test functions do the trick and lead to much better analytic a priori estimates that are also validated in numerical examples.

  • O. Pártl, U. Wilbrandt, J. Mura, A. Caiazzo, Reconstruction of flow domain boundaries from velocity data via multi-step optimization of distributed resistance, Computers & Mathematics with Applications. An International Journal, 129 (2023), pp. 11--33, DOI 10.1016/j.camwa.2022.11.006 .
    Abstract
    We reconstruct the unknown shape of a flow domain using partially available internal velocity measurements. This inverse problem is motivated by applications in cardiovascular imaging where motion-sensitive protocols, such as phase-contrast MRI, can be used to recover three-dimensional velocity fields inside blood vessels. In this context, the information about the domain shape serves to quantify the severity of pathological conditions, such as vessel obstructions. We consider a flow modeled by a linear Brinkman problem with a fictitious resistance accounting for the presence of additional boundaries. To reconstruct these boundaries, we employ a multi-step gradient-based variational method to compute a resistance that minimizes the difference between the computed flow velocity and the available data. Afterward, we apply different post-processing steps to reconstruct the shape of the internal boundaries. To limit the overall computational cost, we use a stabilized equal-order finite element method. We prove the stability and the well-posedness of the considered optimization problem. We validate our method on three-dimensional examples based on synthetic velocity data and using realistic geometries obtained from cardiovascular imaging.

Contributions to Collected Editions

  • L. Ermoneit, B. Schmidt, J. Fuhrmann, Th. Koprucki, L.R. Schreiber, M. Kantner, Simulation of single-electron shuttling for spin-qubit transport in a SiGe quantum bus, in: Book of Abstracts of the International Workshop on Computational Nanotechnology 2023 (IWCN 2023), X. Orios Plaedvall, G. Abadal Berini, X. Cartoixà Soler, A. Cummings, C.F. Destefani, D. Jiménez Jiménez, J. Mart'in Mart'inez, R. Rodr'iguez Mart'inez, A. Benali, eds., pp. 88-89.

  • C. Chainais-Hillairet, R. Eymard, J. Fuhrmann, An approximate two-point Dirichlet flux for quasilinear convection diffusion equations, in: Finite Volumes for Complex Applications X -- Volume 1, Elliptic and Parabolic Problems: FVCA10, Strasbourg, France, October 30, 2023 -- November 03, 2023, Invited Contributions, E. Franck, J. Fuhrmann, V. Michel-Dansac, L. Navoret, eds., 432 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, pp. 225--233, DOI 10.1007/978-3-031-40864-9_17 .

  • R. Finn, M. O'Donovan, P. Farrell, T. Streckenbach, J. Moatti, Th. Koprucki, S. Schulz, Theoretical investigation of carrier transport and recombination processes for deep UV (Al,Ga)N light emitters, in: 23nd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), P. Bardella, A. Tibaldi, eds., IEEE, 2023, pp. 83--84, DOI 10.1109/NUSOD59562.2023.10273485 .
    Abstract
    We present a theoretical study on the impact of alloy disorder on carrier transport and recombination rates in an (Al,Ga)N single quantum well based LED operating in the deep UV spectral range. Our calculations indicate that alloy fluctuations enable percolative pathways which can result in improved carrier injection into the well, but may also increase carrier leakage from the well. Additionally, we find that alloy disorder induces carrier localization effects, a feature particularly noticeable for the holes. These localization effects can lead to locally increased carrier densities when compared to a virtual crystal approximation which neglects alloy disorder. We observe that both radiative and non-radiative recombination rates are increased. Our calculations also indicate that Auger--Meitner recombination increases faster than the radiative rate, based on a comparison with a virtual crystal approximation.

  • S. Matera, Ch. Merdon, D. Runge, Reduced basis approach for convection-diffusion equations with non-linear boundary reaction conditions, in: Finite Volumes for Complex Applications X -- Volume 1, Elliptic and Parabolic Problems: FVCA10, Strasbourg, France, October 30, 2023 -- November 03, 2023, Invited Contributions, E. Franck, J. Fuhrmann, V. Michel-Dansac, L. Navoret, eds., 432 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, pp. 335--343, DOI 10.1007/978-3-031-40864-9_28 .
    Abstract
    This paper aims at an efficient strategy to solve drift-diffusion problems with non-linear boundary conditions as they appear, e.g., in heterogeneous catalysis. Since the non-linearity only involves the degrees of freedom along (a part of) the boundary, a reduced basis ansatz is suggested that computes discrete Green?s-like functions for the present drift-diffusion operator such that the global non-linear problem reduces to a smaller non-linear problem for a boundary method. The computed basis functions are completely independent of the non-linearities. Thus, they can be reused for problems with the same differential operator and geometry. Corresponding scenarios might be inverse problems in heterogeneous catalysis but also modeling the effect of different catalysts in the same reaction chamber. The strategy is explained for a mass-conservative finite volume method and demonstrated on a simple numerical example for catalytic CO oxidation.

  • J. Fuhrmann, B. Gaudeul, Ch. Keller, Two entropic finite volume schemes for a Nernst--Planck--Poisson system with ion volume constraints, in: Finite Volumes for Complex Applications X -- Volume 1, Elliptic and Parabolic Problems: FVCA10, Strasbourg, France, October 30, 2023 -- November 03, 2023, Invited Contributions, E. Franck , J. Fuhrmann, V. Michel-Dansac, L. Navoret, eds., Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, pp. 285--294, DOI 10.1007/978-3-031-40864-9_23 .

Preprints, Reports, Technical Reports

  • L. Ermoneit, B. Schmidt, Th. Koprucki, J. Fuhrmann, T. Breiten, A. Sala, N. Ciroth, R. Xue, L.R. Schreiber, M. Kantner, Optimal control of conveyor-mode spin-qubit shuttling in a Si/SiGe quantum bus in the presence of charged defects, Preprint no. 3082, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3082 .
    PDF (9473 kByte)

  • D. Frerichs-Mihov, M. Zainelabdeen, V. John, On collocation points for physics-informed neural networks applied to convection-dominated convection-diffusion problems, Preprint no. 3074, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3074 .
    Abstract, PDF (623 kByte)
    In recent years physics-informed neural networks (PINNs) for approximating the solution to (initial-)boundary value problems gained a lot of interest. PINNs are trained to minimize several residuals of the problem in collocation points. In this work we tackle convection-dominated convection-diffusion problems, whose solutions usually possess layers, which are small regions where the solution has a steep gradient. Inspired by classical Shishkin meshes, we compare hard- constrained PINNs trained with layer-adapted collocation points with ones trained with equispaced and uniformly randomly chosen points. We observe that layer-adapted points work the best for a problem with an interior layer and the worst for a problem with boundary layers. For both problems at most acceptable solutions can be obtained with PINNs.

  • CH. Keller, J. Fuhrmann, M. Landstorfer, B. Wagner, A model framework for ion channels with selectivity filters based on continuum non-equilibrium thermodynamics, Preprint no. 3072, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3072 .
    Abstract, PDF (7287 kByte)
    A mathematical model framework to describe ion transport in nanopores is presented. The model is based on non-equilibrium thermodynamics and considers finite size effects, solvation phenomena as well as the electrical charges of membrane surfaces and channel proteins. Par- ticular emphasis is placed on the consistent modelling of the selectivity filter in the pore. It is treated as an embedded domain in which the constituents can change their chemical properties. The diffusion process through the filter is governed by an independent diffusion coefficient and at the interfaces, de- and resolvation reactions are introduced as Neumann interface conditions. The evolution of the molar densities is described by drift-diffusion equations, where the fluxes depend on the gradient of the chemical potentials and the electric force. The chemical potentials depend on the molar fractions and on the pressure in the electrolyte and accounts for solvation effects. The framework allows the calculation of current-voltage relations for a variety of chan- nel properties and ion concentrations. We compare our model framework to experimental results for calcium-selective ion channels and show the general validity of our approach. Our parameter studies show that calcium and sodium currents are proportional to the surface charge in the se- lectivity filter and to the diffusion coefficients of the ions. Moreover, they show that the negative charges inside the pore have a decisive influence on the selectivity of divalent over monovalent ions.

  • Y. Hadjimichael, Ch. Merdon, M. Liero, P. Farrell, An energy-based finite-strain model for 3D heterostructured materials and its validation by curvature analysis, Preprint no. 3064, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3064 .
    Abstract, PDF (6517 kByte)
    This paper presents a comprehensive study of the intrinsic strain response of 3D het- erostructures arising from lattice mismatch. Combining materials with different lattice constants induces strain, leading to the bending of these heterostructures. We propose a model for nonlinear elastic heterostructures such as bimetallic beams or nanowires that takes into account local prestrain within each distinct material region. The resulting system of partial differential equations (PDEs) in Lagrangian coordinates incorporates a nonlinear strain and a linear stress-strain relationship governed by Hooke?s law. To validate our model, we apply it to bimetallic beams and hexagonal hetero-nanowires and perform numerical simulations using finite element methods (FEM). Our simulations ex- amine how these structures undergo bending under varying material compositions and cross-sectional geometries. In order to assess the fidelity of the model and the accuracy of simulations, we compare the calculated curvature with analytically derived formula- tions. We derive these analytical expressions through an energy-based approach as well as a kinetic framework, adeptly accounting for the lattice constant mismatch present at each compound material of the heterostructures. The outcomes of our study yield valuable insights into the behavior of strained bent heterostructures. This is particularly significant as the strain has the potential to influence the electronic band structure, piezoelectricity, and the dynamics of charge carriers.

  • D. Frerichs-Mihov, L. Henning, V. John, On loss functionals for physics-informed neural networks for convection-dominated convection-diffusion problems, Preprint no. 3063, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3063 .
    Abstract, PDF (1989 kByte)
    In the convection-dominated regime, solutions of convection-diffusion problems usually possesses layers, which are regions where the solution has a steep gradient. It is well known that many classical numerical discretization techniques face difficulties when approximating the solution to these problems. In recent years, physics-informed neural networks (PINNs) for approximating the solution to (initial-)boundary value problems received a lot of interest. In this work, we study various loss functionals for PINNs that are novel in the context of PINNs and are especially designed for convection-dominated convection-diffusion problems. They are numerically compared to the vanilla and a $hp$-variational loss functional from the literature based on two benchmark problems whose solutions possess different types of layers. We observe that the best novel loss functionals reduce the $L^2(Omega)$ error by $17.3%$ for the first and $5.5%$ for the second problem compared to the methods from the literature.

  • M. Demir, V. John, Pressure-robust approximation of the incompressible Navier--Stokes equations in a rotating frame of reference, Preprint no. 3058, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3058 .
    Abstract, PDF (337 kByte)
    A pressure-robust space discretization of the incompressible Navier--Stokes equations in a rotating frame of reference is considered. The discretization employs divergence-free, H1-conforming mixed finite element methods like Scott--Vogelius pairs. An error estimate for the velocity is derived that tracks the dependency of the error bound on the coefficients of the problem, in particular on the angular velocity. Numerical examples illustrate the theoretical results.

  • V. John, Ch. Merdon, M. Zainelabdeen, Augmenting the grad-div stabilization for Taylor--Hood finite elements with a vorticity stabilization, Preprint no. 3055, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3055 .
    Abstract, PDF (2139 kByte)
    The least squares vorticity stabilization (LSVS), proposed in Ahmed et al. for the Scott--Vogelius finite element discretization of the Oseen equations, is studied as an augmentation of the popular grad-div stabilized Taylor--Hood pair of spaces. An error analysis is presented which exploits the situation that the velocity spaces of Scott--Vogelius and Taylor--Hood are identical. Convection-robust error bounds are derived under the assumption that the Scott--Vogelius discretization is well posed on the considered grid. Numerical studies support the analytic results and they show that the LSVS-grad-div method might lead to notable error reductions compared with the standard grad-div method.

Talks, Poster

  • M. Zainelabdeen, Physics-informed neural networks for convection-dominated convection-diffusion problems, 20th Annual Workshop on Numerical Methods for Problems with Layer Phenomena, May 22 - 25, 2024, Department of Mathematics and Statistics, University of Cyprus, Protaras, May 24, 2024.

  • M. Zainelabdeen, Physics-informed neural networks for convection-dominated convection-diffusion problems, International Conference on Boundary and Interior Layers, BAIL 2024, June 10 - 14, 2024, University of A Coruña, Department of Mathematic, Spain, June 11, 2024.

  • C. Cárcamo, Total pressure-based frequency-domain formulation and convergence analysis of Biot's poroelasticity equations with a new finite element stabilization, Minisymposium "Full and reduced-order modeling of multiphysics problems", WONAPDE 2024: Seventh Chilean Workshop on Numerical Analysis of Partial Differential Equations, Concepción, January 15 - 19, 2024, Universidad de Concepción, Barrio Universitario s/n, Region of Bío-Bío, Chile, January 16, 2024.

  • V. John, Finite element methods respecting the discrete maximum principle for convection-diffusion equations, International Conference on 'Latest Advances in Computational and Applied Mathematics' (LACAM) 2024, February 21 - 24, 2024, Indian Institute of Science Education and Research Thiruvananthapuram, Kerala, India, February 21, 2024.

  • V. John, Numerical methods for convection-dominated problems, ALGORITMY 2024, Central-European Conference on Scientific Computing, March 16 - 20, 2024, Department of Mathematics and Descriptive Geometry, Slovak University of Technology in Bratislava, High Tatra Mountains, Slovakia, March 16, 2024.

  • V. John, On two modeling issues in aortic blood flow simulations, Seminar of Dr. Nagaiah Chamakuri, Scientific Computing Group (SCG), School of Mathematics, Indian Institute of Science Education and Research, Thiruvananthapuram, Kerala, India, February 20, 2024.

  • V. John, On using machine learning techniques for the numerical solution of convection-diffusion problems, Seminar-talk, Prof. Sashikumaar Ganesan, Indian Institute of Science Bangalore, Department of Computational and Data Sciences, Bangalore, India, February 13, 2024.

  • S. Katz, Impact of turbulence modeling on the full and reduced simulations of aortic blood flow, 22nd Computational Fluids Conference (CFC 2023), April 25 - 28, 2023, International Association for Computational Mechanics (IACM), Cannes, France, April 28, 2023.

  • L. Ermoneit, B. Schmidt, J. Fuhrmann, Th. Koprucki, M. Kantner, Coherent spin-qubit shuttling in a SiGe quantum bus: Device-scale modeling, simulation and optimal control, Leibniz MMS Days 2023, Potsdam, April 17 - 19, 2023.

  • Y. Hadjimichael, An energy-based finite-strain constitutive model for bent heterostructured materials, GAMM 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics, May 30 - June 2, 2023, Technische Universität Dresden, June 2, 2023.

  • D. Runge, Reduced basis approach for convection-diffusion equations with non-linear boundary reaction conditions, Finite Volumes for Complex Applications 10 (FVCA10), Université de Strasbourg, France, November 2, 2023.

  • D. Runge, Mass-conservative reduced basis approach for convection-diffusion equations with non-linear boundary reaction conditions, Leibniz MMS Days 2023, Leibniz Network "Mathematical Modeling and Simulation", Leibniz Institute for Agricultural Engineering and Bioeconomy Potsdam (ATB), Potsdam, April 18, 2023.

  • M. Demir, Subgrid artificial viscosity modelling based defect-deferred correction method for fluid-fluid interaction, 2023 International CMMSE conference and the Second conference on high performance computing (CHPC), July 3 - 8, 2023, Universidad de Cádiz, Spain, July 6, 2023.

  • M. Demir, Subgrid artificial viscosity modelling based defect-deferred correction method for fluid-fluid interaction, European Conference on Numerical Mathematics and Advanced Applications (ENUMATH), September 4 - 8, 2023, Instituto Superior Técnico, Lisboa, Portugal, September 7, 2023.

  • A. Caiazzo, Multiscale and reduced-order modeling for poroelasticity, European Conference on Numerical Mathematics and Advanced Applications (ENUMATH), September 4 - 8, 2023, Instituto Superior Técnico, Lisboa, Portugal.

  • C. Cárcamo Sanchez, F. Galarce Marín, A. Caiazzo, I. Sack, K. Tabelow, Quantitative tissue pressure imaging via PDE-informed assimilation of MR-data, MATH+ Day, Humboldt-Universität zu Berlin, October 20, 2023.

  • P. Farrell, Charge transport in Perovskites devices: modeling, numerical analysis and simulations, Workshop on Applied Mathematics: Quantum and Classical Models, Università degli Studi di Firenze, Dipartimento di Matematica e Informatica 'Ulisse Dini', Italy, November 29, 2023.

  • P. Farrell, Modeling and numerical simulation of two-dimensional TMDC memristive devices, 10th International Congress on Industrial and Applied Mathematics (ICIAM 2023), Tokyo, Japan, August 20 - 25, 2023.

  • P. Farrell, Modeling and numerical simulation of two-dimensional memristive devices, 22nd European Consortium for Mathematics in Industry (ECMI) Conference on Industrial and Applied Mathematics, June 26 - 30, 2023, Wrocław University of Science and Technology, Poland.

  • P. Farrell, Device physics characterization and interpretation in perovskite and organic materials (DEPERO), October 3 - 5, 2023, Eidgenössische Technische Hochshcule Zürich, nanoGe, Switzerland.

  • D. Frerichs-Mihov, On deep learning techniques for solving convection-dominated convection-diffusion equations, 10th International Congress on Industrial and Applied Mathematics (ICIAM), Minisymposium 00747 ''Analysis and Numerics on Deep Learning Based Methods for Solving PDEs'', August 20 - 25, 2023, Waseda University, Tokyo, Japan, August 23, 2023.

  • D. Frerichs-Mihov, Using deep learning techniques for solving convection-dominated convection-diffusion equations, 22nd Computational Fluids Conference (CFC 2023), April 25 - 28, 2023, International Association for Computational Mechanics (IACM), Cannes, France, April 28, 2023.

  • J. Fuhrmann, Ch. Keller, M. Landstorfer, B. Wagner, Development of an ion-channel model-framework for in-vitro assisted interpretation of current voltage relations, MATH+ Day, Humboldt-Universität zu Berlin, October 20, 2023.

  • J. Fuhrmann, Thermodynamically consistent finite volume schemes for electrolyte simulations, 10th International Congress on Industrial and Applied Mathematics (ICIAM 2023), August 20 - 25, 2023, Waseda University, Tokyo, Japan, August 22, 2023.

  • J. Fuhrmann, Two entropic finite volume schemes for a Nernst--Planck--Poisson system with ion volume constraints, Finite Volumes for Complex Applications 10 (FVCA10), Université de Strasbourg, France, November 1, 2023.

  • J. Fuhrmann, VORONOIFVM.JL -- A multiphysics finite volume solver for elliptic and parabolic systems, SIAM Conference on Computational Science and Engineering (CSE23), Minisymposium MS67 ``Research Software Engineering with Julia'', February 26 - March 3, 2023, Amsterdam, Netherlands, February 27, 2023.

  • J. Fuhrmann, Voronoi finite volume methods for complex applications in Julia, International Conference on Numerical Analysis of Partial Differential Equations (ANEDP 2023), October 16 - 18, 2023, Moulay Ismail University, Faculty of Sciences, Meknes, Morocco.

  • V. John, A SUPG-stabilized POD-ROM method for convection-diffusion-reaction problems (online talk), Numerical Analysis of Galerkin ROMs seminar series (Online Event), February 28, 2023.

  • V. John, On slope limiters in discontinuous Galerkin discretizations of convection-diffusion problems, European Conference on Numerical Mathematics and Advanced Applications (ENUMATH), Minisymposium MS06 ' 'Theoretical and computational aspects of the discontinuous Galerkin method' ', September 4 - 8, 2023, Instituto Superior Técnico, Lisboa, Portugal, September 5, 2023.

  • V. John, Finite element methods respecting the discrete maximum principle for convection-diffusion equations I, 19th Workshop on Numerical Methods for Problems with Layer Phenomena, Charles University, Faculty of Mathematics and Physics, Department of Numerical Mathematics, Prague, Czech Republic, May 26, 2023.

  • V. John, On recent topics in the finite element analysis of convection-diffusion problems (online talk), Numerical Analysis Seminar (Hybrid Event), University of Waterloo, Applied Mathematics, Canada, April 11, 2023.

  • M. Kantner, L. Ermoneit, B. Schmidt, J. Fuhrmann, A. Sala, L.R. Schreiber, Th. Koprucki, Optimal control of a SiGe-quantum bus for coherent electron shuttling in scalable quantum computing architectures, Silicon Quantum Electronics Workshop 2023, Kyoto, Japan, October 31 - November 2, 2023.

  • CH. Merdon, Gradient-robust hybrid discontinuous Galerkin discretizations for the compressible Stokes equations, Forschungsseminar von Prof. Carsten Carstensen, Humboldt-Universität zu Berlin, Institut fuer Mathematik, October 17, 2023.

  • CH. Merdon, Raviart--Thomas enriched Scott--Vogelius finite element methods for the Navier--Stokes equations (online talk), City University of Hong Kong, Department of Mathematics, Hong Kong, January 18, 2023.

  • CH. Merdon, Raviart--Thomas enriched Scott--Vogelius FEM for the Navier--Stokes equations, Capita Selecta Seminar, SACS - Systems, Analysis and Computational Sciences, Department of Mathematics University of Twente (DAMUT), Enschede, Netherlands, May 10, 2023.

  • CH. Merdon, Raviart--Thomas enriched Scott--Vogelius finite element methods for the Navier--Stokes equations, GAMM 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics, May 30 - June 2, 2023, Technische Universität Dresden, June 2, 2023.

  • CH. Merdon, Raviart--Thomas-enriched Scott--Vogelius finite element methods for the Stokes equations on general meshes, The 29th Biennial Numerical Analysis Conference 2023, June 27 - 30, 2023, University of Strathclyde, Department of Mathematics and Statistics, Glasgow, UK, June 27, 2023.

  • O. Pártl, A computational framework for sustainable geothermal energy production in fracture-controlled reservoir based on well placement optimization, Leibniz MMS Days 2023, Potsdam, April 17, 2023.

  • O. Pártl, Finite element methods respecting the discrete maximum principle for convection-diffusion equations III, 19th Workshop on Numerical Methods for Problems with Layer Phenomena, Charles University, Faculty of Mathematics and Physics, Department of Numerical Mathematics, Prague, Czech Republic, May 26, 2023.

  • O. Pártl, Reconstruction of flow domain boundaries from velocity data via multi-step optimization of distributed resistance, European Conference on Numerical Mathematics and Advanced Applications (ENUMATH), September 4 - 8, 2023, Instituto Superior Técnico, Lisboa, Portugal.

  • J.P. Thiele, Competencies and responsibilities of an RSE and how to acquire them (in Germany), FG RSE 2023: Fachgruppentreffen, Gesellschaft für Informatik, October 10 - 11, 2023, Leibniz Universität Hannover, October 10, 2023.

External Preprints

  • R. Araya, A. Caiazzo, F. Chouly, Stokes problem with slip boundary conditions using stabilized finite elements combined with Nitsche, Preprint no. hal-04077986, Hyper Articles en Ligne (HAL), 2023.
    Abstract
    We discuss how slip conditions for the Stokes equation can be handled using Nitsche method, for a stabilized finite element discretization. Emphasis is made on the interplay between stabilization and Nitsche terms. Well-posedness of the discrete problem and optimal convergence rates are established, and illustrated with various numerical experiments.

  • R. Finn, M. O'Donovan, P. Farrell, J. Moatti, T. Streckenbach, Th. Koprucki, S. Schulz, Theoretical study of the impact of alloy disorder on carrier transport and recombination processes in deep UV (Al, Ga)N light emitters, Preprint no. hal-04037215, Hyper Articles en Ligne (HAL), 2023.
    Abstract
    Aluminium gallium nitride ((Al,Ga)N) has gained significant attention in recent years due to its potential for highly efficient light emitters operating in the deep ultra-violet (UV) range (< 280 nm). However, given that current devices exhibit extremely low efficiencies, understanding the fundamental properties of (Al,Ga)N-based systems is of key importance. Here, using a multi-scale simulation framework, we study the impact of alloy disorder on carrier transport, radiative and non-radiative recombination processes in a c-plane Al0.7Ga0.3N/Al0.8Ga0.2N quantum well embedded in a p-i-n junction. Our calculations reveal that alloy fluctuations can open "percolative" pathways that promote transport for the electrons and holes into the quantum well region. Such an effect is neglected in conventional, and widely used transport simulations. Moreover, we find also that the resulting increased carrier density and alloy induced carrier localization effects significantly increase non-radiative Auger-Meitner recombination in comparison to the radiative process. Thus, to avoid such non-radiative process and potentially related material degradation, a careful design (wider well, multi quantum wells) of the active region is required to improve the efficiency of deep UV light emitters.

  • F. Galarce, J. Mura, A. Caiazzo, Bias and multiscale correction methods for variational state estimation algorithms, Preprint no. arXiv:2311.14031, Cornell University, 2023, DOI 10.48550/arXiv.2311.14031 .
    Abstract
    The integration of experimental data into mathematical and computational models is crucial for enhancing their predictive power in real-world scenarios. However, the performance of data assimilation algorithms can be significantly degraded when measurements are corrupted by biased noise, altering the signal magnitude, or when the system dynamics lack smoothness, such as in the presence of fast oscillations or discontinuities. This paper focuses on variational state estimation using the so-called 'Parameterized Background Data Weak' method, which relies on a parameterized background by a set of constraints, enabling state estimation by solving a minimization problem on a reduced-order background model, subject to constraints imposed by the input measurements. To address biased noise in observations, a modified formulation is proposed, incorporating a correction mechanism to handle rapid oscillations by treating them as slow-decaying modes based on a two-scale splitting of the classical reconstruction algorithm. The effectiveness of the proposed algorithms is demonstrated through various examples, including discontinuous signals and simulated Doppler ultrasound data.

  • B. García-Archilla, V. John, S. Katz, J. Novo, POD-ROMs for incompressible flows including snapshots of the temporal derivative of the full order solution: Error bounds for the pressure, Preprint no. arXiv:2304.08313, Cornell University, 2023, DOI 10.48550/arXiv.2304.08313 .
    Abstract
    Reduced order methods (ROMs) for the incompressible Navier--Stokes equations, based on proper orthogonal decomposition (POD), are studied that include snapshots which approach the temporal derivative of the velocity from a full order mixed finite element method (FOM). In addition, the set of snapshots contains the mean velocity of the FOM. Both the FOM and the POD-ROM are equipped with a grad-div stabilization. A velocity error analysis for this method can be found already in the literature. The present paper studies two different procedures to compute approximations to the pressure and proves error bounds for the pressure that are independent of inverse powers of the viscosity. Numerical studies support the analytic results and compare both methods.

  • F. Goth, R. Alves, M. Braun, L.J. Castro, G. Chourdakis, S. Christ, J. Cohen, F. Erxleben, J.-N. Grad, M. Hagdorn, T. Hodges, G. Juckeland, D. Kempf, A.-L. Lamprecht, J. Linxweiler, M. Schwarzmeier, H. Seibold, J.P. Thiele, H. VON Waldow, S. Wittke, Foundational competencies and responsibilities of a research software engineer, Preprint no. arXiv:2311.11457, Cornell University, 2023, DOI 10.48550/arXiv.2311.11457 .
    Abstract
    The term Research Software Engineer, or RSE, emerged a little over 10 years ago as a way to represent individuals working in the research community but focusing on software development. The term has been widely adopted and there are a number of high-level definitions of what an RSE is. However, the roles of RSEs vary depending on the institutional context they work in. At one end of the spectrum, RSE roles may look similar to a traditional research role. At the other extreme, they resemble that of a software engineer in industry. Most RSE roles inhabit the space between these two extremes. Therefore, providing a straightforward, comprehensive definition of what an RSE does and what experience, skills and competencies are required to become one is challenging. In this community paper we define the broad notion of what an RSE is, explore the different types of work they undertake, and define a list of fundamental competencies as well as values that define the general profile of an RSE. On this basis, we elaborate on the progression of these skills along different dimensions, looking at specific types of RSE roles, proposing recommendations for organisations, and giving examples of future specialisations. An appendix details how existing curricula fit into this framework.

  • P.L. Lederer, Ch. Merdon, Gradient-robust hybrid DG discretizations for the compressible Stokes equations, Preprint no. arXiv.2311.06098, Cornell University, 2023, DOI 10.48550/arXiv.2311.06098 .
    Abstract
    This paper studies two hybrid discontinuous Galerkin (HDG) discretizations for the velocity-density formulation of the compressible Stokes equations with respect to several desired structural properties, namely provable convergence, the preservation of non-negativity and mass constraints for the density, and gradient-robustness. The later property dramatically enhances the accuracy in well-balanced situations, such as the hydrostatic balance where the pressure gradient balances the gravity force. One of the studied schemes employs an H(div)-conforming velocity ansatz space which ensures all mentioned properties, while a fully discontinuous method is shown to satisfy all properties but the gradient-robustness. Also higher-order schemes for both variants are presented and compared in three numerical benchmark problems. The final example shows the importance also for non-hydrostatic well-balanced states for the compressible Navier-Stokes equations.

  • B. Spetzler, D. Abdel, F. Schwierz, M. Ziegler, P. Farrell, The role of mobile point defects in two-dimensional memristive devices, Preprint no. arXiv:2304.06527, Cornell University, 2023, DOI 10.48550/arXiv.2304.06527 .
    Abstract
    Two-dimensional (2D) layered transition metal dichalcogenides (TMDCs) are promising memristive materials for neuromorphic computing systems as they could solve the problem of the excessively high energy consumption of conventional von Neumann computer architectures. Despite extensive experimental work, the underlying switching mechanisms are still not understood, impeding progress in material and device functionality. This study reveals the dominant role of mobile defects in the switching dynamics of 2D TMDC materials. The switching process is governed by the formation and annihilation dynamics of a local vacancy depletion zone. Moreover, minor changes in the interface potential barriers cause fundamentally different device behavior previously thought to originate from multiple mechanisms. The key mechanisms are identified with a charge transport model for electrons, holes, and ionic point defects, including image-charge-induced Schottky barrier lowering (SBL). The model is validated by comparing simulations to measurements for various 2D MoS2-based devices, strongly corroborating the relevance of vacancies in TMDC devices and offering a new perspective on the switching mechanisms. The insights gained from this study can be used to extend the functional behavior of 2D TMDC memristive devices in future neuromorphic computing applications.

  • P. Farrell, J. Moatti, M. O'Donovan, S. Schulz, Th. Koprucki, Importance of satisfying thermodynamic consistency in light emitting diode simulations, Preprint no. hal-04012467, Hyper Articles en Ligne (HAL), 2023.
    Abstract
    We show the importance of using a thermodynamically consistent flux discretization when describing drift-diffusion processes within light emitting diode simulations. Using the classical Scharfetter-Gummel scheme with Fermi-Dirac statistics is an example of such an inconsistent scheme. In this case, for an (In,Ga)N multi quantum well device, the Fermi levels show steep gradients on one side of the quantum wells which are not to be expected. This result originates from neglecting diffusion enhancement associated with Fermi-Dirac statistics in the numerical flux approximation. For a thermodynamically consistent scheme, such as the SEDAN scheme, the spikes in the Fermi levels disappear. We will show that thermodynamic inconsistency has far reaching implications on the current-voltage curves and recombination rates.