Search results for "Computational Mathematic"

Scheduled Relaxation Jacobi method: improvements and applications

2016

Elliptic partial differential equations (ePDEs) appear in a wide variety of areas of mathematics, physics and engineering. Typically, ePDEs must be solved numerically, which sets an ever growing demand for efficient and highly parallel algorithms to tackle their computational solution. The Scheduled Relaxation Jacobi (SRJ) is a promising class of methods, atypical for combining simplicity and efficiency, that has been recently introduced for solving linear Poisson-like ePDEs. The SRJ methodology relies on computing the appropriate parameters of a multilevel approach with the goal of minimizing the number of iterations needed to cut down the residuals below specified tolerances. The efficien…

A 1D coupled Schrödinger drift-diffusion model including collisions

2005

We consider a one-dimensional coupled stationary Schroedinger drift-diffusion model for quantum semiconductor device simulations. The device domain is decomposed into a part with large quantum effects (quantum zone) and a part where quantum effects are negligible (classical zone). We give boundary conditions at the classic-quantum interface which are current preserving. Collisions within the quantum zone are introduced via a Pauli master equation. To illustrate the validity we apply the model to three resonant tunneling diodes.

A partially reflecting random walk on spheres algorithm for electrical impedance tomography

2015

In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias…

Anharmonic effects on the dynamic behavior’s of Klein Gordon model’s

2021

Abstract This work completes and extends the Ref. Tchakoutio Nguetcho et al. (2017), in which we have focused our attention only on the dynamic behavior of gap soliton solutions of the anharmonic Klein-Gordon model immersed in a parameterized on-site substrate potential. We expand our work now inside the permissible frequency band. These considerations have crucial effects on the response of nonlinear excitations that can propagate along this model. Moreover, working in the allowed frequency band is not only interesting from a physical point of view, it also provides an extraordinary mathematical model, a new class of differential equations possessing vital parameters and vertical singular …

Collision orbits in the oblate planet problem

1984

Some of the properties of the oblate planet problem are derived. We use the technique of blowing up the singularity to study the collision orbits. We define some families of them in terms of their asymptotic behavior.

Classes of orbits in the main problem of satellite theory

1986

We consider the main problem in satellite theory restricted to the polar plane. For suitable values of the energy the system has two unstable periodic orbits. We classify the trajectories in terms of their ultimate behavior with respect these periodic orbits in: oscillating, asymptotic and capture orbits. We study the energy level set and the existence and properties of the mentioned types of motion.

High-order methods for the simulation of hydromagnetic instabilities in core-collapse supernovae

2011

AbstractWe present an assessment of the accuracy of a recently developed MHD code used to study hydromagnetic flows in supernovae and related events. The code, based on the constrained transport formulation, incorporates unprecedented ultra-high-order methods (up to 9th order) for the reconstruction and the most accurate approximate Riemann solvers. We estimate the numerical resistivity of these schemes in tearing instability simulations.

Existence of global weak solutions to the kinetic Peterlin model

2018

Abstract We consider a class of kinetic models for polymeric fluids motivated by the Peterlin dumbbell theories for dilute polymer solutions with a nonlinear spring law for an infinitely extensible spring. The polymer molecules are suspended in an incompressible viscous Newtonian fluid confined to a bounded domain in two or three space dimensions. The unsteady motion of the solvent is described by the incompressible Navier–Stokes equations with the elastic extra stress tensor appearing as a forcing term in the momentum equation. The elastic stress tensor is defined by Kramer’s expression through the probability density function that satisfies the corresponding Fokker–Planck equation. In thi…

MAGPACK1A package to calculate the energy levels, bulk magnetic properties, and inelastic neutron scattering spectra of high nuclearity spin clusters

2001

Second generation Car-Parrinello molecular dynamics

2014

Computer simulation methods, such as Monte Carlo or molecular dynamics, are very powerful theoretical techniques to provide detailed and essentially exact informations on rather complex classical many-body problems. With the advent of ab initio molecular dynamics (AIMD), where finite-temperature dynamical trajectories are generated using interatomic forces which are calculated on the fly using accurate electronic structure calculations, the scope of computational research has been greatly extended. This review is intended to outline the basic principles as well as being a survey of the field. Beginning with the derivation of Born–Oppenheimer molecular dynamics, the Car–Parrinello method and…