Researcher profile

J. W. Banks

J. W. Banks contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Direct verification of the kinetic description of wave turbulence for finite-size systems dominated by interactions among groups of 6 waves

The present work considers systems whose dynamics are governed by the nonlinear interactions among groups of 6 nonlinear waves, such as those described by the unforced quintic nonlinear Schrödinger equation. Specific parameter regimes in which ensemble-averaged dynamics of such systems with finite size are accurately described by a wave kinetic equation, as used in wave turbulence theory, are theoretically predicted. In addition, the underlying reasons that the wave kinetic equation may be a poor predictor of wave dynamics outside these regimes are also discussed. These theoretical predictions are directly verified by comparing ensemble averages of solutions to the dynamical equation with corresponding solutions of the wave kinetic equation.

preprint2016arXiv

Vlasov Simulations of Electron-Ion Collision Effects on Damping of Electron Plasma Waves

Collisional effects can play an essential role in the dynamics of plasma waves by setting a minimum damping rate and by interfering with wave-particle resonances. Kinetic simulations of the effects of electron-ion pitch angle scattering on Electron Plasma Waves (EPWs) are presented here. In particular, the effects of such collisions on the frequency and damping of small-amplitude EPWs for a range of collision rates and wave phase velocities are computed and compared with theory. Both the Vlasov simulations and linear kinetic theory find the direct contribution of electron-ion collisions to wave damping is about a factor of two smaller than is obtained from linearized fluid theory. To our knowledge, this simple result has not been published before. Simulations have been carried out using a grid-based (Vlasov) approach, based on a high-order conservative finite difference method for discretizing the Fokker-Planck equation describing the evolution of the electron distribution function. Details of the implementation of the collision operator within this framework are presented. Such a grid-based approach, which is not subject to numerical noise, is of particular interest for the accurate measurements of the wave damping rates.

preprint2015arXiv

A stable partitioned FSI algorithm for incompressible flow and deforming beams

An added-mass partitioned (AMP) algorithm is described for solving fluid-structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier-Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler-Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet-Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.

preprint2013arXiv

A stable FSI algorithm for light rigid bodies in compressible flow

In this article we describe a stable partitioned algorithm that overcomes the added mass instability arising in fluid-structure interactions of light rigid bodies and inviscid compressible flow. The new algorithm is stable even for bodies with zero mass and zero moments of inertia. The approach is based on a local characteristic projection of the force on the rigid body and is a natural extension of the recently developed algorithm for coupling compressible flow and deformable bodies. Normal mode analysis is used to prove the stability of the approximation for a one-dimensional model problem and numerical computations confirm these results. In multiple space dimensions the approach naturally reveals the form of the added mass tensors in the equations governing the motion of the rigid body. These tensors, which depend on certain surface integrals of the fluid impedance, couple the translational and angular velocities of the body. Numerical results in two space dimensions, based on the use of moving overlapping grids and adaptive mesh refinement, demonstrate the behavior and efficacy of the new scheme. These results include the simulation of the difficult problem of a shock impacting an ellipse of zero mass.

preprint2013arXiv

A stable partitioned FSI algorithm for incompressible flow and elastic solids

A stable partitioned algorithm for coupling incompressible flows with compressible elastic solids is described. This added-mass partitioned (AMP) scheme requires no sub-iterations, can be made fully second- or higher-order accurate, and remains stable even in the presence of strong added-mass effects. The approach extends the scheme of Banks et al. [1,2] for compressible flow, and uses Robin (mixed) boundary conditions with the fluid and solid solvers at the interface. The AMP Robin conditions are derived from a local characteristic decomposition in the solid at the interface. Two forms of the Robin conditions are derived depending on whether the fluid equations are advanced with a fractional-step method or not. A normal mode analysis for a discretization of an FSI model problem is performed to show that the new AMP algorithm is stable for any ratio of the solid and fluid densities, including the case of very light solids when the added-mass effects are large. In contrast, it is shown that a traditional partitioned algorithm involving a Dirichlet-Neumann coupling for the same FSI problem is formally unconditionally unstable for any ratio of densities. Exact traveling wave solutions are derived for three FSI model problems of increasing complexity, and these solutions are used to verify the stability and accuracy of the corresponding numerical results obtained from the AMP algorithm for the cases of light, medium and heavy solids.

preprint2013arXiv

A stable partitioned FSI algorithm for incompressible flow and structural shells

A stable partitioned algorithm for fluid-structure interaction (FSI) problems that couple viscous incompressible flow with structural shells or beams is described. This added-mass partitioned (AMP) scheme uses Robin (mixed) interface conditions for the pressure and velocity in the fluid that are derived directly from the governing equations. The AMP scheme is stable even for very light structures, requires no sub-iterations, and can be made fully second-order, or higher-order, accurate. The new scheme is shown to be stable through the analysis of a model problem. Exact traveling wave solutions for three FSI model problems are derived. Numerical results for a linearized FSI problem in two-dimensions, using these exact solutions, demonstrate that the scheme is stable and accurate, even for very light structures.

preprint2013arXiv

Block-Structured Adaptive Mesh Refinement Algorithms for Vlasov Simulation

Direct discretization of continuum kinetic equations, like the Vlasov equation, are under-utilized because the distribution function generally exists in a high-dimensional (>3D) space and computational cost increases geometrically with dimension. We propose to use high-order finite-volume techniques with block-structured adaptive mesh refinement (AMR) to reduce the computational cost. The primary complication comes from a solution state comprised of variables of different dimensions. We develop the algorithms required to extend standard single-dimension block structured AMR to the multi-dimension case. Specifically, algorithms for reduction and injection operations that transfer data between mesh hierarchies of different dimensions are explained in detail. In addition, modifications to the basic AMR algorithm that enable the use of high-order spatial and temporal discretizations are discussed. Preliminary results for a standard 1D+1V Vlasov-Poisson test problem are presented. Results indicate that there is potential for significant savings for some classes of Vlasov problems.

preprint2013arXiv

Richardson Extrapolation for Linearly Degenerate Discontinuities

In this paper we investigate the use of Richardson extrapolation to estimate the convergence rates for numerical solutions to advection problems involving discontinuities. We use modified equation analysis to describe the expectation of the approach. In general, the results do not agree with a-priori estimates of the convergence rates. However, we identify one particular use case where Richardson extrapolation does yield the proper result. We then demonstrate this result using a number of numerical examples.

preprint2012arXiv

A Note on the Convergence of the Godunov Method for Impact Problems

This paper identifies a new pathology that can be found for numerical simulations of nonlinear conservation law systems. Many of the difficulties already identified in the literature (rarefaction shocks, carbuncle phenomena, slowly moving shocks, wall heating, etc) can be traced to insufficient numerical dissipation, and the current case is no different. However, the details of the case we study here are somewhat unique in that the solution which is found by the numerics is very weak and can fail to have a derivative anywhere in the post-shock region.

preprint2012arXiv

Threshold for Electron Trapping Nonlinearity in Langmuir Waves

We assess when electron trapping nonlinearity is expected to be important in Langmuir waves. The basic criterion is that the inverse of the detrapping rate nu_d of electrons in the trapping region of velocity space must exceed the bounce period of deeply-trapped electrons, tau_B = (n_e/delta n)^{1/2} 2pi/omega_pe. A unitless figure of merit, the "bounce number" N_B = 1/(nu_d tau_B), encapsulates this condition and defines a trapping threshold amplitude for which N_B=1. The detrapping rate is found for convective loss (transverse and longitudinal) out of a spatially finite Langmuir wave. Simulations of driven waves with a finite transverse profile, using the 2D-2V Vlasov code Loki, show trapping nonlinearity increases continuously with N_B for transverse loss, and is significant for N_B ~ 1. The detrapping rate due to Coulomb collisions (both electron-electron and electron-ion) is also found, with pitch-angle scattering and parallel drag and diffusion treated in a unified manner. A simple way to combine convective and collisional detrapping is given. Application to underdense plasma conditions in inertial confinement fusion targets is presented. The results show that convective transverse loss is usually the most potent detrapping process in a single f/8 laser speckle. For typical plasma and laser conditions on the inner laser cones of the National Ignition Facility, local reflectivities ~3% are estimated to produce significant trapping effects.