Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space charge field”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Efficient approach to kinetic simulation in the inner magnetically insulated transmission line on Z

This project explores the idea of performing kinetic numerical simulations in the Z inner magnetically insulated transmission line (inner MITL) by reduced physics models such as a guiding center drift kinetic approximation for particles and electrostatic and magnetostatic approximation for the fields. The basic problem explored herein is the generation, formation, and evolution of vortices by electron space charge limited (SCL) emission. The results indicate that for relevant to Z values of peak current and pulse length, these approximations are excellent, while also providing tens to hundreds of times reduction in the computational load. The benefits could be enormous: Implementation of these reduced physics models in present particle-in-cell (PIC) codes could enable them to be routinely used for experimental design while still capturing essential non-thermal (kinetic) physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A method to determine the electric field of liquid argon time projection chambers using a UV laser system and its application in MicroBooNE

Liquid argon time projection chambers (LArTPCs) are now a standard detector technology for making accelerator neutrino measurements, due to their high material density, precise tracking, and calorimetric capabilities. An electric field (E-field) is required in such detectors to 2020 JINST 15 P07010 drift ionization electrons to the anode where they are collected. The E-field of a TPC is often approximated to be uniform between the anode and the cathode planes. However, significant distortions can appear from effects such as mechanical deformations, electrode failures, or the accumulation of space charge generated by cosmic rays. The latter effect is particularly relevant for detectors placed near the Earth’s surface and with large drift distances and long drift time. To determine the E-field in situ, an ultraviolet (UV) laser system is installed in the MicroBooNE experiment at Fermi National Accelerator Laboratory. The purpose of this system is to provide precise measurements of the E-field, and to make it possible to correct for 3D spatial distortions due to E-field non-uniformities. Here we describe the methodology developed for deriving spatial distortions, the drift velocity and the E-field from UV-laser measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms

Non-thermal advanced fuel fusion trades the requirement of a large amount of recirculating tritium in the system for that of large recirculating power. Phase space engineering technologies utilizing externally injected electromagnetic fields can be applied to meet the challenge of maintaining non-thermal particle distributions at a reasonable cost. The physical processes of the phase space engineering are studied from a theoretical and algorithmic perspective. It is emphasized that the operational space of phase space engineering is limited by the underpinning symplectic dynamics of charged particles. The phase space incompressibility according to the Liouville theorem is just one of many constraints, and Gromov's non-squeezing theorem determines the minimum footprint of the charged particles on every conjugate phase space plane. In this sense and level of sophistication, the mathematical abstraction of phase space engineering is symplectic topology. To simulate the processes of phase space engineering, such as the Maxwell demon and electromagnetic energy extraction, and to accurately calculate the minimum footprints of charged particles, recently developed structure-preserving geometric algorithms can be used. The family of algorithms conserves exactly, on discretized spacetime, symplecticity and thus incompressibility, non-squeezability, and symplectic capacities. The algorithms apply to the dynamics of charged particles under the influence of external electromagnetic fields as well as the charged particle–electromagnetic field system governed by the Vlasov–Maxwell equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Designing a boron nitride polyethylene composite for shielding neutrons

Neutrons are encountered in many different fields, including condensed matter physics, space exploration, nuclear power, and healthcare. Neutrons interacting with a biological target produce secondary charged particles that are damaging to human health. The most effective way to shield neutrons is to slow them to thermal energies and then capture the thermalized neutrons. These factors lead us to consider potential materials solutions for neutron shields that maximize the protection of humans while minimizing the shield mass and adapt well to modern additive manufacturing techniques. Using hexagonal boron nitride (hBN) as a capture medium and high-density polyethylene (HDPE) as a thermalization medium, we aim to design the optimal internal structure of h 10 BN/HDPE composites by minimizing the effective dose, which is a measure of the estimated radiation damage exposure for a human. Through Monte Carlo simulations in Geant4, we find that the optimal structure reduces the effective dose up to a factor of 72 over aluminum (Al) and up to a factor of 4 over HDPE; this is a significant improvement in shielding effectiveness that could dramatically reduce the radiation exposure of occupational workers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High magnetic field response of superconductivity dome in quantum artificial high−⁢𝑇 𝐶 superlattices with variable geometry

It is known that cuprate artificial high-𝑇 𝐶 superlattices (AHTS) with period 𝑑, composed of quantum wells confining interface space charge in stoichiometric Mott insulator layers (𝑆), with thickness 𝐿, at the interface with overdoped normal metallic cuprate layers (𝑁) show a superconducting dome by tuning the geometric 𝐿 over 𝑑 ratio of the SNSN superlattice with the top predicted by quantum material design engineering quantum size effects. Here we report high-field magnetotransport measurements up to 41 Tesla of AHTS across the entire superconducting dome. The results show the universal upward-concave behavior of the temperature-dependent upper critical magnetic field in low-𝑇 𝐶 samples at the rising edge and drop edge of the dome, providing strong evidence consistent with two-band superconductivity in agreement with multigap theory used for quantum design of the SNSN superlattices. The measured superconducting coherence length demonstrates that atomic-scale engineering controls not only the critical temperature but also the intrinsic pair size at Fano-Feshbach resonances physics paving the way toward next-generation quantum devices and shedding light on unconventional superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Implications of Collisions on the Spatial Profile of Electric Potential and the Space-Charge-Limited Current

The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D) 4/3 , where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D) 3/2 . Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D 2 ν 2 at the transition from a vacuum to a collisional diode for any device condition.

Anodes↗
DOE OSTI2021

Design of a compact, cryogen-free superconducting solenoid for the electron lens of the Fermilab Integrable Optics Test Accelerator (IOTA)

The proposed electron lens for Fermilab's Integrable Optics Test Accelerator (IOTA) will broaden its capabilities by enabling new research in nonlinear integrable optics, space-charge compensation, proton beam cooling, and more. The electron lens is based on a 5–10 keV, 1–2 A electron beam, shaped using a 0.7 m long, 0.8 T solenoidal magnetic field. A cryogen-free superconducting solenoid has been designed to provide this solenoidal field, taking into consideration the constraints on space, utilities, and infrastructure in the IOTA experimental hall. The solenoid is made of copper stabilized niobium-titanium conductor, conduction-cooled using 4 K closed-cycle cryocoolers. This paper describes the overall design of the solenoid encompassing its mechanical construction, current leads optimization, cryogenic thermal modeling that provides estimates of cooldown time and static/dynamic heat loads, and quench analysis.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Charge Transport and Space-Charge Formation in Cd 1 - x Zn x Te 1 - y Se y Radiation Detectors

The electron- and hole-transport properties in cadmium zinc telluride selenide (CZTS) crystals are studied using a laser-induced transient-current technique with pulsed and dc bias. The internal electric field profile and velocity of surface recombination are determined by Monte Carlo simulations of electron and hole transient currents combined with a numerical solution of the drift-diffusion equation coupled with Poisson’s equation. Electron and hole drift mobilities of μe = 830 cm 2 /Vs and μh = 40 cm 2 /Vs, respectively, are determined. We also develop a simple technique for evaluating surface recombination directly from measured current waveforms without the need for numerical simulation. The good quality of the prepared detector at pulsed bias, with electron- and hole-mobility-lifetime products of (μτ)e = 1.9 × 10 -3 cm 2 /V and (μτ)h = 1.4 × 10 -4 cm 2 /V, respectively, are observed. The formation of a positive space charge, originating from hole injection combined with a recombination level, is found. We observe a significant position dependence of the lifetime of electrons and holes in dc bias due to hole injection. The experiment is successfully fitted by a simple model dominated by a single deep recombination level with an energy of E t =E C -0.73eV; concentration of 7.3 × 10 11 cm -3 ; and electron- and hole-capture cross sections of 3.5 × 10 -14 cm 2 and 6.5 × 10 -14 cm 2 , respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Mapping valence electron distributions with multipole density formalism using 4D-STEM

Recent advancement in aberration correction and detector technology opened a door to various applications using 4D-STEM, which yields a diffraction pattern for each scanning position within a crystal unit-cell in scanning transmission electron microscopy (STEM) and generates incredible amounts of data in momentum space. Currently 4D-STEM analysis relies on the center-of-mass of the diffraction patterns in electric field and charge density mapping. It only derives the total projected charge density and is limited to phase objects, e.g. extremely thin samples. Here, we propose a new analytical method to accurately map aspherical valence electron distributions with atom-centered multipolar functions formalism using the whole 4D-STEM dataset. We demonstrate that, with the full dynamical calculations for various sample thicknesses, the method is sensitive not only to the miniscule charge transfer, but also to the atomic site symmetry and aspherical electron orbitals. The process of the refinement is much more robust and reliable than quantitative convergent beam electron diffraction.

4D-STEM↗
DOE OSTI2024

Connecting flux vacua through scalar field excursions

We show how flux vacua that differ from each other in flux quanta can be seen as different vacua in a single scalar potential of an enlarged field space, which resolves the separation by thin domain walls. This observation, which is motivated by the anti–de Sitter distance conjecture, allows one to compute distances between different vacua using the usual field-space metric. We verify for explicit examples such as scale-separated IIA flux vacua and the IIB Freund-Rubin vacua that the distance conjecture (for scalar fields) is satisfied and that the asymptotic directions in the enlarged field space are indeed hyperbolic. This enlarged field space contains the tachyon fields on the unstable D ˜ p -branes of type II string theory, which can induce the brane charges of the stable D-branes. We suggest that requiring continuous interpolations refines the cobordism conjecture and postdicts the existence of unstable D ˜ p -branes. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Low-Temperature Annealing of CdZnTeSe under Bias

We performed a gradual low-temperature annealing up to 360 K on a CdZnTeSe radiation detector equipped with gold and indium electrodes under bias at both polarities. We observed significant changes in the detector’s resistance and space-charge accumulation. This could potentially lead to the control and improvement of the electronic properties of the detector because the changes are accompanied with the reduction in the bulk dark current and surface leakage current. In this article, we present the results of a detailed study of the internal electric field and conductivity changes in CdZnTeSe detector for various annealing steps under bias taking into account different polarities during annealing and subsequent characterization. We observed that low-temperature annealing results in an increase in the barrier height at the contacts that, in general, reduces the dark current and decreases the positive space charge present in the sample compared to the pre-annealed condition.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Deep operator network surrogate for phase-field modeling of metal grain growth during solidification

A deep operator network (DeepONet) has been constructed that generates accurate representations of phase-field model simulations for evolving two dimensional metal grain morphology growing from melt. These representations serve as lower resolution, computationally efficient stand-ins for quick parameter space exploration of solutions to the the Allen-Cahn equations that dictate the phase-field model simulations. The experimental target for the phase-field model is a uranium casting system cooling a 434 g uranium charge from a maximum temperature of 1400° C at an average rate of 30° C / min , traversing the crystallographic phases of the pure metal. Experimental parameters inform the phase-field model, whose higher resolution computational model solutions are used to train the DeepONet in a given parameter space with the aim of developing a faster, more efficient method for predicting the solidifying metal's microstructure at different potential experimental values. The final DeepONet generates high accuracy, lower resolution predictions with cumulative relative approximation error over all timesteps of less than 0.5%, while ensuring solutions remain within physically feasible ranges. Further, these relative error values are comparable with other state-of-the-art DeepONet models for microstructure evolution, while significantly reducing the amount of training data required. Training a convolutional neural network simultaneously with the DeepONet, enforcing realistic values at the complex metal grain boundaries, and mathematically encoding boundary conditions into the structure of the DeepONet improved prediction accuracy and computational efficiency over a standard DeepONet model.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Generalized symmetries of the graviton

We find the set of generalized symmetries associated with the free graviton theory in four dimensions. These are generated by gauge invariant topological operators that violate Haag duality in ring-like regions. As expected from general QFT grounds, we find a set of “electric” and a dual set of “magnetic” topological operators and compute their algebra. To do so, we describe the theory using phase space gauge-invariant electric and magnetic dual variables constructed out of the curvature tensor. Electric and magnetic fields satisfy a set of constraints equivalent to the ones of a stress tensor of a 3d CFT. The constraints give place to a group $\mathbb{R}^{20}$ of topological operators that are charged under space-time symmetries. Finally, we discuss similarities and differences between linearized gravity and tensor gauge theories that have been introduced recently in the context of fractonic systems in condensed matter physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗