
Programs and Algorithms of Numerical Mathematics 21 : Jablonec nad Nisou, June 1924, 2022 : Proceedings of Seminar
Chleboun, J. ; Kůs, Pavel ; Papež, Jan ; Rozložník, Miroslav ; Segeth, Karel ; Šístek, Jakub
These proceedings contain peerreviewed papers that are based on the invited lectures, short communications, and posters presented at the 21st seminar Programs and Algorithms of Numerical Mathematics (PANM) held in Merkur Hotel, Jablonec nad Nisou, Czech Republic, June 1924, 2022.\nThe seminar was organized by the Institute of Mathematics of the Czech Academy of Sciences under the auspices of EUMATHSIN.CZ, Czech Network for Mathematics in Industry, and with the nancial support provided by the RSJ Foundation. It continued the previous seminars on mathematical software and numerical methods held (biennially, with only one exception) in Alšovice, Bratříkov, Janov nad Nisou, Kořenov, Lázně Libverda, Dolní Maxov, Prague, and Hejnice in the period 19832020. The objective of this series of seminars is to provide a forum for presenting and discussing advanced topics in numerical analysis, computer implementation of numerical algorithms, new approaches to mathematical modeling, and single or multiprocessor applications of computational methods.


Interpolation with restrictions  role of the boundary conditions and individual restrictions
Valášek, Jan ; Sváček, P.
The contribution deals with the remeshing procedure between two computational finite element meshes. The remeshing represented by the interpolation of an approximate solution onto a new mesh is needed in many applications like e.g. in aeroacoustics, here we are particularly interested in the numerical flow simulation of a gradual channel collapse connected with a~severe deterioration of the computational mesh quality. Since the classical Lagrangian projection from one mesh to another is a dissipative method not respecting conservation laws, a conservative interpolation method introducing constraints is described. The constraints have form of Lagrange multipliers enforcing conservation of desired flow quantities, like e.g. total fluid mass, flow kinetic energy or flow potential energy. Then the interpolation problem turns into an error minimization problem, such that the resulting quantities of proposed interpolation satisfy these physical properties while staying as close as possible to the results of Lagrangian interpolation in the L2 norm. The proposed interpolation scheme does not impose any restrictions on mesh generation process and it has a relatively low computational cost. The implementation details are discussed and test cases are shown.


Identification of quasiperiodic processes in the vicinity of the resonance
Fischer, Cyril ; Náprstek, Jiří
In nonlinear dynamical systems, strong quasiperiodic beating effects appear due to combination of selfexcited and forced vibration. The presence of symmetric or asymmetric beatings indicates an exchange of energy between individual degrees of freedom of the model or by multiple close dominant frequencies. This effect is illustrated by the case of the van der Pol equation in the vicinity of resonance. The approximate analysis of these nonlinear effects uses the harmonic balance method and the multiple scale method.


Estimation of EDZ zones in great depths by elasticplastic models
Sysala, Stanislav
This contribution is devoted to modeling damage zones caused by the excavation of tunnels and boreholes (EDZ zones) in connection with the issue of deep storage of spent nuclear fuel in crystalline rocks. In particular, elasticplastic models with MohrCoulomb or HoekBrown yield criteria are considered. Selected details of the numerical solution to the corresponding problems are mentioned. Possibilities of elastic and elasticplastic approaches are illustrated by a numerical example.


Semantic segmentation using support vector machine classifier
Pecha, Marek ; Langford, Z. ; Horák, David ; Tran Mills, R.
This paper deals with wildfire identification in the Alaska regions as a semantic segmentation task using support vector machine classifiers. Instead of colour information represented by means of BGR channels, we proceed with a normalized reflectance over 152 days so that such time series is assigned to each pixel. We compare models associated with $\mathcal{l}1$loss and $\mathcal{l}2$loss functions and stopping criteria based on a projected gradient and duality gap in the presented benchmarks.


Determination of the initial stress tensor from deformation of underground openingtheoretical background and applications
Malík, Josef ; Kolcun, Alexej
In this paper a method for the detection of initial stress tensor is proposed. The method is based on measuring distances between some pairs of points located on the wall of underground opening in the excavation process. This methods is based on the solution of eighteen auxiliary problems in the theory of elasticity with force boundary conditions. The optimal location of the pairs of points on the wall of underground work is studied. The pairs must be located so that the condition number of a certain matrix has the minimal value, which guarantees a reliable estimation of initial stress tensor.


Numerical realization of the Bayesian inversion accelerated using surrogate models
Bérešová, Simona
The Bayesian inversion is a natural approach to the solution of inverse problems based on uncertain observed data. The result of such an inverse problem is the posterior distribution of unknown parameters. This paper deals with the numerical realization of the Bayesian inversion focusing on problems governed by computationally expensive forward models such as numerical solutions of partial differential equations. Samples from the posterior distribution are generated using the Markov chain Monte Carlo (MCMC) methods accelerated with surrogate models. A surrogate model is understood as an approximation of the forward model which should be computationally much cheaper. The target distribution is not fully replaced by its approximation. Therefore, samples from the exact posterior distribution are provided. In addition, nonintrusive surrogate models can be updated during the sampling process resulting in an adaptive MCMC method. The use of the surrogate models significantly reduces the number of evaluations of the forward model needed for a reliable description of the posterior distribution. Described sampling procedures are implemented in the form of a Python package.


Reduced basis solver for stochastic Galerkin formulation of Darcy flow with uncertain material parameters
Béreš, Michal
In this contribution, we present a solution to the stochastic Galerkin (SG) matrix equations coming from the Darcy flow problem with uncertain material coefficients in the separable form. The SG system of equations is kept in the compressed tensor form and its solution is a very challenging task. Here, we present the reduced basis (RB) method as a solver which looks for a lowrank representation of the solution. The construction of the RB consists of iterative expanding of the basis using Monte Carlo sampling. We discuss the setting of the sampling procedure and an efficient solution of multiple similar systems emerging during the sampling procedure using deflation. We conclude with a demonstration of the use of SG solution for forward uncertainty quantification.

 

Spherical basis function approximation with particular trend functions
Segeth, Karel
The paper is concerned with the measurement of scalar physical quantities at nodes on the $(d1)$dimensional unit sphere surface in the hbox{$d$dimensional} Euclidean space and the spherical RBF interpolation of the data obtained. In particular, we consider $d=3$. We employ an inverse multiquadric as the radial basis function and the corresponding trend is a polynomial of degree 2 defined in Cartesian coordinates. We prove the existence of the interpolation formula of the type considered. The formula can be useful in the interpretation of many physical measurements. We show an example concerned with the measurement of anisotropy of magnetic susceptibility having extensive applications in geosciences and present numerical difficulties connected with the high condition number of the matrix of the system defining the interpolation.
