Using a combination of Kawashima- and Goodman-type energy estimates, we establish spectral stability of general small-amplitude relaxation shocks of symmetric dissipative systems. This extends previous results obtained by Plaza and Zumbrun [8] by singular perturbation techniques under an additional technical assumption, namely, that the background equation be noncharacteristic with respect to the shock.
Asymptotic stability; Hyperbolic systems with relaxation; Shock profiles; Spectral proper ties
In this paper we introduce a simulation algorithm based on fluid dynamic models to reproduce the behavior of traffic in a portion of the urban network in Rome. Numerical results, obtained comparing experimental data with numerical solutions, show the effectiveness of our approximation. (c) 2009 Elsevier Inc. All rights reserved.
Neural current imaging aims at analyzing the functionality of the human brain through the localization of those regions where the neural current flows. The reconstruction of an electric current distribution from its magnetic field measured in the outer space, gives rise to a highly ill-posed and ill-conditioned inverse problem. We use a joint sparsity constraint as a regularization term and we propose an efficient iterative thresholding algorithm to recover the current distribution. Some numerical tests are also displayed.
Electric current imaging
Magnetoencephalograpy
Inverse problem
Sparsity constraint
Iterative thresholding
Multiscale basis
Neuronal current imaging aims at analyzing the functionality of the
human brain through the localization of those regions where the neural
current flows. The reconstruction of an electric current distribution from
its magnetic field measured by sophisticated superconducting devices in a
noninvasive way, gives rise to a highly ill-posed and ill-conditioned inverse
problem.
Assuming that each component of the current density vector possesses the
same sparse representation with respect to a preassigned multiscale basis,
allows us to apply new regularization techniques to the magnetic inverse
problem. In particular, we use a joint sparsity constraint as a regulariza-
tion term and we propose an efficient iterative thresholding algorithm to
reconstruct the current distribution. Some bidimensional experiments are
presented in order to show the algorithm properties.
Magnetoencephalograpy
Inverse problem
Sparsity constraint
It- erative thresholding
Multiscale basis.
Direct numerical approximation of a continuous-time infinite horizon control problem, requires to recast the model as a discrete-time, finite-horizon control model. The quality of the optimization results can be heavily degraded if the discretization process does not take into account features of the original model to be preserved. Restricting their attention to optimal growh problems with a steady state, Mercenier and Michel in [1] and [2], studied the conditions to be imposed for ensuring that discrete first-order approximation models have the same steady states as the infinite-horizon continuous-times counterpart. Here we show that Mercenier and Michel scheme is a first order partitioned Runge-Kutta method applied to the state-costate differential system which arises from the Pontryagin maximum principle. The main consequence is that it is possible to consider high order schemes which generalize that algorithm by preserving the steady-growth invariance of the solutions with respect to the discretization process. Numerical examples show the efficiency and accuracy of the proposed methods when applied to the classical Ramsey growth model.
An improvement of a mathematical model of the galvanic iron corrosion, previously presented by one of the authors, is here proposed. The iron(III)-hydroxide formation is, now, considered in addition to the redox reaction. The PDE system, assembled on the basis of the fundamental holding electro-chemistry laws, is numerically solved by a locally refined FD method. For verification purpose we have assembled an experimental galvanic cell; in the present work, we report two tests cases, with acidic and neutral electrolitical solution, where the computed electric potential compares well with the measured experimental one.
convection-diffusion equations
electrochemistry
iron