Award Abstract # 2404521
Development of Adaptive Sparse Grid Discontinuous Galerkin Methods for Multiscale Kinetic Simulations in Plasmas

NSF Org: DMS
Division Of Mathematical Sciences
Recipient: VIRGINIA POLYTECHNIC INSTITUTE & STATE UNIVERSITY
Initial Amendment Date: December 19, 2023
Latest Amendment Date: December 19, 2023
Award Number: 2404521
Award Instrument: Standard Grant
Program Manager: Yuliya Gorb
ygorb@nsf.gov
 (703)292-2113
DMS
 Division Of Mathematical Sciences
MPS
 Directorate for Mathematical and Physical Sciences
Start Date: December 15, 2023
End Date: July 31, 2024 (Estimated)
Total Intended Award Amount: $200,000.00
Total Awarded Amount to Date: $55,191.00
Funds Obligated to Date: FY 2020 = $55,191.00
History of Investigator:
  • Yingda Cheng (Principal Investigator)
    yingda@vt.edu
Recipient Sponsored Research Office: Virginia Polytechnic Institute and State University
300 TURNER ST NW
BLACKSBURG
VA  US  24060-3359
(540)231-5281
Sponsor Congressional District: 09
Primary Place of Performance: Virginia Polytechnic Institute and State University
300 TURNER ST NW
BLACKSBURG
VA  US  24060-3359
Primary Place of Performance
Congressional District:
09
Unique Entity Identifier (UEI): QDE5UHE5XD16
Parent UEI: X6KEFGLHSJX7
NSF Program(s): COMPUTATIONAL MATHEMATICS
Primary Program Source: 01002021DB NSF RESEARCH & RELATED ACTIVIT
Program Reference Code(s): 9263
Program Element Code(s): 127100
Award Agency Code: 4900
Fund Agency Code: 4900
Assistance Listing Number(s): 47.049

ABSTRACT

This project aims at designing efficient numerical schemes for simulating complex plasma phenomena. Plasma is a state of matter similar to gas in which a certain portion of the particles is ionized. Understanding the complex behaviors of plasmas has led to important advances in areas ranging from space physics, fusion energy, to high-power microwave generation and large scale particle accelerators. There is strong need for laying out mathematical and algorithmic foundations for the design of efficient numerical methods so that we can advance basic research in plasma simulations. The algorithms developed in this project have the potential to provide high fidelity simulations in plasma physics with manageable computational cost and will have applications and impacts in multiscale simulations in fusion devices. The principal investigator (PI) will organize special events at professional meetings and workshops to promote the participation of female researchers. This project provides research training opportunities for graduate students.

The objective of the project is to make significant advances on the design and analysis of a class of numerical methods called adaptive sparse grid (aSG) discontinuous Galerkin (DG) methods. The methods incorporate high order accurate DG solver that excels at transport simulations and the dimension reduction technique by aSG approach. The aim of this proposal is to advance the algorithmic foundations of the schemes for time-dependent PDEs, and push them onto the broader arena of multiscale simulations and applications for fusion science. The PI will investigate several fundamental issues including the analysis of CFL conditions, development of multiscale time stepping, postprocessing and hybrid aSG schemes. For a class of multiscale kinetic problems bridging kinetic and fluid models, by utilizing the multiresolution offered in the aSG-DG framework, the research will take advantage of both multiscale simulation tools and multiresolution on hierarchically defined meshes to achieve acceleration in computations. The schemes will be applied to simulations of runaway electrons in tokamak devices.

This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.

PUBLICATIONS PRODUCED AS A RESULT OF THIS RESEARCH

Note:  When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

Peng, Zhichao and Chen, Yanlai and Cheng, Yingda and Li, Fengyan "A Micro-Macro Decomposed Reduced Basis Method for the Time-Dependent Radiative Transfer Equation" Multiscale Modeling & Simulation , v.22 , 2024 https://doi.org/10.1137/22M1533487 Citation Details

PROJECT OUTCOMES REPORT

Disclaimer

This Project Outcomes Report for the General Public is displayed verbatim as submitted by the Principal Investigator (PI) for this award. Any opinions, findings, and conclusions or recommendations expressed in this Report are those of the PI and do not necessarily reflect the views of the National Science Foundation; NSF has not approved or endorsed its content.

In many science and engineering applications, the underlying models are posed as partial differential equations in high dimensions. Those include Boltzmann equations in statistical physics, Schrodinger equations in quantum mechanics. They go beyond the physical dimensions in 3D, and are challenging to simulate due to limited computational resources. 

This project develops a class of highly accurate numerical methods to tackle the challenge of curse of dimensions. We developed the adaptive sparse grid discontinous Galerkin methods, which can use adaptive choice of multiwavelet bases to efficiently compute such problems. The method can compute highly challenging nonlinear equations by a novel adaptive sparse grid by interpolatory multiwavelet bases. Beyond the computational efficiency, the method can preserve many physical quantities of interest due to the weak formulation underlying the discontinuous Galerkin framework. 

The building blocks of the method includes the following: (1) a weak formulation from discontinuous Galerkin finite element method (2) tools from signal processing: multiwavelets of various kinds (3) fast linear algebra routines exploring hierarchical structures of the grid. We applied the schemes to several applications: including Vlasov/Boltzmann equations in plasmas, nonlinear dispersive equations, Hamilton-Jacobi equations in control. We publish a publicly-available C++ software on Github providing benchmark code and algorithm implementation.  The package is capable of treating a large class of high dimensional linear and nonlinear PDEs. We establish theoretical results for stability and convergence for linear equations.

The project has trained one Ph.D. students and involves two other postdoctoral scholars. They have continued their academic careers in postdoc and tenure track positions.

 

 

 

 


Last Modified: 09/10/2024
Modified by: Yingda Cheng

Please report errors in award information by writing to: awardsearch@nsf.gov.

Print this page

Back to Top of page