Search results for " numerical analysis."
showing 10 items of 103 documents
Stochastic Galerkin method for cloud simulation
2018
AbstractWe develop a stochastic Galerkin method for a coupled Navier-Stokes-cloud system that models dynamics of warm clouds. Our goal is to explicitly describe the evolution of uncertainties that arise due to unknown input data, such as model parameters and initial or boundary conditions. The developed stochastic Galerkin method combines the space-time approximation obtained by a suitable finite volume method with a spectral-type approximation based on the generalized polynomial chaos expansion in the stochastic space. The resulting numerical scheme yields a second-order accurate approximation in both space and time and exponential convergence in the stochastic space. Our numerical results…
Multivariate morphometric analysis of the Stipa turkestanica group (Poaceae: Stipa sect. Stipa)
2015
Based on numerical analyses of macromorphological characters (cluster analysis, principal coordinate analysis and principal component analysis), scanning electron microscopy observation of lemma and lamina micromorphology, as well as field observations, five taxa belonging to the Stipa turkestanica group have been recognized in the mountain area of Central Asia. They are S. turkestanica subsp. turkestanica, S. turkestanica subsp. trichoides, S. macroglossa subsp. macroglossa, S. macroglossa subsp. kazachstanica and S. kirghisorum. As a result of this study, we propose one new combination, S. macroglossa var. pubescens, and designate lectotypes for S. turkestanica subsp. trichoides and S. ma…
A posteriori modelling-discretization error estimate for elliptic problems with L ∞-Coefficients
2017
We consider elliptic problems with complicated, discontinuous diffusion tensor A0. One of the standard approaches to numerically treat such problems is to simplify the coefficient by some approximation, say Aϵ, and to use standard finite elements. In [19] a combined modelling-discretization strategy has been proposed which estimates the discretization and modelling errors by a posteriori estimates of functional type. This strategy allows to balance these two errors in a problem adapted way. However, the estimate of the modelling error was derived under the assumption that the difference A0 - Aϵ becomes small with respect to the L∞-norm. This implies in particular that interfaces/discontinui…
Better numerical approximation by Durrmeyer type operators
2018
The main object of this paper is to construct new Durrmeyer type operators which have better features than the classical one. Some results concerning the rate of convergence and asymptotic formulas of the new operator are given. Finally, the theoretical results are analyzed by numerical examples.
Estimates for the differences of positive linear operators and their derivatives
2019
The present paper deals with the estimate of the differences of certain positive linear operators and their derivatives. Oxur approach involves operators defined on bounded intervals, as Bernstein operators, Kantorovich operators, genuine Bernstein-Durrmeyer operators, and Durrmeyer operators with Jacobi weights. The estimates in quantitative form are given in terms of the first modulus of continuity. In order to analyze the theoretical results in the last section, we consider some numerical examples.
Rough McKean-Vlasov dynamics for robust ensemble Kalman filtering
2021
Motivated by the challenge of incorporating data into misspecified and multiscale dynamical models, we study a McKean-Vlasov equation that contains the data stream as a common driving rough path. This setting allows us to prove well-posedness as well as continuity with respect to the driver in an appropriate rough-path topology. The latter property is key in our subsequent development of a robust data assimilation methodology: We establish propagation of chaos for the associated interacting particle system, which in turn is suggestive of a numerical scheme that can be viewed as an extension of the ensemble Kalman filter to a rough-path framework. Finally, we discuss a data-driven method bas…
Guaranteed error control bounds for the stabilised space-time IgA approximations to parabolic problems
2017
The paper is concerned with space-time IgA approximations of parabolic initial-boundary value problems. We deduce guaranteed and fully computable error bounds adapted to special features of IgA approximations and investigate their applicability. The derivation method is based on the analysis of respective integral identities and purely functional arguments. Therefore, the estimates do not contain mesh-dependent constants and are valid for any approximation from the admissible (energy) class. In particular, they provide computable error bounds for norms associated with stabilised space-time IgA approximations as well as imply efficient error indicators enhancing the performance of fully adap…
Existence of dynamical low-rank approximations to parabolic problems
2021
The existence and uniqueness of weak solutions to dynamical low-rank evolution problems for parabolic partial differential equations in two spatial dimensions is shown, covering also non-diagonal diffusion in the elliptic part. The proof is based on a variational time-stepping scheme on the low-rank manifold. Moreover, this scheme is shown to be closely related to practical methods for computing such low-rank evolutions.
A degenerating convection-diffusion system modelling froth flotation with drainage
2022
Abstract Froth flotation is a common unit operation used in mineral processing. It serves to separate valuable mineral particles from worthless gangue particles in finely ground ores. The valuable mineral particles are hydrophobic and attach to bubbles of air injected into the pulp. This creates bubble-particle aggregates that rise to the top of the flotation column where they accumulate to a froth or foam layer that is removed through a launder for further processing. At the same time, the hydrophilic gangue particles settle and are removed continuously. The drainage of liquid due to capillarity is essential for the formation of a stable froth layer. This effect is included into a previous…
A generalized Newton iteration for computing the solution of the inverse Henderson problem
2020
We develop a generalized Newton scheme IHNC for the construction of effective pair potentials for systems of interacting point-like particles.The construction is made in such a way that the distribution of the particles matches a given radial distribution function. The IHNC iteration uses the hypernetted-chain integral equation for an approximate evaluation of the inverse of the Jacobian of the forward operator. In contrast to the full Newton method realized in the Inverse Monte Carlo (IMC) scheme, the IHNC algorithm requires only a single molecular dynamics computation of the radial distribution function per iteration step, and no further expensive cross-correlations. Numerical experiments…