Search results for "modeling"

Corrigendum to “Fractional differential equations solved by using Mellin transform” [Commun Nonlinear Sci Numer Simul 19(7) (2014) 2220–2227]

2015

Implicit analytic solutions for a nonlinear fractional partial differential beam equation

2020

Abstract Analytic solutions in implicit form are derived for a nonlinear partial differential equation (PDE) with fractional derivative elements, which can model the dynamics of a deterministically excited Euler-Bernoulli beam resting on a viscoelastic foundation. Specifically, the initial-boundary value problem for the corresponding PDE is reduced to an initial value problem for a nonlinear ordinary differential equation in a Hilbert space. Next, by employing the cosine and sine families of operators, a variation of parameters representation of the solution map is introduced. Due to the presence of a nonlinear term, a local fixed point theorem is employed to prove the local existence and u…

Controllability method for the Helmholtz equation with higher-order discretizations

2007

We consider a controllability technique for the numerical solution of the Helmholtz equation. The original time-harmonic equation is represented as an exact controllability problem for the time-dependent wave equation. This problem is then formulated as a least-squares optimization problem, which is solved by the conjugate gradient method. Such an approach was first suggested and developed in the 1990s by French researchers and we introduce some improvements to its practical realization. We use higher-order spectral elements for spatial discretization, which leads to high accuracy and lumped mass matrices. Higher-order approximation reduces the pollution effect associated with finite elemen…

Monotonic solution of heterogeneous anisotropic diffusion problems

2013

Anisotropic problems arise in various areas of science and engineering, for example groundwater transport and petroleum reservoir simulations. The pure diffusive anisotropic time-dependent transport problem is solved on a finite number of nodes, that are selected inside and on the boundary of the given domain, along with possible internal boundaries connecting some of the nodes. An unstructured triangular mesh, that attains the Generalized Anisotropic Delaunay condition for all the triangle sides, is automatically generated by properly connecting all the nodes, starting from an arbitrary initial one. The control volume of each node is the closed polygon given by the union of the midpoint of…

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…