08 Fakultät Mathematik und Physik
Permanent URI for this collectionhttps://elib.uni-stuttgart.de/handle/11682/9
Browse
32 results
Search Results
Item Open Access On a stochastic Camassa-Holm type equation with higher order nonlinearities(2020) Rohde, Christian; Tang, HaoThe subject of this paper is a generalized Camassa-Holm equation under random perturbation. We first establish local existence and uniqueness results as well as blow-up criteria for pathwise solutions in the Sobolev spaces Hs with s>3/2. Then we analyze how noise affects the dependence of solutions on initial data. Even though the noise has some already known regularization effects, much less is known concerning the dependence on initial data. As a new concept we introduce the notion of stability of exiting times and construct an example showing that multiplicative noise (in Itô sense) cannot improve the stability of the exiting time, and simultaneously improve the continuity of the dependence on initial data. Finally, we obtain global existence theorems and estimate associated probabilities.Item Open Access General interface problems. 1(1994) Nicaise, Serge; Sändig, Anna-MargareteWe study transmission problems for elliptic operators of order 2m with general boundary and interface conditions, introducing new covering conditions. This allows to prove solvability, regularity and asymptotics of solutions in weighted Sobolev spaces. We give some numerical examples for the location of the singular exponents.Item Open Access Analysis of target data-dependent greedy kernel algorithms : convergence rates for f-, f· P- and f/P-greedy(2022) Wenzel, Tizian; Santin, Gabriele; Haasdonk, BernardData-dependent greedy algorithms in kernel spaces are known to provide fast converging interpolants, while being extremely easy to implement and efficient to run. Despite this experimental evidence, no detailed theory has yet been presented. This situation is unsatisfactory, especially when compared to the case of the data-independent P-greedy algorithm, for which optimal convergence rates are available, despite its performances being usually inferior to the ones of target data-dependent algorithms. In this work, we fill this gap by first defining a new scale of greedy algorithms for interpolation that comprises all the existing ones in a unique analysis, where the degree of dependency of the selection criterion on the functional data is quantified by a real parameter. We then prove new convergence rates where this degree is taken into account, and we show that, possibly up to a logarithmic factor, target data-dependent selection strategies provide faster convergence. In particular, for the first time we obtain convergence rates for target data adaptive interpolation that are faster than the ones given by uniform points, without the need of any special assumption on the target function. These results are made possible by refining an earlier analysis of greedy algorithms in general Hilbert spaces. The rates are confirmed by a number of numerical examples.Item Open Access Resilience and fault tolerance in high-performance computing for numerical weather and climate prediction(2021) Benacchio, Tommaso; Bonaventura, Luca; Altenbernd, Mirco; Cantwell, Chris D.; Düben, Peter D.; Gillard, Mike; Giraud, Luc; Göddeke, Dominik; Raffin, Erwan; Teranishi, Keita; Wedi, NilsProgress in numerical weather and climate prediction accuracy greatly depends on the growth of the available computing power. As the number of cores in top computing facilities pushes into the millions, increased average frequency of hardware and software failures forces users to review their algorithms and systems in order to protect simulations from breakdown. This report surveys hardware, application-level and algorithm-level resilience approaches of particular relevance to time-critical numerical weather and climate prediction systems. A selection of applicable existing strategies is analysed, featuring interpolation-restart and compressed checkpointing for the numerical schemes, in-memory checkpointing, user-level failure mitigation and backup-based methods for the systems. Numerical examples showcase the performance of the techniques in addressing faults, with particular emphasis on iterative solvers for linear systems, a staple of atmospheric fluid flow solvers. The potential impact of these strategies is discussed in relation to current development of numerical weather prediction algorithms and systems towards the exascale. Trade-offs between performance, efficiency and effectiveness of resiliency strategies are analysed and some recommendations outlined for future developments.Item Open Access Permeability estimation of regular porous structures : a benchmark for comparison of methods(2021) Wagner, Arndt; Eggenweiler, Elissa; Weinhardt, Felix; Trivedi, Zubin; Krach, David; Lohrmann, Christoph; Jain, Kartik; Karadimitriou, Nikolaos; Bringedal, Carina; Voland, Paul; Holm, Christian; Class, Holger; Steeb, Holger; Rybak, IrynaThe intrinsic permeability is a crucial parameter to characterise and quantify fluid flow through porous media. However, this parameter is typically uncertain, even if the geometry of the pore structure is available. In this paper, we perform a comparative study of experimental, semi-analytical and numerical methods to calculate the permeability of a regular porous structure. In particular, we use the Kozeny-Carman relation, different homogenisation approaches (3D, 2D, very thin porous media and pseudo 2D/3D), pore-scale simulations (lattice Boltzmann method, Smoothed Particle Hydrodynamics and finite-element method) and pore-scale experiments (microfluidics). A conceptual design of a periodic porous structure with regularly positioned solid cylinders is set up as a benchmark problem and treated with all considered methods. The results are discussed with regard to the individual strengths and limitations of the used methods. The applicable homogenisation approaches as well as all considered pore-scale models prove their ability to predict the permeability of the benchmark problem. The underestimation obtained by the microfluidic experiments is analysed in detail using the lattice Boltzmann method, which makes it possible to quantify the influence of experimental setup restrictions.Item Open Access Coupled simulations and parameter inversion for neural system and electrophysiological muscle models(2024) Homs‐Pons, Carme; Lautenschlager, Robin; Schmid, Laura; Ernst, Jennifer; Göddeke, Dominik; Röhrle, Oliver; Schulte, MiriamThe functioning of the neuromuscular system is an important factor for quality of life. With the aim of restoring neuromuscular function after limb amputation, novel clinical techniques such as the agonist‐antagonist myoneural interface (AMI) are being developed. In this technique, the residual muscles of an agonist‐antagonist pair are (re‐)connected via a tendon in order to restore their mechanical and neural interaction. Due to the complexity of the system, the AMI can substantially profit from in silico analysis, in particular to determine the prestretch of the residual muscles that is applied during the procedure and determines the range of motion of the residual muscle pair. We present our computational approach to facilitate this. We extend a detailed multi‐X model for single muscles to the AMI setup, that is, a two‐muscle‐one‐tendon system. The model considers subcellular processes as well as 3D muscle and tendon mechanics and is prepared for neural process simulation. It is solved on high performance computing systems. We present simulation results that show (i) the performance of our numerical coupling between muscles and tendon and (ii) a qualitatively correct dependence of the range of motion of muscles on their prestretch. Simultaneously, we pursue a Bayesian parameter inference approach to invert for parameters of interest. Our approach is independent of the underlying muscle model and represents a first step toward parameter optimization, for instance, finding the prestretch, to be applied during surgery, that maximizes the resulting range of motion. Since our multi‐X fine‐grained model is computationally expensive, we present inversion results for reduced Hill‐type models. Our numerical results for cases with known ground truth show the convergence and robustness of our approach.Item Open Access On some distributional properties of subordinated Gaussian random fields(2022) Merkle, Robin; Barth, AndreaMotivated by the subordinated Brownian motion, we define a new class of (in general discontinuous) random fields on higher-dimensional parameter domains: the subordinated Gaussian random field. We investigate the pointwise marginal distribution of the constructed random fields, derive a Lévy-Khinchin-type formula and semi-explicit formulas for the covariance function. Further, we study the pointwise stochastic regularity and present various numerical examples.Item Open Access On the construction of the Stokes flow in a domain with cylindrical ends(2024) Wendland, Wolfgang L.Based on existence results for the Stokes operator and its solution properties in manifolds with cylindrical ends by Große et al. and Kohr et al., the Stokes flow in a three-dimensional compact domain Ω+ with circular openings Σ𝑗 ( 𝑗 = 1, 2) through which the fluid enters and leaves Ω+ through unbounded cylindrical pipes the Stokes flow is modeled as a mixed boundary value problemΩ+ whereas in the cylindrical ends, the velocities and pressures are constant on every straight line in the cylindrical directions with initial values from the openings Σ𝑗 of Ω+. These values equal the velocities and pressures which are obtained from the mixed boundary values' solution in Ω+ at the openings Σ𝑗.Item Open Access Molecular mechanics of disordered solids(2023) Bamer, Franz; Ebrahem, Firaz; Markert, Bernd; Stamm, BenjaminDisordered solids are ubiquitous in engineering and everyday use. Although research has made considerable progress in the last decades, our understanding of the mechanics of these materials is, at best, in an embryonic state. Since the nature of disorder complicates the realization of physically meaningful continuum-mechanical models, particle-based molecular descriptions provide a powerful alternative. This paper reviews the numerical realization of classical molecular dynamics from an engineer’s perspective, starting with selecting potential functions, boundary conditions, time integration, and thermodynamic ensembles. Then, we discuss the concept of the potential energy landscape and the computational realization of the most suitable minimization methods. Subsequently, we discuss the algorithms necessary to numerically generate disordered materials, considering their thermodynamic properties and structural identification. We comprehensively and critically review computational methods and strategies available to mimic disordered materials on a molecular level and discuss some intriguing phenomena that are, to date, mostly ignored when applying models based on continuum-mechanical frameworks. We present the crucial difference between the shear response of a crystalline and a disordered structure. In this context, we elaborate on why it is beneficial to use an overdamped, athermal description to disentangle the complex deformation mechanics of disordered solids and comprehensively discuss the theory of the mechanics of disordered materials, including the problems of prediction and reversibility. Furthermore, we examine the fracture process on the nanoscale and investigate the response behavior to more complex deformation protocols. Finally, we provide critical conclusions, including challenges and future perspectives for engineers.Item Open Access Validation and calibration of coupled porous-medium and free-flow problems using pore-scale resolved models(2020) Rybak, Iryna; Schwarzmeier, Christoph; Eggenweiler, Elissa; Rüde, UlrichThe correct choice of interface conditions and effective parameters for coupled macroscale free-flow and porous-medium models is crucial for a complete mathematical description of the problem under consideration and for accurate numerical simulation of applications. We consider single-fluid-phase systems described by the Stokes–Darcy model. Different sets of coupling conditions for this model are available. However, the choice of these conditions and effective model parameters is often arbitrary. We use large-scale lattice Boltzmann simulations to validate coupling conditions by comparison of the macroscale simulations against pore-scale resolved models. We analyse three settings (lid-driven cavity over a porous bed, infiltration problem and general filtration problem) with different geometrical configurations (channelised and staggered distributions of solid grains) and different sets of interface conditions. Effective parameters for the macroscale models (permeability tensor, boundary layer constants) are computed numerically for each geometrical configuration. Numerical simulation results demonstrate the sensitivity of the coupled Stokes–Darcy problem to the location of the sharp fluid–porous interface, the effective model parameters and the interface conditions.