Paulina Sepulveda

Juliusz Wasieleski, Tomasz Służalec, Maciej Woźniak et al.

We develop a wildfire simulation model that evolves the temperature scalar field using an energy balance equation accounting for heat generation, transport, and loss. For these equations, we develop quasi-implicit time integration schemes using direction splitting of the differential operators. We use the Peaceman-Rachford and Strang splitting methods, including the Crank-Nicolson method. Based on these discretizations, we derive variational formulations and explore the Kronecker product structure of the matrices. In the wildfire model, there are some non-linear terms that we treat explicitly. We perform a detailed analysis of how treating these terms affects the stability of the time integration scheme. Namely, we show that a quasi-implicit time integration scheme achieves 10 times higher simulation accuracy. We present two wildfire simulations. The first is a simulation of the 2024 wildfire disaster in the Valparaíso region of Chile. The second one is a simulation of the 2019 wildfire disaster in Las Palmas de Gran Canaria, Spain. We discuss the numerical results and compare them against satellite images and measurement records. We also present a numerical experiment for comparison with the state-of-the-art wildfire simulation model FARSITE. Our sequential code has a linear computational cost of ${\cal O}(N)$. We also present the parallel scalability of the WILDFIRE-IGA-ADS code to illustrate the possibility of running the code on a local workstation.

Jay Gopalakrishnan, Paulina Sepulveda

A spacetime discontinuous Petrov-Galerkin (DPG) method for the linear wave equation is presented. This method is based on a weak formulation that uses a broken graph space. The wellposedness of this formulation is established using a previously presented abstract framework. One of the main tasks in the verification of the conditions of this framework is proving a density result. This is done in detail for a simple domain in arbitrary dimensions. The DPG method based on the weak formulation is then studied theoretically and numerically. Error estimates and numerical results are presented for triangular, rectangular, tetrahedral, and hexahedral meshes of the spacetime domain. The potential for using the built-in error estimator of the DPG method for an adaptivity mesh refinement strategy in two and three dimensions is also presented.

Leszek Demkowicz, Jay Gopalakrishnan, Sriram Nagaraj et al.

A spacetime Discontinuous Petrov Galerkin (DPG) method for the linear time-dependent Schrodinger equation is proposed. The spacetime approach is particularly attractive for capturing irregular solutions. Motivated by the fact that some irregular Schrodinger solutions cannot be solutions of certain first order reformulations, the proposed spacetime method uses the second order Schrodinger operator. Two variational formulations are proved to be well-posed: a strong formulation (with no relaxation of the original equation) and a weak formulation (also called the ultraweak formulation, that transfers all derivatives onto test functions). The convergence of the DPG method based on the ultraweak formulation is investigated using an interpolation operator. A standalone appendix analyzes the ultraweak formulation for general differential operators. Reports of numerical experiments motivated by pulse propagation in dispersive optical fibers are also included.

Jay Gopalakrishnan, Peter Monk, Paulina Sepulveda

Certain Friedrichs systems can be posed on Hilbert spaces normed with a graph norm. Functions in such spaces arising from advective problems are found to have traces with a weak continuity property at points where the inflow and outflow boundaries meet. Motivated by this continuity property, an explicit space-time finite element scheme of the tent pitching type, with spaces that conform to the continuity property, is designed. Numerical results for a model one-dimensional wave propagation problem are presented.