Motivated by numerical modeling of ultrasound waves, we investigate robust conforming finite element discretizations of quasilinear and possibly nonlocal equations of Westervelt type. These wave equations involve either a strong dissipation or damping of fractional-derivative type and we unify them into one class by introducing a memory kernel that satisfies non-restrictive regularity and positivity assumptions. As the involved damping parameter is relatively small and can become negligible in certain (inviscid) media, it is important to develop methods that remain stable as the said parameter vanishes. To this end, the contributions of this work are twofold. First, we determine sufficient conditions under which conforming finite element discretizations of (non)local Westervelt equations can be made robust with respect to the dissipation parameter. Secondly, we establish the rate of convergence of the semi-discrete solutions in the singular vanishing dissipation limit. The analysis hinges upon devising appropriate energy functionals for the semi-discrete solutions that remain uniformly bounded with respect to the damping parameter.
相關內容
Group testing enables to identify infected individuals in a population using a smaller number of tests than individual testing. To achieve this, group testing algorithms commonly assume knowledge of the number of infected individuals; nonadaptive and several adaptive algorithms fall in this category. Some adaptive algorithms, like binary splitting, operate without this assumption, but require a number of stages that may scale linearly with the size of the population. In this paper we contribute a new algorithm that enables a balance between the number of tests and the number of stages used, and which we term diagonal group testing. Diagonal group testing, like binary splitting, does not require knowledge of the number of infected individuals, yet unlike binary splitting, is order-optimal w.r.t. the expected number of tests it requires and is guaranteed to succeed in a small number of stages that scales at most logarithmically with the size of the population. Numerical evaluations, for diagonal group testing and a hybrid approach we propose, support our theoretical findings.
We present a rigorous theoretical analysis of the convergence rate of the deep mixed residual method (MIM) when applied to a linear elliptic equation with various types of boundary conditions. The MIM method has been proposed as a more effective numerical approximation method compared to the deep Galerkin method (DGM) and deep Ritz method (DRM) in various cases. Our analysis shows that MIM outperforms DRM and deep Galerkin method for weak solution (DGMW) in the Dirichlet case due to its ability to enforce the boundary condition. However, for the Neumann and Robin cases, MIM demonstrates similar performance to the other methods. Our results provides valuable insights into the strengths of MIM and its comparative performance in solving linear elliptic equations with different boundary conditions.
We are interested in the discretisation of a drift-diffusion system in the framework of hybrid finite volume (HFV) methods on general polygonal/polyhedral meshes. The system under study is composed of two anisotropic and nonlinear convection-diffusion equations with nonsymmetric tensors, coupled with a Poisson equation and describes in particular semiconductor devices immersed in a magnetic field. We introduce a new scheme based on an entropy-dissipation relation and prove that the scheme admits solutions with values in admissible sets - especially, the computed densities remain positive. Moreover, we show that the discrete solutions to the scheme converge exponentially fast in time towards the associated discrete thermal equilibrium. Several numerical tests confirm our theoretical results. Up to our knowledge, this scheme is the first one able to discretise anisotropic drift-diffusion systems while preserving the bounds on the densities.
In this work, two Crank-Nicolson schemes without corrections are developed for sub-diffusion equations. First, we propose a Crank-Nicolson scheme without correction for problems with regularity assumptions only on the source term. Second, since the existing Crank-Nicolson schemes have a severe reduction of convergence order for solving sub-diffusion equations with singular source terms in time, we then extend our scheme and propose a new Crank-Nicolson scheme for problems with singular source terms in time. Second-order error estimates for both the two Crank-Nicolson schemes are rigorously established by a Laplace transform technique, which are numerically verified by some numerical examples.
High-order finite-difference entropy stable schemes for two-fluid relativistic plasma flow equations
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux components. The first two components describe the ion and electron flows, which are modeled using the relativistic hydrodynamics equation. The third component is Maxwell's equations, which are linear systems. The coupling of the ion and electron flows, and electromagnetic fields is via source terms only. Furthermore, we also show that the source terms do not affect the entropy evolution. To design semi-discrete entropy stable schemes, we extend the RHD entropy stable schemes in Bhoriya et al. to three dimensions. This is then coupled with entropy stable discretization of the Maxwell's equations. Finally, we use SSP-RK schemes to discretize in time. We also propose ARK-IMEX schemes to treat the stiff source terms; the resulting nonlinear set of algebraic equations is local (at each discretization point). These equations are solved using the Newton's Method, which results in an efficient method. The proposed schemes are then tested using various test problems to demonstrate their stability, accuracy and efficiency.
In this paper, we are interested in constructing a scheme solving compressible Navier--Stokes equations, with desired properties including high order spatial accuracy, conservation, and positivity-preserving of density and internal energy under a standard hyperbolic type CFL constraint on the time step size, e.g., $\Delta t=\mathcal O(\Delta x)$. Strang splitting is used to approximate convection and diffusion operators separately. For the convection part, i.e., the compressible Euler equation, the high order accurate postivity-preserving Runge--Kutta discontinuous Galerkin method can be used. For the diffusion part, the equation of internal energy instead of the total energy is considered, and a first order semi-implicit time discretization is used for the ease of achieving positivity. A suitable interior penalty discontinuous Galerkin method for the stress tensor can ensure the conservation of momentum and total energy for any high order polynomial basis. In particular, positivity can be proven with $\Delta t=\mathcal{O}(\Delta x)$ if the Laplacian operator of internal energy is approximated by the $\mathbb{Q}^k$ spectral element method with $k=1,2,3$. So the full scheme with $\mathbb{Q}^k$ ($k=1,2,3$) basis is conservative and positivity-preserving with $\Delta t=\mathcal{O}(\Delta x)$, which is robust for demanding problems such as solutions with low density and low pressure induced by high-speed shock diffraction. Even though the full scheme is only first order accurate in time, numerical tests indicate that higher order polynomial basis produces much better numerical solutions, e.g., better resolution for capturing the roll-ups during shock reflection.
In this work, following the Discrete de Rham (DDR) paradigm, we develop an arbitrary-order discrete divdiv complex on general polyhedral meshes. The construction rests 1) on discrete spaces that are spanned by vectors of polynomials whose components are attached to mesh entities and 2) on discrete operators obtained mimicking integration by parts formulas. We provide an in-depth study of the algebraic properties of the local complex, showing that it is exact on mesh elements with trivial topology. The new DDR complex is used to design a numerical scheme for the approximation of biharmonic problems, for which we provide detailed stability and convergence analyses.
This is Part II of our paper in which we prove finite time blowup of the 2D Boussinesq and 3D axisymmetric Euler equations with smooth initial data of finite energy and boundary. In Part I of our paper [ChenHou2023a], we establish an analytic framework to prove stability of an approximate self-similar blowup profile by a combination of a weighted $L^\infty$ norm and a weighted $C^{1/2}$ norm. Under the assumption that the stability constants, which depend on the approximate steady state, satisfy certain inequalities stated in our stability lemma, we prove stable nearly self-similar blowup of the 2D Boussinesq and 3D Euler equations with smooth initial data and boundary. In Part II of our paper, we provide sharp stability estimates of the linearized operator by constructing space-time solutions with rigorous error control. We also obtain sharp estimates of the velocity in the regular case using computer assistance. These results enable us to verify that the stability constants obtained in Part I [ChenHou2023a] indeed satisfy the inequalities in our stability lemma. This completes the analysis of the finite time singularity of the axisymmetric Euler equations with smooth initial data and boundary.
We introduce a priori Sobolev-space error estimates for the solution of nonlinear, and possibly parametric, PDEs using Gaussian process and kernel based methods. The primary assumptions are: (1) a continuous embedding of the reproducing kernel Hilbert space of the kernel into a Sobolev space of sufficient regularity; and (2) the stability of the differential operator and the solution map of the PDE between corresponding Sobolev spaces. The proof is articulated around Sobolev norm error estimates for kernel interpolants and relies on the minimizing norm property of the solution. The error estimates demonstrate dimension-benign convergence rates if the solution space of the PDE is smooth enough. We illustrate these points with applications to high-dimensional nonlinear elliptic PDEs and parametric PDEs. Although some recent machine learning methods have been presented as breaking the curse of dimensionality in solving high-dimensional PDEs, our analysis suggests a more nuanced picture: there is a trade-off between the regularity of the solution and the presence of the curse of dimensionality. Therefore, our results are in line with the understanding that the curse is absent when the solution is regular enough.
Knowledge graph (KG) embeddings learn low-dimensional representations of entities and relations to predict missing facts. KGs often exhibit hierarchical and logical patterns which must be preserved in the embedding space. For hierarchical data, hyperbolic embedding methods have shown promise for high-fidelity and parsimonious representations. However, existing hyperbolic embedding methods do not account for the rich logical patterns in KGs. In this work, we introduce a class of hyperbolic KG embedding models that simultaneously capture hierarchical and logical patterns. Our approach combines hyperbolic reflections and rotations with attention to model complex relational patterns. Experimental results on standard KG benchmarks show that our method improves over previous Euclidean- and hyperbolic-based efforts by up to 6.1% in mean reciprocal rank (MRR) in low dimensions. Furthermore, we observe that different geometric transformations capture different types of relations while attention-based transformations generalize to multiple relations. In high dimensions, our approach yields new state-of-the-art MRRs of 49.6% on WN18RR and 57.7% on YAGO3-10.
小貼士
登錄享
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191