Abstract

Nested numerical simulations of ionospheric plasma density structures associated with nonlinear evolution of the Rayleigh-Taylor (RT) instability in equatorial spread F (ESF) are presented. The numerical implementation of the nested model uses a spatial discretization with a C grid staggering configuration where normal velocities of ions and electrons are staggered one-half grid length from the density of charged particles. The advection of charged particles is computed with a fifth order accurate in space weighted essentially nonoscillatory (WENO) scheme. The continuity equation is integrated using a third-order Runge-Kutta (RK) time integration scheme. The equation for the electric potential is solved at each time step with a multigrid method. For the limited area and nested simulations, the lateral boundary conditions are treated via implicit relaxation applied in buffer zones where the density of charged particles for each nest is relaxed to that obtained from the parent domain. The high resolution in targeted regions offered by the nested model was able to resolve secondary RT instabilities, and to improve the resolution of the primary RT bubble compared to the coarser large domain model. The computational results are validated by conducting a large domain simulation where the resolution is increased everywhere.

Original languageEnglish (US)
Article number061202
JournalJournal of Fluids Engineering, Transactions of the ASME
Volume136
Issue number6
DOIs
StatePublished - 2014

Fingerprint

Charged particles
Plasma density
Advection
Boundary conditions
Electrons
Computer simulation
Ions
Electric potential

ASJC Scopus subject areas

  • Mechanical Engineering

Cite this

@article{2efa2ec12ffc46579b6a752a8152a4d6,
title = "Multiscale nested simulations of Rayleigh-Taylor instabilities in ionospheric flows",
abstract = "Nested numerical simulations of ionospheric plasma density structures associated with nonlinear evolution of the Rayleigh-Taylor (RT) instability in equatorial spread F (ESF) are presented. The numerical implementation of the nested model uses a spatial discretization with a C grid staggering configuration where normal velocities of ions and electrons are staggered one-half grid length from the density of charged particles. The advection of charged particles is computed with a fifth order accurate in space weighted essentially nonoscillatory (WENO) scheme. The continuity equation is integrated using a third-order Runge-Kutta (RK) time integration scheme. The equation for the electric potential is solved at each time step with a multigrid method. For the limited area and nested simulations, the lateral boundary conditions are treated via implicit relaxation applied in buffer zones where the density of charged particles for each nest is relaxed to that obtained from the parent domain. The high resolution in targeted regions offered by the nested model was able to resolve secondary RT instabilities, and to improve the resolution of the primary RT bubble compared to the coarser large domain model. The computational results are validated by conducting a large domain simulation where the resolution is increased everywhere.",
author = "Alex Mahalov and Mohamed Moustaoui",
year = "2014",
doi = "10.1115/1.4025657",
language = "English (US)",
volume = "136",
journal = "Journal of Fluids Engineering, Transactions of the ASME",
issn = "0098-2202",
publisher = "American Society of Mechanical Engineers(ASME)",
number = "6",

}

TY - JOUR

T1 - Multiscale nested simulations of Rayleigh-Taylor instabilities in ionospheric flows

AU - Mahalov, Alex

AU - Moustaoui, Mohamed

PY - 2014

Y1 - 2014

N2 - Nested numerical simulations of ionospheric plasma density structures associated with nonlinear evolution of the Rayleigh-Taylor (RT) instability in equatorial spread F (ESF) are presented. The numerical implementation of the nested model uses a spatial discretization with a C grid staggering configuration where normal velocities of ions and electrons are staggered one-half grid length from the density of charged particles. The advection of charged particles is computed with a fifth order accurate in space weighted essentially nonoscillatory (WENO) scheme. The continuity equation is integrated using a third-order Runge-Kutta (RK) time integration scheme. The equation for the electric potential is solved at each time step with a multigrid method. For the limited area and nested simulations, the lateral boundary conditions are treated via implicit relaxation applied in buffer zones where the density of charged particles for each nest is relaxed to that obtained from the parent domain. The high resolution in targeted regions offered by the nested model was able to resolve secondary RT instabilities, and to improve the resolution of the primary RT bubble compared to the coarser large domain model. The computational results are validated by conducting a large domain simulation where the resolution is increased everywhere.

AB - Nested numerical simulations of ionospheric plasma density structures associated with nonlinear evolution of the Rayleigh-Taylor (RT) instability in equatorial spread F (ESF) are presented. The numerical implementation of the nested model uses a spatial discretization with a C grid staggering configuration where normal velocities of ions and electrons are staggered one-half grid length from the density of charged particles. The advection of charged particles is computed with a fifth order accurate in space weighted essentially nonoscillatory (WENO) scheme. The continuity equation is integrated using a third-order Runge-Kutta (RK) time integration scheme. The equation for the electric potential is solved at each time step with a multigrid method. For the limited area and nested simulations, the lateral boundary conditions are treated via implicit relaxation applied in buffer zones where the density of charged particles for each nest is relaxed to that obtained from the parent domain. The high resolution in targeted regions offered by the nested model was able to resolve secondary RT instabilities, and to improve the resolution of the primary RT bubble compared to the coarser large domain model. The computational results are validated by conducting a large domain simulation where the resolution is increased everywhere.

UR - http://www.scopus.com/inward/record.url?scp=84899645900&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84899645900&partnerID=8YFLogxK

U2 - 10.1115/1.4025657

DO - 10.1115/1.4025657

M3 - Article

AN - SCOPUS:84899645900

VL - 136

JO - Journal of Fluids Engineering, Transactions of the ASME

JF - Journal of Fluids Engineering, Transactions of the ASME

SN - 0098-2202

IS - 6

M1 - 061202

ER -