[{"id":889,"date":"2026-10-01T15:15:00","date_gmt":"2026-10-01T19:15:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=889"},"modified":"2026-10-02T09:13:06","modified_gmt":"2026-10-02T13:13:06","slug":"achieving-generalization-in-scientific-machine-through-topological-structure","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/achieving-generalization-in-scientific-machine-through-topological-structure\/","title":{"rendered":"Achieving Generalization in Scientific Machine Through Topological Structure"},"content":{"rendered":"\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<p class=\"has-text-align-center has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Presented by: Dr. Nat Trask<\/strong><\/p>\n\n\n\n<p class=\"has-text-align-center has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>University of Pennsylvania<\/strong><\/p>\n\n\n\n<p class=\"has-text-align-center has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Hosted by:<\/p>\n\n\n\n<p class=\"has-text-align-center has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Mathematics in Computation Talk Series<\/strong><\/p>\n\n\n\n<p class=\"has-text-align-center has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract: <\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Within the last year, the artificial intellI community has turned a renewed focus toward engineering and physical sciences, both within the U.S. Department of Energy (DOE) and industry. &nbsp;Many of these efforts hinge upon an assumption that scale alone will be able to achieve predictive scientific models. &nbsp;To date, however, works show that neural physics models rapidly degrade as they are used at the periphery of their training data distribution; this presents a major risk to applying them to design or discovery tasks which, by definition, aim to operate specifically at this frontier. &nbsp;This talk will present recent work in which topological structure in both the representation of geometry and of physics can be used to construct models that are able to extrapolate: both in the sense of predicting physics on geometries or forecasting to future times unseen during training. &nbsp;This structure can be accessed by developing hybrid neural architectures which wrap state-of-the-art transformers in a \u201ctopological scaffold\u201d provided by mixed finite element machinery. &nbsp;By exposing variational structure, provable properties are established regarding stability and accuracy for learned models which is intractable for contemporary \u201cblack-box\u201d architectures. &nbsp;Several examples from both community benchmarks and DOE digital twinning efforts are provided where these architectures are able to achieve state of the art performance, often with orders of magnitude smaller models.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Trask is an associate professor in the Mechanical Engineering and Applied Mechanics department at the University of Pennsylvania. &nbsp;He maintains a joint appointment with Sandia National Laboratories, and is co-director of Scalable, Efficient and Accelerated Causal Reasoning Operators, Graphs and Spikes for Earth and Embedded Systems, a DOE Mathematical Multifaceted Integrated Capability Center developing next generation scientific machine learning architectures. &nbsp;He has worked on a number of application areas spanning multi-scale\/multi-physics in fluids, solids, electromagnetism, materials, and solid-state physics.<strong>&nbsp;<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Host<\/strong>:&nbsp;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Viktor Reshniak,\u00a0<a href=\"mailto:reshniakv@ornl.gov\">reshniakv@ornl.gov<\/a>,\u00a0865-341-4394<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Presented by: Dr. Nat Trask University of Pennsylvania Hosted by: Mathematics in Computation Talk Series Abstract: Within the last year, the artificial intellI community has turned a renewed focus toward engineering and physical sciences, both within the U.S. Department of Energy (DOE) and industry. &nbsp;Many of these efforts hinge upon an assumption that scale alone<\/p>\n","protected":false},"author":110,"featured_media":585,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-889","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/889","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=889"}],"version-history":[{"count":1,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/889\/revisions"}],"predecessor-version":[{"id":890,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/889\/revisions\/890"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/585"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=889"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=889"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=889"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":567,"date":"2026-09-25T10:01:00","date_gmt":"2026-09-25T14:01:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=567"},"modified":"2026-10-08T16:38:37","modified_gmt":"2026-10-08T20:38:37","slug":"fast-high-dimensional-kronecker-tensor-operations-for-sparse-grids","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/fast-high-dimensional-kronecker-tensor-operations-for-sparse-grids\/","title":{"rendered":"Fast High-Dimensional Kronecker Tensor Operations for Sparse Grids"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Miroslav Stoyanov, ORNL<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Sparse Grids (SG) is a common tool for applications that require multidimensional function approximation, such as constructing surrogate models in parameter space or solving high-dimensional partial differential equations.&nbsp; From uncertainty quantification to Kinetic models for plasma physics, the key advantage of SG is the mitigation of the computational cost associated with high-dimensional problems, also known as the curse of dimesionality.&nbsp; The grids are constructed from hierarchical superposition of multidimensional tensors of single dimensional basis functions, and the resulting complex structure allows for fine tuning and adapting to the behavior of a specific function with minimal number of degrees of freedom.&nbsp; The linear algebra operations associated with SG resemble algebraic operations using classical Kronecker matrices and Kronecker theory is in turn a powerful tool for tensor linear algebra.&nbsp; The Kronecker structure allows combinations of like-terms and reduction of the overall computational cost by a factor exponential in the number of dimensions.&nbsp; However, the same like-term relationship is lost in a general sparse context.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Several sufficient conditions are presented that enable recovery of the Kronecker cost reduction in a hierarchical sparse-grid context.&nbsp; The analysis considers hierarchy completeness, the direction of information transfer (upward or downward through the hierarchy), and minimal padding strategies. &nbsp;Numerical results demonstrate speedups of up to 100\u00d7, even for moderate dimensions, when computing surplus coefficients for sparse-grid interpolation. &nbsp;Greater speedups are obtained when sparse grids are combined with the Discontinuous Galerkin (DG) method for the solution of high-dimensional Partial Differential Equations (PDEs).<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Miroslav Stoyanov obtained his Ph.D. in Mathematics from Virginia Tech in 2009. &nbsp;After a postdoctoral assignment at Florida State University, Dr. Stoyanov joined Oak Ridge National Laboratory in 2012, where he has continued his research. His research interests include uncertainty quantification, high-dimensional approximation theory and high-performance computing and the development of mathematical libraries.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Host<\/strong>:<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Pablo Seleson,&nbsp;<a href=\"mailto:selesonpd@ornl.gov\">selesonpd@ornl.gov<\/a>, 865-323-2859 <\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Miroslav Stoyanov, ORNL Abstract: Sparse Grids (SG) is a common tool for applications that require multidimensional function approximation, such as constructing surrogate models in parameter space or solving high-dimensional partial differential equations.&nbsp; From uncertainty quantification to Kinetic models for plasma physics, the key advantage of SG is the mitigation of the computational cost associated<\/p>\n","protected":false},"author":110,"featured_media":565,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-567","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/567","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=567"}],"version-history":[{"count":2,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/567\/revisions"}],"predecessor-version":[{"id":1066,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/567\/revisions\/1066"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/565"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=567"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=567"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=567"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":560,"date":"2026-09-22T13:22:49","date_gmt":"2026-09-22T17:22:49","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=560"},"modified":"2026-10-08T16:39:18","modified_gmt":"2026-10-08T20:39:18","slug":"closed-form-logical-error-rate-approximations-for-surface-codes","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/closed-form-logical-error-rate-approximations-for-surface-codes\/","title":{"rendered":"Closed Form Logical Error Rate Approximations for Surface Codes"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Shaked Regev, ORNL<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">The authors propose a novel method to calculate logical error rates in surface codes, assuming independent and identically distributed physical errors. &nbsp;These results fit well known scaling laws for studied physical error rates and code distances but break down in some other regimes. &nbsp;They show how to use their method to analyze hypothetical quantum computers with various configurations and select designs with lower error rates. &nbsp;Currently, this requires expensive classical simulations of quantum decoders for various distances and physical error rates or inaccurate extrapolation from minimal experimental data. &nbsp;Instead, they use the symmetry of the problem to count the configurations that result in a logical error with their novel software. &nbsp;Given a physical error rate, the method can deduce the probability of a logical error to provably good accuracy.&nbsp; The authors include an analysis of measurement errors to allow for a more complete comparison of different surface code implementations. &nbsp;They extend this counting method to a particular locally correlated error model.&nbsp;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong><strong>Speaker\u2019s Bio:<\/strong><\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Shaked Regev is a Research Scientist at Oak Ridge National Laboratory whose research interests include Quantum Error Correction, Numerical Linear Algebra, Mathematical Optimization, GPU Programming, and High-Performance Computing. &nbsp;He received his Ph.D. in Computational and Mathematical Engineering from Stanford University for his thesis on &#8220;Preconditioning Techniques for Sparse Linear Systems<strong>.<\/strong>&#8220;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Host<\/strong>:<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Elaine Wong, <a href=\"mailto:wongey@ornl.gov\">wongey@ornl.gov<\/a>, 865-341-2325<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Shaked Regev, ORNL Abstract: The authors propose a novel method to calculate logical error rates in surface codes, assuming independent and identically distributed physical errors. &nbsp;These results fit well known scaling laws for studied physical error rates and code distances but break down in some other regimes. &nbsp;They show how to use their method<\/p>\n","protected":false},"author":110,"featured_media":498,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-560","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/560","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=560"}],"version-history":[{"count":2,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/560\/revisions"}],"predecessor-version":[{"id":1067,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/560\/revisions\/1067"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/498"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=560"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=560"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=560"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":562,"date":"2026-09-17T15:15:00","date_gmt":"2026-09-17T19:15:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=562"},"modified":"2026-10-08T16:39:55","modified_gmt":"2026-10-08T20:39:55","slug":"adaptive-sparse-grid-discontinuous-galerkin-approximations-to-the-bhatnagar-gross-krook-model","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/adaptive-sparse-grid-discontinuous-galerkin-approximations-to-the-bhatnagar-gross-krook-model\/","title":{"rendered":"Adaptive Sparse-Grid Discontinuous Galerkin Approximations to the Bhatnagar\u2013Gross\u2013Krook Model"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Stefan Schnake, ORNL<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">This talk presents adaptive sparse-grid discontinuous Galerkin (DG) discretizations for the Bhatnagar\u2013Gross\u2013Krook (BGK) model, a kinetic equation posed in four- and six-dimensional phase-space.&nbsp; Standard DG methods are rendered impractical for the BGK model by the curse of dimensionality, motivating compressed representations that adapt to the solution in time.&nbsp; Using the adaptive sparse-grid DG method, we quantify accuracy and compression by comparing the adaptive degrees of freedom to full-grid DG methods and by assessing the resulting kinetic and fluid quantities in both fluid and rarefied regimes.&nbsp; Test cases include a relaxation problem, a multidimensional Sod shock tube, and shear\/expansion flows used in prior low-rank BGK studies.&nbsp; To build an efficient Maxwellian evaluation without violating conservation, a central obstacle for structure-persevering BGK simulations, we introduce a hybrid interpolation strategy that exploits velocity separability to recover the correct discrete collision invariants and prove conservation of the resulting discrete collision operator on adaptive sparse grids.&nbsp; The results show that the adaptive sparse-grid strategy can recover accurate and physically relevant solutions with sharp gradients, and the method reduces the active degrees of freedom by factors ranging from several-fold to several orders of magnitude, with the largest reductions occurring in the six-dimensional examples.&nbsp; All computations are performed with the open-source ASGarD adaptive sparse-grid DG library.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Stefan&nbsp;Schnake is a research scientist in the Multiscale Methods and Dynamics Group.&nbsp; He received his Ph.D. in 2017 from the University of Tennessee and joined Oak Ridge National Laboratory in 2020 after a postdoctoral appointment at the University of Oklahoma.&nbsp; His research interests include low-rank and sparse-grid methods for compressing tensor representations of solutions to dynamical systems.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Host<\/strong>:<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Pablo Seleson,&nbsp;<a href=\"mailto:selesonpd@ornl.gov\">selesonpd@ornl.gov<\/a>, 865-323-2859 &nbsp;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Stefan Schnake, ORNL Abstract: This talk presents adaptive sparse-grid discontinuous Galerkin (DG) discretizations for the Bhatnagar\u2013Gross\u2013Krook (BGK) model, a kinetic equation posed in four- and six-dimensional phase-space.&nbsp; Standard DG methods are rendered impractical for the BGK model by the curse of dimensionality, motivating compressed representations that adapt to the solution in time.&nbsp; Using the<\/p>\n","protected":false},"author":110,"featured_media":553,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-562","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/562","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=562"}],"version-history":[{"count":2,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/562\/revisions"}],"predecessor-version":[{"id":1068,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/562\/revisions\/1068"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/553"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=562"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=562"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=562"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":543,"date":"2026-08-27T15:15:00","date_gmt":"2026-08-27T19:15:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=543"},"modified":"2026-10-02T16:49:23","modified_gmt":"2026-10-02T20:49:23","slug":"exact-limits-of-random-projections-for-preserving-geometry-distance-recovery-nearest-neighbor-rankings-and-covariance-shape-in-gaussian-models","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/exact-limits-of-random-projections-for-preserving-geometry-distance-recovery-nearest-neighbor-rankings-and-covariance-shape-in-gaussian-models\/","title":{"rendered":"Exact Limits of Random Projections for Preserving Geometry: Distance Recovery, Nearest-Neighbor Rankings, and Covariance Shape in Gaussian Models"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Piyush Sao , ORNL<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>The random-projection guarantee and what it leaves unsaid.<\/strong>&nbsp; The Johnson\u2013Lindenstrauss (JL) lemma is a cornerstone of high-dimensional dimensionality reduction, ensuring that rank-m linear projections preserve pairwise distances within tight relative error bounds.&nbsp; However, in modern high-dimensional statistics, critical geometric decisions &#8211; such as nearest-neighbor ordering, statistical inference, and detecting covariance anisotropy &#8211; depend on centered fluctuations that are substantially smaller than the baseline distance itself.&nbsp; Standard JL guarantees maintain relative accuracy but say nothing about these local fluctuations, so algorithms that rely on fine structure can fail unexpectedly even when classical JL bounds are satisfied.&nbsp; This talk asks: how much of this fine-grained geometric structure actually survives low-rank linear sketching?<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Decomposing Distance Variance, Ranking Failure, and Spectral Contraction.<\/strong>&nbsp; This talk analyzes rank-$m$ linear projections of $d$-dimensional data across both isotropic and anisotropic regimes.&nbsp; For isotropic data, squared distances decompose into independent retained and discarded components, establishing an exact nonlinear recovery ceiling of $m\/d$ for distance-feature variance (derived via Laguerre polynomial expansions).&nbsp; In anisotropic settings, distance recovery depends primarily on the squared-eigenvalue mass captured by top-$m$ directions; however, nonlinear features can exceed this linear share when projections retain eigenspaces unevenly.&nbsp; Furthermore, covariance-shape information contracts at rate $(m\/d)^2$ &#8211; faster than the $m\/d$ rate for mean and scale &#8211; reflecting a &#8220;sphericalization&#8221; effect toward isotropy.&nbsp; Strikingly, even when pairwise distances meet standard JL tolerances, low-dimensional projections rank distances no better than chance when $\\log n &lt;&lt;m &lt;&lt; d$, and fixed-candidate nearest-neighbor overlaps degrade to random chance as $m\/d \\to 0$.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Implications for High-Dimensional Algorithm Design.<\/strong>&nbsp; These results establish that preserving relative distances is insufficient to preserve local structure.&nbsp; By defining explicit recovery ceilings and contraction rates, the findings demonstrate that tasks such as ranking data and inferring shape require mechanisms beyond a single oblivious projection, including reranking, metric correction, or multiple independent sketches.&nbsp; For general linear maps, three spectral summaries govern these downstream tasks: rank controls optimal nonlinear inference, spectral spread governs uncorrected norm estimation, and small singular values determine stability under additive noise.&nbsp; Together, these principles expose the limits of random projections and guide the design of sketching operators tailored to target applications.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Piyush K. Sao is a Research Scientist in the Computer Science and Mathematics Division at Oak Ridge National Laboratory, where he develops scalable numerical and combinatorial algorithms for high-performance computing platforms.&nbsp; His research spans theoretically optimal, energy-efficient, and fault-tolerant algorithms for sparse and dense linear algebra, graph analytics, and discrete algorithms; his recent work establishes information-theoretic limits of mixed-precision computation, fixed-point emulation, and randomized algorithms for scientific computing and machine learning.&nbsp; Algorithms he co-developed have been incorporated into widely used libraries, including SuperLU DIST and ArborX.&nbsp; He received the 2022 Best Paper Prize of the Society for Industrial and Applied Mathematics Activity Group on Supercomputing, and his work on large-scale knowledge-graph analytics was an Association for Computing Machinery Gordon Bell Prize finalist at SC20 and SC22.&nbsp; He has served on the program committees of the Supercomputing Conference and International Parallel and Distributed Processing Symposium since 2019.&nbsp; He received his Ph.D. in Computational Science and Engineering from Georgia Tech and an undergraduate degree in Electrical Engineering from the Institute of Technology, Madras.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Host<\/strong>:<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Elaine Wong, <a href=\"mailto:wongey@ornl.gov\">wongey@ornl.gov<\/a>, 865-341-2325<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Piyush Sao , ORNL Abstract: The random-projection guarantee and what it leaves unsaid.&nbsp; The Johnson\u2013Lindenstrauss (JL) lemma is a cornerstone of high-dimensional dimensionality reduction, ensuring that rank-m linear projections preserve pairwise distances within tight relative error bounds.&nbsp; However, in modern high-dimensional statistics, critical geometric decisions &#8211; such as nearest-neighbor ordering, statistical inference, and detecting<\/p>\n","protected":false},"author":110,"featured_media":467,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-543","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/543","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=543"}],"version-history":[{"count":3,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/543\/revisions"}],"predecessor-version":[{"id":1054,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/543\/revisions\/1054"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/467"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=543"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=543"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=543"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":540,"date":"2026-08-20T15:00:00","date_gmt":"2026-08-20T19:00:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=540"},"modified":"2026-09-09T15:28:14","modified_gmt":"2026-09-09T19:28:14","slug":"multifidelity-training-methods-for-machine-learning","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/multifidelity-training-methods-for-machine-learning\/","title":{"rendered":"Multifidelity Training Methods for Machine Learning"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Elizabeth Qian<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Machine learning (ML) methods have garnered significant interest as potential approaches for constructing surrogate models of complex engineering systems for which traditional simulation is computationally expensive. &nbsp;However, in many scientific and engineering applications, training data are scarce because of the high cost of generating data from traditional high-fidelity simulations. &nbsp;ML models trained on limited data often yield unreliable predictions outside the regime represented by the training data.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">This talk presents multifidelity training approaches for ML that exploit scientific settings in which data of varying fidelities and computational costs are available. &nbsp;For example, high-fidelity data may be generated by an expensive, fully resolved physics simulation, whereas lower-fidelity data may be obtained from a less computationally expensive model based on simplifying assumptions. &nbsp;Two recent approaches developed in the speaker&#8217;s research group for multifidelity training are highlighted: one based on linear regression and another based on nonlinear autoregressive Gaussian process regression. &nbsp;The performance of these approaches is demonstrated on problems arising in the modeling of reacting flows.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Elizabeth Qian holds a joint appointment at the Georgia Institute of Technology as an Assistant Professor in the School of Aerospace Engineering and the School of Computational Science and Engineering. &nbsp;Her interdisciplinary research focuses on developing computational methods that enable engineering design and decision-making for complex systems. &nbsp;Her research specialties include the development of efficient surrogate models through model reduction and scientific ML, as well as multifidelity approaches for accelerating expensive computations in uncertainty quantification, optimization, and control.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Qian previously held a postdoctoral appointment as a von K\u00e1rm\u00e1n Instructor in the Department of Computing + Mathematical Sciences at the California Institute of Technology. &nbsp;She earned her Doctor of Philosophy, Master of Science, and Bachelor of Science degrees from the Department of Aeronautics and Astronautics at the Massachusetts Institute of Technology. &nbsp;Her awards and honors include a Caltech-wide teaching award selected by the undergraduate student body, the 2020 Society for Industrial and Applied Mathematics Student Paper Prize, a Fannie and John Hertz Foundation Fellowship, a National Science Foundation Graduate Research Fellowship, and a United States Fulbright Student grant. &nbsp;&nbsp;She currently holds a visiting faculty appointment as a Hans Fischer Fellow at the Technical University of Munich<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Elizabeth Qian Abstract: Machine learning (ML) methods have garnered significant interest as potential approaches for constructing surrogate models of complex engineering systems for which traditional simulation is computationally expensive. &nbsp;However, in many scientific and engineering applications, training data are scarce because of the high cost of generating data from traditional high-fidelity simulations. &nbsp;ML models<\/p>\n","protected":false},"author":110,"featured_media":541,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-540","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/540","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=540"}],"version-history":[{"count":1,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/540\/revisions"}],"predecessor-version":[{"id":542,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/540\/revisions\/542"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/541"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=540"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=540"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=540"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":537,"date":"2026-04-09T15:15:00","date_gmt":"2026-04-09T19:15:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=537"},"modified":"2026-09-09T15:24:48","modified_gmt":"2026-09-09T19:24:48","slug":"accelerated-and-communication-efficient-federated-learning-algorithms","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/accelerated-and-communication-efficient-federated-learning-algorithms\/","title":{"rendered":"Accelerated and Communication Efficient Federated Learning Algorithms"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Paul Laiu , Oak Ridge National Laboratory<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Federated Learning (FL) is a distributed Machine Learning (ML) paradigm that enables multiple local devices, that is, clients, and a central server to collaboratively train a ML model using data stored locally on the clients without transferring the data. &nbsp;Common challenges in FL include data privacy preservation, communication efficiency, and convergence guarantees.&nbsp; This presentation introduces two new FL algorithms: Accelerating Federated Learning with One-Step Anderson Acceleration (FedOSAA) and Federated Dynamical Low-Rank Training (FeDLRT). &nbsp;FedOSAA leverages a widely used acceleration scheme, Anderson acceleration, to improve the convergence rate of standard FL schemes while incurring minimal additional communication cost.&nbsp; FeDLRT adopts a dynamical low-rank training scheme that constrains model training to low-rank factors of the model parameters, resulting in reductions in both communication and memory costs during the training process.&nbsp; By incorporating a variance-reduction scheme, convergence guarantees are established for both algorithms, even when the data distribution is heterogeneous, that is, non\u2013independent and identically distributed, among the local clients.&nbsp; The advantages of these two FL algorithms are demonstrated through numerical experiments on ML problems ranging from logistic regression to the training of transformer models in a FL setting.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Paul Laiu is&nbsp;a Staff Mathematician in the&nbsp;Multiscale Methods and Dynamics Group&nbsp;at the Oak Ridge National Laboratory. &nbsp;He&nbsp;received his Ph.D. degree in Electrical and Computer Engineering from&nbsp;University of Maryland&nbsp;College Park in 2016. &nbsp;Paul\u2019s research interest includes scientific machine learning, surrogate modeling, iterative solvers, and numerical schemes for various partial differential equations. &nbsp;His work focuses on the design, development, and analysis of mathematical tools that accelerate the simulation and learning of multiscale systems, with applications in astrophysics, computational fluid dynamics, cybersecurity, fusion, and rarefied gas.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Paul Laiu , Oak Ridge National Laboratory Abstract: Federated Learning (FL) is a distributed Machine Learning (ML) paradigm that enables multiple local devices, that is, clients, and a central server to collaboratively train a ML model using data stored locally on the clients without transferring the data. &nbsp;Common challenges in FL include data privacy<\/p>\n","protected":false},"author":110,"featured_media":538,"comment_status":"closed","ping_status":"","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-537","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/537","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=537"}],"version-history":[{"count":1,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/537\/revisions"}],"predecessor-version":[{"id":539,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/537\/revisions\/539"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/538"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=537"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=537"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=537"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":11,"date":"2026-04-02T03:15:22","date_gmt":"2026-04-02T07:15:22","guid":{"rendered":"http:\/\/localhost\/golden-wp\/?p=11"},"modified":"2026-09-09T15:33:44","modified_gmt":"2026-09-09T19:33:44","slug":"toward-consistent-data-assimilation","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/toward-consistent-data-assimilation\/","title":{"rendered":"Toward Consistent Data Assimilation with Structured Generative Models"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Ricardo Baptista , University of Toronto<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Accurate state estimation, also known as data assimilation, is essential for geophysical forecasts, ranging from numerical weather prediction to long-term climate studies.&nbsp; While ensemble Kalman methods are widely adopted for this task in high-dimensional settings, these methods are inconsistent in capturing true uncertainty in non-Gaussian settings.&nbsp; This presentation introduces a framework for consistent data assimilation.&nbsp; First, inference methods based on conditional generative models are shown to generalize ensemble Kalman methods and to better characterize both the state and its associated uncertainty in nonlinear filtering problems.&nbsp; Second, a dimension-reduction approach for limited-data settings is presented by identifying and encoding low-dimensional structure in generative models, with guarantees on approximation error.&nbsp; The benefits of this framework are demonstrated through applications in fluid mechanics involving chaotic dynamics, where classical methods exhibit instability in small-sample regimes.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Ricardo Baptista is an Assistant Professor in Statistical Sciences at the&nbsp;University of Toronto&nbsp;and a Faculty Affiliate at the&nbsp;Vector Institute. &nbsp;His research focuses on establishing the&nbsp;mathematical foundations and theoretical guarantees of probabilistic machine&nbsp;learning models, with broad applications in science and engineering. &nbsp;Before&nbsp;joining the University of Toronto, he served as an Instructor in Computing and Mathematical Sciences at California Institute of Technology,&nbsp;and as a Postdoctoral Scientist at Amazon. &nbsp;He holds a Ph.D.&nbsp;in Computational&nbsp;Science and Engineering from the Massachusetts Institute of Technology, where he specialized in uncertainty&nbsp;quantification, and a Bachelor of Applied Science in Engineering Science from the University of&nbsp;Toronto.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">April 02<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">3:15pm &#8211; 4:15pm<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">J302 5600<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Ricardo Baptista , University of Toronto Abstract: Accurate state estimation, also known as data assimilation, is essential for geophysical forecasts, ranging from numerical weather prediction to long-term climate studies.&nbsp; While ensemble Kalman methods are widely adopted for this task in high-dimensional settings, these methods are inconsistent in capturing true uncertainty in non-Gaussian settings.&nbsp; This<\/p>\n","protected":false},"author":110,"featured_media":506,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-11","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/11","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=11"}],"version-history":[{"count":1,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/11\/revisions"}],"predecessor-version":[{"id":507,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/11\/revisions\/507"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/506"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=11"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=11"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=11"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":511,"date":"2026-03-05T15:00:00","date_gmt":"2026-03-05T20:00:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=511"},"modified":"2026-09-25T13:40:00","modified_gmt":"2026-09-25T17:40:00","slug":"on-the-choice-of-optimization-norm-for-anderson-acceleration-of-the-picard-iteration-for-navier-stokes-equations","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/on-the-choice-of-optimization-norm-for-anderson-acceleration-of-the-picard-iteration-for-navier-stokes-equations\/","title":{"rendered":"On the Choice of Optimization Norm for Anderson Acceleration of the Picard Iteration for Navier-Stokes Equations"},"content":{"rendered":"\n<h3 class=\"wp-block-heading has-dark-color has-text-color\">Dr. Leo Rebholz , Clemson University<\/h3>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Recent convergence theory for Anderson Acceleration (AA) requires that the AA optimization norm coincide with the Hilbert space norm associated with the fixed-point operator.&nbsp; In practical implementations, however, the $\\ell^2$ norm is among the most commonly selected norms.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Limited research has addressed this discrepancy or identified conditions under which the $\\ell^2$ norm is theoretically justified.&nbsp; To investigate this issue, AA is applied to the Picard iteration for the Navier-Stokes Equations (NSE) with varying choices of optimization norm.&nbsp; A sharpened and generalized convergence estimate is first established for depth $m$ Anderson acceleration\u2013Picard applied to the NSE when the $H^1_0$ optimization norm is employed.&nbsp; The analysis is problem-specific and incorporates a refined treatment of the nonlinear terms relative to prior AA-Picard convergence studies.&nbsp; The small-data assumption present in earlier analyses is removed.&nbsp; New Anderson acceleration term identities are developed to strengthen the nonlinear term estimates.&nbsp; A convergence result is then proved for the case in which the $L^2$ optimization norm is used.&nbsp; This estimate is shown to be closely aligned with the corresponding $H^1_0$ result.&nbsp; Although an analogous theoretical framework does not appear attainable for the $\\ell^2$ norm, several numerical experiments are conducted to compare AA-Picard convergence under different optimization norm selections.&nbsp; The results indicate that convergence behavior is consistently similar for the $L^2$ and $H^1_0$ norms and generally, though not universally, similar for the $\\ell^2$ norm.&nbsp; In a benchmark problem involving channel flow past a cylinder on relatively coarse meshes, AA-Picard employing the $\\ell^2$ norm exhibits substantially poorer convergence than when the $L^2$ or $H^1_0$ norms are used.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Leo Rebholz is a Dean\u2019s Distinguished Professor in the School of Mathematical and Statistical Sciences at Clemson University.&nbsp; His scholarly work encompasses a broad spectrum of computational fluid dynamics, including the numerical analysis and development of stable, high-fidelity discretization techniques for incompressible flow models. In recent years, his research has emphasized data assimilation methodologies and the acceleration of nonlinear solvers, with particular focus on enhancing robustness, convergence properties, and computational efficiency for large-scale, multiphysics simulations.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dr. Rebholz has authored more than 130 peer-reviewed journal articles and five scholarly monographs.&nbsp; He has supervised and mentored an extensive cohort of graduate and undergraduate researchers, including 13 doctoral dissertations, 14 master\u2019s theses, and 15 undergraduate research projects, contributing significantly to the advancement of computational and applied mathematics research training.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">March 05<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">3:15pm &#8211; 4:15pm<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">H308 5600<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Dr. Leo Rebholz , Clemson University Abstract: Recent convergence theory for Anderson Acceleration (AA) requires that the AA optimization norm coincide with the Hilbert space norm associated with the fixed-point operator.&nbsp; In practical implementations, however, the $\\ell^2$ norm is among the most commonly selected norms. Limited research has addressed this discrepancy or identified conditions under<\/p>\n","protected":false},"author":110,"featured_media":510,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-511","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/511","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=511"}],"version-history":[{"count":2,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/511\/revisions"}],"predecessor-version":[{"id":514,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/511\/revisions\/514"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/510"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=511"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=511"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=511"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}},{"id":515,"date":"2026-02-26T03:15:00","date_gmt":"2026-02-26T08:15:00","guid":{"rendered":"https:\/\/events.ornl.gov\/mictalks\/?p=515"},"modified":"2026-09-09T15:17:29","modified_gmt":"2026-09-09T19:17:29","slug":"computational-aspects-of-radiative-transfer-asymptotics-to-uncertainty-quantification","status":"publish","type":"post","link":"https:\/\/events.ornl.gov\/mictalks\/computational-aspects-of-radiative-transfer-asymptotics-to-uncertainty-quantification\/","title":{"rendered":"Computational Aspects of Radiative Transfer: Asymptotics to Uncertainty Quantification"},"content":{"rendered":"\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Abstract:<\/strong>&nbsp;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Thermal radiative transfer (TRT) governs phenomena ranging from supernovas in astrophysics to laser-driven fusion experiments in plasma physics.&nbsp; The interaction of radiation and matter is characterized by prohibitively small-time scales, nonlinear coupling, and high-dimensional particle dynamics, rendering conventional numerical methods computationally expensive.&nbsp; Moreover, physical systems are rarely deterministic; uncertainties arising from modelling assumptions, measurement errors, and device imperfections must be incorporated to obtain realistic descriptions of the underlying phenomena.&nbsp; Quantifying these uncertainties is therefore essential for robustness and reliability in applications.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Dynamical Low-Rank Approximation (DLRA) provides a promising approach for addressing both computational complexity and uncertainty quantification.&nbsp; The first part of this talk presents strategies for constructing low-rank integrators that remain robust in the presence of prohibitively small-time scales, commonly referred to as asymptotic preserving schemes.&nbsp; Challenges encountered in real-world applications are discussed along with potential mitigation strategies.&nbsp; The second part introduces a Low-Rank Multilevel Monte Carlo framework for uncertainty quantifications in TRT equations based on geometric spatial refinement. The efficacy of these methods is demonstrated through a range of test cases and benchmark problems.&nbsp;<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><strong>Speaker\u2019s Bio:<\/strong><\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\">Chinmay Patwardhan&nbsp;is a Ph.D. student in the Scientific Computing and Mathematics Department at the Karlsruhe Institute of Technology.&nbsp; His research focuses on model order reduction techniques and uncertainty quantification methods, with an emphasis on developing cost-effective computational approaches based on DLRA.&nbsp; These methods&nbsp;aim to reduce the computational cost of simulating kinetic equations and addressing high-dimensional problems.&nbsp; His additional research interests include scientific machine learning and the development of simulation frameworks. &nbsp;Chinmay\u2019s&nbsp;work lies&nbsp;at the intersection of applied mathematics and computational science, with the goal of advancing efficient numerical methods for complex systems.<\/p>\n\n\n\n<p class=\"has-dark-color has-text-color has-lg-font-size wp-block-paragraph\"><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Abstract:&nbsp; Thermal radiative transfer (TRT) governs phenomena ranging from supernovas in astrophysics to laser-driven fusion experiments in plasma physics.&nbsp; The interaction of radiation and matter is characterized by prohibitively small-time scales, nonlinear coupling, and high-dimensional particle dynamics, rendering conventional numerical methods computationally expensive.&nbsp; Moreover, physical systems are rarely deterministic; uncertainties arising from modelling assumptions, measurement<\/p>\n","protected":false},"author":110,"featured_media":516,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"_acf_changed":false,"footnotes":""},"categories":[1],"tags":[],"class_list":["post-515","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-uncategorized"],"acf":[],"_links":{"self":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/515","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/users\/110"}],"replies":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/comments?post=515"}],"version-history":[{"count":3,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/515\/revisions"}],"predecessor-version":[{"id":532,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/posts\/515\/revisions\/532"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media\/516"}],"wp:attachment":[{"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/media?parent=515"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/categories?post=515"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/events.ornl.gov\/mictalks\/wp-json\/wp\/v2\/tags?post=515"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}]