Stability of microstructure for tetragonal to monoclinic martensitic transformations
We give an analysis of the stability and uniqueness of the simply laminated microstructure for all three tetragonal to monoclinic martensitic transformations. The energy density for tetragonal to monoclinic transformations has four rotationally invariant wells since the transformation has four variants. One of these tetragonal to monoclinic martensitic transformations corresponds to the shearing of the rectangular side, one corresponds to the shearing of the square base, and one corresponds to...
The magnetization of a ferromagnetic sample solves a non-convex variational problem, where its relaxation by convexifying the energy density resolves relevant macroscopic information. The numerical analysis of the relaxed model has to deal with a constrained convex but degenerated, nonlocal energy functional in mixed formulation for magnetic potential and magnetization . In [C. Carstensen and A. Prohl, Numer. Math. 90 (2001) 65–99], the conforming -element in spatial dimensions is shown to...
The magnetization of a ferromagnetic sample solves a non-convex variational problem, where its relaxation by convexifying the energy density resolves relevant macroscopic information. The numerical analysis of the relaxed model has to deal with a constrained convex but degenerated, nonlocal energy functional in mixed formulation for magnetic potential u and magnetization m. In [C. Carstensen and A. Prohl, Numer. Math.90 (2001) 65–99], the conforming P1 - (P0)d-element in d=2,3 spatial dimensions...
In dimension one it is proved that the solution to a total variation-regularized least-squares problem is always a function which is "constant almost everywhere" , provided that the data are in a certain sense outside the range of the operator to be inverted. A similar, but weaker result is derived in dimension two.
In this paper, we discuss the numerical simulation for a class of constrained optimal control problems governed by integral equations. The Galerkin method is used for the approximation of the problem. A priori error estimates and a superconvergence analysis for the approximation scheme are presented. Based on the results of the superconvergence analysis, a recovery type a posteriori error estimator is provided, which can be used for adaptive mesh refinement.
In this paper we are concerned with finite element approximations to the evaluation of American options. First, following W. Allegretto etc., SIAM J. Numer. Anal. 39 (2001), 834–857, we introduce a novel practical approach to the discussed problem, which involves the exact reformulation of the original problem and the implementation of the numerical solution over a very small region so that this algorithm is very rapid and highly accurate. Secondly by means of a superapproximation and interpolation...
The powerful Hamilton-Jacobi theory is used for constructing regularizations and error estimates for optimal design problems. The constructed Pontryagin method is a simple and general method for optimal design and reconstruction: the first, analytical, step is to regularize the hamiltonian; next the solution to its stationary hamiltonian system, a nonlinear partial differential equation, is computed with the Newton method. The method is efficient for designs where the hamiltonian function can be...
The powerful Hamilton-Jacobi theory is used for constructing regularizations and error estimates for optimal design problems. The constructed Pontryagin method is a simple and general method for optimal design and reconstruction: the first, analytical, step is to regularize the Hamiltonian; next the solution to its stationary Hamiltonian system, a nonlinear partial differential equation, is computed with the Newton method. The method is efficient for designs where the Hamiltonian function...
Some dynamical systems are characterized by more than one time-scale, e.g. two well separated time-scales are typical for quasiperiodic systems. The aim of this paper is to show how singular perturbation methods based on the slow-fast decomposition can serve for an enhanced parameter estimation when the slowly changing features are rigorously treated. Although the ultimate goal is to reduce the standard error for the estimated parameters, here we test two methods for numerical approximations of...
For the symmetric Pareto Eigenvalue Complementarity Problem (EiCP), by reformulating it as a constrained optimization problem on a differentiable Rayleigh quotient function, we present a class of descent methods and prove their convergence. The main features include: using nonlinear complementarity functions (NCP functions) and Rayleigh quotient gradient as the descent direction, and determining the step size with exact linear search. In addition, these algorithms are further extended to solve the...
Searching for the optimal partitioning of a domain leads to the use of the adjoint method in topological asymptotic expansions to know the influence of a domain perturbation on a cost function. Our approach works by restricting to local subproblems containing the perturbation and outperforms the adjoint method by providing approximations of higher order. It is a universal tool, easily adapted to different kinds of real problems and does not need the fundamental solution of the problem; furthermore...
Searching for the optimal partitioning of a domain leads to the use of the adjoint method in topological asymptotic expansions to know the influence of a domain perturbation on a cost function. Our approach works by restricting to local subproblems containing the perturbation and outperforms the adjoint method by providing approximations of higher order. It is a universal tool, easily adapted to different kinds of real problems and does not need...