More Like this

1. Multidimensional WENO-AO Reconstructions Using a Simplified Smoothness Indicator and Applications to Conservation Laws

   https://doi.org/10.1007/s10915-023-02319-x  

   Huang, Chieh-Sen ; Arbogast, Todd ; Tian, Chenyu ( August 2023 , Journal of Scientific Computing)

   Finite volume, weighted essentially non-oscillatory (WENO) schemes require the computation of a smoothness indicator. This can be expensive, especially in multiple space dimensions. We consider the use of the simple smoothness indicator$$\sigma ^{\textrm{S}}= \frac{1}{N_{\textrm{S}}-1}\sum _{j} ({\bar{u}}_{j} - {\bar{u}}_{m})^2$$${\sigma }^{\text{S}}=\frac{1}{N_{\text{S}}-1}{\sum }_j{({\overline{u}}_j-{\overline{u}}_m)}^2$, where$$N_{\textrm{S}}$$$N_{\text{S}}$is the number of mesh elements in the stencil,$${\bar{u}}_j$$${\overline{u}}_j$is the local function average over mesh elementj, and indexmgives the target element. Reconstructions utilizing standard WENO weighting fail with this smoothness indicator. We develop a modification of WENO-Z weighting that gives a reliable and accurate reconstruction of adaptive order, which we denote as SWENOZ-AO. We prove that it attains the order of accuracy of the large stencil polynomial approximation when the solution is smooth, and drops to the order of the small stencil polynomial approximations when there is a jump discontinuity in the solution. Numerical examples in one and two space dimensions on general meshes verify the approximation properties of the reconstruction. They also show it to be about 10 times faster in two space dimensions than reconstructions using the classic smoothness indicator. The new reconstruction is applied to define finite volume schemes to approximate the solution of hyperbolic conservation laws. Numerical tests show results of the same quality as standard WENO schemes using the classic smoothness indicator, but with an overall speedup in the computation time of about 3.5–5 times in 2D tests. Moreover, the computational efficiency (CPU time versus error) is noticeably improved.

    

   more » « less
2. Asymptotic Simplification of Aggregation-Diffusion Equations Towards the Heat kernel

   https://doi.org/10.1007/s00205-022-01838-5  

   Carrillo, José A. ; Gómez-Castro, David ; Yao, Yao ; Zeng, Chongchun ( January 2023 , Archive for Rational Mechanics and Analysis)

   We give sharp conditions for the large time asymptotic simplification of aggregation-diffusion equations with linear diffusion. As soon as the interaction potential is bounded and its first and second derivatives decay fast enough at infinity, then the linear diffusion overcomes its effect, either attractive or repulsive, for large times independently of the initial data, and solutions behave like the fundamental solution of the heat equation with some rate. The potential$$W(x) \sim \log |x|$$$W\left(x\right)\sim log\left|x\right|$for$$|x| \gg 1$$$\left|x\right|\gg 1$appears as the natural limiting case when the intermediate asymptotics change. In order to obtain such a result, we produce uniform-in-time estimates in a suitable rescaled change of variables for the entropy, the second moment, Sobolev norms and the$$C^\alpha $$$C^{\alpha }$regularity with a novel approach for this family of equations using modulus of continuity techniques.

    

   more » « less
3. Identifying a stochastic clock network with light entrainment for single cells of Neurospora crassa

   https://doi.org/10.1038/s41598-020-72213-1  

   Caranica, C. ; Al-Omari, A. ; Schüttler, H. -B. ; Arnold, J. ( September 2020 , Scientific Reports)

   Stochastic networks for the clock were identified by ensemble methods using genetic algorithms that captured the amplitude and period variation in single cell oscillators ofNeurosporacrassa. The genetic algorithms were at least an order of magnitude faster than ensemble methods using parallel tempering and appeared to provide a globally optimum solution from a random start in the initial guess of model parameters (i.e., rate constants and initial counts of molecules in a cell). The resulting goodness of fit$${x}^{2}$$$x^2$was roughly halved versus solutions produced by ensemble methods using parallel tempering, and the resulting$${x}^{2}$$$x^2$per data point was only$${\chi }^{2}/n$$${\chi }^2/n$= 2,708.05/953 = 2.84. The fitted model ensemble was robust to variation in proxies for “cell size”. The fitted neutral models without cellular communication between single cells isolated by microfluidics provided evidence for onlyoneStochastic Resonance at one common level of stochastic intracellular noise across days from 6 to 36 h of light/dark (L/D) or in a D/D experiment. When the light-driven phase synchronization was strong as measured by the Kuramoto (K), there was degradation in the single cell oscillations away from the stochastic resonance. The rate constants for the stochastic clock network are consistent with those determined on a macroscopic scale of 10⁷cells.

    

   more » « less
4. Global Well-Posedness for $$H^{-1}(\mathbb {R})$$ Perturbations of KdV with Exotic Spatial Asymptotics

   https://doi.org/10.1007/s00220-022-04522-7  

   Laurens, Thierry ( November 2022 , Communications in Mathematical Physics)

   Given a suitable solutionV(t, x) to the Korteweg–de Vries equation on the real line, we prove global well-posedness for initial data$$u(0,x) \in V(0,x) + H^{-1}(\mathbb {R})$$$u(0,x)\in V(0,x)+H^{-1}\left(R\right)$. Our conditions onVdo include regularity but do not impose any assumptions on spatial asymptotics. We show that periodic profiles$$V(0,x)\in H^5(\mathbb {R}/\mathbb {Z})$$$V(0,x)\in H^5(R/Z)$satisfy our hypotheses. In particular, we can treat localized perturbations of the much-studied periodic traveling wave solutions (cnoidal waves) of KdV. In the companion paper Laurens (Nonlinearity. 35(1):343–387, 2022.https://doi.org/10.1088/1361-6544/ac37f5) we show that smooth step-like initial data also satisfy our hypotheses. We employ the method of commuting flows introduced in Killip and Vişan (Ann. Math. (2) 190(1):249–305, 2019.https://doi.org/10.4007/annals.2019.190.1.4) where$$V\equiv 0$$$V\equiv 0$. In that setting, it is known that$$H^{-1}(\mathbb {R})$$$H^{-1}\left(R\right)$is sharp in the class of$$H^s(\mathbb {R})$$$H^s\left(R\right)$spaces.

    

   more » « less
5. Spherical convex hull of random points on a wedge

   https://doi.org/10.1007/s00208-023-02704-9  

   Besau, Florian ; Gusakova, Anna ; Reitzner, Matthias ; Schütt, Carsten ; Thäle, Christoph ; Werner, Elisabeth M. ( August 2023 , Mathematische Annalen)

   Consider two half-spaces$$H_1^+$$$H_1^{+}$and$$H_2^+$$$H_2^{+}$in$${\mathbb {R}}^{d+1}$$$R^{d+1}$whose bounding hyperplanes$$H_1$$$H_1$and$$H_2$$$H_2$are orthogonal and pass through the origin. The intersection$${\mathbb {S}}_{2,+}^d:={\mathbb {S}}^d\cap H_1^+\cap H_2^+$$$S_{2,+}^d:=S^d\cap H_1^{+}\cap H_2^{+}$is a spherical convex subset of thed-dimensional unit sphere$${\mathbb {S}}^d$$$S^d$, which contains a great subsphere of dimension$$d-2$$$d-2$and is called a spherical wedge. Choosenindependent random points uniformly at random on$${\mathbb {S}}_{2,+}^d$$$S_{2,+}^d$and consider the expected facet number of the spherical convex hull of these points. It is shown that, up to terms of lower order, this expectation grows like a constant multiple of$$\log n$$$logn$. A similar behaviour is obtained for the expected facet number of a homogeneous Poisson point process on$${\mathbb {S}}_{2,+}^d$$$S_{2,+}^d$. The result is compared to the corresponding behaviour of classical Euclidean random polytopes and of spherical random polytopes on a half-sphere.

    

   more » « less