Search results for "numerical analysis"

Numerical approximation of the viscous quantum hydrodynamic model for semiconductors

2006

The viscous quantum hydrodynamic equations for semiconductors with constant temperature are numerically studied. The model consists of the one-dimensional Euler equations for the electron density and current density, including a quantum correction and viscous terms, coupled to the Poisson equation for the electrostatic potential. The equations can be derived formally from a Wigner-Fokker-Planck model by a moment method. Two different numerical techniques are used: a hyperbolic relaxation scheme and a central finite-difference method. By simulating a ballistic diode and a resonant tunneling diode, it is shown that numerical or physical viscosity changes significantly the behavior of the solu…

Stochastic 0-dimensional Biogeochemical Flux Model: Effect of temperature fluctuations on the dynamics of the biogeochemical properties in a marine e…

2021

Abstract We present a new stochastic model, based on a 0-dimensional version of the well known biogeochemical flux model (BFM), which allows to take into account the temperature random fluctuations present in natural systems and therefore to describe more realistically the dynamics of real marine ecosystems. The study presents a detailed analysis of the effects of randomly varying temperature on the lower trophic levels of the food web and ocean biogeochemical processes. More in detail, the temperature is described as a stochastic process driven by an additive self-correlated Gaussian noise. Varying both correlation time and intensity of the noise source, the predominance of different plank…

Free surfaces: shape sensitivity analysis and numerical methods

1999

Matrices A such that A^{s+1}R = RA* with R^k = I

2018

[EN] We study matrices A is an element of C-n x n such that A(s+1)R = RA* where R-k = I-n, and s, k are nonnegative integers with k >= 2; such matrices are called {R, s+1, k, *}-potent matrices. The s = 0 case corresponds to matrices such that A = RA* R-1 with R-k = I-n, and is studied using spectral properties of the matrix R. For s >= 1, various characterizations of the class of {R, s + 1, k, *}-potent matrices and relationships between these matrices and other classes of matrices are presented. (C) 2018 Elsevier Inc. All rights reserved.

The local boundedness of solutions for a class of degenerate nonlinear elliptic higher-order equations withL1-data

2008

We prove local boundedness of solutions for a class of degenerate nonlinear elliptic higher-order equations with L(1)-data.

Unified Analysis of Periodization-Based Sampling Methods for Matérn Covariances

2020

The periodization of a stationary Gaussian random field on a sufficiently large torus comprising the spatial domain of interest is the basis of various efficient computational methods, such as the ...

On the implementation of weno schemes for a class of polydisperse sedimentation models

2011

The sedimentation of a polydisperse suspension of small rigid spheres of the same density, but which belong to a finite number of species (size classes), can be described by a spatially one-dimensional system of first-order, nonlinear, strongly coupled conservation laws. The unknowns are the volume fractions (concentrations) of each species as functions of depth and time. Typical solutions, e.g. for batch settling in a column, include discontinuities (kinematic shocks) separating areas of different composition. The accurate numerical approximation of these solutions is a challenge since closed-form eigenvalues and eigenvectors of the flux Jacobian are usually not available, and the characte…

Capturing Shock Reflections: An Improved Flux Formula

1996

Godunov type schemes, based on exact or approximate solutions to the Riemann problem, have proven to be an excellent tool to compute approximate solutions to hyperbolic systems of conservation laws. However, there are many instances in which a particular scheme produces inappropriate results. In this paper we consider several situations in which Roe's scheme gives incorrect results (or blows up all together) and we propose an alternative flux formula that produces numerical approximations in which the pathological behavior is either eliminated or reduced to computationally acceptable levels.

Power ENO methods: a fifth-order accurate Weighted Power ENO method

2004

In this paper we introduce a new class of ENO reconstruction procedures, the Power ENO methods, to design high-order accurate shock capturing methods for hyperbolic conservation laws, based on an extended class of limiters, improving the behavior near discontinuities with respect to the classical ENO methods. Power ENO methods are defined as a correction of classical ENO methods [J. Comput. Phys. 71 (1987) 231], by applying the new limiters on second-order differences or higher. The new class of limiters includes as a particular case the minmod limiter and the harmonic limiter used for the design of the PHM methods [see SIAM J. Sci. Comput. 15 (1994) 892]. The main features of these new ENO…

A Flux-Split Algorithm Applied to Relativistic Flows

1998

The equations of RFD can be written as a hyperbolic system of conservation laws by choosing an appropriate vector of unknowns. We give an explicit formulation of the full spectral decomposition of the Jacobian matrices associated with the fluxes in each spatial direction, which is the essential ingredient of the techniques we propose in this paper. These techniques are based on the recently derived flux formula of Marquina, a new way to compute the numerical flux at a cell interface which leads to a conservative, upwind numerical scheme. Using the spectral decompositions in a fundamental way, we construct high order versions of the basic first-order scheme described by R. Donat and A. Marqu…