Computational Science Laboratory
Browse by
The mission of the Computational Science Laboratory (CSL) is to develop innovative computational solutions for complex realworld problems, and to foster a productive research and education environment emphasizing collaboration and innovation.
Collections in this community

Scholarly Works, Computational Science Laboratory [19]
Research articles, presentations, and other scholarship
Recent Submissions

Linearly Implicit General Linear Methods
(20211201)Linearly implicit Runge–Kutta methods provide a fitting balance of implicit treat ment of stiff systems and computational cost. In this paper we extend the class of linearly implicit Runge–Kutta methods to include multistage ... 
Physicsinformed neural networks for PDEconstrained optimization and control
(20220506)A fundamental problem of science is designing optimal control policies that manipulate a given environment into producing the desired outcome. Control PhysicsInformed Neural Networks simultaneously solve a given system ... 
Linearly implicit GARK schemes
(Elsevier, 20210301)Systems driven by multiple physical processes are central to many areas of science and engineering. Time discretization of multiphysics systems is challenging, since different processes have different levels of stiffness ... 
Machine learning based algorithms for uncertainty quantification in numerical weather prediction models
(Elsevier, 20210301)Complex numerical weather prediction models incorporate a variety of physical processes, each described by multiple alternative physical schemes with specific parameters. The selection of the physical schemes and the choice ... 
Partitioned exponential methods for coupled multiphysics systems
(Elsevier, 20210301)Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. ... 
Symplectic GARK methods for Hamiltonian systems
(20210306)Generalized Additive RungeKutta schemes have shown to be a suitable tool for solving ordinary differential equations with additively partitioned righthand sides. This work generalizes these GARK schemes to symplectic ... 
Subspace adaptivity in RosenbrockKrylov methods for the time integration of initial value problems
(Elsevier, 20210315)The Rosenbrock–Krylov family of time integration schemes is an extension of RosenbrockW methods that employs a specific Krylov based approximation of the linear system solutions arising within each stage of the integrator. ... 
Alternating directions implicit integration in a general linear method framework
(Elsevier, 20210515)Alternating Directions Implicit (ADI) integration is an operator splitting approach to solve parabolic and elliptic partial differential equations in multiple dimensions based on solving sequentially a set of related ... 
Multirate implicit Euler schemes for a class of differentialalgebraic equations of index1
(Elsevier, 20210515)Systems of differential equations which consist of subsystems with widely different dynamical behaviour can be integrated by multirate time integration schemes to increase the efficiency. These schemes allow the usage of ... 
Investigation of Nonlinear Model Order Reduction of the Quasigeostrophic Equations through a PhysicsInformed Convolutional Autoencoder
(20210827)Reduced order modeling (ROM) is a field of techniques that approximates complex physicsbased models of realworld processes by inexpensive surrogates that capture important dynamical characteristics with a smaller number ... 
A Stochastic Covariance Shrinkage Approach to Particle Rejuvenation in the Ensemble Transform Particle Filter
(20210920)Rejuvenation in particle filters is necessary to prevent the collapse of the weights when the number of particles is insufficient to sample the high probability regions of the state space. Rejuvenation is often implemented ... 
AdjointMatching Neural Network Surrogates for Fast 4DVar Data Assimilation
(20211116)The data assimilation procedures used in many operational numerical weather forecasting systems are based around variants of the 4DVar algorithm. The cost of solving the 4DVar problem is dominated by the cost of forward ... 
An Ensemble Variational FokkerPlanck Method for Data Assimilation
(20211127)Particle flow filters that aim to smoothly transform particles from samples from a prior distribution to samples from a posterior are a major topic of active research. In this work we introduce a generalized framework which ... 
Multirate linearlyimplicit GARK schemes
(Springer, 20211228)Many complex applications require the solution of initialvalue problems where some components change fast, while others vary slowly. Multirate schemes apply different step sizes to resolve different components of the ... 
Eliminating Order Reduction on Linear, TimeDependent ODEs with GARK Methods
(20220119)When applied to stiff, linear differential equations with timedependent forcing, RungeKutta methods can exhibit convergence rates lower than predicted by the classical order condition theory. Commonly, this order reduction ... 
A unified formulation of splittingbased implicit time integration schemes
(ACADEMIC PRESS INC ELSEVIER SCIENCE, 20220101)Splittingbased time integration approaches such as fractional step, alternating direction implicit, operator splitting, and locally one dimensional methods partition the system of interest into components, and solve ... 
Robust data assimilation using L1 and Huber norms
(20151106)Data assimilation is the process to fuse information from priors, observations of nature, and numerical models, in order to obtain best estimates of the parameters or state of a physical system of interest. Presence of ... 
Efficient Construction of Local Parametric Reduced Order Models Using Machine Learning Techniques
(20151111)Reduced order models are computationally inexpensive approximations that capture the important dynamical characteristics of large, highfidelity computer models of physical systems. This paper applies machine learning ... 
The reducedorder hybrid Monte Carlo sampling smoother
(WileyBlackwell, 20170110)Hybrid MonteCarlo (HMC) sampling smoother is a fully nonGaussian fourdimensional data assimilation algorithm that works by directly sampling the posterior distribution formulated in the Bayesian framework. The smoother ...