Related papers: Solving Hodgkin-Huxley equations using the compact…
We develop and analyze the first hybridizable discontinuous Galerkin (HDG) method for solving fifth-order Korteweg-de Vries (KdV) type equations. We show that the semi-discrete scheme is stable with proper choices of the stabilization…
In this paper we suggest the use of light for performing useful computations. Namely, we propose a special device which uses light rays for solving the Hamiltonian path problem on a directed graph. The device has a graph-like representation…
The Hodgkin-Huxley (HH) model is a powerful model to explain different aspects of spike generation in excitable cells. However, the HH model was proposed in 1952 when the real structure of the ion channel was unknown. It is now common…
We describe a rapidly converging algorithm for solving the Kohn--Sham equations and equations of similar structure that appear frequently in calculations of the structure of inhomogeneous electronic many--body systems. The algorithm has its…
This paper studies a new class of integration schemes for the numerical solution of semi-explicit differential-algebraic equations of differentiation index 2 in Hessenberg form. Our schemes provide the flexibility to choose different…
Synchronizations of processing elements (PEs) in massively parallel simulations, which arise due to communication or load imbalances between PEs, significantly affect the scalability of scientific applications. We have recently proposed a…
In a previous work, the author and D.C. Dobson proposed a numerical method for solving the complex Helmholtz equation based on the minimization variational principles developed by Milton, Seppecher, and Bouchitte. This method results in a…
We leverage commutative hypercomplex analysis to find closed-form solutions of some systems of stochastic differential equations. Specifically, we obtain necessary and sufficient conditions under which a system of stochastic differential…
The density functional approach in the Kohn-Sham approximation is widely used to study properties of many-electron systems. Due to the nonlinearity of the Kohn-Sham equations, the general self-consistence searching method involves…
In this paper, a class of high-order methods to numerically solve Functional Differential Equations with Piecewise Continuous Arguments (FDEPCAs) is discussed. The framework stems from the expansion of the vector field associated with the…
Complex geometric optics solutions to a system of d-bar equations appearing in the context of electrical impedance tomography and the scattering theory of the integrable Davey-Stewartson II equations are studied for large values of the…
A Godunov-type finite volume scheme on unstructured triangular grids is proposed to numerically solve the Savage-Hutter equations in curvilinear coordinate. We show the direct observation that the model is a not Galilean invariant system.…
The purpose of this work is to introduce and analyze a numerical scheme to efficiently solve boundary value problems involving the spectral fractional Laplacian. The approach is based on a reformulation of the problem posed on a…
In this paper a hybridized weak Galerkin (HWG) finite element method for solving the Stokes equations in the primary velocity-pressure formulation is introduced. The WG method uses weak functions and their weak derivatives which are defined…
The two-fluid plasma model has a wide range of timescales which must all be numerically resolved regardless of the timescale on which plasma dynamics occurs. The answer to solving numerically stiff systems is generally to utilize…
A novel implicit integration scheme for the Discrete Element Method (DEM) based on the variational integrator approach is presented. The numerical solver provides a fully dynamical description that, notably, reduces to an energy…
The dynamics of a many-particle system are often modeled by mapping the Hamiltonian onto a Schr\"odinger equation. An alternative approach is to solve the Hamiltonian equations directly in a model space of many-body configurations. In a…
We propose a multilevel approach for trace systems resulting from hybridized discontinuous Galerkin (HDG) methods. The key is to blend ideas from nested dissection, domain decomposition, and high-order characteristic of HDG discretizations.…
To describe highly heterogeneous systems using the Cahn-Hilliard equation, the standard form of the thermodynamic potential with a constant coefficient in the gradient term and a polynomial of the fourth degree may not be sufficient. The…
Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of…