Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Order”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study

The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Switching interlayer magnetic order in bilayer CrI 3 by stacking reversal

CrI 3 , a hot two-dimensional (2D) magnet, exhibits complex magnetism depending on the number of layers and interlayer stacking patterns. For bilayer CrI 3 , the interlayer magnetism can be tuned between ferromagnetic (FM) and antiferromagnetic (AFM) order by manipulating the stacking order. However, the stacking is mostly modified through translation between the layers, while the effect of rotation between the layers on the interlayer magnetic order has not yet been fully investigated. Here, we considered three energetically stable stacking patterns R$\bar3$ with combining macron], C2/m and AA in bilayer CrI 3 , and their reversed counterparts R$\bar3$ with combining macron]-r, C2/m-r and AA-r through rotating one layer by 180° with respect to the other layer. Our first-principles calculations suggest that the interlayer magnetic ground state can be switched from AFM to FM (or FM to AFM) by reversing the stacking pattern. A detailed microscopic analysis was carried out by magnetic force theory calculations on C2/m stacking which favors AFM and C2/m-r stacking which favors FM. The interlayer magnetic interactions and the origin of the magnetic order change were revealed through specific orbital analysis. Our work demonstrates that stacking rotation can also tune the interlayer magnetism of CrI 3 and provides insight into its interlayer magnetic properties at the microscopic level.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Non-covalent planarizing interactions yield highly ordered and thermotropic liquid crystalline conjugated polymers

Controlling the multi-level assembly and morphological properties of conjugated polymers through structural manipulation has contributed significantly to the advancement of organic electronics. In this work, a redox active conjugated polymer, TPT–TT, composed of alternating 1,4-(2-thienyl)-2,5-dialkoxyphenylene (TPT) and thienothiophene (TT) units is reported with non-covalent intramolecular S⋯O and S⋯H–C interactions that induce controlled main-chain planarity and solid-state order. As confirmed by density functional theory (DFT) calculations, these intramolecular interactions influence the main chain conformation, promoting backbone planarization, while still allowing dihedral rotations at higher kinetic energies (higher temperature), and give rise to temperature-dependent aggregation properties. Thermotropic liquid crystalline (LC) behavior is confirmed by cross-polarized optical microscopy (CPOM) and closely correlated with multiple thermal transitions observed by differential scanning calorimetry (DSC). This LC behavior allows us to develop and utilize a thermal annealing treatment that results in thin films with notable long-range order, as shown by grazing-incidence X-ray diffraction (GIXD). Specifically, we identified a first LC phase, ranging from 218 °C to 107 °C, as a nematic phase featuring preferential face-on π–π stacking and edge-on lamellar stacking exhibiting a large extent of disorder and broad orientation distribution. A second LC phase is observed from 107 °C to 48 °C, as a smectic A phase featuring sharp, highly ordered out-of-plane lamellar stacking features and sharp tilted backbone stacking peaks, while the structure of a third LC phase with a transition at 48 °C remains unclear, but resembles that of the solid state at ambient temperature. Furthermore, the significance of thermal annealing is evident in the ∼3-fold enhancement of the electrical conductivity of ferric tosylate-doped annealed films reaching 55 S cm −1 . More importantly, thermally annealed TPT–TT films exhibit both a narrow distribution of charge-carrier mobilities (1.4 ± 0.1) × 10 −2 cm 2 V −1 s −1 along with a remarkable device yield of 100% in an organic field-effect transistor (OFET) configuration. This molecular design approach to obtain highly ordered conjugated polymers in the solid state affords a deeper understanding of how intramolecular interactions and repeat-unit symmetry impact liquid crystallinity, solution aggregation, solution to solid-state transformation, solid-state morphology, and ultimately device applications.

Chemistry↗
DOE OSTI2023

Direct observation of magnetic ordering induced via systematic lattice disorder in artificial rhombus spin ices

It is critical to understand the effect of lattice geometry on the order parameter of a condensed matter system, as it controls phase transitions in such systems. Artificial spin ices (ASIs) are two-dimensional lattices of Ising-like nanomagnets that provide an opportunity to explore such phenomena by lithographically controlling the lattice geometry to observe its influence on magnetic ordering and frustration effects. Here we report a systematic approach to studying the effects of disorder in rhombus ASIs generated from combinations of five vertex motifs. We investigate four geometries characterized by a geometric order parameter, with symmetries ranging from periodic to quasiperiodic to random. Lorentz transmission electron microscopy data indicates magnetic domain behavior depends on chains of strongly-coupled islands in the periodic and sixfold-twinned lattices, while the behavior of the disordered lattice is dominated by vertex motifs with large configurational degeneracy. Utilizing micromagnetic simulations, a quantitative analysis of the lattice energetics showed that the experimental rotationally-demagnetized state of the disordered ASI was closer in energy to the idealized ground state compared to other periodic and twinned ASIs. Finally, our work provides a unique pathway for using degeneracy, magnetic frustration, and order to control the magnetization behavior of designer disordered systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Uncharacteristic second order martensitic transformation in metals via epitaxial stress fields

While most phase transformations, e.g., ferroelectric or ferromagnetic, can be first or second order depending on external applied fields, martensitic transformations in metallic alloys are nearly universally first order. We demonstrate that epitaxial stress originating from the incorporation of a tailored second phase can modify the free energy landscape that governs the phase transition and change its order from first to second. High-fidelity molecular dynamics simulations show a remarkable change in the character of the martensitic transformation in Ni–Al alloys near the critical point. We observe the continuous evolution of the transformation order parameter and scaling with power-law exponents comparable to those in other ferroic transitions exhibiting critical behavior. Our theoretical work provides a foundation for recent experimental and computational results on martensites near critical points.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning for nonintrusive model order reduction of the parametric inviscid transonic flow past an airfoil

Fluid flow in the transonic regime finds relevance in aerospace engineering, particularly in the design of commercial air transportation vehicles. Computational fluid dynamics models of transonic flow for aerospace applications are computationally expensive to solve because of the high degrees of freedom as well as the coupled nature of the conservation laws. While these issues pose a bottleneck for the use of such models in aerospace design, computational costs can be significantly minimized by constructing special, structure-preserving surrogate models called reduced-order models. In this work, we propose a machine learning method to construct reduced-order models via deep neural networks and we demonstrate its ability to preserve accuracy with a significantly lower computational cost. In addition, our machine learning methodology is physics-informed and constrained through the utilization of an interpretable encoding by way of proper orthogonal decomposition. Application to the inviscid transonic flow past the RAE2822 airfoil under varying freestream Mach numbers and angles of attack, as well as airfoil shape parameters with a deforming mesh, shows that the proposed approach adapts to high-dimensional parameter variation well. Notably, the proposed framework precludes the knowledge of numerical operators utilized in the data generation phase, thereby demonstrating its potential utility in the fast exploration of design space for diverse engineering applications. Comparison against a projection-based nonintrusive model order reduction method demonstrates that the proposed approach produces comparable accuracy and yet is orders of magnitude computationally cheap to evaluate, despite being agnostic to the physics of the problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

First-principles modeling of nonlinear optical spectra in the condensed phase is highly challenging because both environment and vibronic interactions can play a large role in determining spectral shapes and excited state dynamics. In this work, we compute two dimensional electronic spectroscopy (2DES) signals based on a cumulant expansion of the energy gap fluctuation operator, with specific focus on analyzing mode mixing effects introduced by the Duschinsky rotation and the role of the third order term in the cumulant expansion for both model and realistic condensed phase systems. We show that for a harmonic model system, the third order cumulant correction captures effects introduced by a mismatch in curvatures of ground and excited state potential energy surfaces, as well as effects of mode mixing. We also demonstrate that 2DES signals can be accurately reconstructed from purely classical correlation functions using quantum correction factors. We then compute nonlinear optical spectra for the Nile red and methylene blue chromophores in solution, assessing the third order cumulant contribution for realistic systems. We show that the third order cumulant correction is strongly dependent on the treatment of the solvent environment, revealing the interplay between environmental polarization and the electronic-vibrational coupling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

2nd and 3rd order spectral energy corrections with penumbral de-blurring methodology for opacity platform used on the National Ignition Facility

The Opacity Spectrometer (OpSpec) used in the National Ignition Facility’s opacity experiments measures x-ray spectra from 0.9 to 2.1 keV from the different experimental regions: the backlight source, emission source, and the absorption region with the transmission calculated from these regions. The OpSpec designs have gone through several iterations to help improve the signal-to-noise ratio, remove alternate crystal plane reflections, and improve spectral resolution, which helps to increase the validity of the opacity measurements. However, the source spans well outside the current working spectral range, and higher-order reflections are intrinsic to the crystal, which increases the overall signal seen in the data regions. The recorded data are the convolution of 1st order transmission, higher-order reflections, and the penumbra blurring. In conclusion, this work represents the details for deconvolving the 2nd and 3rd order spectral energy corrections with a penumbral de-blurring to correct the relative measurement of x-ray intensity of different spectral energies and further analysis of datasets relevant to the opacity experiments.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree–Fock reference is also), and orbital invariant. Here, our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H 2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Room temperature colossal superparamagnetic order in aminoferrocene–graphene molecular magnets

Intensive studies are published for graphene-based molecular magnets due to their remarkable electric, thermal, and mechanical properties. However, to date, most of all produced molecular magnets are ligand based and subject to challenges regarding the stability of the ligand(s). The lack of long-range coupling limits high operating temperature and leads to a short-range magnetic order. Herein, we introduce an aminoferrocene-based graphene system with room temperature superparamagnetic behavior in the long-range magnetic order that exhibits colossal magnetocrystalline anisotropy of 8 × 10 5 and 3 × 10 7 J/m 3 in aminoferrocene and graphene-based aminoferrocene, respectively. These values are comparable to and even two orders of magnitude larger than pure iron metal. Aminoferrocene [C 10 H 11 FeN] + is synthesized by an electrophilic substitution reaction. It was then reacted with graphene oxide that was prepared by the modified Hammers method. The phase structure and functionalization of surface groups were characterized and confirmed by XRD, FT-IR, and Raman spectroscopy. To model the behavior of the aminoferrocene between two sheets of hydroxylated graphene, we have used density functional theory by placing the aminoferrocene molecule between two highly ordered hydroxylated sheets and allowing the structure to relax. Here, the strong bowing of the isolated graphene sheets suggests that the charge transfer and resulting magnetization could be strongly influenced by pressure effects. In contrast to strategies based on ligands surface attachment, our present work that uses interlayer intercalated aminoferrocene opens routes for future molecular magnets as well as the design of qubit arrays and quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Symmetry-specific characterization of bond orientation order in DNA-assembled nanoparticle lattices

Bond-orientational order in DNA-assembled nanoparticles lattices is explored with the help of recently introduced Symmetry-specific Bond Order Parameters (SymBOPs). Here, this approach provides a more sensitive analysis of local order than traditional scalar BOPs, facilitating the identification of coherent domains at the single bond level. The present study expands the method initially developed for assemblies of anisotropic particles to the isotropic ones or cases where particle orientation information is unavailable. The SymBOP analysis was applied to experiments on DNA-frame-based assembly of nanoparticle lattices. It proved highly sensitive in identifying coherent crystalline domains with different orientations, as well as detecting topological defects, such as dislocations. Furthermore, the analysis distinguishes individual sublattices within a single crystalline domain, such as pair of interpenetrating FCC lattices within a cubic diamond. The results underscore the versatility and robustness of SymBOPs in characterizing ordering phenomena, making them valuable tools for investigating structural properties in various systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Learning physics-based reduced-order models from data using nonlinear manifolds

Here we present a novel method for learning reduced-order models of dynamical systems using nonlinear manifolds. First, we learn the manifold by identifying nonlinear structure in the data through a general representation learning problem. The proposed approach is driven by embeddings of low-order polynomial form. A projection onto the nonlinear manifold reveals the algebraic structure of the reduced-space system that governs the problem of interest. The matrix operators of the reduced-order model are then inferred from the data using operator inference. Numerical experiments on a number of nonlinear problems demonstrate the generalizability of the methodology and the increase in accuracy that can be obtained over reduced-order modeling methods that employ a linear subspace approximation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter ($α$) was set to recover the exact dissociation limit of H 2 in a minimal basis set. Alternatively $α$ can be viewed as a regularization parameter, where each value of $α$ represents a valid variant of BW-s2, which we denote as BW-s2($α$). In this work, we semi-empirically optimize $α$ for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal $α$ parameter ($α$ = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2($α$) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2($α$) is iterative $\mathcal{O}$($N^5$), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Direct neutron-diffraction-based measurement of magnetic order in brownmillerite SrCoO 2.5 and La 0.5 Sr 0.5 CoO 2.5 thin films

Epitaxial cobaltites have emerged as exemplary materials for electrochemical gating, in large part due to their topotactic perovskite (P) ↔ brownmillerite (BM) transformations. SrCoO 3-δ , for example, can be cycled between metallic ferromagnetic P SrCoO 3 and insulating BM SrCoO 2.5 , realizing exceptional modulation of electronic, thermal, and optical properties. It is often presumed that such cycling also generates ferromagnetic–antiferromagnetic (F-AF) modulation due to the G-type AF order in bulk SrCoO 2.5 . Little is understood about magnetism in thin-film BM SrCoO 2.5 , however, meaning that the true magnetic property modulation is unclear. We address this here through a neutron diffraction study of BM La 1-x Sr x CoO 2.5 films at x = 0.5 and 1.0. Lightly compressively strained SrCoO 2.5 films are shown to retain G-type AF order, albeit with suppressed Néel temperature (~340 K). Of high interest for AF spintronics, room-temperature F–AF cycling is thus possible across the SrCoO 3-δ P ↔ BM transformation. At x = 0.5, however, BM La 0.5 Sr 0.5 CoO 2.5 films are found to exhibit no detectable G-type AF order but instead weak F order (Curie temperature ~115 K), unveiling a La 0.5 Sr 0.5 CoO 3-δ phase diagram with two distinct F phases. These results thus uncover new, unanticipated magnetic phase behavior in these materials, in addition to being directly relevant to cobaltite-based magnetoionics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of improved higher-order correction for the NIF opacity spectrometer

X-ray opacity measurements on the National Ignition Facility (NIF) are in the process of reproducing earlier measurements from the Sandia Z Facility, in particular for oxygen and iron plasmas. These measurements have the potential to revise our understanding of the “solar problem” and of the hot degenerate Q class white dwarf structure by probing plasma conditions near the base of their convection zones. Accurate opacity measurements using soft x-ray Bragg crystal spectrometers require correction for higher-order diffraction effects. Extending prior work in this area [Dutra et al., Review of Scientific Instruments 93, 113527 (2022)], we have developed a new method to remove higher-order spectral components from NIF opacity spectrometer data. Furthermore, by modeling absorption and backlighting continuum spectra and subtracting the second- and third-order components from the measured data, we are able to perform this correction while avoiding imprinting first-order model line features onto the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗