Sökning: onr:"swepub:oai:DiVA.org:liu-76798" >
Monotonicity recove...
Monotonicity recovering and accuracy preserving optimization methods for postprocessing finite element solutions
-
- Burdakov, Oleg (författare)
- Linköpings universitet,Optimeringslära,Tekniska högskolan
-
- Kapyrin, Ivan (författare)
- Russian Academy of Science
-
- Vassilevski, Yuri (författare)
- Russian Academy of Science
-
(creator_code:org_t)
- Elsevier, 2012
- 2012
- Engelska.
-
Ingår i: Journal of Computational Physics. - : Elsevier. - 0021-9991 .- 1090-2716. ; 231:8, s. 3126-3142
- Relaterad länk:
-
https://liu.diva-por... (primary) (Raw object)
-
visa fler...
-
http://liu.diva-port...
-
https://urn.kb.se/re...
-
https://doi.org/10.1...
-
visa färre...
Abstract
Ämnesord
Stäng
- We suggest here a least-change correction to available finite element (FE) solution. This postprocessing procedure is aimed at recovering the monotonicity and some other important properties that may not be exhibited by the FE solution. Although our approach is presented for FEs, it admits natural extension to other numerical schemes, such as finite differences and finite volumes. For the postprocessing, a priori information about the monotonicity is assumed to be available, either for the whole domain or for a subdomain where the lost monotonicity is to be recovered. The obvious requirement is that such information is to be obtained without involving the exact solution, e.g. from expected symmetries of this solution. less thanbrgreater than less thanbrgreater thanThe postprocessing is based on solving a monotonic regression problem with some extra constraints. One of them is a linear equality-type constraint that models the conservativity requirement. The other ones are box-type constraints, and they originate from the discrete maximum principle. The resulting postprocessing problem is a large scale quadratic optimization problem. It is proved that the postprocessed FE solution preserves the accuracy of the discrete FE approximation. less thanbrgreater than less thanbrgreater thanWe introduce an algorithm for solving the postprocessing problem. It can be viewed as a dual ascent method based on the Lagrangian relaxation of the equality constraint. We justify theoretically its correctness. Its efficiency is demonstrated by the presented results of numerical experiments.
Nyckelord
- Finite element solution
- Accuracy analysis
- Constrained monotonic regression
- Large scale quadratic optimization
- Lagrangian relaxation
- Dual ascent method
- MATHEMATICS
- MATEMATIK
Publikations- och innehållstyp
- ref (ämneskategori)
- art (ämneskategori)
Hitta via bibliotek
Till lärosätets databas