This study proposes a novel spatial discretization procedure for the compressible Euler equations which guarantees entropy conservation at a discrete level when an arbitrary equation of state is assumed. The proposed method, based on a locally-conservative discretization, guarantees also the spatial conservation of mass, momentum, and total energy and is kinetic energy-preserving. In order to achieve the entropy-conservation property for an arbitrary non-ideal gas, a general strategy is adopted based on the manipulation of discrete balance equations through the imposition of global entropy conservation and the use of a summation by parts rule. The procedure, which is extended to an arbitrary order of accuracy, conducts to a general form of the internal-energy numerical flux which results in a nonlinear function of thermodynamic and dynamic variables and still admits the mass flux as a residual degree of freedom. The effectiveness of the novel entropy-conservative formulation is demonstrated through numerical tests making use of some of the most popular cubic equations of state.
相關內容
The main purpose of this paper is to design a local discontinuous Galerkin (LDG) method for the Benjamin-Ono equation. We analyze the stability and error estimates for the semi-discrete LDG scheme. We prove that the scheme is $L^2$-stable and it converges at a rate $\mathcal{O}(h^{k+1/2})$ for general nonlinear flux. Furthermore, we develop a fully discrete LDG scheme using the four-stage fourth order Runge-Kutta method and ensure the devised scheme is strongly stable in case of linear flux using two-step and three-step stability approach under an appropriate time step constraint. Numerical examples are provided to validate the efficiency and accuracy of the method.
We consider non-linear Bayesian inverse problems of determining the parameter $f$. For the posterior distribution with a class of Gaussian process priors, we study the statistical performance of variational Bayesian inference to the posterior with variational sets consisting of Gaussian measures or a mean-field family. We propose certain conditions on the forward map $\mathcal{G}$, the variational set $\mathcal{Q}$ and the prior such that, as the number $N$ of measurements increases, the resulting variational posterior distributions contract to the ground truth $f_0$ generating the data, and derive a convergence rate with polynomial order or logarithmic order. As specific examples, we consider a collection of non-linear inverse problems, including the Darcy flow problem, the inverse potential problem for a subdiffusion equation, and the inverse medium scattering problem. Besides, we show that our convergence rates are minimax optimal for these inverse problems.
Predicting the response of nonlinear dynamical systems subject to random, broadband excitation is important across a range of scientific disciplines, such as structural dynamics and neuroscience. Building data-driven models requires experimental measurements of the system input and output, but it can be difficult to determine whether inaccuracies in the model stem from modelling errors or noise. This paper presents a novel method to identify the causal component of the input-output data from measurements of a system in the presence of output noise, as a function of frequency, without needing a high fidelity model. An output prediction, calculated using an available model, is optimally combined with noisy measurements of the output to predict the input to the system. The parameters of the algorithm balance the two output signals and are utilised to calculate a nonlinear coherence metric as a measure of causality. This method is applicable to a broad class of nonlinear dynamical systems. There are currently no solutions to this problem in the absence of a complete benchmark model.
To solve the Cahn-Hilliard equation numerically, a new time integration algorithm is proposed, which is based on a combination of the Eyre splitting and the local iteration modified (LIM) scheme. The latter is employed to tackle the implicit system arising each time integration step. The proposed method is gradient-stable and allows to use large time steps, whereas, regarding its computational structure, it is an explicit time integration scheme. Numerical tests are presented to demonstrate abilities of the new method and to compare it with other time integration methods for Cahn-Hilliard equation.
Approximating solutions of ordinary and partial differential equations constitutes a significant challenge. Based on functional expressions that inherently depend on neural networks, neural forms are specifically designed to precisely satisfy the prescribed initial or boundary conditions of the problem, while providing the approximate solutions in closed form. Departing from the important class of ordinary differential equations, the present work aims to refine and validate the neural forms methodology, paving the ground for further developments in more challenging fields. The main contributions are as follows. First, it introduces a formalism for systematically crafting proper neural forms with adaptable boundary matches that are amenable to optimization. Second, it describes a novel technique for converting problems with Neumann or Robin conditions into equivalent problems with parametric Dirichlet conditions. Third, it outlines a method for determining an upper bound on the absolute deviation from the exact solution. The proposed augmented neural forms approach was tested on a set of diverse problems, encompassing first- and second-order ordinary differential equations, as well as first-order systems. Stiff differential equations have been considered as well. The resulting solutions were subjected to assessment against existing exact solutions, solutions derived through the common penalized neural method, and solutions obtained via contemporary numerical analysis methods. The reported results demonstrate that the augmented neural forms not only satisfy the boundary and initial conditions exactly, but also provide closed-form solutions that facilitate high-quality interpolation and controllable overall precision. These attributes are essential for expanding the application field of neural forms to more challenging problems that are described by partial differential equations.
Error estimates of cubic interpolated pseudo-particle scheme (CIP scheme) for the one-dimensional advection equation with periodic boundary conditions are presented. The CIP scheme is a semi-Lagrangian method involving the piecewise cubic Hermite interpolation. Although it is numerically known that the space-time accuracy of the scheme is third order, its rigorous proof remains an open problem. In this paper, denoting the spatial and temporal mesh sizes by $ h $ and $ \Delta t $ respectively, we prove an error estimate $ O(\Delta t^3 + \frac{h^4}{\Delta t}) $ in $ L^2 $ norm theoretically, which justifies the above-mentioned prediction if $ h = O(\Delta t) $. The proof is based on properties of the interpolation operator; the most important one is that it behaves as the $ L^2 $ projection for the second-order derivatives. We remark that the same strategy perfectly works as well to address an error estimate for the semi-Lagrangian method with the cubic spline interpolation.
We study the decidability and complexity of non-cooperative rational synthesis problem (abbreviated as NCRSP) for some classes of probabilistic strategies. We show that NCRSP for stationary strategies and Muller objectives is in 3-EXPTIME, and if we restrict the strategies of environment players to be positional, NCRSP becomes NEXPSPACE solvable. On the other hand, NCRSP_>, which is a variant of NCRSP, is shown to be undecidable even for pure finite-state strategies and terminal reachability objectives. Finally, we show that NCRSP becomes EXPTIME solvable if we restrict the memory of a strategy to be the most recently visited t vertices where t is linear in the size of the game.
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The $L2$-$1_\sigma$ scheme is used to discretize the time-fractional operator, whereas a modified version of the composite trapezoidal approximation is employed to discretize the Volterra operator. Subsequently, it helps to convert the proposed model into a second-order BVP in a semi-discrete form. The multi-dimensional Haar wavelets are then used for grid adaptation and efficient computations for the 2D problem, whereas the standard second-order approximations are employed to approximate the spatial derivatives for the 1D case. The stability analysis is carried out on an adaptive mesh in time. The convergence analysis leads to $O(N^{-2}+M^{-2})$ accurate solution in the space-time domain for the 1D problem having time singularity based on the $L^\infty$ norm for a suitable choice of the grading parameter. Furthermore, it provides $O(N^{-2}+M^{-3})$ accurate solution for the 2D problem having unbounded time derivative at $t=0$. The analysis also highlights a higher order accuracy for a sufficiently smooth solution resides in $C^3(\overline{\Omega}_t)$ even if the mesh is discretized uniformly. The truncation error estimates for the time-fractional operator, integral operator, and spatial derivatives are presented. Numerous tests are performed on several examples in support of the theoretical analysis. The advancement of the proposed methodology is demonstrated through the application of the time-fractional Fokker-Planck equation and the fractional-order viscoelastic dynamics having weakly singular kernels. It also confirms the superiority of the proposed method compared with existing approaches available in the literature.
We consider the configuration space of points on the two-dimensional sphere that satisfy a specific system of quadratic equations. We construct periodic orbits in this configuration space using elliptic theta functions and show that they satisfy semi-discrete analogues of mKdV and sine-Gordon equations. The configuration space we investigate corresponds to the state space of a linkage mechanism known as the Kaleidocycle, and the constructed orbits describe the characteristic motion of the Kaleidocycle. Our approach is founded on the relationship between the deformation of spatial curves and integrable systems, offering an intriguing example where an integrable system generates an orbit in the space of real solutions to polynomial equations defined by geometric constraints.
The inverse problems about fractional Calder\'on problem and fractional Schr\"odinger equations are of interest in the study of mathematics. In this paper, we propose the inverse problem to simultaneously reconstruct potentials and sources for fractional Schr\"odinger equations with internal source terms. We show the uniqueness for reconstructing the two terms under measurements from two different nonhomogeneous boundary conditions. By introducing the variational Tikhonov regularization functional, numerical method based on conjugate gradient method(CGM) is provided to realize this inverse problem. Numerical experiments are given to gauge the performance of the numerical method.
小貼士
登錄享
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191