Related papers: Fast computation of optimal damping parameters for…
We perform a linear stability analysis of three-layer radial porous media and Hele-Shaw flows with variable viscosity in the middle layer. A nonlinear change of variables results in an eigenvalue problem that has time-dependent coefficients…
We consider computing the $k$-th eigenvalue and its corresponding eigenvector of a generalized Hermitian eigenvalue problem of $n\times n$ large sparse matrices. In electronic structure calculations, several properties of materials, such as…
We introduce a randomized algorithm for computing the minimal-norm solution to an underdetermined system of linear equations. Given an arbitrary full-rank m x n matrix A with m<n, any m x 1 vector b, and any positive real number epsilon…
We consider the bi-Laplacian eigenvalue problem for the modes of vibration of a thin elastic plate with a discrete set of clamped points. A high-order boundary integral equation method is developed for efficient numerical determination of…
We study the large time behavior of solutions to two-dimensional Euler and Navier-Stokes equations linearized about shear flows of the mixing layer type in the unbounded channel $\mathbb{T} \times \mathbb{R}$. Under a simple spectral…
Lipid vesicles appear ubiquitously in biological systems. Understanding how the mechanical and intermolecular interations deform vesicle membrane is a fundamental question in biophysics. In this article we developed a fast algorithm to…
We study the Cauchy problem for the equation of the form $$ \ddot{u}(t) + (\aa A + B)\dot{u}(t) + (A+G)u(t) = 0,\tag* $$ where $A$, $B$, and $G$ are \o s in a Hilbert space $\Cal H$ with $A$ selfadjoint, $\sigma(A)=[0,\infty)$, $B\ge0$…
We provide faster algorithms and improved sample complexities for approximating the top eigenvector of a matrix. Offline Setting: Given an $n \times d$ matrix $A$, we show how to compute an $\epsilon$ approximate top eigenvector in time…
In this paper we propose a numerical method to solve the Cauchy problem based on the viscous shallow water equations in an horizontally moving domain. More precisely, we are interested in a flooding and drying model, used to modelize the…
We apply the approximate dynamics derived from the Gaussian time-dependent variational principle to the Hamiltonian $ \hat H= {1/2}(\hat p_x ^2+ \hat p_y ^2)+ {1/2}\hat x^2\hat y^2$, which is strongly chaotic in the classical limit. We are…
The largest eigenvalue of a matrix is always larger or equal than its largest diagonal entry. We show that for a large class of random Laplacian matrices, this bound is essentially tight: the largest eigenvalue is, up to lower order terms,…
We devise a new time-stepping algorithm for two-dimensional nonlinear unsteady surface and interfacial waves. The algorithm uses Cauchy's integral formula, which only requires information on the interface, to solve Laplace equation by using…
We introduce a high-order finite element method for approximating the Vlasov-Poisson equations. This approach employs continuous Lagrange polynomials in space and explicit Runge-Kutta schemes for time discretization. To stabilize the…
In this paper, the Cauchy problem for the three-dimensional (3-D) isentropic compressible Navier-Stokes equations with degenerate viscosities is considered. By introducing some new variables and making use of the "quasi-symmetric…
In this paper, by means of the Riesz basis approach, we study the stability of a weakly damped system of two second order evolution equations coupled through the velocities. If the fractional order damping becomes viscous and the waves…
The computation of damping rates of an oscillating fluid with a free surface in which viscosity is small and surface tension high is numerically challenging. A typical application requiring such computation is drop-on-demand (DoD)…
This work concerns the minimization of the pseudospectral abscissa of a matrix-valued function dependent on parameters analytically. The problem is motivated by robust stability and transient behavior considerations for a linear control…
We consider the problem of computing the optimal solution and objective of a linear program under linearly changing linear constraints. The problem studied is given by $\min c^t x \text{ s.t } Ax + \lambda Dx \leq b$ where $\lambda$ belongs…
We present an O(N) algorithm to study the vibrational properties of amorphous silicon within the framework of tight-binding approach. The dynamical matrix elements have been evaluated numerically in the harmonic approximation exploiting the…
In this article, we consider the frame hydrodynamics of biaxial nematic phases, a coupled system between the evolution of the orthonormal frame and the Navier--Stokes equation, which is derived from a molecular-theory-based dynamical tensor…