2D Burgers equation with large Reynolds number using POD/DEIM and calibration

Model order reduction of the two‐dimensional Burgers equation is investigated. The mathematical formulation of POD/discrete empirical interpolation method (DEIM)‐reduced order model (ROM) is derived based on the Galerkin projection and DEIM from the existing high fidelity‐implicit finite‐difference full model. For validation, we numerically compared the POD ROM, POD/DEIM, and the full model in two cases of Re = 100 and Re = 1000, respectively. We found that the POD/DEIM ROM leads to a speed‐up of CPU time by a factor of O(10). The computational stability of POD/DEIM ROM is maintained by means of a careful selection of POD modes and the DEIM interpolation points. The solution of POD/DEIM in the case of Re = 1000 has an accuracy with error O(10^?3) versus O(10^?4) in the case of Re = 100 when compared with the high fidelity model. For this turbulent flow, a closure model consisting of a Tikhonov regularization is carried out in order to recover the missing information and is developed to account for the small‐scale dissipation effect of the truncated POD modes. It is shown that the computational results of this calibrated ROM exhibit considerable agreement with the high fidelity model, which implies the efficiency of the closure model used. Copyright © 2016 John Wiley & Sons, Ltd.     

Symplectic Self-adjointness of Infinite Dimensional Hamiltonian Operators

Symplectic self-adjointness of infinite dimensional Hamiltonian operators is studied, the necessary and sufficient conditions are given. Using the relatively bounded perturbation, the sufficient conditions about symplectic self-adjointness are shown.     

Harnack’s inequality for a class of non-divergent equations in the Heisenberg group

We prove an invariant Harnack’s inequality for operators in non-divergence form structured on Heisenberg vector fields when the coe?cient matrix is uniformly positive definite, continuous, and symplectic. The method consists in constructing appropriate barriers to obtain pointwise-to-measure estimates for supersolutions in small balls, and then invoking the axiomatic approach developed by Di Fazio, Gutiérrez, and Lanconelli to obtain Harnack’s inequality.     

Automorphisms of Hilbert schemes of points on a generic projective K3 surface

We study automorphisms of the Hilbert scheme of n points on a generic projective K3 surface S, for any . We show that is either trivial or generated by a non‐symplectic involution and we determine numerical and divisorial conditions which allow us to distinguish between the two cases. As an application of these results we prove that, for any , there exist infinitely many admissible degrees for the polarization of the K3 surface S such that admits a non‐natural involution. This provides a generalization of the results of [7] for .     

A Flexible Framework for Cubic Regularization Algorithms for Nonconvex Optimization in Function Space

We propose a cubic regularization algorithm that is constructed to deal with nonconvex minimization problems in function space. It allows for a flexible choice of the regularization term and thus accounts for the fact that in such problems one often has to deal with more than one norm. Global and local convergence results are established in a general framework.     

A coefficient identification problem for a mathematical model related to ductal carcinoma in situ

We consider a coefficient identification problem for a mathematical model with free boundary related to ductal carcinoma in situ (DCIS). This inverse problem aims to determine the nutrient consumption rate from additional measurement data at a boundary point. We first obtain a global‐in‐time uniqueness of our inverse problem. Then based on the optimization method, we present a regularization algorithm to recover the nutrient consumption rate. Finally, our numerical experiment shows the effectiveness of the proposed numerical method.     

Regularized inversion of the Laplace transform for series of experiments

Not only in low-field nuclear magnetic resonance, Laplace inversion is a relevant and challenging topic. Considerable conceptual and technical progress has been made, especially for the inversion of data encoding two decay dimensions. Distortion of spectra by overfitting of even moderate noise is counteracted requiring a priori smooth spectra. In this contribution, we treat the case of simple and fast one-dimensional decay experiments that are repeated many times in a series in order to study the evolution of a sample or process. Incorporating the a priori knowledge that also in the series dimension evolution should be smooth, peak position can be stabilized and resolution improved in the decay dimension. It is explained how the standard one-dimensional regularized Laplace inversion can be extended quite simply in order to include regularization in the series dimension. Obvious improvements compared with series of one-dimensional inversions are presented for simulated as well as experimental data. For the latter, comparison with multiexponential fitting is performed.     

Contingent derivatives and regularization for noncoercive inverse problems

Abstract

We study the inverse problem of parameter identification in noncoercive variational problems that commonly appear in applied models. We examine the differentiability of the set-valued parameter-to-solution map using the first-order and the second-order contingent derivatives. We explore the inverse problem using the output least-squares and the modified output least-squares objectives. By regularizing the noncoercive variational problem, we obtain a single-valued regularized parameter-to-solution map and investigate its smoothness and boundedness. We also consider optimization problems using the output least-squares and the modified output least-squares objectives for the regularized variational problem. We give a complete convergence analysis showing that for the output least-squares and the modified output least-squares, the regularized minimization problems approximate the original optimization problems suitably. We also provide the first-order and the second-order adjoint method for the computation of the first-order and the second-order derivatives of the output least-squares objective. We provide discrete formulas for the gradient and the Hessian calculation and present numerical results.