Search results for "numerical analysis"
showing 10 items of 883 documents
Efficient numerical method for simulating static and dynamic properties of superfluid helium
2004
Density functional theory (DFT) offers computationally affordable way of describing static and dynamic properties of superfluid 4He. In general, the DFT models yield single particle-like Schrodinger equations with a nonlinear potential term that accounts for all the many-body interactions. The resulting equations can be solved for small amplitude plane wave excitations in the bulk whereas fully numerical solution must be sought in more complicated cases. In this paper we propose a numerical method that can be used in solving the time-dependent nonlinear Schrodinger equation in both real and imaginary times. The method is based on operator splitting technique where each component operator is…
The MAST-edge centred lumped scheme for the flow simulation in variably saturated heterogeneous porous media
2012
A novel methodology is proposed for the solution of the flow equation in a variably saturated heterogeneous porous medium. The computational domain is descretized using triangular meshes and the governing PDEs are discretized using a lumped in the edge centres numerical technique. The dependent unknown variable of the problem is the piezometric head. A fractional time step methodology is applied for the solution of the original system, solving consecutively a prediction and a correction problem. A scalar potential of the flow field exists and in the prediction step a MArching in Space and Time (MAST) formulation is applied for the sequential solution of the Ordinary Differential Equation of…
Time-harmonic elasticity with controllability and higher-order discretization methods
2008
The time-harmonic solution of the linear elastic wave equation is needed for a variety of applications. The typical procedure for solving the time-harmonic elastic wave equation leads to difficulties solving large-scale indefinite linear systems. To avoid these difficulties, we consider the original time dependent equation with a method based on an exact controllability formulation. The main idea of this approach is to find initial conditions such that after one time-period, the solution and its time derivative coincide with the initial conditions.The wave equation is discretized in the space domain with spectral elements. The degrees of freedom associated with the basis functions are situa…
Solution of time-independent Schrödinger equation by the imaginary time propagation method
2007
Numerical solution of eigenvalues and eigenvectors of large matrices originating from discretization of linear and non-linear Schrodinger equations using the imaginary time propagation (ITP) method is described. Convergence properties and accuracy of 2nd and 4th order operator-splitting methods for the ITP method are studied using numerical examples. The natural convergence of the method is further accelerated with a new dynamic time step adjustment method. The results show that the ITP method has better scaling with respect to matrix size as compared to the implicitly restarted Lanczos method. An efficient parallel implementation of the ITP method for shared memory computers is also demons…
Fictitious Domain Methods for the Numerical Solution of Two-Dimensional Scattering Problems
1998
Fictitious domain methods for the numerical solution of two-dimensional scattering problems are considered. The original exterior boundary value problem is approximated by truncating the unbounded domain and by imposing a nonreflecting boundary condition on the artificial boundary. First-order, second-order, and exact nonreflecting boundary conditions are tested on rectangular and circular boundaries. The finite element discretizations of the corresponding approximate boundary value problems are performed using locally fitted meshes, and the discrete equations are solved with fictitious domain methods. A special finite element method using nonmatching meshes is considered. This method uses …
Implementation Aspects of 3D Lattice-BGK: Boundaries, Accuracy, and a New Fast Relaxation Method
1999
In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to present a comprehensive overview of some pitfalls and shortcomings of the lattice-BGK method and to introduce some new ideas useful in practical simulations. We begin with an evaluation of the widely used bounce-back boundary condition in staircase geometries by simulating flow in an inclined tube. It is shown that the bounce-back scheme is first-order accurate in space when the location of the non…
A novel method to simulate the 3D chlorophyll distribution in marine oligotrophic waters
2021
Abstract A 3D advection-diffusion-reaction model is proposed to investigate the abundance of phytoplankton in a difficult-to-access ecosystem such as the Gulf of Sirte (southern Mediterranean Sea) characterized by oligotrophic waters. The model exploits experimentally measured environmental variables to reproduce the dynamics of four populations that dominate phytoplankton community in the studied area: Synechococcus, Prochlorococcus HL, Prochlorococcus LL and picoeukaryotes. The theoretical results obtained for phytoplankton abundances are converted into chl-a and Dvchl-a concentrations, and the simulated vertical chlorophyll profiles are compared to the corresponding experimentally acquir…
Dynamics of a minimal consumer network with bi -directional influence
2018
Abstract We study the dynamics of a model of interdependent consumer behavior defined by a family of two-dimensional noninvertible maps. This family belongs to a class of coupled logistic maps with different nonlinearity parameters and coupling terms that depend on one variable only. In our companion paper we considered the case of independent consumers as well as the case of uni-directionally connected consumers. The present paper aims at describing the dynamics in the case of a bi-directional connection. In particular, we investigate the bifurcation structure of the parameter plane associated with the strength of coupling between the consumers, focusing on the mechanisms of qualitative tr…
Standard polynomials and matrices with superinvolutions
2016
Abstract Let M n ( F ) be the algebra of n × n matrices over a field F of characteristic zero. The superinvolutions ⁎ on M n ( F ) were classified by Racine in [12] . They are of two types, the transpose and the orthosymplectic superinvolution. This paper is devoted to the study of ⁎-polynomial identities satisfied by M n ( F ) . The goal is twofold. On one hand, we determine the minimal degree of a standard polynomial vanishing on suitable subsets of symmetric or skew-symmetric matrices for both types of superinvolutions. On the other, in case of M 2 ( F ) , we find generators of the ideal of ⁎-identities and we compute the corresponding sequences of cocharacters and codimensions.
On the reconstruction of discontinuous functions using multiquadric RBF–WENO local interpolation techniques
2020
Abstract We discuss several approaches involving the reconstruction of discontinuous one-dimensional functions using parameter-dependent multiquadric radial basis function (MQ-RBF) local interpolants combined with weighted essentially non-oscillatory (WENO) techniques, both in the computation of the locally optimized shape parameter and in the combination of RBF interpolants. We examine the accuracy of the proposed reconstruction techniques in smooth regions and their ability to avoid Gibbs phenomena close to discontinuities. In this paper, we propose a true MQ-RBF–WENO method that does not revert to the classical polynomial WENO approximation near discontinuities, as opposed to what was pr…