Risultati di regolarità, stabilità, instabilità e biforcazione per un problema a frontiera libera
In this work we derive a posteriori error estimates based on equations residuals for the heat equation with discontinuous diffusivity coefficients. The estimates are based on a fully discrete scheme based on conforming finite elements in each time slab and on the A-stable -scheme with . Following remarks of [Picasso, Comput. Methods Appl. Mech. Engrg. 167 (1998) 223–237; Verfürth, Calcolo 40 (2003) 195–212] it is easy to identify a time-discretization error-estimator and a space-discretization...
In this work we derive a posteriori error estimates based on equations residuals for the heat equation with discontinuous diffusivity coefficients. The estimates are based on a fully discrete scheme based on conforming finite elements in each time slab and on the A-stable θ-scheme with 1/2 ≤ θ ≤ 1. Following remarks of [Picasso, Comput. Methods Appl. Mech. Engrg. 167 (1998) 223–237; Verfürth, Calcolo40 (2003) 195–212] it is easy to identify a time-discretization error-estimator and a space-discretization...
This article discusses the numerical approximation of time dependent Ginzburg-Landau equations. Optimal error estimates which are robust with respect to a large Ginzburg-Landau parameter are established for a semi-discrete in time and a fully discrete approximation scheme. The proofs rely on an asymptotic expansion of the exact solution and a stability result for degree-one Ginzburg-Landau vortices. The error bounds prove that degree-one vortices can be approximated robustly while unstable higher...
This article discusses the numerical approximation of time dependent Ginzburg-Landau equations. Optimal error estimates which are robust with respect to a large Ginzburg-Landau parameter are established for a semi-discrete in time and a fully discrete approximation scheme. The proofs rely on an asymptotic expansion of the exact solution and a stability result for degree-one Ginzburg-Landau vortices. The error bounds prove that degree-one vortices can be approximated robustly while unstable higher...
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations. The constant in the corresponding abstract...
An abstract framework for constructing stable decompositions of the spaces corresponding to general symmetric positive definite problems into “local” subspaces and a global “coarse” space is developed. Particular applications of this abstract framework include practically important problems in porous media applications such as: the scalar elliptic (pressure) equation and the stream function formulation of its mixed form, Stokes’ and Brinkman’s equations....
Motivated by downside risk minimization on the wealth process in an incomplete market model, we have studied in the recent work the asymptotic behavior as time horizon T → ∞ of the minimizing probability that the empirical mean of a controlled semi-martingale falls below a certain level on the time horizon T. This asymptotic behavior relates to a risk-sensitive stochastic control problem in the risk-averse case. Indeed, we obtained an expression of the decay rate of the probability by the Legendre...
Singularly perturbed problems often yield solutions with strong directional features, e.g. with boundary layers. Such anisotropic solutions lend themselves to adapted, anisotropic discretizations. The quality of the corresponding numerical solution is a key issue in any computational simulation. To this end we present a new robust error estimator for a singularly perturbed reaction–diffusion problem. In contrast to conventional estimators, our proposal is suitable for anisotropic finite element...
Singularly perturbed problems often yield solutions with strong directional features, e.g. with boundary layers. Such anisotropic solutions lend themselves to adapted, anisotropic discretizations. The quality of the corresponding numerical solution is a key issue in any computational simulation. To this end we present a new robust error estimator for a singularly perturbed reaction-diffusion problem. In contrast to conventional estimators, our proposal is suitable for anisotropic finite element...
We discuss a family of discontinuous Petrov–Galerkin (DPG) schemes for quite general partial differential operators. The starting point of our analysis is the DPG method introduced by [Demkowicz et al., SIAM J. Numer. Anal. 49 (2011) 1788–1809; Zitelli et al., J. Comput. Phys. 230 (2011) 2406–2432]. This discretization results in a sparse positive definite linear algebraic system which can be obtained from a saddle point problem by an element-wise Schur complement reduction applied to the test space....
While a large amount of papers are dealing with robust multilevel methods and algorithms for linear FEM elliptic systems, the related higher order FEM problems are much less studied. Moreover, we know that the standard hierarchical basis two-level splittings deteriorate for strongly anisotropic problems. A first robust multilevel preconditioner for higher order FEM systems obtained after discretizations of elliptic problems with an anisotropic diffusion tensor is presented in this paper. We study...
In this article we derive a macroscopic model for the time evolution of root density, starting from a discrete mesh of roots, using homogenization techniques. In the microscopic model each root grows vertically according to an ordinary differential equation. The roots growth rates depend on the spatial distribution of nutrient in the soil, which also evolves in time, leading to a fully coupled non-linear problem. We derive an effective partial differential equation for the root tip surface and for...
In this article we derive a macroscopic model for the time evolution of root density, starting from a discrete mesh of roots, using homogenization techniques. In the microscopic model each root grows vertically according to an ordinary differential equation. The roots growth rates depend on the spatial distribution of nutrient in the soil, which also evolves in time, leading to a fully coupled non-linear problem. We derive an effective partial differential equation for the root tip surface and for...