High-order finite-volume modeling of drift waves

Dorf, M.; Dorr, M.; Hittinger, J.; Lee, W.; Ghosh, D.

doi:10.1016/j.jcp.2018.07.009

Title: High-order finite-volume modeling of drift waves

Full Record
Other Related Research

Abstract

The study discusses high-order finite-volume numerical modeling of drift waves, which is an ubiquitous phenomenon in magnetized plasmas. It is found that some standard discretization methods applied to the conservative form of the governing equations can lead to a numerical instability. A method to stabilize high-order discretization is proposed and demonstrated to work in numerical simulations performed with the fourth-order finite-volume code COGENT. As practical examples, a stable drift-wave solution with adiabatic electrons and the collisionless (universal) drift-wave instability driven by electron kinetic effects are considered. Finally, application of the present analysis to a broader range of computational fluid dynamics systems is discussed.

Authors:

Dorf, M. ^[1]; Dorr, M. ^[1]; Hittinger, J. ^[1]; Lee, W. ^[2]; Ghosh, D. ^[1]

Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Univ. of California, San Diego, CA (United States)

Publication Date:: Tue Jul 10 00:00:00 EDT 2018

Research Org.:: Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

Sponsoring Org.:: USDOE Office of Science (SC), Fusion Energy Sciences (FES)

OSTI Identifier:: 1476214

Alternate Identifier(s):: OSTI ID: 1564510

Report Number(s):: LLNL-JRNL-733873
Journal ID: ISSN 0021-9991; 885907

Grant/Contract Number:: AC52-07NA27344

Resource Type:: Accepted Manuscript

Journal Name:: Journal of Computational Physics

Additional Journal Information:: Journal Volume: 373; Journal ID: ISSN 0021-9991

Publisher:: Elsevier

Country of Publication:: United States

Language:: English

Subject:: 97 MATHEMATICS AND COMPUTING; 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; finite-volume method; high-order discretization; incompressible flow; gyrokinetic simulation; drift waves

Citation Formats


                    Dorf, M., Dorr, M., Hittinger, J., Lee, W., and Ghosh, D. High-order finite-volume modeling of drift waves.  United States: N. p., 2018. 
Web.  doi:10.1016/j.jcp.2018.07.009.

Copy to clipboard


                    Dorf, M., Dorr, M., Hittinger, J., Lee, W., & Ghosh, D. High-order finite-volume modeling of drift waves.  United States.  https://doi.org/10.1016/j.jcp.2018.07.009

Copy to clipboard


                    Dorf, M., Dorr, M., Hittinger, J., Lee, W., and Ghosh, D. Tue .  
"High-order finite-volume modeling of drift waves".  United States.  https://doi.org/10.1016/j.jcp.2018.07.009.  https://www.osti.gov/servlets/purl/1476214.

Copy to clipboard


                    
@article{osti_1476214,

  title        = {High-order finite-volume modeling of drift waves},

  author       = {Dorf, M. and Dorr, M. and Hittinger, J. and Lee, W. and Ghosh, D.},

  abstractNote = {The study discusses high-order finite-volume numerical modeling of drift waves, which is an ubiquitous phenomenon in magnetized plasmas. It is found that some standard discretization methods applied to the conservative form of the governing equations can lead to a numerical instability. A method to stabilize high-order discretization is proposed and demonstrated to work in numerical simulations performed with the fourth-order finite-volume code COGENT. As practical examples, a stable drift-wave solution with adiabatic electrons and the collisionless (universal) drift-wave instability driven by electron kinetic effects are considered. Finally, application of the present analysis to a broader range of computational fluid dynamics systems is discussed.},

  doi          = {10.1016/j.jcp.2018.07.009},

  journal      = {Journal of Computational Physics},

  number       = ,

  volume       = 373,

  place        = {United States},

  year         = {Tue Jul 10 00:00:00 EDT 2018},

  month        = {Tue Jul 10 00:00:00 EDT 2018}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (Publisher)

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.1016/j.jcp.2018.07.009

Other availability

Search WorldCat to find libraries that may hold this journal

Citation Metrics:

Cited by: 2 works

Citation information provided by
Web of Science

Save / Share:

Export Metadata

Save to My Library

Similar Records in DOE PAGES and OSTI.GOV collections:

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

Journal Article Zhang, Lu ; Wang, Siyang ; Petersson, N. Anders - SIAM Journal on Scientific Computing

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinementmore »« less
https://doi.org/10.1137/20m1339702

Full Text Available
Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r,v_r,v_θ) phase space coordinates

Journal Article Vogman, G. V. ; Shumlak, U. ; Colella, P. - Journal of Computational Physics

© 2018 Elsevier Inc. A fourth-order finite-volume Vlasov–Poisson algorithm is developed for simulating axisymmetric plasma configurations in (r,vr,vθ) phase space coordinates. The Vlasov equation for cylindrical phase space coordinates is cast into conservation-law form and is discretized on a structured grid. The conservative finite-volume discretization is based on fifth-order upwind reconstructions of the distribution function and a fourth-order quadrature rule that accounts for transverse variations of fluxes along control volume surfaces. High-order specular reflection boundary conditions enable high-fidelity treatment of plasma distribution functions at the axis and at wall boundaries. The numerical method is applied to simulate a confined uniformmore »« less
Cited by 13
https://doi.org/10.1016/j.jcp.2018.07.029

Full Text Available
A high-performance finite-volume algorithm for solving partial differential equations governing compressible viscous flows on structured grids

Journal Article Guzik, S. M. ; Gao, X. ; Olschanowsky, C. - Computers and Mathematics with Applications (Oxford)

Here this work focuses on the development of a high-performance fourth-order finite-volume method to solve the nonlinear partial differential equations governing the compressible Navier–Stokes equations on a Cartesian grid with adaptive mesh refinement. The novelty of the present study is to introduce the loop chaining concept to this complex fourth-order fluid dynamics algorithm for significant improvement in code performance on parallel machines. Specific operations involved in the algorithm include the finite-volume formulation of fourth-order spatial discretization stencils and optimal inter-loop parallelization strategies. Numerical fluxes of the Navier–Stokes equations comprise the hyperbolic (inviscid) and elliptic (viscous) components. The hyperbolic flux ismore »« less
Cited by 8
https://doi.org/10.1016/j.camwa.2016.08.004
A stable finite difference method for the elastic wave equation on complex geometries with free surfaces

Journal Article Appelo, D ; Petersson, N A - Communications in Computational Physics, vol. 5, no. 1, January 1, 2009, pp. 84-107

The isotropic elastic wave equation governs the propagation of seismic waves caused by earthquakes and other seismic events. It also governs the propagation of waves in solid material structures and devices, such as gas pipes, wave guides, railroad rails and disc brakes. In the vast majority of wave propagation problems arising in seismology and solid mechanics there are free surfaces. These free surfaces have, in general, complicated shapes and are rarely flat. Another feature, characterizing problems arising in these areas, is the strong heterogeneity of the media, in which the problems are posed. For example, on the characteristic length scalesmore »« less
Full Text Available
Field theory and a structure-preserving geometric particle-in-cell algorithm for drift wave instability and turbulence

Journal Article Xiao, Jianyuan ; Qin, Hong - Nuclear Fusion

A field theory and the associated structure-preserving geometric particle-in-cell (PIC) algorithm are developed to study low frequency electrostatic perturbations with fully kinetic ions and adiabatic electrons in magnetized plasmas. The algorithm is constructed by geometrically discretizing the field theory using discrete exterior calculus, high-order Whitney interpolation forms, and non-canonical Hamiltonian splitting method. The discretization preserves the non-canonical symplectic structure of the particle-field system, as well as the electromagnetic gauge symmetry. As a result, the algorithm is charge-conserving and possesses long-term conservation properties. Because drift wave turbulence and anomalous transport intrinsically involve multi time-scales, simulation studies using fully kinetic particles demandmore »« less
Cited by 10
https://doi.org/10.1088/1741-4326/ab38dc

Full Text Available

Similar Records