For the Maxwell equations in time-dependent media only finite difference schemes with time-dependent conductivity are known. In this paper we present a numerical scheme based on the Magnus expansion and operator splitting that can handle time-dependent permeability and permittivity too. We demonstrate our results with numerical tests.
Intrinsic equations represent promising approach for the description of rotor blade dynamics. They are the system of non-linear partial differential equations. Stability of numeric solution by the finite difference method is described. The stability is studied for various numerical schemes with different methods for the computation of spatial derivatives from time level (i.e., mean values of old and new time step) to (i.e., only from new time step). Stable solution was obtained only for schemes...
The topic of this paper is the numerical analysis of time periodic solution for electro-magnetic phenomena. The Limit Absorption Method (LAM) which forms the basis of our study is presented. Theoretical results have been proved in the linear finite dimensional case. This method is applied to scattering problems and transport of charged particles.
The topic of this paper is the numerical analysis of time periodic solution for electro-magnetic phenomena. The Limit Absorption Method (LAM) which forms the basis of our study is presented. Theoretical results have been proved in the linear finite dimensional case. This method is applied to scattering problems and transport of charged particles.
In this article we implement different numerical schemes to simulate the Schrödinger-Debye equations that occur in nonlinear optics. Since the existence of blow-up solutions is an open problem, we try to compute such solutions. The convergence of the methods is proved and simulations seem indeed to show that for at least small delays self-focusing solutions may exist.
In this article we implement different numerical schemes to simulate the Schrödinger-Debye equations that occur in nonlinear optics. Since the existence of blow-up solutions is an open problem, we try to compute such solutions. The convergence of the methods is proved and simulations seem indeed to show that for at least small delays self-focusing solutions may exist.
We deal with numerical analysis and simulations of the Davey-Stewartson equations which model, for example, the evolution of water surface waves. This time dependent PDE system is particularly interesting as a generalization of the 1-d integrable NLS to 2 space dimensions. We use a time splitting spectral method where we give a convergence analysis for the semi-discrete version of the scheme. Numerical results are presented for various blow-up phenomena of the equation, including blowup of defocusing,...
We deal with numerical analysis and simulations of the Davey-Stewartson equations which model, for example, the evolution of water surface waves. This time dependent PDE system is particularly interesting as a generalization of the 1-d integrable NLS to 2 space dimensions. We use a time splitting spectral method where we give a convergence analysis for the semi-discrete version of the scheme. Numerical results are presented for various blow-up phenomena of the equation, including blowup of defocusing,...
In this paper, we consider the blow-up solutions for a quasilinear parabolic partial differential equation . We numerically investigate the blow-up rates of these solutions by using a numerical method which is recently proposed by the authors [3].
The time-dependent intensity of a UV-photon source, located inside an interstellar cloud, is determined by formulating and solving an inverse problem for the integro-differential transport equation of photons in a one-dimensional slab. Starting from a discretizazion of the forward problem, an iterative procedure is used to compute the values of the source intensity at increasing values of the time.
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
The goal of this article is to analyze the observability properties for a space semi-discrete approximation scheme derived from a mixed finite element method of the 1d wave equation on nonuniform meshes. More precisely, we prove that observability properties hold uniformly with respect to the mesh-size under some assumptions, which, roughly, measures the lack of uniformity of the meshes, thus extending the work [Castro and Micu, Numer. Math.102 (2006) 413–462] to nonuniform meshes. Our results...
A 1-D model of a slab of glass of a small thickness is considered. The governing equations are those of the classical 1-D linear viscoelasticity. A load due to the temperature gradients is assumed. The aim is to model the process called annealing. It is shown that an additional load due to structural strain is crucial for the success of the model. Algorithms of a numerical solution of the governing equations are proposed. Numerical results are presented and commented.
Experimental evidence collected over the years shows that blood exhibits non-Newtonian characteristics such as shear-thinning, viscoelasticity, yield stress and thixotropic behaviour. Under certain conditions these characteristics become relevant and must be taken into consideration when modelling blood flow. In this work we deal with incompressible generalized Newtonian fluids, that account for the non-constant viscosity of blood, and present a new numerical method to handle fluid-rigid body interaction...