Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multispecies”

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 19 records

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

A Scaling Study for Incompressible Multispecies Solver in Vertex-CFD

Multispecies incompressible flows occur widely in engineering and environmental applications, such as chemical reactors, fuel cells, ocean mixing, and biomedical systems. However, accurately resolving the complex transport and mixing phenomena associated with multiple interacting species remains computationally challenging, especially for large-scale problems. In this study, we present a robust, high-performance computing--enabled multispecies incompressible Navier–Stokes solver integrated within the Vertex-CFD framework. Our solver employs a fully coupled, implicit, finite element--based formulation that accurately captures the advection, diffusion, and interaction of multiple species in incompressible flows by leveraging the Kokkos library for parallel computing to achieve high computational efficiency. For pressure coupling, the entropically damped artificial compressibility method is utilized. We validated the solver against canonical test cases, including multispecies advection, diffusion, and Bateman systems; the results demonstrate second- and third-order spatial accuracy and consistent convergence. Additionally, we demonstrated the strong and weak scaling study results obtained on the leadership-class high-performance computing system, Frontier at Oak Ridge National Laboratory.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗

Multispecies Ion Acceleration in 3D Magnetic Reconnection with Hybrid-Kinetic Simulations

Magnetic reconnection drives multispecies particle acceleration broadly in space and astrophysics. Here we perform the first 3D hybrid simulations (fluid electrons, kinetic ions) that contain sufficient scale separation to produce nonthermal heavy-ion acceleration, with fragmented flux ropes critical for accelerating all species. We demonstrate the acceleration of all ion species (up to Fe) into power-law spectra with similar indices, by a common Fermi acceleration mechanism. The upstream ion velocities influence the first Fermi reflection for injection. The subsequent onsets of Fermi acceleration are delayed for ions with lower charge-mass ratios (Q/M), until growing flux ropes magnetize them. This leads to a species-dependent maximum energy/nucleon ∝(Q/M) α . These findings are consistent with in situ observations in reconnection regions, suggesting Fermi acceleration as the dominant multispecies ion acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct Evidence of Multispecies Hydrodynamics in Ignition-Scale Hohlraums

A targeted experiment at the National Ignition Facility (NIF) confirms the presence of multispecies hydrodynamics in inertial confinement fusion hohlraums relevant to ignition. Here, the effects are identified by filling the gold hohlraum with a deuterium-tritium (DT) gas mixture instead of helium. As the hohlraum is heated by the NIF lasers, it implodes inward, compressing and heating the DT, which leads to fusion. The resulting DT-fusion neutrons are measured in space, time, yield, angle, and energy. A distinct, peaked, triangular shape in the radial neutron emission profile provides evidence of a “leaky piston” effect caused by the interpenetration of DT into the expanding gold. This process reduces the reversibility of DT compression and decreases neutron generation on the temporal trailing edge compared to the leading edge. These results are well described by multispecies hydrodynamics simulations, which reproduce the observed spatial and temporal features, as well as the total neutron yield and angularly resolved energy spectra. In contrast, conventional simulations that use only single-species physics fail to match the experimental data. They overpredict the neutron yield, produce a flatter-than-observed spatial profile, and show excessive emission at late times compared to the experimental measurements.

Higginson, Drew P. [Lawrence Livermore National La↗

Global weak solutions for a nonlocal multispecies Fokker–Planck–Landau system

The global-in-time existence of weak solutions to a spatially homogeneous multispecies Fokker–Planck–Landau system for plasmas in the three-dimensional whole space is shown. The Fokker–Planck–Landau system is a simplification of the Landau equations assuming a linearized, velocity-independent, and isotropic kernel. The resulting equations depend nonlocally and nonlinearly on the moments of the distribution functions via the multispecies local Maxwellians. Furthermore, the existence proof is based on a three-level approximation scheme, energy and entropy estimates, as well as compactness results, and it holds for both soft and hard potentials.

97 MATHEMATICS AND COMPUTING↗

The effects of multispecies Hohlraum walls on stimulated Brillouin scattering, Hohlraum dynamics, and beam propagation

We report experiments and simulations have been conducted to investigate the efficacy of Ta 2 O 5 -lined Hohlraum walls at reducing stimulated Brillouin backscattering (SBS) as well as any subsequent effects on the Hohlraum dynamics and capsule implosions in indirect drive experiments at the National Ignition Facility. Using a 1.1 MJ 400 TW, 351 nm, shaped laser pulse, we measure a 5× reduction in SBS power in the peak of the pulse from the wall on the outer 50° cone beams. The SBS spectrum indicates a reduction in the high-Z spectral signature when using multispecies wall materials. Detailed hydrodynamic simulations were performed using different heat conduction models with flux limiters. Additional simulations were performed on the plasma maps using the 3D parallel paraxial code pF3D to compare backscatter powers between the pure Au and Ta 2 O 5 -lined Hohlraums. Further analysis, using hydrodynamically equivalent plasmas, shows that the SBS reduction is clearly a result of the added ion Landau damping caused by the oxygen ions and not from differences in plasma conditions. The experimental and simulation results also show an increase in the wall plasma expansion when using the Ta 2 O 5 liner leading to a 70% more oblate implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of a new multispecies tumor growth model coupling 3D phase-fields with a 1D vascular network

In this study, we present and analyze a mathematical model for tumor growth incorporating ECM erosion, interstitial flow, and the effect of vascular flow and nutrient transport. The model is of phase-field or diffused-interface type in which multiple phases of cell species and other constituents are separated by smooth evolving interfaces. The model involves a mesoscale version of Darcy’s law to capture the flow mechanism in the tissue matrix. Modeling flow and transport processes in the vasculature supplying the healthy and cancerous tissue, one-dimensional (1D) equations are considered. Since the models governing the transport and flow processes are defined together with cell species models on a three-dimensional (3D) domain, we obtain a 3D–1D coupled model.

3D–1D coupled blood flow models↗

A multispecies turbulence model for the mixing and de-mixing of miscible fluids

A de-mix capable model for turbulence in compressible, variable density flows is proposed. The model is based on the Besnard-Harlow-Rauenzahn (BHR) family of models (Besnard D, Harlow F, Rauenzahn R, et al. Turbulence transport equations for variable-density turbulence and their relationship to two-field models. NM (United States): Los Alamos National Laboratory; 1992 (Technical Report LA-12303-MS), but is extended to track the evolution of the turbulent fluxes and fluctuations in the material mass fractions for each species present. The new evolution equations are introduced without requiring additional closures or new empirically tuned coefficients relative to previous BHR models, and are shown to improve the model’s ability to reproduce the behaviour of simulations containing mixing layers that are constrained by a stabilising force. As a result, the model is tested in a range of canonical flows including Rayleigh-Taylor driven, shock driven, and shear driven turbulence, and is shown to produce reasonable agreement with simulations and experiments in these scenarios.

42 ENGINEERING↗

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗

Relaxation time approximation for a multispecies relativistic gas

We generalize a recent prescription for the relaxation time approximation for the relativistic Boltzmann equation for systems with multiple particle species at finite temperature. This is performed by adding counter-terms to the traditional Anderson-Witting ansatz for each particle species. Our approach allows for the use of momentum-dependent relaxation times and the obedience of local conservation laws regardless of the definition of the local equilibrium state. As an application, we derive the first order Chapman-Enskog corrections to the equilibrium distribution and display results for the hadron-resonance gas. We also demonstrate that our collision term ansatz obeys the second law of thermodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Design and performance of a novel low energy multispecies beamline for an antihydrogen experiment

The ALPHA Collaboration, based at the CERN Antiproton Decelerator, has recently implemented a novel beamline for low energy (≲100 eV) positron and antiproton transport between cylindrical Penning traps that have strong axial magnetic fields. Here, we describe how a combination of semianalytical and numerical calculations was used to optimize the layout and design of this beamline. Using experimental measurements taken during the initial commissioning of the instrument, we evaluate its performance and validate the models used for its development. By combining data from a range of sources, we show that the beamline has a high transfer efficiency and estimate that the percentage of particles captured in the experiments from each bunch is (78 ± 3)% for up to 10 5 antiprotons and (71 ± 5)% for bunches of up to 10 7 positrons.

43 PARTICLE ACCELERATORS↗

Structural Basis of Ligand Selectivity by a Bacterial Adhesin Lectin Involved in Multispecies Biofilm Formation

Carbohydrate recognition by lectins governs critical host-microbe interactions. MpPA14 (Marinomonas primoryensis PA14 domain) lectin is a domain of a 1.5-MDa adhesin responsible for a symbiotic bacterium-diatom interaction in Antarctica. Here, we show that MpPA14 binds various monosaccharides, with l-fucose and N-acetylglucosamine being the strongest ligands (dissociation constant [K d ], ~150 μM). High-resolution structures of MpPA14 with 15 different sugars bound elucidated the molecular basis for the lectin’s apparent binding promiscuity but underlying selectivity. MpPA14 mediates strong Ca 2+ -dependent interactions with the 3,4-diols of l-fucopyranose and glucopyranoses, and it binds other sugars via their specific minor isomers. Thus, MpPA14 only binds polysaccharides like branched glucans and fucoidans with these free end groups. Consistent with our findings, adhesion of MpPA14 to diatom cells was selectively blocked by l-fucose, but not by N-acetyl galactosamine. The MpPA14 lectin homolog present in a Vibrio cholerae adhesin was produced and was shown to have the same sugar binding preferences as MpPA14. The pathogen’s lectin was unable to effectively bind the diatom in the presence of fucose, thus demonstrating the antiadhesion strategy of blocking infection via ligand-based antagonists.

59 BASIC BIOLOGICAL SCIENCES↗

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Abstract Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

Microbiology↗