Related papers: The energy technique for the six-step BDF method
Consider a multidimensional SDE of the form $X_t = x+\int_{0}^{t} b(X_{s-})ds+\int{0}^{t} f(X_{s-})dZ_s$ where $(Z_s)_{s\ge 0}$ is a symmetric stable process. Under suitable assumptions on the coefficients the unique strong solution of the…
A filtration system, comprising a Biot poroelastic solid coupled to an incompressible Stokes free-flow, is considered in 3D. Across the flat 2D interface, the Beavers-Joseph-Saffman coupling conditions are taken. In the inertial, linear,…
We present an effective adaptive procedure for the numerical approximation of the steady-state Gross-Pitaevskii equation. Our approach is solely based on energy minimization, and consists of a combination of gradient flow iterations and…
We propose and implement a third-order accurate numerical scheme for the Landau-Lifshitz-Gilbert equation, which describes magnetization dynamics in ferromagnetic materials under large damping parameters. This method offers two key…
Consider the Schr\"odinger--Bopp--Podolsky system \[ \begin{cases} -\epsilon^2\Delta u+(V+K\phi)u=u|u|^{p-1};\newline \Delta^2\phi-\Delta\phi=4\pi K u^2 \end{cases} ~\text{in}~\mathbb{R}^3 \] for sufficiently small $\epsilon>0$, where…
A modified primal-dual weak Galerkin (M-PDWG) finite element method is designed for the second order elliptic equation in non-divergence form. Compared with the existing PDWG methods proposed in \cite{wwnondiv}, the system of equations…
We extend the Vanderbilt ultrasoft pseudopotential scheme by adding kinetic energy density terms, in order to use meta-GGA exchange potentials, such as the Becke-Johnson or Tran-Blaha potentials, in the planewave-pseudopotential…
Recently, energetic variational approach was employed to derive models for non-isothermal electrokinetics by Liu et. al \cite{Liu-Wu-Liu-CMS2018}. In particular, the Poisson-Nernst-Planck-Fourier (PNPF) system for the dynamics of $N$-ionic…
In this paper we investigate the existence of solutions to the following Schr\"{o}dinger system in the critical case \begin{equation*} -\Delta u_{i}+\lambda_{i}u_{i}=u_{i}\sum_{j = 1}^{d}\beta_{ij}u_{j}^{2} \text{ in } \Omega, \quad u_i=0…
We consider numerical methods for the Poisson-Nernst-Planck-Cahn-Hilliard (PNPCH) equations with steric interactions. We propose a novel energy stable numerical scheme that respects mass conservation and positivity at the discrete level.…
Based on previous works, analytical calculational procedures for dealing with the strongly interacting fermions ground state are further developed through a medium dependent potential in terms of the Bethe-Peierls contact interaction model.…
We study the existence of positive solutions on $\R^{N+1}$ to semilinear elliptic equation $-\Delta u+u=f(u)$ where $N\geq 1$ and $f$ is modeled on the power case $f(u)=|u|^{p-1}u$. Denoting with $c$ the mountain pass level of $\f(u)=\tfrac…
The stability of a Bose-Einstein condensed state of trapped ultra-cold atoms is investigated under the assumption of an attractive two-body and a repulsive three-body interaction. The Ginzburg-Pitaevskii-Gross (GPG) nonlinear Schr\"odinger…
The numerical matrix Numerov algorithm is used to solve the stationary Schr\"odinger equation for central Coulomb potentials. An efficient approximation for accelerating the convergence is proposed. The Numerov method is error-prone if the…
We report the computational discovery of complex, topologically charged, and spectrally stable states in three-dimensional multi-component nonlinear wave systems of nonlinear Schr{\"o}dinger type. While our computations relate to…
The multiplicative Newton-like method developed by the author et al. is extended to the situation where the dynamics is restricted to the orthogonal group. A general framework is constructed without specifying the cost function. Though the…
For the Landau--Lifshitz--Gilbert (LLG) equation of micromagnetics we study linearly implicit backward difference formula (BDF) time discretizations up to order $5$ combined with higher-order non-conforming finite element space…
We consider an elastic/viscoelastic transmission problem for the Bresse system with fully Dirichlet or Dirichlet-Neumann-Neumann boundary conditions. The physical model consists of three wave equations coupled in certain pattern. The system…
The Poisson-Nernst-Planck (PNP) system is a widely accepted model for simulation of ionic channels. In this paper, we design, analyze, and numerically validate a second order unconditional positivity-preserving scheme for solving a reduced…
We prove a Br\'ezis--Oswald type existence theorem for positive solutions of semilinear equations in an abstract setting in which the underlying linear operator has a compact positivity-improving resolvent. The assumptions imposed on the…