Temperature decay in an aluminium plate is observed using Galerkin finite element method for 2D transient heat conduction equation. Taking ∆t of 0.001 hr , temperature variation is studied using unconditionally stable first order and second order accurate schemes backward Euler and modified CrankNicholson respectively. A comparative study has been made taking different combinations of meshes an...

The linearized Cahn-Hilliard-Cook equation is discretized in the spatial variables by a standard finite element method. Strong convergence estimates are proved under suitable assumptions on the covariance operator of the Wiener process, which is driving the equation. The backward Euler time stepping is also studied. The analysis is set in a framework based on analytic semigroups. The main part ...

The semi-implicit Euler discretization method is studied for abstract evolution equations in a Hilbert space H, like _ u = f(t; u; u) ; t 2 (0; T ] ; u(0) = u0, where f(t; ; v) is one-sided Lipschitz and R(I hf(t; ; v)) = H for h > 0 su ciently small, and f(t; u; ) is Lipschitz-continuous. Extension to Banach spaces is then pointed out. Ordinary and partial , di erential and integro-di erential...

We explain an interface for the implementation of rate-independent elastoplasticity which separates the pointwise evaluation of the elastoplastic material law and the global solution of the momentum balance equation. The elastoplastic problem is discretized in time by the implicit Euler method and every time step is solved with a Newton iteration. For the discretization in space the material pa...

In population balance equations, the distribution of the entities depends not only on space and time but also on their own properties referred to as internal coordinates. The operator splitting method is used to transform the whole time-dependent problem into two unsteady subproblems of a smaller complexity. The first subproblem is a timedependent convection-diffusion problem while the second o...

We consider the numerical solution of hyperbolic systems of conservation laws with relaxation using a shock capturing nite diierence scheme on a xed, uniform spatial grid. We conjecture that certain a priori criteria insure that the numerical method does not produce spurious solutions as the relaxation time vanishes. One criterion is that the limits of vanishing relaxation time and vanishing vi...

We analyze the spatially semidiscrete piecewise linear finite element method for a nonlocal parabolic equation resulting from thermistor problem. Our approach is based on the properties of the elliptic projection defined by the bilinear form associated with the variational formulation of the finite element method. We assume minimal regularity of the exact solution that yields optimal order erro...

The classical two-phase Stefan problem in level set formulation is considered. The implementation of a solver on triangular grids is described. Extended finite elements (X-FEM) in space and an implicit Euler method in time are used to approximate the temperature. For the level set equation, a discontinuous Galerkin (DG) and a strong stability preserving (SSP) RungeKutta scheme are employed. Pol...

We propose an efficient and spectrally accurate numerical method for computing vortex lattice structures in fast rotating Bose-Einstein condensates (BECs) with strongly repulsive interactions. The key ingredients of the method is to discretize the normalized gradient flow under rotational frame by Fourier spectral method in space and by backward Euler method in time. Different vortex lattice st...

In this paper, an adaptive approach based on a time discontinuous Galerkin method (dG) is proposed. For large time steps and in a linear non-autonomous case, it reduces to the asymptotically second-order accurate Lobatto IIIC type implicit Runge-Kutta method, which is equal to the Padé-(0,2) approximation of the exponential function. In the asymptotic range it will result in the asymptotically ...