Search NASA⌕ Search

SEARCH · Search NASA

Results for “charged particle channeling”

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 181 records · Page 10

Transport of time-varying plasma currents by whistler wave packets

The relationship between pulsed currents and electromagnetic waves is examined in a regime characterized by electron MHD. Pulsed currents are generated by (1) collection/emission of charged particles by/from biased electrodes and (2) induction of currents by time-varying and moving magnetic fields. Pulsed currents are observed to propagate at the speed of whistler wave packets. Their field structure forms ropelike configurations which are electromagnetically force-free. Moving sources induce 'eddy' currents which excite waves and form Cerenkov-like whistler 'wings'. The radiation patterns of moving magnetic antennas and electrodynamic tethers are investigated. Nonlinear effects of large-amplitude, antenna-launched whistler pulses are observed. These involve a new modulational instability in which a channel of high conductivity which permits the wave/currents to penetrate deeply into a collisional plasma is formed.

Stenzel, R. L.↗

Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes

Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. Here, to address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.

organic↗

Imaging viscous flow of the Dirac fluid in graphene

The electron–hole plasma in charge-neutral graphene is predicted to realize a quantum critical system in which electrical transport features a universal hydrodynamic description, even at room temperature. This quantum critical ‘Dirac fluid’ is expected to have a shear viscosity close to a minimum bound, with an interparticle scattering rate saturating1 at the Planckian time, the shortest possible timescale for particles to relax. Although electrical transport measurements at finite carrier density are consistent with hydrodynamic electron flow in graphene a clear demonstration of viscous flow at the charge-neutrality point remains elusive. In this work, we directly image viscous Dirac fluid flow in graphene at room temperature by measuring the associated stray magnetic field. Nanoscale magnetic imaging is performed using quantum spin magnetometers realized with nitrogen vacancy centres in diamond. Scanning single-spin and wide-field magnetometry reveal a parabolic Poiseuille profile for electron flow in a high-mobility graphene channel near the charge-neutrality point, establishing the viscous transport of the Dirac fluid. This measurement is in contrast to the conventional uniform flow profile imaged in a metallic conductor and also in a low-mobility graphene channel. Furthermore, via combined imaging and transport measurements, we obtain viscosity and scattering rates, and observe that these quantities are comparable to the universal values expected at quantum criticality. This finding establishes a nearly ideal electron fluid in charge-neutral, high-mobility graphene at room temperature4. Our results will enable the study of hydrodynamic transport in quantum critical fluids relevant to strongly correlated electrons in high-temperature superconductors9. This work also highlights the capability of quantum spin magnetometers to probe correlated electronic phenomena at the nanoscale.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Beam particle identification and tagging of incompletely stripped heavy beams with HEIST

A challenge preventing successful inverse kinematics measurements with heavy nuclei that are not fully stripped is identifying and tagging the beam particles. For this purpose, the HEavy ISotope Tagger (HEIST) has been developed. HEIST utilizes two micro-channel plate timing detectors to measure the time-of-flight, a multi-sampling ion chamber to measure energy loss, and a high-purity germanium detector to identify isomer decays and calibrate the isotope identification system. HEIST has successfully identified 198 Pb and other nearby nuclei at energies of about 75 MeV/A. In the experiment discussed, a typical cut containing 89% of all 198 Pb 80+ in the beam had a purity of 86%. We examine the issues of charge state contamination. Here, the observed charge state populations of these ions are presented and, using an adjusted beam energy, are well described by the charge state model GLOBAL.

47 OTHER INSTRUMENTATION↗

Search for charge non-conservation and Pauli exclusion principle violation with the M AJORANA D EMONSTRATOR

Charge conservation and the Pauli exclusion principle result from fundamental symmetries in the standard model of particle physics, and are typically taken as axiomatic. High-precision tests for small violations of these symmetries could point to new physics. Here we consider three models for violation of these processes, which would produce detectable ionization in the high-purity germanium detectors of the M AJORANA D EMONSTRATOR experiment. Using a 37.5 kg yr exposure, we report a lower limit on the electron mean lifetime, improving the previous best limit for the e → v e $\overline{v_e}$ v e decay channel by more than an order of magnitude. Here, we also present searches for two types of violation of the Pauli exclusion principle, setting limits on the probability of an electron to be found in a symmetric quantum state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Review of the High Performance Antiproton Trap (HiPAT) Experiment at the Marshall Space Flight Center

The significant energy density of matter-antimatter annihilation is attractive to the designers of future space propulsion systems, with the potential to offer a highly compact source of power. Many propulsion concepts exist that could take advantage of matter-antimatter reactions, and current antiproton production rates are sufficient to support basic proof-of-principle evaluation of technology associated with antimatter- derived propulsion. One enabling technology for such experiments is portable storage of low energy antiprotons, allowing antiprotons to be trapped, stored, and transported for use at an experimental facility. To address this need, the Marshall Space Flight Center's Propulsion Research Center is developing a storage system referred to as the High Performance Antiproton Trap (HiPAT) with a design goal of containing 10(exp 12) particles for up to 18 days. The HiPAT makes use of an electromagnetic system (Penning- Malmberg design) consisting of a 4 Telsa superconductor, high voltage electrode structure, radio frequency (RF) network, and ultra high vacuum system. To evaluate the system normal matter sources (both electron guns and ion sources) are used to generate charged particles. The electron beams ionize gas within the trapping region producing ions in situ, whereas the ion sources produce the particles external to the trapping region and required dynamic capture. A wide range of experiments has been performed examining factors such as ion storage lifetimes, effect of RF energy on storage lifetime, and ability to routinely perform dynamic ion capture. Current efforts have been focused on improving the FW rotating wall system to permit longer storage times and non-destructive diagnostics of stored ions. Typical particle detection is performed by extracting trapped ions from HiPAT and destructively colliding them with a micro-channel plate detector (providing number and energy information). This improved RF system has been used to detect various plasma modes for both electron and ion plasmas in the two traps at MSFC, including axial, cyclotron, and diocotron modes. New diagnostics are also being added to HiPAT to measure the axial density distribution of the trapped cloud to match measured RF plasma modes to plasma conditions.

Pearson, J. B.↗

Selecting Charged Current Muon Neutrino Interactions on Argon with the Short-Baseline Near Detector

The uncertainties associated with neutrino-nucleus cross section modelling need to be reduced in order to measure neutrino oscillation parameters, the mass hierarchy and CP violation to a high degree of precision. SBND, a 112 ton liquid argon time projection chamber, will be the near detector of the short-baseline neutrino program at Fermilab. Once data taking begins in 2022, it will provide flux constraints for sterile neutrino searches and produce world leading neutrino-argon cross sections with up to seven million neutrino events in three years. This thesis describes the simulation, calibration and reconstruction tools developed to measure charged current muon neutrino interactions.Backgrounds from cosmic ray muons, neutrino interactions outside of the detector and different neutrino interaction channels were identified and removed as part of this process.The signal events were selected with purity of 90% and an efficiency of 67%, providing a strong starting point for meas urements of exclusive channels. The expected statistical and systematic uncertainties on rate predictions and cross section measurements at the full 6.6×10$^{20}$ proton on target exposure were evaluated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Jet fragmentation function and groomed substructure of bottom quark jets in proton-proton collisions at 5.02 TeV

A measurement of the substructure of bottom quark jets (b jets) in proton-proton (pp) collisions is presented. The measurement uses data collected in pp collisions at $\sqrt{s}=5.02$ TeV, with a low number of simultaneous interactions per bunch crossing, recorded by the CMS experiment in 2017, corresponding to an integrated luminosity of 301 pb −1 . An algorithm to identify and cluster the charged decay daughters of b hadrons is developed for this analysis, which facilitates the exposure of the gluon radiation pattern of b jets using iterative Cambridge-Aachen declustering. The soft-drop-groomed jet radius, R g , and momentum balance, z g , of b quark jets are presented. These observables can be used to test perturbative quantum chromodynamics predictions that account for mass effects. Because the b hadron is partially reconstructed from its charged decay daughters, only charged particles are used for the jet substructure studies. In addition, a jet fragmentation function, z b,ch , is measured, which is defined as the distribution of the ratio of the transverse momentum (p T ) of the partially reconstructed b hadron with respect to the charged-particle component of the jet p T . The substructure variable distributions are unfolded to the charged-particle level. The b jet substructure is compared to the substructure of jets in an inclusive jet sample that is dominated by light-quark and gluon jets in order to assess the role of the b quark mass. A strong suppression of emissions at small R g values is observed for b jets when compared to inclusive jets, consistent with the dead-cone effect. The measurement is also compared with theoretical predictions from Monte Carlo event generators. This is the first substructure measurement of b jets that clusters together the b hadron decay daughters independent of the b hadron species and decay channel.

boosted jets↗

Superconductivity in a strange metal

In a gas of charged particles with the density n, mass m, and charge e, the electrical conductivity σ is given by the Drude formula: σ = ne 2 /mΓ, where Γ is the scattering rate. Standard metals are well described by Landau's Fermi Liquid (FL) theory, in which the electric current is carried by quasi-particles, low-energy excitations of the FL that resemble electrons with some effective mass m*. The scattering rate Γ= 1/τ= v F /l, where t is the (momentum) relaxation) time, l is the mean-free path, and v F is the Fermi velocity, can be expressed ("Matthiessen's Rule") as a sum of contributions from various scattering channels: Γ = Γ 0 + Γ el-el + Γ el-ph + …, where Γ 0 describes scattering on lattice imperfections, Γ el-el the electron-electron scattering, Γ el-ph the electron-phonon scattering, etc. Of these, Γ 0 = v F /l 0 , where l 0 is the average distance between the defects, is temperature-independent. Γ el-el should scale as T 2 because of Fermi statistics; for two electrons to scatter on one another, both must come from the "Debye shell" of the width k B T/E F , where k B is the Boltzmann constant and E F is the Fermi energy. Γ el-ph typically grows as T 5 , so we expect this to overwhelm the other channels at a high enough T. However, since l cannot be shorter than the distance between the atoms, the total Γ saturates at Mott-Ioffe-Regel (MIR) limit, roughly v F /a 0 , where v F is the Fermi velocity and a 0 is the lattice constant. The resistivity should also saturate at low T, at ρ 0 = m*v F /ne 2 l 0 T. FL theory also describes other electronic properties; e.g., it predicts that in the magnetic field B, the resistivity of the metal should increase with B 2 , because σ(B) = σ(B=0)/(1 + (ω c /Γ) 2 ), where ω c = eB/m*.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Simple Mechanism for Fractionating Oxygen Isotopes in the Solar Nebula

Lightning in the Solar Nebula is caused by the tribo-electric charging of dust grains carried by massive turbulent flows and driven by the accretion energy in the disk: it has long been one agent assumed responsible for the formation of chondrules. The degree to which charge separation can occur is dependent upon a number of factors, including the concentration of radioactive sources and the total level of ionization in the nebula, and these factors determine the maximum energy likely to be released by a single bolt. While chondrule formation requires a massive discharge, even a small lightning bolt can vaporize grains in the ionized discharge channel. Experimental studies have shown that silica, iron silicate and iron oxide grains formed from a high voltage discharge in hydrogen rich gas containing some oxygen produces solids that are enriched in O-17 and O-18 relative to the composition of the starting gas. Vaporization of silicates produces SiO, metal and free oxygen atoms in each discharge and these species will immediately begin to recondense from the hot plasma. Freshly condensed grains are incrementally enriched in heavy oxygen while the gas is enriched in O-16. Repeated evaporation and condensation of silicates in continuously occurring lightning discharges will monotonically increase the fractionation of oxygen isotopes between the O-17 and O-18 rich dust and the O-16 rich gas. The first mass independently fractionated refractory oxide particles were produced in the lab following the condensation of a flowing gas mixture containing variable amounts of hydrogen, silane, pentacarbonyl iron and oxygen that passed through a high voltage discharge powered by a Tesla coil. While the exact chemical pathway is still uncertain, the most probable reaction mechanisms involve oxidation of the growing refractory clusters by O3, OH or O atoms. This model has some interesting consequences for chemical processes in the early solar nebula. Chemical fractionation of recondensed dust evaporated via lightning discharges should be strongly time dependent. At earlier times, the accretion rate is maximal, thus driving strong turbulence, energetic grain-grain collisions, tribo-electric charging and charge separation, leading to frequent, powerful lightning discharges. As the accretion rate diminishes, turbulence decreases and lightning discharges will become both less powerful and less frequent, thus decreasing the rate of dust-gas fractionation. The most rapid increase in the formation of O-16 poor dust will occur early in nebular history. Generation of fractionated dust should be distributed throughout the inner disk. Once condensed, grain dispersal would average out any significant isotopic anomalies within the inner disk.

Nuth, Joseph A., III↗

Measurement of the polarisation of single top quarks and antiquarks produced in the $t$-channel at $ \sqrt{s} $ = 13 TeV and bounds on the $tWb$ dipole operator from the ATLAS experiment

A simultaneous measurement of the three components of the top-quark and top-antiquark polarisation vectors in t-channel single-top-quark production is presented. This analysis is based on data from proton–proton collisions at a centre-of-mass energy of 13 TeV corresponding to an integrated luminosity of 139 fb –1 , collected with the ATLAS detector at the LHC. Selected events contain exactly one isolated electron or muon, large missing transverse momentum and exactly two jets, one being b-tagged. Stringent selection requirements are applied to discriminate t-channel single-top-quark events from the background contributions. The top-quark and top-antiquark polarisation vectors are measured from the distributions of the direction cosines of the charged-lepton momentum in the top-quark rest frame. The three components of the polarisation vector for the selected top-quark event sample are P x' = 0.01 ± 0.18, P y' = –0.029 ± 0.027, P z' = 0.91 ± 0.10 and for the top-antiquark event sample they are P x' = –0.02 ± 0.20, P y' = –0.007 ± 0.051, P z' = 0.79 ± 0.16. Normalised differential cross-sections corrected to a fiducial region at the stable-particle level are presented as a function of the charged-lepton angles for top-quark and top-antiquark events inclusively and separately. These measurements are in agreement with Standard Model predictions. The angular differential cross-sections are used to derive bounds on the complex Wilson coefficient of the dimension-six O tW operator in the framework of an effective field theory. The obtained bounds are C tW ϵ [–0.9, 1.4] and C itW ϵ [–0.8, 0.2], both at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Standard versus non-standard CP phases in neutrino oscillation in matter with non-unitarity

We formulate a perturbative framework for the flavor transformation of the standard active three neutrinos but with a non-unitary flavor mixing matrix, a system which may be relevant for the leptonic unitarity test. We use the $\alpha$ parametrization of the non-unitary matrix and take its elements $\alpha_{\beta \gamma}$ ($\beta,\gamma = e,\mu,\tau$) and the ratio $\epsilon \simeq \Delta m^2_{21} / \Delta m^2_{31}$ as the small expansion parameters. Two qualitatively new features that hold in all the oscillation channels are uncovered in the probability formula obtained to first order in the expansion: (1) The phases of the complex $\alpha$ elements always come into the observable in the particular combination with the $\nu$SM CP phase $\delta$ in the form $[e^{- i \delta } \bar{\alpha}_{\mu e}, ~e^{ - i \delta} \bar{\alpha}_{\tau e}, ~\bar{\alpha}_{\tau \mu}]$ under the Particle Data Group convention of a unitary $\nu$SM mixing matrix. (2) The diagonal $\alpha$ parameters appear in particular combinations $\left( a/b - 1 \right) \alpha_{ee} + \alpha_{\mu \mu}$ and $\alpha_{\mu \mu} - \alpha_{\tau \tau}$, where $a$ and $b$ denote, respectively, the matter potential due to charged current and neutral current reactions. This property holds only in the unitary evolution part of the probability, and there is no such feature in the genuine non-unitary part, while the $\delta$–$\alpha$ parameter phase correlation exists for both. The reason for such remarkable stability of the phase correlation is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields

Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.

36 MATERIALS SCIENCE↗

Lightning and 85-GHz MCSs in the Global Tropics

Numerous observations of tropical convection show that tropical continental mesoscale convective systems (MCSs) are much more prolific lightning producers than their oceanic counterparts. Satellite-based climatologies using 85-GHz passive microwave ice-scattering signatures from the Special Sensor Microwave/Imager (SSM/I) indicate that MCSs of various size and intensity are found throughout the global tropics. In contrast, global lightning distributions show a strong land bias with an order of magnitude difference between land and ocean lightning. This is somewhat puzzling, since 85-GHz ice-scattering and the charge separation processes that lead to lightning are both thought to depend upon the existence of large graupel particles. The fact that low 85-GHz brightness temperatures are observed in tropical oceanic MCSs containing virtually no lightning leads to the postulate that tropical oceanic and tropical continental MCSs have fundamentally different hydrometeor profiles through the mixed phase region of the cloud (0 C <= T <= 20 C). Until recently, validation of this postulate has not been practicable on a global scale. Recent deployment of the Tropical Rainfall Measuring Mission (TRMM) satellite presents a unique opportunity for MCS studies. The multi-sensor instrument ensemble aboard TRMM, including a multi-channel microwave radiometer, the Lightning Imaging Sensor (LIS), and the first space-borne radar, facilitates high-resolution case studies of MCS structure throughout the global tropics. An important precursor, however, is to better understand the distribution of MCSs and lightning in the tropics. With that objective in mind, this research undertakes a systematic comparison of 85-GHz-defined MCSs and lightning over the global tropics for a full year, as an initial step toward quantifying differences between land and ocean convective systems.

Toracinta, E. Richard↗

LHC probes of the 10 TeV scale

The usual range of new particle masses, up to a few TeV, searched for at the LHC may be substantially extended if ultraheavy diquark particles exist. A diquark scalar, $S_{uu}$, that interacts perturbatively with two up quarks may be as heavy as 10 TeV and would still produce tens of spectacular events at the 14 TeV LHC. It is shown here that an ultraheavy $S_{uu}$ could be discovered through final states of very high energy in various channels, especially if the diquark can decay into other new heavy particles. Examples include cascade decays of $S_{uu}$ via a second scalar produced in pairs, which leads to two dijet resonances, or to more exotic signals with top quarks, Higgs bosons, electroweak bosons, and high-$p_T$ jets. Another possibility is that the diquark decays into a vectorlike quark of multi-TeV mass and a top or up quark. Signal events include one or two highly boosted top quarks and a Higgs boson or a $Z$, without counterparts containing top antiquarks. Similarly, direct decays of the diquark into $tj$ or $tt$ with leptonic top decays involve only positively charged leptons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Investigating high-energy proton-induced reactions on spherical nuclei: Implications for the preequilibrium exciton model

A number of accelerator-based isotope production facilities utilize $100-200$ MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed between the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer ($E_p=200$ MeV) and the Los Alamos Isotope Production Facility ($E_p=100$ MeV) to measure $^{93}$Nb(p,x) cross sections between $50-200$ MeV. The results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. The default code predictions largely failed to reproduce the measurements. Therefore, we developed a standardized procedure, which determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. Overall, the primary focus of the procedure was to determine the best parameterization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend towards a relative decrease for internal transition rates at intermediate proton energies ($E_p=20-60$ MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.

43 PARTICLE ACCELERATORS↗

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗