Search results for "Jacobian matrix"
showing 10 items of 47 documents
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…
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…
Relativistic positioning: errors due to uncertainties in the satellite world lines
2014
Global navigation satellite systems use appropriate satellite constellations to get the coordinates of an user -close to Earth- in an almost inertial reference system. We have simulated both GPS and GALILEO constellations. Uncertainties in the satellite world lines lead to dominant positioning errors. In this paper, a detailed analysis of these errors is developed inside a great region surrounding Earth. This analysis is performed in the framework of the so-called relativistic positioning systems. Our study is based on the Jacobian, J, of the transformation giving the emission coordinates in terms of the inertial ones. Around points of vanishing J, positioning errors are too large. We show …
The Numerical Simulation of Relativistic Fluid Flow with Strong Shocks
2001
In this review we present and analyze the performance of a Go-dunov type method applied to relativistic fluid flow. Our model equations are the corresponding Euler equations for special relativistic hydrodynamics. By choosing an appropriate vector of unknowns, the equations of special relativistic fluid dynamics (RFD) can be written as a hyperbolic system of conservation laws. We give a complete description of the spectral decomposition of the Jacobian matrices associated to the fluxes in each spatial direction, (see (Donat et al., 1998), for details), which is the essential ingredient of the Godunov-type numerical method we propose in this paper. We also review a numerical flux formula tha…
Computing Strong Shocks in Ultrarelativistic Flows: A Robust Alternative
1999
In recent years, shock capturing methods have started to be used in numerical simulations in Relativistic Fluid Dynamics (RFD). These techniques lead to explicit numerical codes that are able to successfully simulate the extreme conditions of the ultrarelativistic regime. After [2], an explicit, ready-to-use description of the full spectral decomposition of the Jacobian matrices of the RFD system is available, and this allows us to implement Marquina’s scheme [3] in RFD. The scheme is seen to maintain the good behavior shown in [3] with respect to certain numerical pathologies.
Jacobian-free approximate solvers for hyperbolic systems: Application to relativistic magnetohydrodynamics
2017
Abstract We present recent advances in PVM (Polynomial Viscosity Matrix) methods based on internal approximations to the absolute value function, and compare them with Chebyshev-based PVM solvers. These solvers only require a bound on the maximum wave speed, so no spectral decomposition is needed. Another important feature of the proposed methods is that they are suitable to be written in Jacobian-free form, in which only evaluations of the physical flux are used. This is particularly interesting when considering systems for which the Jacobians involve complex expressions, e.g., the relativistic magnetohydrodynamics (RMHD) equations. On the other hand, the proposed Jacobian-free solvers hav…
Jacobian-Free Incomplete Riemann Solvers
2018
The purpose of this work is to present some recent developments about incomplete Riemann solvers for general hyperbolic systems. Polynomial Viscosity Matrix (PVM) methods based on internal approximations to the absolute value function are introduced, and they are compared with Chebyshev-based PVM solvers. These solvers only require a bound on the maximum wave speed, so no spectral decomposition is needed. Moreover, they can be written in Jacobian-free form, in which only evaluations of the physical flux are used. This is particularly interesting when considering systems for which the Jacobians involve complex expressions. Some numerical experiments involving the relativistic magnetohydrodyn…
Jacobian of weak limits of Sobolev homeomorphisms
2016
Abstract Let Ω be a domain in ℝ n {\mathbb{R}^{n}} , where n = 2 , 3 {n=2,3} . Suppose that a sequence of Sobolev homeomorphisms f k : Ω → ℝ n {f_{k}\colon\Omega\to\mathbb{R}^{n}} with positive Jacobian determinants, J ( x , f k ) > 0 {J(x,f_{k})>0} , converges weakly in W 1 , p ( Ω , ℝ n ) {W^{1,p}(\Omega,\mathbb{R}^{n})} , for some p ⩾ 1 {p\geqslant 1} , to a mapping f. We show that J ( x , f ) ⩾ 0 {J(x,f)\geqslant 0} a.e. in Ω. Generalizations to higher dimensions are also given.
On the accurate determination of nonisolated solutions of nonlinear equations
1981
A simple but efficient method to obtain accurate solutions of a system of nonlinear equations with a singular Jacobian at the solution is presented. This is achieved by enlarging the system to a higher dimensional one whose solution in question is isolated. Thus it can be computed e. g. by Newton's method, which is locally at least quadratically convergent and selfcorrecting, so that high accuracy is attainable.
Planar Mappings of Finite Distortion
2010
We review recent results on planar mappings of finite distortion. This class of mappings contains all analytic functions and quasiconformal mappings.