数值分析
Support for arithmetic in multiple precisions and number formats is becoming increasingly common in emerging high-performance architectures. From a computational scientist's perspective, our goal is to determine how and where we can safely…
This paper is part of a series developing $C^0$ finite element methods for fourth-order elliptic equations on polygonal domains. Here, we investigate how boundary conditions influence the design of effective $C^0$ schemes, specifically…
We revisit the inverse problem of reconstructing a spatially varying diffusion coefficient in stationary elliptic equations from boundary Cauchy data. From a theoretical perspective, we introduce a gradient-weighted modification of the…
We introduce a framework for repurposing error estimators for source problems to compute an estimator for the gap between eigenspaces and their discretizations. Of interest are eigenspaces of finite clusters of eigenvalues of unbounded…
This article is concerned with the efficient computation of modular matrix multiplication C=AB mod p, a key kernel in computer algebra. We focus on floating-point arithmetic, which allows for using efficient matrix multiplication libraries.…
Accurately simulating wave propagation is crucial in fields such as acoustics, electromagnetism, and seismic analysis. Traditional numerical methods, like finite difference and finite element approaches, are widely used to solve governing…
This paper investigates a class of non-autonomous highly oscillatory ordinary differential equations characterized by a linear component inversely proportional to a small parameter $\varepsilon$, with purely imaginary eigenvalues, and an…
This work presents a mathematical model to enable rapid prediction of airborne contaminant transport based on scarce sensor measurements. The method is designed for applications in critical infrastructure protection (CIP), such as…
This paper concerns the numerical valuation of swing options with discrete action times under a linear two-factor mean-reverting model with jumps. The resulting sequence of two-dimensional partial integro-differential equations (PIDEs) are…
This paper is interested in the computation of stresses within jammed packings of rigid polygonal cells. The cells are considered to follow a Tresca friction law. First, a constrained minimization problem is introduced where the friction…
Rapid and reliable solvers for parametric partial differential equations (PDEs) are needed in many scientific and engineering disciplines. For example, there is a growing demand for composites and architected materials with heterogeneous…
We propose and analyze a structure-preserving approximation of the non-isothermal Cahn-Hilliard equation using conforming finite elements for the spatial discretization and a problem-specific mixed explicit-implicit approach for the…
Background and Objective: Transcranial temporal interference stimulation (tTIS) is a promising non-invasive brain stimulation technique in which interference between electrical current fields extends the possibilities of electrical brain…
We consider a nonlinear Klein--Gordon equation in the nonrelativistic limit regime with highly oscillatory initial data in the form of a modulated plane wave. In this regime, the solution exhibits rapid oscillations in both time and space,…
A conforming finite element scheme with mixed explicit-implicit time discretization for quasi-incompressible Navier-Stokes-Maxwell-Stefan systems in a bounded domain with periodic boundary conditions is presented. The system consists of the…
We analyze a coupled Cahn-Hilliard-Forchheimer system featuring concentration-dependent mobility, mass source and convective transport. The velocity field is governed by a generalized quasi-incompressible Forchheimer equation with…
The two-phase Navier-Stokes Cahn-Hilliard (NSCH) mixture model is a key framework for simulating multiphase flows with non-matching densities. Developing fully discrete, energy-stable schemes for this model remains challenging, due to the…
This paper proposes an adaptive hyper-reduction method to reduce the computational cost associated with the simulation of parametric particle-based kinetic plasma models, specifically focusing on the Vlasov-Poisson equation. Conventional…
The Ohta-Kawasaki equation models the mesoscopic phase separation of immiscible polymer chains that form diblock copolymers, with applications in directed self-assembly for lithography. We perform a mathematical analysis of this model under…
In view of the advantages of simplicity and effectiveness of the Kaczmarz method, which was originally employed to solve the large-scale system of linear equations $Ax=b$, we study the greedy randomized block Kaczmarz method (ME-GRBK) and…