Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron Fluid Model”

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 595 records · Page 33

Coupling a Computational Fluid Dynamics Model to a Spacecraft Thermal System Model for the DraMS Instrument Thermal Analysis

The Dragonfly Mass Spectrometer (DraMS) is an instrument on the Dragonfly mission, which will spend 7 years in deep space cruise before landing and operating on the surface of Titan. Vacuum thermal analyses are required for deep space cruise, and convection analyses are required for the Titan surface operations. Model exchanges across multiple thermal teams are needed for all phases of the mission. For DraMS, Thermal Desktop® (TD) has been the main thermal analytical tool of choice due to its capability in modeling complex thermal systems with relatively low computational power and for its availability across thermal teams. However, TD does not have computational fluid dynamics (CFD) capability and struggles to accurately capture complex convective behavior. DraMS has fans operating in tandem and gas flow behaviors are not easily predicted due to its complex flow paths. CFD software, such as Fluent, can model and predict such complex flow behaviors, but CFD models are computationally expensive, and its workflow processes are not tailored towards simulating large and complex systems. Therefore, a coupled modeling approach was chosen for DraMS: A TD model was used for simulating all the conductive, radiative, and source terms, while a Fluent CFD model was added on, as needed, to the TD model to provide the convective boundary conditions using the System Coupling software. The coupling software allows the TD and Fluent models to communicate data and arrive at a co-solved and co-converged solution. Furthermore, Thermal Iso-value Exchange (TIE) method was developed to facilitate and improve the TD-Fluent data exchange process. This paper will discuss the analytical studies that were done to verify the accuracy and usability of the coupled approach and the challenges associated, which lead to the development of the TIE approach. DraMS thermal design and co-solved analysis results will also be discussed.

Heat transfer↗

HTMT-class Latency Tolerant Parallel Architecture for Petaflops Scale Computation

Computational Aero Sciences and other numeric intensive computation disciplines demand computing throughputs substantially greater than the Teraflops scale systems only now becoming available. The related fields of fluids, structures, thermal, combustion, and dynamic controls are among the interdisciplinary areas that in combination with sufficient resolution and advanced adaptive techniques may force performance requirements towards Petaflops. This will be especially true for compute intensive models such as Navier-Stokes are or when such system models are only part of a larger design optimization computation involving many design points. Yet recent experience with conventional MPP configurations comprising commodity processing and memory components has shown that larger scale frequently results in higher programming difficulty and lower system efficiency. While important advances in system software and algorithms techniques have had some impact on efficiency and programmability for certain classes of problems, in general it is unlikely that software alone will resolve the challenges to higher scalability. As in the past, future generations of high-end computers may require a combination of hardware architecture and system software advances to enable efficient operation at a Petaflops level. The NASA led HTMT project has engaged the talents of a broad interdisciplinary team to develop a new strategy in high-end system architecture to deliver petaflops scale computing in the 2004/5 timeframe. The Hybrid-Technology, MultiThreaded parallel computer architecture incorporates several advanced technologies in combination with an innovative dynamic adaptive scheduling mechanism to provide unprecedented performance and efficiency within practical constraints of cost, complexity, and power consumption. The emerging superconductor Rapid Single Flux Quantum electronics can operate at 100 GHz (the record is 770 GHz) and one percent of the power required by convention semiconductor logic. Wave Division Multiplexing optical communications can approach a peak per fiber bandwidth of 1 Tbps and the new Data Vortex network topology employing this technology can connect tens of thousands of ports providing a bi-section bandwidth on the order of a Petabyte per second with latencies well below 100 nanoseconds, even under heavy loads. Processor-in-Memory (PIM) technology combines logic and memory on the same chip exposing the internal bandwidth of the memory row buffers at low latency. And holographic storage photorefractive storage technologies provide high-density memory with access a thousand times faster than conventional disk technologies. Together these technologies enable a new class of shared memory system architecture with a peak performance in the range of a Petaflops but size and power requirements comparable to today's largest Teraflops scale systems. To achieve high-sustained performance, HTMT combines an advanced multithreading processor architecture with a memory-driven coarse-grained latency management strategy called "percolation", yielding high efficiency while reducing the much of the parallel programming burden. This paper will present the basic system architecture characteristics made possible through this series of advanced technologies and then give a detailed description of the new percolation approach to runtime latency management.

Sterling, Thomas↗

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

Capability in Theory, Modeling, and Validation for a Range of Innovative Fusion Concepts using High-Fidelity Moment-Kinetic Models

A computational modeling capability is created and available to the fusion community to understand and design lower-cost and innovative fusion concepts. The approach uses high- fidelity kinetic, moment-kinetic, and moment models and includes sophisticated plasma- boundary interactions. A majority of fusion-relevant simulations are performed with magnetohydrodynamic models and hybrid particle-in-cell codes, with limited-fidelity electron and kinetic physics. However, in fusion configurations like Z-pinches, field-reversed- configurations, plasma jet magneto-inertial fusion, spinning mirrors, and others, kinetic effects (both electron and ions) are critical to understand the physics and design scaling into the highly kinetic regime of a burning fusion plasma. Furthermore, as present fusion machines move towards a burning plasma regime, liquid-metal blankets are needed to handle first-wall heat- flux, reduce erosion, and eventually for energy conversion and fuel breeding. The work performed under this ARPA-E BETHE Capability Team advances the state-of-the-art in modeling and understanding plasma dynamics in fusion devices and its coupling with liquid-metal dynamics. These are critical areas of research for fusion energy to become realizable. To address these complex problems, we have leveraged and extended computational capabilities through the code, Gkeyll (developed jointly with Princeton Plasma Physics Laboratory and academic partners), for kinetic and moment modeling of fusion plasmas. The Concept Teams supported by this Capability Team include the Wisconsin High-field Axisymmetric Mirror (WHAM), Centrifugal Mirror Experiment (CFME), Plasma-Jet Magneto- Inertial Fusion (PJMIF), and solid and liquid wall plasma-material interaction studies relevant to a number of fusion concepts including Zap Energy’s Z-pinch. This software is open-source and available to the fusion community as a high-fidelity tool for the design of lower-cost fusion experiments. 3D gyrokinetic simulations of WHAM are now possible for long enough time scales to understand the evolution of interchange instabilities. 3D multi-fluid simulations of CMFE at higher Mach numbers are now possible for detailed design iterations with the goal of stability. The state-of-the-art in understanding shock formation and shock mitigation regimes in merging liners for PJMIF have been furthered by our kinetic simulations. Our novel models and frameworks studying plasma-material interaction by incorporating wall emission for various solid wall materials of relevance to pulsed and steady fusion concepts have advanced the state-of-the-art in our understanding of particle fluxes, heat fluxes, and other quantities at cathodes and anodes. The results from this work may explain discrepancies between experimental and theoretical predictions of achieved current densities in pulsed concepts such as Z-pinches. Another significant contribution of this Capability Team is the development and deployment of a novel experimental platform, LEX (Liquid Electrode eXperiment), at Virginia Tech to understand liquid metal free-surface response to electromagnetic pulses. The novel experiments along with model validation quantified the effect of different materials and sizes of liquid metal droplets on the radiative power balance of fusion plasmas for pulsed concepts. Furthermore, these experiments provided mitigation strategies for violent liquid metal response for high current pulses as would be expected in fusion regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗

Development, construction and qualification tests of the mechanical structures of the electromagnetic calorimeter of the Mu2e experiment at Fermilab

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the Charged Lepton Flavour Violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of physics beyond the Standard Model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. The calorimeter provides excellent electron identification, complementary information to aid pattern recognition and track reconstruction, and a fast calorimetric online trigger. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1340 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements.The hostile Mu2e operational conditions, in terms of radiation levels (total ionizing dose of 12 krad and a neutron fluence of 5x1010 n/cm2 @ 1 MeVeq (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on the design of the detector mechanical structures and materials choice. The support structure of the two 670 crystal matrices employs two aluminum hollow rings and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and service front-end electronics for each crystal are assembled in a unique mechanical unit inserted in a machined copper holder. The 670 units are supported by a machined plate made of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics is hosted in aluminum custom crates positioned on the external lateral surface of the two disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system.The constraints on the calorimeter mechanical structures design, the development from the conceptual design to the specifications of all the structural components, the status of components production and the components quality assurance tests are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Centrifugal acceleration of the polar wind

The effect of parallel ion acceleration associated with convection was first applied to energization of test particle polar ions by Cladis (1986). However, this effect is typically neglected in 'self-consistent' models of polar plasma outflow, apart from the fluid simulation by Swift (1990). Here we include approximations for this acceleration, which we broadly characterize as centrifugal in nature, in our time-dependent, semikinetic model of polar plasma outflow and describe the effects on the bulk parameter profiles and distribution functions of H+ and O+. For meridional convection across the pole the approximate parallel force along a polar magnetic field line may be written as F(sub cent, pole) = 1.5m(E(sub i))/B(sub i))squared (r(squared)/r(sup 3)(sub i)) where m is ion mass, r is geometric distance; and E(sub i), B(sub i) and r(sub i) refer to the electric and magnetic field magnitudes and geocentric distance at the ionosphere, respectively. For purely longitudinal convection along a constant L shell the parallel force is F(cent. long) = F(sub cent, pole)(1 - (r/(r(sub i)L))(sup 3/2)/(1 - 3r/(4 r(sub i)L))(sup 5/2). For high latitudes the difference between these two cases is relatively unimportant below approximately 5 R(sub E). We find that the steady state O+ bulk velocities and parallel temperatures strongly increase and decrease, respectively, with convection strength. In particular, the bulk velocities increase from near 0 km/s at 4000 km altitude to approximately 10 km/s at 5 R(sub E) geocentric distance for 50-mV/m ionospheric convection electric field. However, the centrifugal effect on the steady O+ density profiles depends on the exobase ion and electron temperatures: for low-base temperatures (T(sub i) = T(sub e) = 3000 K) the O+ density at high altitudes increases greatly with convection, while for higher base temperatures (T(sub i) = 5000 K, T(sub e) = 9000 K), the high-altitude O+ density decreases somewhat as convection is enhanced. The centrifugal force further has a pronounced effect on the escaping O+ flux, especially for cool exobase conditions; as referenced to the 4000-km altitude, the steady state O+ flux increases from 10(exp 5) ions/sq cm/s when the ionospheric convection field E(sub i) = 0 mV/m to approximately 10(exp 7) ions/sq cm/s when E(sub i) = 100 mV/m. The centrifugal effect also decreases the time scale for approach to steady-state. For example, in the plasma expansion for T(sub i) = T(sub e) = 3000 K, the O+ density at 7 R(sub E) reaches only 10(exp -7) of it final value approximately 1.5 hours after expansion onset for E(sub i) = 0. For meridional convection driven by E(sub i) = 50 mV/m, the density at the same time after initial injection is 30-50% of its asymptotic level. The centrifugal acceleration described here is a possible explanation for the large (up to approximately 10 km/s or more) o+ outflow velocities observed in the midlatitude polar magnetosphere with the Dynamics Explorer 1 and Akebono spacecraft.

Horwitz, J. L.↗

Acquisition of Long-Duration, Low-Gravity Slosh Data Utilizing Existing ISS Equipment (SPHERES) for Calibration of CFD Models of Coupled Fluid-Vehicle Behavior

Accurate prediction of coupled fluid slosh and launch vehicle or spacecraft dynamics (e.g., nutation/precessional movement about various axes, attitude changes, ect.) requires Computational Fluid Dynamics (CFD) models calibrated with low-gravity, long duration slosh data. Recently completed investigations of reduced gravity slosh behavior have demonstrated the limitations of utilizing parabolic flights on specialized aircraft with respect to the specific objectives of the experiments. Although valuable data was collected, the benefits of longer duration low-gravity environments were clearly established. The proposed research provides the first data set from long duration tests in zero gravity that can be directly used to benchmark CFD models, including the interaction between the sloshing fluid and the tank/vehicle dynamics. To explore the coupling of liquid slosh with the motion of an unconstrained tank in microgravity, NASA's Kennedy Space Center, Launch Services Program has teamed up with the Florida Institute of Technology (FIT), Massachusetts Institute of Technology (MIT) and the Office of the Chief Technologist (OCT) to perform a series of slosh dynamics experiments on the International Space Station using the SPHERES platform. The Synchronized Position Hold Engage Reorient Experimental Satellites (SPHERES) testbed provides a unique, free-floating instrumented platform on ISS that can be utilized in a manner that would solve many of the limitations of the current knowledge related to propellant slosh dynamics on launch vehicle and spacecraft fuel tanks. The six degree of freedom (6-DOF) motion of the SPHERES free-flyer is controlled by an array of cold-flow C02 thrusters, supplied from a built-in liquid C02 tank. These SPHERES can independently navigate and re-orient themselves within the ISS. The intent of this project is to design an externally mounted tank to be driven inside the ISS by a set of two SPHERES devices (Figure 1 ). The tank geometry simulates a launch vehicle upper stage propellant tank and the maneuvers replicate those of real vehicles. The design includes inertial sensors, data acquisition, image capture and data storage interfaces to the SPHERES VERTIGO computer system on board the flight article assembly. The design also includes mechanical and electronic interfaces to the existing SPHERES hardware, which include self-contained packages that can operate in conjunction with the existing SPHERES electronics.

Schallhorn, Paul↗

Acquisition of Long-Duration, Low-Gravity Slosh Data Utilizing Existing ISS Equipment (SPHERES) for Calibration of CFD Models of Coupled Fluid-Vehicle Behavior

Accurate prediction of coupled fluid slosh and launch vehicle or spacecraft dynamics (e.g., nutation/precessional movement about various axes, attitude changes, ect.) requires Computational Fluid Dynamics (CFD) models calibrated with low-gravity, long duration slosh data. Recently completed investigations of reduced gravity slosh behavior have demonstrated the limitations of utilizing parabolic flights on specialized aircraft with respect to the specific objectives of the experiments. Although valuable data was collected, the benefits of longer duration low-gravity environments were clearly established. The proposed research provides the first data set from long duration tests in zero gravity that can be directly used to benchmark CFD models, including the interaction between the sloshing fluid and the tank/vehicle dynamics. To explore the coupling of liquid slosh with the motion of an unconstrained tank in microgravity, NASA's Kennedy Space Center, Launch Services Program has teamed up with the Florida Institute of Technology (FIT), Massachusetts Institute of Technology (MIT) and the NASA Game Changing Development Program (GCD) to perform a series of slosh dynamics experiments on the International Space Station using the SPHERES platform. The Synchronized Position Hold Engage Reorient Experimental Satellites (SPHERES) testbed provides a unique, free-floating instrumented platform on ISS that can be utilized in a manner that would solve many of the limitations of the current knowledge related to propellant slosh dynamics on launch vehicle and spacecraft fuel tanks. The six degree of freedom (6-DOF) motion of the SPHERES free-flyer is controlled by an array of cold-flow C02 thrusters, supplied from a built-in liquid C02 tank. These SPHERES can independently navigate and re-orient themselves within the ISS. The intent of this project is to design an externally mounted tank to be driven inside the ISS by a set of two SPHERES devices (Figure 1). The tank geometry simulates a launch vehicle upper stage propellant tank and the maneuvers replicate those of real vehicles. The design includes inertial sensors, data acquisition, image capture and data storage interfaces to the SPHERES VERTIGO computer system on board the flight article assembly. The design also includes mechanical and electronic interfaces to the existing SPHERES hardware, which include self-contained packages that can operate in conjunction with the existing SPHERES electronics

Schallhorn, Paul↗

A Synthesis of Fluid Dynamics and Quantum Chemistry for the Design of Nanoelectronics

In 1959, during a famous lecture entitled "There's Plenty of Room at the Bottom", Richard Feynman focused on the startling technical possibilities that would exist at the limit of miniaturization, that being atomically precise devices with dimensions in the nanometer range. A nanometer is both a convenient unit of length for medium to large sized molecules, and the root of the name of the new interdisciplinary field of "nanotechnology". Essentially, "nanoelectronics" denotes the goal of shrinking electronic devices, such as diodes and transistors, as well as integrated circuits of such devices that can perform logical operations, down to dimensions in the range of 100 nanometers. The thirty-year hiatus in the development of nanotechnology can figuratively be seen as a period of waiting for the bottom-up and atomically precise construction skills of synthetic chemistry to meet the top-down reductionist aspirations of device physics. The sub-nanometer domain of nineteenth-century classical chemistry has steadily grown, and state-of-the-art supramolecular chemistry can achieve atomic precision in non-repeating molecular assemblies of the size desired for nanotechnology. For nanoelectronics in particular, a basic understanding of the electron transport properties of molecules must also be developed. Quantum chemistry provides powerful computational methods that can accurately predict the properties of small to medium sized molecules on a desktop workstation, and those of large molecules if one has access to a supercomputer. Of the many properties of a molecule that quantum chemistry routinely predicts, the ability to carry a current is one that had not even been considered until recently. "Currently", there is a controversy over just how to define this key property. Reminiscent of the situation in high-Tc superconductivity, much of the difficulty arises from the different models that are used to simplify the complex electronic structure of real materials. A model-independent approach has been proposed, that sacrifices the plentiful molecular orbitals and Bloch functions of conventional approaches, for a single three-dimensional observable quantity, the electron momentum density Pi(sub rho). This quantity is simply the probability of any electron having momentum rho, multiplied by the total number of electrons in the system (the position of the electron is uncertain). We have explored the utility of this new approach in providing a fundamental understanding of the electron transport properties of molecules that have provi been nominated as candidates for components in the design of nanoelectronics; phenylene-ethynylene oligomers. Some of the molecular systems that have been studied are sketched below.

MacDougall, Preston J.↗

Development, construction and tests of the Mu2e electromagnetic calorimeter mechanical structures

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the charged lepton flavour violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of the existence of physics beyond the standard model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. In particular, the calorimeter provides excellent electron identification, a fast calorimetric online trigger, and complementary information to aid pattern recognition and track reconstruction. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1348 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements. The hostile Mu2e operational conditions, in terms of radiation levels (total expected ionizing dose of 12 krad and a neutron fluence of 5 × 10$^{10}$ n/cm$^{2}$ @ 1 MeV$_{eq}$ (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on scintillating materials, sensors, electronics and on the design of the detector mechanical structures and material choice. The support structure of each 674 crystal matrix is composed of an aluminum hollow ring and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and front-end electronics for the readout of each crystal are inserted in a machined copper holder and make a unique mechanical unit. The resulting 674 mechanical units are supported by a machined plate of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics are hosted in aluminum custom crates positioned on the external lateral surface of the disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system. In this paper we report on the calorimeter mechanical structure design, the mechanical and thermal simulations that have determined the design technological choices, and the status of component production, quality assurance tests and plans for assembly at Fermilab.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Simulation studies of vestibular macular afferent-discharge patterns using a new, quasi-3-D finite volume method

A quasi-three-dimensional finite-volume numerical simulator was developed to study passive voltage spread in vestibular macular afferents. The method, borrowed from computational fluid dynamics, discretizes events transpiring in small volumes over time. The afferent simulated had three calyces with processes. The number of processes and synapses, and direction and timing of synapse activation, were varied. Simultaneous synapse activation resulted in shortest latency, while directional activation (proximal to distal and distal to proximal) yielded most regular discharges. Color-coded visualizations showed that the simulator discretized events and demonstrated that discharge produced a distal spread of voltage from the spike initiator into the ending. The simulations indicate that directional input, morphology, and timing of synapse activation can affect discharge properties, as must also distal spread of voltage from the spike initiator. The finite volume method has generality and can be applied to more complex neurons to explore discrete synaptic effects in four dimensions.

NASA Center ARC↗

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spin transport of a doped Mott insulator in moiré heterostructures

Moiré superlattices of semiconducting transition metal dichalcogenide heterobilayers are model systems for investigating strongly correlated electronic phenomena. Specifically, WSe 2 /WS 2 moiré superlattices have emerged as a quantum simulator for the two-dimensional extended Hubbard model. Experimental studies of charge transport have revealed correlated Mott insulator and generalized Wigner crystal states, but spin transport of the moiré heterostructure has not yet been sufficiently explored. Here, we use spatially and temporally resolved circular dichroism spectroscopy to directly image the spin transport as a function of carrier doping and temperature in WSe 2 /WS 2 moiré heterostructures. We observe diffusive spin transport at all hole concentrations at 11 Kelvin — including the Mott insulator at one hole per moiré unit cell — where charge transport is strongly suppressed. At elevated temperatures the spin diffusion constant remains unchanged in the Mott insulator state, but it increases significantly at finite doping away from the Mott state. The doping- and temperature-dependent spin transport can be qualitatively understood using a t–J model, where spins can move via the hopping of spin-carrying charges and via the exchange interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Corrosion-Resistant Non-Carbon Electrocatalyst Supports for PEFCS

Insufficient durability of commercial catalysts is a persistent issue for PEFC applications.1 Herein, Pt on Nb-doped-TiO 2 is reported as a highly durable electrocatalyst on which the direct 4-electron reduction of oxygen to water is more facile compared to Pt/C. All performance metrics are reported comparing Pt/Nb-TiO2 with a 15% Pt loading, against a commercial Pt/C catalyst (46.5% Pt loading, Tanaka, K. K.). Nb-doped-TiO 2 with high surface area and high electronic conductivity was synthesized using the supercritical fluid method. Initially, the durability of the catalyst was characterized using accelerated stability tests (ASTs) involving 10,000 high potential cycles (DOE/FCCJ protocol) and the Pt/Nb-TiO 2 was found to retain 78% of its initial electrochemically active surface area (ECSA) compared to the 57.6 % retained by Pt/C. These observations were in excellent agreement with previous reports that the Pt particle size of Pt/C grew from 2nm to 8nm during the AST protocol along with severe corrosion and amorphization of the carbon surface.2, 3 TEM and XPS studies of the Pt-Nb-doped-TiO 2 catalyst showed that the Pt particle size grew from 3nm to 6nm and the Nb(IV) and Ti(III) in the support were oxidized to Nb(V) and Ti(IV) after the durability test. Thus, the improvement in Pt/Nb-TiO 2 ECSA retention was attributed to the lower extent of particle growth and lack of oxidative support loss upon oxidation as compared to Pt/C. The oxygen reduction reaction (ORR) performance was characterized by linear polarization using a rotating disk electrode (RDE). The electrochemical surface areas of Pt/Nb-TiO 2 and Pt/C were found to be 48m2·g -1 and 83m2·g -1 respectively, and the mass activity for the ORR at 0.9V vs. RHE were found to be 150 mA·mg-1Pt and 124 mA·mg-1Pt respectively. The improved mass activity on Pt/Nb-TiO 2 was attributed to strong metal support interaction (SMSI) between the Nb-TiO 2 support and the Pt catalyst based on the 625 meV decrease in the binding energy of the Pt4f x-ray photo-electron spectroscopy (XPS) peaks of Pt/Nb-TiO 2 compared to Pt/C. To quantify the impact of the SMSI, a kinetic model was applied to calculate the elementary reaction rate constants for the various steps of the ORR on both catalysts. The reaction rate constant (k1) for the direct 4-electron transfer pathway to produce H 2 O was significantly larger in Pt/Nb-TiO 2 as compared to Pt/C. Thus, the reduction in the electron binding energy as observed in the XPS was found to aid the facile filling of the higher energy 2p orbitals of O 2 thereby leading to improved 4-electron transfer kinetics and improved overall activity.

08 HYDROGEN↗

Quantitative Species Measurements in Microgravity Combustion Flames using Near-Infrared Diode Lasers

Understanding the physical phenomena controlling the ignition and spread of flames in microgravity has importance for space safety as well as for characterizing dynamical and chemical combustion processes which are normally masked by buoyancy and other gravity-related effects. Unfortunately, combustion is highly complicated by fluid mechanical and chemical kinetic processes, requiring the use of numerical modeling to compare with carefully designed experiments. More sophisticated diagnostic methods are needed to provide the kind of quantitative data necessary to characterize the properties of microgravity combustion as well as provide accurate feedback to improve the predictive capabilities of the models. Diode lasers are a natural choice for use under the severe conditions of low gravity experiments. Reliable, simple solid state operation at low power satisfies the operational restrictions imposed by drop towers, aircraft and space-based studies. Modulation wavelength absorption spectroscopy (WMS) provides a means to make highly sensitive and quantitative measurements of local gas concentration and, in certain cases, temperature. With near-infrared diode lasers, detection of virtually all major combustion species with extremely rapid response time is possible in an inexpensive package. Advancements in near-infrared diode laser fabrication technology and concurrent development of optical fibers for these lasers led to their use in drop towers. Since near-infrared absorption line strengths for overtone and combination vibrational transitions are weaker than the mid-infrared fundamental bands, WMS techniques are applied to increase detection sensitivity and allow measurement of the major combustion gases. In the first microgravity species measurement, Silver et al. mounted a fiber-coupled laser at the top of the NASA 2.2-sec drop tower and piped the light through a single-mode fiber to the drop rig. A fiber splitter divided the light into eight channels that directed the laser beam across a methane or propane diffusion jet flame. The light beams were recaptured by a set of gradient index lenses, coupled back into separate fiber optic lines, and transmitted back to detectors and electronics in the instrument package. In these experiments a 6-mm od fiber cable (containing the nine optical fibers) fell with the drop rig. Using separate detection and demodulation channels, spatial and temporal (up to 20 Hz) maps of water vapor and methane concentrations were obtained at differing heights in the flames. While this apparatus was useful from a demonstration standpoint, several drawbacks needed attention before useful scientific measurements could be obtained. First, eight lines of sight are somewhat insufficient for detailing the spatial profiles of the gas. Second, multiple detection channels operating in parallel are both expensive and present a challenge for accurate calibration. As a result, a newer scanning system was developed in our first contract under this program. The primary characteristic of this system is that it contains a single detection channel and achieves "continuous" spatial resolution by scanning the laser beam across the flame region, then directing this beam onto a single detector. Thus spatial measurements are converted to a temporal series of data. The true spatial resolution is limited only by the beam diameter and width of the sweep. In these experiments the beam is focused to about 1-mm diameter and scans across a region up to 4-cm wide.

Silver, Joel A.↗

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗