{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,12,31]],"date-time":"2025-12-31T00:41:10Z","timestamp":1767141670096,"version":"build-2238731810"},"reference-count":37,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2021,9,12]],"date-time":"2021-09-12T00:00:00Z","timestamp":1631404800000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2021,9,12]],"date-time":"2021-09-12T00:00:00Z","timestamp":1631404800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"DOI":"10.13039\/501100003958","name":"Stichting voor de Technische Wetenschappen","doi-asserted-by":"publisher","award":["Perspectief P15-36"],"award-info":[{"award-number":["Perspectief P15-36"]}],"id":[{"id":"10.13039\/501100003958","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["J Sci Comput"],"published-print":{"date-parts":[[2021,10]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>\n                    Liouville\u2019s equation on phase space in geometrical optics describes the evolution of an energy distribution through an optical system, which is discontinuous across optical interfaces. The discontinuous Galerkin spectral element method is conservative and can achieve higher order of convergence locally, making it a suitable method for this equation. When dealing with optical interfaces in phase space, non-local boundary conditions arise. Besides being a difficulty in itself, these non-local boundary conditions must also satisfy energy conservation constraints. To this end, we introduce an energy conservative treatment of optical interfaces. Numerical experiments are performed to prove that the method obeys energy conservation. Furthermore, the method is compared to the industry standard ray tracing. The numerical experiments show that the discontinuous Galerkin spectral element method outperforms ray tracing by reducing the computation time by up to three orders of magnitude for an error of\n                    <jats:inline-formula>\n                      <jats:alternatives>\n                        <jats:tex-math>$$10^{-6}$$<\/jats:tex-math>\n                        <mml:math xmlns:mml=\"http:\/\/www.w3.org\/1998\/Math\/MathML\">\n                          <mml:msup>\n                            <mml:mn>10<\/mml:mn>\n                            <mml:mrow>\n                              <mml:mo>-<\/mml:mo>\n                              <mml:mn>6<\/mml:mn>\n                            <\/mml:mrow>\n                          <\/mml:msup>\n                        <\/mml:math>\n                      <\/jats:alternatives>\n                    <\/jats:inline-formula>\n                    .\n                  <\/jats:p>","DOI":"10.1007\/s10915-021-01612-x","type":"journal-article","created":{"date-parts":[[2021,9,12]],"date-time":"2021-09-12T10:04:01Z","timestamp":1631441041000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":2,"title":["An Energy Conservative hp-method for Liouville\u2019s Equation of Geometrical Optics"],"prefix":"10.1007","volume":"89","author":[{"ORCID":"https:\/\/orcid.org\/0000-0003-0769-0502","authenticated-orcid":false,"given":"R. A. M.","family":"van Gestel","sequence":"first","affiliation":[]},{"given":"M. J. H.","family":"Anthonissen","sequence":"additional","affiliation":[]},{"given":"J. H. M. ten Thije","family":"Boonkkamp","sequence":"additional","affiliation":[]},{"given":"W. L.","family":"IJzerman","sequence":"additional","affiliation":[]}],"member":"297","published-online":{"date-parts":[[2021,9,12]]},"reference":[{"key":"1612_CR1","volume-title":"Mathematical Methods of Classical Mechanics","author":"VI Arnold","year":"2013","unstructured":"Arnold, V.I.: Mathematical Methods of Classical Mechanics, vol. 60. Springer Science & Business Media, Berlin (2013)"},{"issue":"23","key":"1612_CR2","doi-asserted-by":"publisher","first-page":"27797","DOI":"10.1364\/OE.22.027797","volume":"22","author":"M Bahrami","year":"2014","unstructured":"Bahrami, M., Goncharov, A.V.: Geometry-invariant GRIN lens: finite ray tracing. Opt. Express 22(23), 27797\u201327810 (2014)","journal-title":"Opt. Express"},{"issue":"1\u20132","key":"1612_CR3","doi-asserted-by":"publisher","first-page":"65","DOI":"10.1016\/0168-9274(89)90055-X","volume":"6","author":"C Bernardi","year":"1989","unstructured":"Bernardi, C., Maday, Y., Landriani, G.S.: Noncorming matching conditions for coupling spectral and finite element methods. Appl. Numer. Math. 6(1\u20132), 65\u201384 (1989)","journal-title":"Appl. Numer. Math."},{"key":"1612_CR4","doi-asserted-by":"crossref","unstructured":"Bernardi, C., Maday, Y., Patera, A. T.: Domain Decomposition by the Mortar Element Method. In: Asymptotic and Numerical Methods for Partial Differential Equations with Critical Parameters, pages 269\u2013286. Springer (1993)","DOI":"10.1007\/978-94-011-1810-1_17"},{"issue":"3","key":"1612_CR5","doi-asserted-by":"publisher","first-page":"1801","DOI":"10.1137\/110828010","volume":"50","author":"T Bui-Thanh","year":"2012","unstructured":"Bui-Thanh, T., Ghattas, O.: Analysis of an hp-nonconforming discontinuous Galerkin spectral element method for wave propagation. SIAM J. Numer. Anal. 50(3), 1801\u20131826 (2012)","journal-title":"SIAM J. Numer. Anal."},{"key":"1612_CR6","unstructured":"Carpenter, M. H., Kennedy, C. A.: Fourth-order 2N-storage Runge-Kutta schemes. NASA TM 109112 (1994)"},{"key":"1612_CR7","doi-asserted-by":"publisher","first-page":"10905","DOI":"10.1016\/j.jcp.2019.109095","volume":"403","author":"N Chalmers","year":"2020","unstructured":"Chalmers, N., Krivodonova, L.: A robust CFL condition for the discontinuous Galerkin method on triangular meshes. J. Comput. Phys. 403, 10905 (2020)","journal-title":"J. Comput. Phys."},{"issue":"6","key":"1612_CR8","doi-asserted-by":"publisher","first-page":"A2395","DOI":"10.1137\/16M1089198","volume":"39","author":"J Chan","year":"2017","unstructured":"Chan, J., Hewett, R.J., Warburton, T.: Weight-adjusted discontinuous Galerkin methods: curvilinear meshes. SIAM J. Sci. Comput. 39(6), A2395\u2013A2421 (2017)","journal-title":"SIAM J. Sci. Comput."},{"key":"1612_CR9","doi-asserted-by":"publisher","DOI":"10.1201\/b18785","volume-title":"Introduction to Nonimaging Optics","author":"J Chaves","year":"2017","unstructured":"Chaves, J.: Introduction to Nonimaging Optics. CRC Press, Boca Raton (2017)"},{"issue":"26","key":"1612_CR10","doi-asserted-by":"publisher","first-page":"13014","DOI":"10.1364\/OE.14.013014","volume":"14","author":"A Cvetkovic","year":"2006","unstructured":"Cvetkovic, A., Dross, O., Chaves, J., Benitez, P., Mi\u00f1ano, J.C., Mohedano, R.: Etendue-preserving mixing and projection optics for high-luminance LEDs, applied to automotive headlamps. Opt. Express 14(26), 13014 (2006)","journal-title":"Opt. Express"},{"key":"1612_CR11","unstructured":"Filosa, C.: Phase Space Ray Tracing for Illumination Optics. PhD thesis, Eindhoven University of Technology (2018)"},{"issue":"13","key":"1612_CR12","doi-asserted-by":"publisher","first-page":"3599","DOI":"10.1364\/AO.55.003599","volume":"55","author":"C Filosa","year":"2016","unstructured":"Filosa, C., ten Thije Boonkkamp, J.H.M., Ijzerman, W.L.: Ray tracing method in phase space for two-dimensional optical systems. Appl. Op. 55(13), 3599\u20133606 (2016)","journal-title":"Appl. Op."},{"issue":"4","key":"1612_CR13","doi-asserted-by":"publisher","first-page":"135","DOI":"10.13189\/ms.2017.050401","volume":"5","author":"C Filosa","year":"2017","unstructured":"Filosa, C., ten Thije Boonkkamp, J.H.M., Ijzerman, W.L.: Phase space ray tracing for a two-dimensional parabolic reflector. Math. Stat. 5(4), 135\u2013142 (2017)","journal-title":"Math. Stat."},{"key":"1612_CR14","volume-title":"An Introduction to Ray Tracing","author":"AS Glassner","year":"1989","unstructured":"Glassner, A.S.: An Introduction to Ray Tracing. Elsevier, Amsterdam (1989)"},{"key":"1612_CR15","doi-asserted-by":"crossref","unstructured":"Griffiths, D. J.: Introduction to Electrodynamics (2005)","DOI":"10.1016\/B978-1-85573-953-6.50026-X"},{"key":"1612_CR16","volume-title":"Optics","author":"E Hecht","year":"2002","unstructured":"Hecht, E., et al.: Optics, vol. 4. Addison Wesley, San Francisco (2002)"},{"issue":"3","key":"1612_CR17","doi-asserted-by":"publisher","first-page":"031304","DOI":"10.1117\/1.OE.53.3.031304","volume":"53","author":"AM Herkommer","year":"2013","unstructured":"Herkommer, A.M.: Phase space optics: an alternate approach to freeform optical systems. Opt. Eng. 53(3), 031304 (2013)","journal-title":"Opt. Eng."},{"key":"1612_CR18","volume-title":"Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications","author":"JS Hesthaven","year":"2007","unstructured":"Hesthaven, J.S., Warburton, T.: Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications. Springer Science & Business Media, Berlin (2007)"},{"issue":"5","key":"1612_CR19","doi-asserted-by":"publisher","first-page":"1763","DOI":"10.1016\/j.jcp.2009.11.006","volume":"229","author":"DI Ketcheson","year":"2010","unstructured":"Ketcheson, D.I.: Runge-Kutta methods with minimum storage implementations. J. Comput. Phys. 229(5), 1763\u20131773 (2010)","journal-title":"J. Comput. Phys."},{"issue":"2","key":"1612_CR20","doi-asserted-by":"publisher","first-page":"475","DOI":"10.1006\/jcph.1996.0225","volume":"128","author":"DA Kopriva","year":"1996","unstructured":"Kopriva, D.A.: A conservative staggered-grid Chebyshev multidomain method for compressible flows. II. A semi-structured method. J. Comput. Phys. 128(2), 475\u2013488 (1996)","journal-title":"J. Comput. Phys."},{"key":"1612_CR21","doi-asserted-by":"publisher","DOI":"10.1007\/978-90-481-2261-5","volume-title":"Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers","author":"DA Kopriva","year":"2009","unstructured":"Kopriva, D.A.: Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers. Springer Science & Business Media, Berlin (2009)"},{"key":"1612_CR22","first-page":"274","volume":"272","author":"DA Kopriva","year":"2016","unstructured":"Kopriva, D.A., Gassner, G.J.: Geometry effects in nodal discontinuous Galerkin methods on curved elements that are provably stable. Appl. Math. Comput. 272, 274\u2013290 (2016)","journal-title":"Appl. Math. Comput."},{"issue":"1","key":"1612_CR23","doi-asserted-by":"publisher","first-page":"105","DOI":"10.1002\/nme.394","volume":"53","author":"DA Kopriva","year":"2002","unstructured":"Kopriva, D.A., Woodruff, S.L., Hussaini, M.Y.: Computation of electromagnetic scattering with a non-conforming discontinuous spectral element method. Int. J. Numer. Meth. Eng. 53(1), 105\u2013122 (2002)","journal-title":"Int. J. Numer. Meth. Eng."},{"key":"1612_CR24","doi-asserted-by":"publisher","DOI":"10.1007\/978-3-319-03425-6","volume-title":"Introduction to Quasi-Monte Carlo Integration and Applications","author":"G Leobacher","year":"2014","unstructured":"Leobacher, G., Pillichshammer, F.: Introduction to Quasi-Monte Carlo Integration and Applications. Springer, Berlin (2014)"},{"key":"1612_CR25","volume-title":"Introduction to Radiometry and Photometry","author":"WR McCluney","year":"2014","unstructured":"McCluney, W.R.: Introduction to Radiometry and Photometry. Artech House, New York (2014)"},{"issue":"5","key":"1612_CR26","doi-asserted-by":"publisher","first-page":"368","DOI":"10.1119\/1.1969512","volume":"31","author":"FE Nicodemus","year":"1963","unstructured":"Nicodemus, F.E.: Radiance. Am. J. Phys. 31(5), 368\u2013377 (1963)","journal-title":"Am. J. Phys."},{"issue":"3","key":"1612_CR27","doi-asserted-by":"publisher","first-page":"411","DOI":"10.1364\/JOSAA.378829","volume":"37","author":"H Ohno","year":"2020","unstructured":"Ohno, H.: Symplectic ray tracing based on Hamiltonian optics in gradient-index media. JOSA A 37(3), 411\u2013416 (2020)","journal-title":"JOSA A"},{"key":"1612_CR28","first-page":"15","volume":"5186","author":"DG Pelka","year":"2003","unstructured":"Pelka, D.G., Patel, K.: An overview of LED applications for general illumination. Des. Eff. Illum. Syst. 5186, 15\u201326 (2003)","journal-title":"Des. Eff. Illum. Syst."},{"issue":"6","key":"1612_CR29","doi-asserted-by":"publisher","first-page":"065103","DOI":"10.1117\/1.OE.56.6.065103","volume":"56","author":"D Rausch","year":"2017","unstructured":"Rausch, D., Rommel, M., Herkommer, A.M., Talpur, T.: Illumination design for extended sources based on phase space mapping. Op. Eng. 56(6), 065103 (2017)","journal-title":"Op. Eng."},{"issue":"12","key":"1612_CR30","doi-asserted-by":"publisher","first-page":"4657","DOI":"10.1016\/j.jcp.2011.02.040","volume":"230","author":"T Toulorge","year":"2011","unstructured":"Toulorge, T., Desmet, W.: CFL conditions for Runge-Kutta discontinuous Galerkin methods on triangular grids. J. Comput. Phys. 230(12), 4657\u20134678 (2011)","journal-title":"J. Comput. Phys."},{"key":"1612_CR31","unstructured":"van Lith, B. S.: Principles of Computational Illumination Optics. PhD thesis, Eindhoven University of Technology (2017)"},{"key":"1612_CR32","first-page":"100030","volume":"3","author":"BS van Lith","year":"2019","unstructured":"van Lith, B.S., ten Thije Boonkkamp, J.H.M., IJzerman, W.L.: Active flux schemes on moving meshes with applications to geometric optics. J. Comput. Phys. X 3, 100030 (2019)","journal-title":"J. Comput. Phys. X"},{"issue":"2","key":"1612_CR33","doi-asserted-by":"publisher","first-page":"739","DOI":"10.1007\/s10915-015-0157-6","volume":"68","author":"BS van Lith","year":"2016","unstructured":"van Lith, B.S., ten Thije Boonkkamp, J.H.M., IJzerman, W.L., Tukker, T.W.: A novel scheme for Liouvilles equation with a discontinuous Hamiltonian and applications to geometrical optics. J. Sci. Comput. 68(2), 739\u2013771 (2016)","journal-title":"J. Sci. Comput."},{"issue":"1","key":"1612_CR34","doi-asserted-by":"publisher","first-page":"48","DOI":"10.1016\/0021-9991(80)90033-9","volume":"35","author":"J Williamson","year":"1980","unstructured":"Williamson, J.: Low-storage Runge-Kutta schemes. J. Comput. Phys. 35(1), 48\u201356 (1980)","journal-title":"J. Comput. Phys."},{"key":"1612_CR35","volume-title":"Geometric Optics on Phase Space","author":"KB Wolf","year":"2004","unstructured":"Wolf, K.B.: Geometric Optics on Phase Space. Springer Science & Business Media, Berlin (2004)"},{"key":"1612_CR36","doi-asserted-by":"publisher","first-page":"262","DOI":"10.1016\/j.optlastec.2012.05.040","volume":"45","author":"X Zhu","year":"2013","unstructured":"Zhu, X., Zhu, Q., Wu, H., Chen, C.: Optical design of LED-based automotive headlamps. Op. Laser Technol. 45, 262\u2013266 (2013)","journal-title":"Op. Laser Technol."},{"issue":"2","key":"1612_CR37","doi-asserted-by":"publisher","first-page":"227","DOI":"10.1016\/S0168-9274(98)00129-9","volume":"31","author":"DW Zingg","year":"1999","unstructured":"Zingg, D.W., Chisholm, T.T.: Runge-Kutta methods for linear ordinary differential equations. Appl. Numer. Math. 31(2), 227\u2013238 (1999)","journal-title":"Appl. Numer. Math."}],"updated-by":[{"DOI":"10.1007\/s10915-021-01754-y","type":"correction","label":"Correction","source":"publisher","updated":{"date-parts":[[2022,1,24]],"date-time":"2022-01-24T00:00:00Z","timestamp":1642982400000}}],"container-title":["Journal of Scientific Computing"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s10915-021-01612-x.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/article\/10.1007\/s10915-021-01612-x\/fulltext.html","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/link.springer.com\/content\/pdf\/10.1007\/s10915-021-01612-x.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,1,24]],"date-time":"2022-01-24T02:17:16Z","timestamp":1642990636000},"score":1,"resource":{"primary":{"URL":"https:\/\/link.springer.com\/10.1007\/s10915-021-01612-x"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,9,12]]},"references-count":37,"journal-issue":{"issue":"1","published-print":{"date-parts":[[2021,10]]}},"alternative-id":["1612"],"URL":"https:\/\/doi.org\/10.1007\/s10915-021-01612-x","relation":{},"ISSN":["0885-7474","1573-7691"],"issn-type":[{"value":"0885-7474","type":"print"},{"value":"1573-7691","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,9,12]]},"assertion":[{"value":"28 May 2020","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"9 July 2021","order":2,"name":"revised","label":"Revised","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"11 July 2021","order":3,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"12 September 2021","order":4,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"24 January 2022","order":5,"name":"change_date","label":"Change Date","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"Correction","order":6,"name":"change_type","label":"Change Type","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"A Correction to this paper has been published:","order":7,"name":"change_details","label":"Change Details","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"https:\/\/doi.org\/10.1007\/s10915-021-01754-y","URL":"https:\/\/doi.org\/10.1007\/s10915-021-01754-y","order":8,"name":"change_details","label":"Change Details","group":{"name":"ArticleHistory","label":"Article History"}},{"order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Declarations"}},{"value":"The authors declare that they have no conflict of interest.","order":2,"name":"Ethics","group":{"name":"EthicsHeading","label":"Conflict of interest"}},{"value":"Custom code was written in MATLAB and will not be made available online.","order":3,"name":"Ethics","group":{"name":"EthicsHeading","label":"Code Availability"}}],"article-number":"27"}}