@article{HijaziFreitagLandwehr2023, author = {Hijazi, Saddam and Freitag, Melina A. and Landwehr, Niels}, title = {POD-Galerkin reduced order models and physics-informed neural networks for solving inverse problems for the Navier-Stokes equations}, series = {Advanced modeling and simulation in engineering sciences : AMSES}, volume = {10}, journal = {Advanced modeling and simulation in engineering sciences : AMSES}, number = {1}, publisher = {SpringerOpen}, address = {Berlin}, issn = {2213-7467}, doi = {10.1186/s40323-023-00242-2}, pages = {38}, year = {2023}, abstract = {We present a Reduced Order Model (ROM) which exploits recent developments in Physics Informed Neural Networks (PINNs) for solving inverse problems for the Navier-Stokes equations (NSE). In the proposed approach, the presence of simulated data for the fluid dynamics fields is assumed. A POD-Galerkin ROM is then constructed by applying POD on the snapshots matrices of the fluid fields and performing a Galerkin projection of the NSE (or the modified equations in case of turbulence modeling) onto the POD reduced basis. A POD-Galerkin PINN ROM is then derived by introducing deep neural networks which approximate the reduced outputs with the input being time and/or parameters of the model. The neural networks incorporate the physical equations (the POD-Galerkin reduced equations) into their structure as part of the loss function. Using this approach, the reduced model is able to approximate unknown parameters such as physical constants or the boundary conditions. A demonstration of the applicability of the proposed ROM is illustrated by three cases which are the steady flow around a backward step, the flow around a circular cylinder and the unsteady turbulent flow around a surface mounted cubic obstacle.}, language = {en} } @article{MicheletBindelliniMelinetal.2023, author = {Michelet, Robin and Bindellini, Davide and Melin, Johanna and Neumann, Uta and Blankenstein, Oliver and Huisinga, Wilhelm and Johnson, Trevor N. and Whitaker, Martin J. and Ross, Richard and Kloft, Charlotte}, title = {Insights in the maturational processes influencing hydrocortisone pharmacokinetics in congenital adrenal hyperplasia patients using a middle-out approach}, series = {Frontiers in Pharmacology}, volume = {13}, journal = {Frontiers in Pharmacology}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {1663-9812}, doi = {10.3389/fphar.2022.1090554}, pages = {14}, year = {2023}, abstract = {Introduction: Hydrocortisone is the standard of care in cortisol replacement therapy for congenital adrenal hyperplasia patients. Challenges in mimicking cortisol circadian rhythm and dosing individualization can be overcome by the support of mathematical modelling. Previously, a non-linear mixed-effects (NLME) model was developed based on clinical hydrocortisone pharmacokinetic (PK) pediatric and adult data. Additionally, a physiologically-based pharmacokinetic (PBPK) model was developed for adults and a pediatric model was obtained using maturation functions for relevant processes. In this work, a middle-out approach was applied. The aim was to investigate whether PBPK-derived maturation functions could provide a better description of hydrocortisone PK inter-individual variability when implemented in the NLME framework, with the goal of providing better individual predictions towards precision dosing at the patient level. Methods: Hydrocortisone PK data from 24 adrenal insufficiency pediatric patients and 30 adult healthy volunteers were used for NLME model development, while the PBPK model and maturation functions of clearance and cortisol binding globulin (CBG) were developed based on previous studies published in the literature. Results: Clearance (CL) estimates from both approaches were similar for children older than 1 year (CL/F increasing from around 150 L/h to 500 L/h), while CBG concentrations differed across the whole age range (CBG(NLME) stable around 0.5 mu M vs. steady increase from 0.35 to 0.8 mu M for CBG (PBPK)). PBPK-derived maturation functions were subsequently included in the NLME model. After inclusion of the maturation functions, none, a part of, or all parameters were re-estimated. However, the inclusion of CL and/or CBG maturation functions in the NLME model did not result in improved model performance for the CL maturation function (\& UDelta;OFV > -15.36) and the re-estimation of parameters using the CBG maturation function most often led to unstable models or individual CL prediction bias. Discussion: Three explanations for the observed discrepancies could be postulated, i) non-considered maturation of processes such as absorption or first-pass effect, ii) lack of patients between 1 and 12 months, iii) lack of correction of PBPK CL maturation functions derived from urinary concentration ratio data for the renal function relative to adults. These should be investigated in the future to determine how NLME and PBPK methods can work towards deriving insights into pediatric hydrocortisone PK.}, language = {en} } @article{GerlachGlueckKunze2023, author = {Gerlach, Moritz and Gl{\"u}ck, Jochen and Kunze, Markus}, title = {Stability of transition semigroups and applications to parabolic equations}, series = {Transactions of the American Mathematical Society}, volume = {376}, journal = {Transactions of the American Mathematical Society}, number = {1}, publisher = {American Mathematical Soc.}, address = {Providence}, issn = {0002-9947}, doi = {10.1090/tran/8620}, pages = {153 -- 180}, year = {2023}, abstract = {This paper deals with the long-term behavior of positive operator semigroups on spaces of bounded functions and of signed measures, which have applications to parabolic equations with unbounded coefficients and to stochas-tic analysis. The main results are a Tauberian type theorem characterizing the convergence to equilibrium of strongly Feller semigroups and a generalization of a classical convergence theorem of Doob. None of these results requires any kind of time regularity of the semigroup.}, language = {en} } @article{FalkenhagenKnoechelKloftetal.2023, author = {Falkenhagen, Undine and Kn{\"o}chel, Jane and Kloft, Charlotte and Huisinga, Wilhelm}, title = {Deriving mechanism-based pharmacodynamic models by reducing quantitative systems pharmacology models}, series = {CPT: Pharmacometrics \& Systems Pharmacology}, volume = {12}, journal = {CPT: Pharmacometrics \& Systems Pharmacology}, number = {4}, publisher = {Wiley}, address = {Hoboken}, issn = {2163-8306}, doi = {10.1002/psp4.12903}, pages = {432 -- 443}, year = {2023}, abstract = {Quantitative systems pharmacology (QSP) models integrate comprehensive qualitative and quantitative knowledge about pharmacologically relevant processes. We previously proposed a first approach to leverage the knowledge in QSP models to derive simpler, mechanism-based pharmacodynamic (PD) models. Their complexity, however, is typically still too large to be used in the population analysis of clinical data. Here, we extend the approach beyond state reduction to also include the simplification of reaction rates, elimination of reactions, and analytic solutions. We additionally ensure that the reduced model maintains a prespecified approximation quality not only for a reference individual but also for a diverse virtual population. We illustrate the extended approach for the warfarin effect on blood coagulation. Using the model-reduction approach, we derive a novel small-scale warfarin/international normalized ratio model and demonstrate its suitability for biomarker identification. Due to the systematic nature of the approach in comparison with empirical model building, the proposed model-reduction algorithm provides an improved rationale to build PD models also from QSP models in other applications.}, language = {en} } @article{StueblerKloftHuisinga2023, author = {St{\"u}bler, Sabine and Kloft, Charlotte and Huisinga, Wilhelm}, title = {Cell-level systems biology model to study inflammatory bowel diseases and their treatment options}, series = {CPT: pharmacometrics \& systems pharmacology}, volume = {12}, journal = {CPT: pharmacometrics \& systems pharmacology}, number = {5}, publisher = {Nature Publ. Group}, address = {London}, issn = {2163-8306}, doi = {10.1002/psp4.12932}, pages = {690 -- 705}, year = {2023}, abstract = {To help understand the complex and therapeutically challenging inflammatory bowel diseases (IBDs), we developed a systems biology model of the intestinal immune system that is able to describe main aspects of IBD and different treatment modalities thereof. The model, including key cell types and processes of the mucosal immune response, compiles a large amount of isolated experimental findings from literature into a larger context and allows for simulations of different inflammation scenarios based on the underlying data and assumptions. In the context of a large and diverse virtual IBD population, we characterized the patients based on their phenotype (in contrast to healthy individuals, they developed persistent inflammation after a trigger event) rather than on a priori assumptions on parameter differences to a healthy individual. This allowed to reproduce the enormous diversity of predispositions known to lead to IBD. Analyzing different treatment effects, the model provides insight into characteristics of individual drug therapy. We illustrate for anti-TNF-alpha therapy, how the model can be used (i) to decide for alternative treatments with best prospects in the case of nonresponse, and (ii) to identify promising combination therapies with other available treatment options.}, language = {en} } @article{SharmaHainzlZoeller2023, author = {Sharma, Shubham and Hainzl, Sebastian and Z{\"o}ller, Gert}, title = {Seismicity parameters dependence on main shock-induced co-seismic stress}, series = {Geophysical journal international}, volume = {235}, journal = {Geophysical journal international}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0956-540X}, doi = {10.1093/gji/ggad201}, pages = {509 -- 517}, year = {2023}, abstract = {The Gutenberg-Richter (GR) and the Omori-Utsu (OU) law describe the earthquakes' energy release and temporal clustering and are thus of great importance for seismic hazard assessment. Motivated by experimental results, which indicate stress-dependent parameters, we consider a combined global data set of 127 main shock-aftershock sequences and perform a systematic study of the relationship between main shock-induced stress changes and associated seismicity patterns. For this purpose, we calculate space-dependent Coulomb Stress (\& UDelta;CFS) and alternative receiver-independent stress metrics in the surrounding of the main shocks. Our results indicate a clear positive correlation between the GR b-value and the induced stress, contrasting expectations from laboratory experiments and suggesting a crucial role of structural heterogeneity and strength variations. Furthermore, we demonstrate that the aftershock productivity increases nonlinearly with stress, while the OU parameters c and p systematically decrease for increasing stress changes. Our partly unexpected findings can have an important impact on future estimations of the aftershock hazard.}, language = {en} } @article{Lewandowski2022, author = {Lewandowski, Max}, title = {Hadamard states for bosonic quantum field theory on globally hyperbolic spacetimes}, series = {Journal of mathematical physics}, volume = {63}, journal = {Journal of mathematical physics}, number = {1}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0022-2488}, doi = {10.1063/5.0055753}, pages = {34}, year = {2022}, abstract = {According to Radzikowski's celebrated results, bisolutions of a wave operator on a globally hyperbolic spacetime are of the Hadamard form iff they are given by a linear combination of distinguished parametrices i2(G˜aF-G˜F+G˜A-G˜R) in the sense of Duistermaat and H{\"o}rmander [Acta Math. 128, 183-269 (1972)] and Radzikowski [Commun. Math. Phys. 179, 529 (1996)]. Inspired by the construction of the corresponding advanced and retarded Green operator GA, GR as done by B{\"a}r, Ginoux, and Pf{\"a}ffle {Wave Equations on Lorentzian Manifolds and Quantization [European Mathematical Society (EMS), Z{\"u}rich, 2007]}, we construct the remaining two Green operators GF, GaF locally in terms of Hadamard series. Afterward, we provide the global construction of i2(G˜aF-G˜F), which relies on new techniques such as a well-posed Cauchy problem for bisolutions and a patching argument using Čech cohomology. This leads to global bisolutions of the Hadamard form, each of which can be chosen to be a Hadamard two-point-function, i.e., the smooth part can be adapted such that, additionally, the symmetry and the positivity condition are exactly satisfied.}, language = {en} } @article{Zoeller2022, author = {Z{\"o}ller, Gert}, title = {A note on the estimation of the maximum possible earthquake magnitude based on extreme value theory for the Groningen Gas Field}, series = {The bulletin of the Seismological Society of America : BSSA}, volume = {112}, journal = {The bulletin of the Seismological Society of America : BSSA}, number = {4}, publisher = {Seismological Society of America}, address = {El Cerito, Calif.}, issn = {0037-1106}, doi = {10.1785/0120210307}, pages = {1825 -- 1831}, year = {2022}, abstract = {Extreme value statistics is a popular and frequently used tool to model the occurrence of large earthquakes. The problem of poor statistics arising from rare events is addressed by taking advantage of the validity of general statistical properties in asymptotic regimes. In this note, I argue that the use of extreme value statistics for the purpose of practically modeling the tail of the frequency-magnitude distribution of earthquakes can produce biased and thus misleading results because it is unknown to what degree the tail of the true distribution is sampled by data. Using synthetic data allows to quantify this bias in detail. The implicit assumption that the true M-max is close to the maximum observed magnitude M-max,M-observed restricts the class of the potential models a priori to those with M-max = M-max,M-observed + Delta M with an increment Delta M approximate to 0.5... 1.2. This corresponds to the simple heuristic method suggested by Wheeler (2009) and labeled :M-max equals M-obs plus an increment." The incomplete consideration of the entire model family for the frequency-magnitude distribution neglects, however, the scenario of a large so far unobserved earthquake.}, language = {en} } @article{KayaFreitag2022, author = {Kaya, Adem and Freitag, Melina A.}, title = {Conditioning analysis for discrete Helmholtz problems}, series = {Computers and mathematics with applications : an international journal}, volume = {118}, journal = {Computers and mathematics with applications : an international journal}, publisher = {Elsevier Science}, address = {Amsterdam}, issn = {0898-1221}, doi = {10.1016/j.camwa.2022.05.016}, pages = {171 -- 182}, year = {2022}, abstract = {In this paper, we examine conditioning of the discretization of the Helmholtz problem. Although the discrete Helmholtz problem has been studied from different perspectives, to the best of our knowledge, there is no conditioning analysis for it. We aim to fill this gap in the literature. We propose a novel method in 1D to observe the near-zero eigenvalues of a symmetric indefinite matrix. Standard classification of ill-conditioning based on the matrix condition number is not true for the discrete Helmholtz problem. We relate the ill-conditioning of the discretization of the Helmholtz problem with the condition number of the matrix. We carry out analytical conditioning analysis in 1D and extend our observations to 2D with numerical observations. We examine several discretizations. We find different regions in which the condition number of the problem shows different characteristics. We also explain the general behavior of the solutions in these regions.}, language = {en} } @article{HoudebertZass2022, author = {Houdebert, Pierre and Zass, Alexander}, title = {An explicit Dobrushin uniqueness region for Gibbs point processes with repulsive interactions}, series = {Journal of applied probability / Applied Probability Trust}, volume = {59}, journal = {Journal of applied probability / Applied Probability Trust}, number = {2}, publisher = {Cambridge Univ. Press}, address = {Cambridge}, issn = {0021-9002}, doi = {10.1017/jpr.2021.70}, pages = {541 -- 555}, year = {2022}, abstract = {We present a uniqueness result for Gibbs point processes with interactions that come from a non-negative pair potential; in particular, we provide an explicit uniqueness region in terms of activity z and inverse temperature beta. The technique used relies on applying to the continuous setting the classical Dobrushin criterion. We also present a comparison to the two other uniqueness methods of cluster expansion and disagreement percolation, which can also be applied for this type of interaction.}, language = {en} }