Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Numerical assessment of mass conservation on a MAC-type method for viscoelastic free surface flows

A numerical study of mass conservation of MAC-type methods is presented, for viscoelastic free-surface flows. We use an implicit formulation which allows for greater time steps, and therefore time marching schemes for advecting the free surface marker particles have to be accurate in order to preserve the good mass conservation properties of this methodology. We then present an improvement by using a Runge-Kutta scheme coupled with a local linear extrapolation on the free surface. A thorough study of the viscoelastic impacting drop problem, for both Oldroyd-B and XPP fluid models, is presented, investigating the influence of timestep, grid spacing and other model parameters to the overall mass conservation of the method. Furthermore, an unsteady fountain flow is also simulated to illustrate the low mass conservation error obtained.

Stability analysis of Crank-Nicolson and Euler schemes for time-dependent diffusion equations

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Variational and perturbative schemes for a spiked harmonic oscillator

A variational analysis of the spiked harmonic oscillator Hamiltonian operator - d2/dx2 + x2 + l(l + 1)/x2 + λ|x| -α, where α is a real positive parameter, is reported in this work. The formalism makes use of the functional space spanned by the solutions of the Schrödinger equation for the linear harmonic oscillator Hamiltonian supplemented by a Dirichlet boundary condition, and a standard procedure for diagonalizing symmetric matrices. The eigenvalues obtained by increasing the dimension of the basis set provide accurate approximations for the ground state energy of the model system, valid for positive and relatively large values of the coupling parameter λ. Additionally, a large coupling perturbative expansion is carried out and the contributions up to fourth-order to the ground state energy are explicitly evaluated. Numerical results are compared for the special case α = 5/2. © 1989 American Institute of Physics.