亚洲男人的天堂2018av,欧美草比,久久久久久免费视频精选,国色天香在线看免费,久久久久亚洲av成人片仓井空

In this contribution we derive and analyze a new numerical method for kinetic equations based on a variable transformation of the moment approximation. Classical minimum-entropy moment closures are a class of reduced models for kinetic equations that conserve many of the fundamental physical properties of solutions. However, their practical use is limited by their high computational cost, as an optimization problem has to be solved for every cell in the space-time grid. In addition, implementation of numerical solvers for these models is hampered by the fact that the optimization problems are only well-defined if the moment vectors stay within the realizable set. For the same reason, further reducing these models by, e.g., reduced-basis methods is not a simple task. Our new method overcomes these disadvantages of classical approaches. The transformation is performed on the semi-discretized level which makes them applicable to a wide range of kinetic schemes and replaces the nonlinear optimization problems by inversion of the positive-definite Hessian matrix. As a result, the new scheme gets rid of the realizability-related problems. Moreover, a discrete entropy law can be enforced by modifying the time stepping scheme. Our numerical experiments demonstrate that our new method is often several times faster than the standard optimization-based scheme.

In this paper, we consider the problem of jointly performing online parameter estimation and optimal sensor placement for a partially observed infinite dimensional linear diffusion process. We present a novel solution to this problem in the form of a continuous-time, two-timescale stochastic gradient descent algorithm, which recursively seeks to maximise the log-likelihood with respect to the unknown model parameters, and to minimise the expected mean squared error of the hidden state estimate with respect to the sensor locations. We also provide extensive numerical results illustrating the performance of the proposed approach in the case that the hidden signal is governed by the two-dimensional stochastic advection-diffusion equation.

We study approximation theorems for the Euler characteristic of the Vietoris-Rips and Cech filtration. The filtration is obtained from a Poisson or binomial sampling scheme in the critical regime. We apply our results to the smooth bootstrap of the Euler characteristic and determine its rate of convergence in the Kantorovich-Wasserstein distance and in the Kolmogorov distance.

We present a new numerical approach for wave induced dynamic fracture. The method is based on a discontinuous Galerkin approximation of the first-order hyperbolic system for elastic waves and a phase-field approximation of brittle fracture driven by the maximum tension. The algorithm is staggered in time and combines an implicit midpoint rule for the wave propagation followed by an implicit Euler step for the phase-field evolution. At fracture, the material is degraded, and the waves are reflected at the diffusive interfaces. Two and three-dimensional examples demonstrate the advantages of the proposed method for the computation of crack growth and spalling initiated by reflected and superposed waves.

The subject of this work is a new stochastic Galerkin method for second-order elliptic partial differential equations with random diffusion coefficients. It combines operator compression in the stochastic variables with tree-based spline wavelet approximation in the spatial variables. Relying on a multilevel expansion of the given random diffusion coefficient, the method is shown to achieve optimal computational complexity up to a logarithmic factor. In contrast to existing results, this holds in particular when the achievable convergence rate is limited by the regularity of the random field, rather than by the spatial approximation order. The convergence and complexity estimates are illustrated by numerical experiments.

In this paper, we propose two novel inertial-like algorithms for solving the split common null point problem (SCNPP) with respect to set-valued maximal operators. The features of the presented algorithm are using new inertial structure (i.e, the design of the new inertial-like method does neither involve computation of the norm of the difference between $x_n$ and $x_{n-1}$ in advance, nor need to consider the special value of the inertial parameter $\theta_n$ to make the condition $\sum_{n=1}^\infty \alpha_n\|x_n-x_{n-1}\|^2<\infty$ valid) and the selection of the step-sizes does not need prior knowledge of operator norms. Numerical experiments are presented to illustrate the performance of the algorithms.

The paper addresses the problem of sampling discretization of integral norms of elements of finite-dimensional subspaces satisfying some conditions. We prove sampling discretization results under two standard kinds of assumptions -- conditions on the entropy numbers and conditions in terms of the Nikol'skii-type inequalities. We prove some upper bounds on the number of sample points sufficient for good discretization and show that these upper bounds are sharp in a certain sense. Then we apply our general conditional results to subspaces with special structures, namely, subspaces with the tensor product structure. We demonstrate that applications of results based on the Nikol'skii-type inequalities provide somewhat better results than applications of results based on the entropy numbers conditions. Finally, we apply discretization results to the problem of sampling recovery.

A convergence theory for the $hp$-FEM applied to a variety of constant-coefficient Helmholtz problems was pioneered in the papers [Melenk-Sauter, 2010], [Melenk-Sauter, 2011], [Esterhazy-Melenk, 2012], [Melenk-Parsania-Sauter, 2013]. This theory shows that, if the solution operator is bounded polynomially in the wavenumber $k$, then the Galerkin method is quasioptimal provided that $hk/p \leq C_1$ and $p\geq C_2 \log k$, where $C_1$ is sufficiently small, $C_2$ is sufficiently large, and both are independent of $k,h,$ and $p$. The significance of this result is that if $hk/p= C_1$ and $p=C_2\log k$, then quasioptimality is achieved with the total number of degrees of freedom proportional to $k^d$; i.e., the $hp$-FEM does not suffer from the pollution effect. This paper proves the analogous quasioptimality result for the heterogeneous (i.e. variable-coefficient) Helmholtz equation, posed in $\mathbb{R}^d$, $d=2,3$, with the Sommerfeld radiation condition at infinity, and $C^\infty$ coefficients. We also prove a bound on the relative error of the Galerkin solution in the particular case of the plane-wave scattering problem. These are the first ever results on the wavenumber-explicit convergence of the $hp$-FEM for the Helmholtz equation with variable coefficients.

Capturing solution near the singular point of any nonlinear SBVPs is challenging because coefficients involved in the differential equation blow up near singularities. In this article, we aim to construct a general method based on orthogonal polynomials as wavelets. We discuss multiresolution analysis for wavelets generated by orthogonal polynomials, e.g., Hermite, Legendre, Chebyshev, Laguerre, and Gegenbauer. Then we use these wavelets for solving nonlinear SBVPs. These wavelets can deal with singularities easily and efficiently. To deal with the nonlinearity, we use both Newton's quasilinearization and the Newton-Raphson method. To show the importance and accuracy of the proposed methods, we solve the Lane-Emden type of problems and compare the computed solutions with the known solutions. As the resolution is increased the computed solutions converge to exact solutions or known solutions. We observe that the proposed technique performs well on a class of Lane-Emden type BVPs. As the paper deals with singularity, non-linearity significantly and different wavelets are used to compare the results.

Trefftz methods are high-order Galerkin schemes in which all discrete functions are elementwise solution of the PDE to be approximated. They are viable only when the PDE is linear and its coefficients are piecewise constant. We introduce a 'quasi-Trefftz' discontinuous Galerkin method for the discretisation of the acoustic wave equation with piecewise-smooth wavespeed: the discrete functions are elementwise approximate PDE solutions. We show that the new discretisation enjoys the same excellent approximation properties as the classical Trefftz one, and prove stability and high-order convergence of the DG scheme. We introduce polynomial basis functions for the new discrete spaces and describe a simple algorithm to compute them. The technique we propose is inspired by the generalised plane waves previously developed for time-harmonic problems with variable coefficients; it turns out that in the case of the time-domain wave equation under consideration the quasi-Trefftz approach allows for polynomial basis functions.