Abstract
We formulate a new class of optimization-based methods for data transfer (remap) of a scalar conserved quantity between two close meshes with the same connectivity. We present the methods in the context of the remap of a mass density field, which preserves global mass (the integral of the density over the computational domain). The key idea is to formulate remap as a global inequality-constrained optimization problem for mass fluxes between neighboring cells. The objective is to minimize the discrepancy between these fluxes and the given high-order target mass fluxes, subject to constraints that enforce physically motivated bounds on the associated primitive variable. In so doing, we separate accuracy considerations, handled by the objective functional, from the enforcement of physical bounds, handled by the constraints. The resulting second-order, conservative, and bound-preserving optimization-based remap (OBR) formulation is applicable to general, unstructured, heterogeneous grids. Under some weak requirements on grid proximity we prove that the OBR algorithm preserves linear fields in one, two and three dimensions. The chapter also examines connections between the OBR and the flux-corrected remap (FCR), which can be interpreted as a modified version of OBR (M-OBR), with the same objective but a smaller feasible set. The feasible set for M-OBR (FCR) is given by simple box constraints derived by using a “worst-case” scenario approach, which may result in loss of linearity preservation and ultimately accuracy for some grid motions. The optimality of the OBR solution means that, given a set of target fluxes and a distance measure, OBR finds the best possible approximations of these fluxes with respect to this measure, which also satisfy the physically motivated bounds. In this sense, OBR can serve as a natural benchmark for evaluating the accuracy of existing and future numerical methods for data transfer with respect to a given class of flux reconstruction methods and flux distance measures. In this context, we perform numerical comparisons between OBR, FCR and iFCR (a version of FCR which utilizes an iterative procedure to enhance the accuracy of FCR numerical fluxes).
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Notes
- 1.
Typically, in a continuous rezone ALE the rezoned grid is close to the Lagrangian but has better geometric quality.
- 2.
In this chapter, we use the set-relational definitions and the corresponding geometric interpretations of ⊂, ⊆, ∪, ∩, ∖ and ∈ interchangeably. Their meaning will be clear from the context. In particular, relations between entities defined on \(\widetilde {K}_{h}({\varOmega})\) and those defined on K h (Ω) only make sense when interpreted geometrically relative to the common domain Ω.
- 3.
In practice, this also means that the integrals in (17) should be approximated by quadratures that are exact for linear functions.
- 4.
The Euclidean distance is used for simplicity. The objective can be defined using any valid distance function (or, equivalently, norm).
- 5.
This guarantees that taking min and max in (24) is well-defined. Otherwise, the correct statement of this result should involve inf and sup.
- 6.
For a complete match with FCR we can set all free variables to 1.
- 7.
Because side nodes can move in different directions swept regions are not simple extrusions of the sides, which can complicate the computation of integrals. Using Green’s theorem, integrals of polynomials over swept regions can be replaced by integrals of higher-degree polynomials over the (lower-dimensional) boundaries of these regions, see Margolin and Shashkov [17], Dukowicz and Kodis [8]. This provides an efficient way to compute the fluxes, regardless of the shape of the swept regions.
References
Bell, J., Berger, M., Saltzman, J., Welcome, M.: Three-dimensional adaptive mesh refinement for hyperbolic conservation laws. SIAM J. Sci. Comput. 15(1), 127–138 (1994)
Berger, M., Murman, S.M., Aftosmis, M.J.: Analysis of slope limiters on irregular grids. In: Proceedings of the 43rd AIAA Aerospace Sciences Meeting, Reno, NV (2005). AIAA2005-0490
Berger, M.J., Colella, P.: Local adaptive mesh refinement for shock hydrodynamics. J. Comput. Phys. 82, 64–84 (1989)
Bochev, P., Day, D.: Analysis and computation of a least-squares method for consistent mesh tying. J. Comp. Appl. Math 218, 21–33 (2008)
Bochev, P., Ridzal, D., Scovazzi, G., Shashkov, M.: Formulation, analysis and numerical study of an optimization-based conservative interpolation (remap) of scalar fields for arbitrary Lagrangian-Eulerian methods. J. Comput. Phys. 230(12), 5199–5225 (2011)
Carey, G., Bicken, G., Carey, V., Berger, C., Sanchez, J.: Locally constrained projections on grids. Int. J. Numer. Methods Eng. 50, 549–577 (2001)
Coleman, T.F., Li, Y.: A reflective newton method for minimizing a quadratic function subject to bounds on some of the variables. SIAM J. Optim. 6(4), 1040–1058 (1996)
Dukowicz, J.K., Kodis, J.W.: Accurate conservative remapping (rezoning) for arbitrary Lagrangian-Eulerian computations. SIAM J. Sci. Stat. Comput. 8, 305–321 (1987)
Hirt, C., Amsden, A., Cook, J.: An arbitrary Lagrangian-Eulerian computing method for all flow speeds. J. Comput. Phys. 14, 227–253 (1974)
Jones, P.W.: First- and second-order conservative remapping schemes for grids in spherical coordinates. Mon. Weather Rev. 127(9), 2204–2210 (1999)
Kucharik, M., Shashkov, M., Wendroff, B.: An efficient linearity-and-bound-preserving remapping method. J. Comput. Phys. 188(2), 462–471 (2003)
Kuzmin, D., Löhner, R., Turek, S. (eds.): Flux-Corrected Transport: Principles, Algorithms and Applications. Springer, Berlin (2005)
Laursen, T., Heinstein, M.: A three dimensional surface-to-surface projection algorithm for non-coincident domains. Commun. Numer. Methods Eng. 19, 421–432 (2003)
Liska, R., Shashkov, M., Váchal, P., Wendroff, B.: Optimization-based synchronized flux-corrected conservative interpolation (remapping) of mass and momentum for arbitrary Lagrangian-Eulerian methods. J. Comput. Phys. 229(5), 1467–1497 (2010)
Loubere, R., Shashkov, M.J.: A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. J. Comput. Phys. 209(1), 105–138 (2005)
Loubere, R., Staley, M., Wendroff, B.: The repair paradigm: New algorithms and applications to compressible flow. J. Comput. Phys. 211(2), 385–404 (2006)
Margolin, L.G., Shashkov, M.: Second-order sign-preserving conservative interpolation (remapping) on general grids. J. Comput. Phys. 184(1), 266–298 (2003)
Margolin, L.G., Shashkov, M.: Remapping, recovery and repair on a staggered grid. Comput. Methods Appl. Mech. Eng. 193(39–41), 4139–4155 (2004)
Miller, D.S., Burton, D.E., Oliviera, J.S.: Efficient second order remapping on arbitrary two dimensional meshes. Technical Report UCRL-ID-123530, Lawrence Livermore National Laboratory (1996)
Rider, W.J., Kothe, D.B.: Constrained minimization for monotonic reconstruction. In: Proceedings of the 13th AIAA Computational Fluid Dynamics Conference, Snowmass, CO, pp. 955–964 (1997). AIAA 97-2036
Rockafellar, R.T.: Convex Analysis. Princeton University Press, Princeton (1970)
Schär, C., Smolarkiewicz, P.K.: A synchronous and iterative flux-correction formalism for coupled transport equations. J. Comput. Phys. 128(1), 101–120 (1996)
Swartz, B.: Good neighborhoods for multidimensional Van Leer limiting. J. Comput. Phys. 154(1), 237–241 (1999)
Acknowledgements
All authors acknowledge funding by the DOE Office of Science Advanced Scientific Computing Research (ASCR) Program. PB, DR and GS also acknowledge funding by the NNSA Climate Modeling and Carbon Measurement Project. DR and MS also acknowledge funding by the Advanced Simulation & Computing (ASC) Program.
Our colleagues Dmitri Kuzmin, Richard Liska, Kara Peterson, John Shadid, Pavel Váchal and Joseph Young provided many comments and valuable insights that helped improve this work.
Sandia National Laboratories is a multi-program laboratory operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000. The work of MS was carried out under the auspices of the National Nuclear Security Administration of the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396.
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2012 Springer Science+Business Media Dordrecht
About this chapter
Cite this chapter
Bochev, P., Ridzal, D., Scovazzi, G., Shashkov, M. (2012). Constrained-Optimization Based Data Transfer. In: Kuzmin, D., Löhner, R., Turek, S. (eds) Flux-Corrected Transport. Scientific Computation. Springer, Dordrecht. https://doi.org/10.1007/978-94-007-4038-9_10
Download citation
DOI: https://doi.org/10.1007/978-94-007-4038-9_10
Publisher Name: Springer, Dordrecht
Print ISBN: 978-94-007-4037-2
Online ISBN: 978-94-007-4038-9
eBook Packages: Physics and AstronomyPhysics and Astronomy (R0)