Search

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\"odinger 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., Comment: 5 pages, 3 figures

From among the waves whose dynamics are governed by the nonlinear Schr\"odinger (NLS) equation, we find a robust, spatiotemporally disordered family, in which waves initialized with increasing amplitudes, on average, over long time scales, effectively evolve as ever more weakly coupled collections of plane waves. In particular, the relative amount of energy contained in their coupling decays to zero with increasing wave amplitude.

This paper describes a novel partitioned algorithm for fluid-structure interaction (FSI) problems that couples the motion of rigid bodies and incompressible flow. This is the first partitioned algorithm that remains stable and second-order accurate, without sub-time-step iterations, for very light, and even zero-mass, bodies in three dimensions. This new added-mass partitioned (AMP) algorithm extends the previous developments in [1, 2] by generalizing the added-damping tensors to account for arbitrary three-dimensional rotations, and by employing a general quadrature for the surface integral over a rigid body to derive the discrete AMP interface condition for the fluid pressure. Stability analyses for two three-dimensional model problems show that the algorithm remains stable for bodies of any mass when applied to the relevant model problems. The resulting AMP algorithm is implemented in parallel using a moving composite grid framework to treat one or more rigid bodies in complex three-dimensional configurations. The new three-dimensional algorithm is verified and validated though several benchmark problems, including the motion of a sphere in a viscous incompressible fluid and the interaction of a bi-leaflet mechanical heart valve and a pulsating fluid. Numerical simulations confirm the predictions of the stability analysis even for complex problems, and show that the AMP algorithm remains stable, without sub-iterations, for light and even zero-mass three-dimensional rigid bodies of general shape. These benchmark problems are further used to examine the parallel performance of the algorithm and to investigate the conditioning of the linear system for the pressure including the newly derived AMP interface conditions., Comment: 47 pages, 25 figures

A stable partitioned algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This {\em added-mass partitioned} (AMP) algorithm remains stable, without sub-iterations, for light and even zero mass rigid bodies when added-mass and viscous added-damping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Added-damping effects due to the viscous shear forces on the body are treated by inclusion of added-damping tensors that are derived through a linearization of the integrals defining the force and torque. Added-damping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this second part of a two-part series, the general formulation of the AMP scheme is presented including the form of the AMP interface conditions and added-damping tensors for general geometries. A fully second-order accurate implementation of the AMP scheme is developed in two dimensions based on a fractional-step method for the incompressible Navier-Stokes equations using finite difference methods and overlapping grids to handle the moving geometry. The numerical scheme is verified on a number of difficult benchmark problems.

A stable partitioned algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This {\em added-mass partitioned} (AMP) algorithm remains stable, without sub-iterations, for light and even zero mass rigid bodies when added-mass and viscous added-damping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added-mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Added-damping effects due to the viscous shear forces on the body are treated by inclusion of added-damping tensors that are derived through a linearization of the integrals defining the force and torque. Added-damping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this first part of a two-part series, the properties of the AMP scheme are motivated and evaluated through the development and analysis of some model problems. The analysis shows when and why the traditional partitioned scheme becomes unstable due to either added-mass or added-damping effects. The analysis also identifies the proper form of the added-damping which depends on the discrete time-step and the grid-spacing normal to the rigid body. The results of the analysis are confirmed with numerical simulations that also demonstrate a second-order accurate implementation of the AMP scheme.

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., Comment: 44 pages, 11 figures

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.

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., Comment: 22 pages

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., Comment: 37 pages

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., Comment: 16 pages, 15 figures, accepted for publication in Phys. Plasmas

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., Comment: 32 pages, 20 figures

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., Comment: 19 pages, 4 figure

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.

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.