Search NASA⌕ Search

SEARCH · Search NASA

Results for “SPACE SIMULATION CHAMBER”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗
DOE OSTI2026

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING↗
DOE OSTI2020

Measurement of space charge effects in the MicroBooNE LArTPC using cosmic muons

Large liquid argon time projection chambers (LArTPCs), especially those operating near the surface, are susceptible to space charge effects. In the context of LArTPCs, the space charge effect is the build-up of slow-moving positive ions in the detector primarily due to ionization from cosmic rays, leading to a distortion of the electric field within the detector. This effect leads to a displacement in the reconstructed position of signal ionization electrons in LArTPC detectors ("spatial distortions"), as well as to variations in the amount of electron-ion recombination experienced by ionization throughout the volume of the TPC. We present techniques that can be used to measure and correct for space charge effects in large LArTPCs by making use of cosmic muons, including the use of track pairs to unambiguously pin down spatial distortions in three dimensions. Finally, the performance of these calibration techniques are studied using both Monte Carlo simulation and MicroBooNE data, utilizing a UV laser system as a means to estimate the systematic bias associated with the calibration methodology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Wire-cell 3D pattern recognition techniques for neutrino event reconstruction in large LArTPCs: algorithm description and quantitative evaluation with MicroBooNE simulation

Wire-Cell is a 3D event reconstruction package for liquid argon time projection chambers. Through geometry, time, and drifted charge from multiple readout wire planes, 3D space points with associated charge are reconstructed prior to the pattern recognition stage. Pattern recognition techniques, including track trajectory and dQ/dx (ionization charge per unit length) fitting, 3D neutrino vertex fitting, track and shower separation, particle-level clustering, and particle identification are then applied on these 3D space points as well as the original 2D projection measurements. A deep neural network is developed to enhance the reconstruction of the neutrino interaction vertex. Compared to traditional algorithms, the deep neural network boosts the vertex efficiency by a relative 30% for charged-current ν e interactions. Therefore, this pattern recognition achieves 80–90% reconstruction efficiencies for primary leptons, after a 65.8% (72.9%) vertex efficiency for charged-current ν e (ν μ ) interactions. Based on the resulting reconstructed particles and their kinematics, we also achieve 15-20% energy reconstruction resolutions for charged-current neutrino interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Spatio-temporal focal spot characterization and modeling of the NIF ARC kilojoule picosecond laser

The advanced radiographic capability (ARC) laser system, part of theNational Ignition Facility (NIF) at Lawrence Livermore National Laboratory, is a short-pulse laser capability integrated into the NIF. The ARC is designed to provide adjustable pulse lengths of ~1–38 ps in four independent beamlets, each with energies up to 1 kJ (depending on pulse duration). A detailed model of the ARC lasers has been developed that predicts the time- and space-resolved focal spots on target for each shot. Measurements made to characterize static and dynamic wavefront characteristics of the ARC are important inputs to the code. Modeling has been validated with measurements of the time-integrated focal spot at the target chamber center (TCC) at low power, and the space-integrated pulse duration at high power, using currently available diagnostics. In this work, these simulations indicate that each of the four ARC beamlets achieves a peak intensity on target of up to a few 10 18 W=cm 2 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Calibration and simulation of ionization signal and electronics noise in the ICARUS liquid argon time projection chamber

The ICARUS liquid argon time projection chamber (LArTPC) neutrino detector has been taking physics data since 2022 as part of the Short-Baseline Neutrino (SBN) Program. This paper details the equalization of the response to charge in the ICARUS time projection chamber (TPC), as well as data-driven tuning of the simulation of ionization charge signals and electronics noise. The equalization procedure removes non-uniformities in the ICARUS TPC response to charge in space and time. This work leverages the copious number of cosmic ray muons available to ICARUS at the surface. The ionization signal shape simulation applies a novel procedure that tunes the simulation to match what is measured in data. The end result of the equalization procedure and simulation tuning allows for a comparison of charge measurements in ICARUS between Monte Carlo simulation and data, showing good performance with minimal residual bias between the two.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

The Dual Nature of Entrainment-Mixing Signatures Revealed through Large-Eddy Simulations of a Convection-Cloud Chamber

Abstract Entrainment of subsaturated air into a cloud can influence its optical and microphysical properties in various ways, depending on the droplet evaporation and turbulent mixing time scales. Previous experiments in the Pi convection-cloud chamber have revealed that, given a fixed entrained air property, the mixing of entrained subsaturated air results in complete evaporation of some cloud droplets, with the rest remaining unchanged. This is a signature of inhomogeneous mixing. While comparing the results of entrainment with varying air properties, the mixing signature appears as if the subsaturated air is well mixed with the cloud to evenly reduce the droplets’ size. In other words, taken together, the experiments appear to have the signature of homogeneous mixing. To explore these results in a greater depth, we conduct large-eddy simulations combined with a bin microphysics scheme. Our results reproduce the similar signatures of inhomogeneous and homogeneous mixing, implying that LES can resolve the inhomogeneous mixing when the grid spacing is smaller than the entrained air parcel. Additionally, we observe that increasing the aerosol injection rate enhances the signature of inhomogeneous mixing, while coarser grid spacing diminishes it. Finally, the change in wall fluxes in response to various entrained air properties confirms that the homogeneous signature seen in the analysis of an ensemble of simulations is the result of various equilibrium states. This further strengthens the suggestion that the homogeneous mixing signature found in aircraft observations near the cloud top may result from combining entrainment events of different intensities, possibly caused by various-sized eddies. Significance Statement Large-eddy simulation and size-resolved microphysics can resolve time scales for turbulent mixing and evaporation and, therefore, are well suited for reproducing, extending, and interpreting the entrainment experiment in the Pi convection-cloud chamber. Our simulation results confirm (i) the inhomogeneous mixing signature for an individual entrainment event and (ii) the appearance of homogeneous mixing in an ensemble of entrainment episodes. Furthermore, we demonstrate that the inhomogeneous mixing signature is more pronounced in a polluted cloud, but coarser grid spacing in simulations may compromise the accuracy of this signature. Last, the homogeneous mixing signature results from various equilibrium states established for different entrainment intensities and adjusted wall fluxes, which are challenging to measure experimentally but can be easily analyzed in the simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

SBND Calibration

The Short-Baseline Near Detector (SBND) is a 112-ton liquid argon time projection chamber (LArTPC) serving as the near detector of the Short-Baseline Neutrino Program at Fermilab. SBND aims to make precise measurements of neutrino-argon interactions, necessitating an extensive calibration program for the LArTPC to ensure precision/accuracy in these measurements. This talk describes the use of cosmogenic, radiogenic, and neutrino interaction activity with both Monte Carlo simulation and data in calibrating the SBND subsystems of the Time Projection Chamber, Photon Detection System, and Cosmic Ray Tagger. We will highlight some of SBND's first calibration results, including electronics calibrations, measurement of the electron lifetime, shower calibration, and studies of space charge effects in the detector.

Oh, Sungbin [Fermilab]↗
DOE OSTI2020

Laser Transport System Vacuum Simulations and LED Atom Tracker

MAGIS-100 is an experiment using atom interferometry on strontium atoms in a 100m vacuum chamber. The laser travels through a laser transport system (LTS) before going down the chamber. Two projects were completed over the course of the internship. One project was conducting vacuum simulations on the LTS to ensure it meets the experimental pressure requirement of 10-11 torr. This was done by characterizing the pressure profile with variations in pump size, pump spacing, and orifice size. The initial design of the LTS included one 30 L/s ion pump in the high vacuum (HV) region, a lens mount with holes of diameter 0.788”, and three 150 L/s ion pumps in the ultra-high vacuum (UHV) region. After adjusting these variables, it was concluded that the lens mount holes could be reduced to 0.200” in diameter and one 150 L/s ion pump in the UHV could be removed and still meet the experimental requirements. The other project was to design an LED atom tracker that shows the position of the strontium atoms in the 100m shaft. This will help with public outreach by allowing visitors to have a visual of the experiment. An Arduino Uno and Arduino Due were programmed using classical physics formulas to direct the RGB LED strips to flash in sync with the falling atoms. The requirements for powering and connecting the hardware were calculated and recommendations for scaling the system to 100m were made.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

SAF Combustion & Soot Processes [Slides]

The atomization, mixing, combustion and emissions characteristics of aviation fuels were measured using a novel approach based on a non-continuous injection scheme called the single-hole atomizer (SHA). High-speed microscopy revealed differences between fuels in terms of evaporation and mixing regimes over conditions relevant to modern and next generation aero-engine combustors. Measurements of liquid and vapor penetration, mixing fields, combustion and emissions metrics (ignition delay, lift-off length, PAH formation, soot mass) highlighted the effects of fuels and combustor conditions. The experimental results are being leveraged to adjust and validate chemical and CFD models. Detailed analysis of sampled soot showed subtle differences in soot morphology between fuels. The results revealed the presence of contaminants potentially affecting surface chemistry and the nucleation propensity of water droplets on particles. Chemical mechanisms for NJFCP A-2, C-1 and C-4 showed good performance over a large parameter space. Spray breakup at relight conditions is vastly different from the atomization observed at high pressure. CFD simulations of the SHA target conditions confirmed the good behavior of the C-1 kinetic mechanism. The simulations support the strong relationship between low and high temperature reactions. New altitude chamber facility to enable detailed characterization of the heterogeneous nucleation process of water on aerosol particles.

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

Measurement of three-dimensional inclusive muon-neutrino charged-current cross sections on argon with the MicroBooNE detector

We report the measurement of the triple-differential cross section d 3 σ/dE vis d cos(θ μ )dP μ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4 x 10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. The mapping from reconstructed kinematics to truth quantities is validated within uncertainties by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after applying a conditional constraint from the muon momentum distribution in data. The success of this validation provides confidence that the energy transfer in the MicroBooNE detector is well-modeled within simulation uncertainties, enabling a reliable unfolding to a triple-differential cross section defined at the nominal neutrino flux over muon momentum, muon scattering angle, and visible neutrino energy. This validation not only supports accurate cross-section extraction, but also establishes a critical foundation for tuning interaction models used in future neutrino oscillation measurements. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Supersaturation Variability from Scalar Mixing: Evaluation of a New Subgrid-Scale Model Using Direct Numerical Simulations of Turbulent Rayleigh–Bénard Convection

We report that supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars - water vapor and temperature - leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Electron Cloud Simulations for the Electron-Ion Collider in Brookhaven National Laboratory

For high-intensity circular accelerators and storage rings, if the secondary electron emission yield (SEY) of the vacuum chamber surfaces is high, an EC (Electron Cloud) could build up with the passage of the circulating beam. The presence of EC can strongly affect the beam quality, such as transverse instabilities, transverse emittance growth, and beam loss. Furthermore, the heat load from EC can exceed the available cryogenic capability. For the RHIC superconducting (SC) arc magnets, to be used for the hadron storage ring of the Electron-Ion Collider (EIC), the dynamic heat load budget is 0.5 W/m to the 4.5 K stainless steel beam pipe. This can limit the maximum beam bunches or intensity of the EIC, hence diminishing the luminosity provided by the EIC. The EC has affected the beam instability and significantly contributed to cryogenic heat load in the Large Hadron Collider (LHC). Positron storage rings for which ECs have been a critical factor in the design and performance include KEKB in Japan and EC buildup remains one of the concerns for future high-intensity accelerators design. EC considerations have driven the SuperKEKB collider design and the positron damping ring for the proposed International Linear Collider (ILC). The LHC luminosity upgrade is contingent on reducing the bunch spacing to 25 ns; at this bunch spacing, severe EC buildup has been observed. The success of the upgrade is likely contingent on limiting EC buildup. To study the EC heat load, we did some EC simulations with PyECLOUD code for the dipole, quadrupole, sextupole magnets, and the warm (drift) section of the EIC hadron storage ring. PyECLOUD is an EC simulation code developed by CERN. The code has been validated and used to study EC in the LHC, SPS, and PS. To eliminate EC buildup, the sources of electrons must be minimized. First, we will reduce the production of electrons due to residual gas ionization by specifying the maximum gas density, which requires a vacuum chamber with low electron-stimulated desorption (ESD) yields and a sufficient (preferably distributed) pumping speed. Second, we should reduce the secondary electrons with a lower secondary electron yield (SEY) material.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Time-Dependent Boundary Modeling to Inform Design of SPARC Diagnostic and Actuators

The purpose of this work was to inform the design of diagnostic and actuator systems for a new experimental facility for Commonwealth Fusion Systems (CFS), called SPARC [Creely2020], using time-dependent plasma boundary simulations with the Scrape Off Layer Plasma Simulator (SOLPS) code [Wiesen2015]. The SPARC tokamak is a device that aims to demonstrate fusion energy production and high-field plasma operating scenarios. The developed technologies will be incorporated into a net-electricity pilot that will demonstrate approximately 200MW of electric power. A major difference from other private industry fusion companies is the use of a standard aspect ratio, but at a high magnetic field strength, which provides a well-established physics basis and a modest extrapolation to reactor operation. Time-dependent simulations were used to inform the design of main-chamber plasma facing components, which are protected by feedback control of the heat and particle loads during high power operation. The simulations provide the timescales and magnitudes of plasma and neutral particle response useful for determining the position and operation of diagnostics and actuators that will enable control. The plasma response and phase-space diagrams also provide input into advanced model-based control schemes, which are envisioned for later SPARC operation. Model-based control can reduce the risk of component failure, which would require expensive in-situ repairs and the associated delays, by predicting the system response to actuators over a short time-horizon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2021

Quasi-steady testing approach for high-power Hall thrusters

Hall effect thrusters operating at power levels in excess of several hundreds of kilowatts have been identified as enabling technologies for applications such as lunar tugs, large satellite orbital transfer vehicles, and solar system exploration. These large thrusters introduce significant testing challenges due to the propellant flow rate exceeding the pumping speed available in most laboratories. Even with proposed upgrades in mind, the likelihood that multiple vacuum facilities will exist in the near future to allow long duration testing of high-power Hall thrusters operating at power levels in excess of 100 kW remains extremely low. In this article, we numerically explore the feasibility of testing Hall thrusters in a quasi-steady mode defined by pulsing the mass flow rate between a nominal and a low value. Here, our simulations indicate that sub-second durations available before the chamber reaches critical pressure are sufficiently long to achieve the steady-state current and flow field distributions, allowing us to characterize thruster performance and the near plume region.

42 ENGINEERING↗
DOE OSTI2022

The Dynamics of Ions on Phased Radio-frequency Carpets in High Pressure Gases and Application for Barium Tagging in Xenon Gas Time Projection Chambers

Radio-frequency (RF) carpets with ultra-fine pitches are examined for ion transport in gases at atmospheric pressures and above. We develop new analytic and computational methods for modeling RF ion transport at densities where dynamics are strongly influenced by buffer gas collisions. An analytic description of levitating and sweeping forces from phased arrays is obtained, then thermodynamic and kinetic principles are used to calculate ion loss rates in the presence of collisions. This methodology is validated against detailed microscopic SIMION simulations. We then explore a parameter space of special interest for neutrinoless double beta decay experiments: transport of barium ions in xenon at pressures from 1 to 10 bar. Our computations account for molecular ion formation and pressure dependent mobility as well as finite temperature effects. We discuss the challenges associated with achieving suitable operating conditions, which lie beyond the capabilities of existing devices, using presently available or near-future manufacturing techniques.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Simulations and analysis tools for charge-exchange (d, 2 He) reactions in inverse kinematics with the AT-TPC

Charge-exchange (d, 2 He) reactions in inverse kinematics at intermediate energies are a very promising method to investigate the Gamow–Teller transition strength in unstable nuclei. A simulation and analysis software based on the attpcroot package was developed to study this type of reactions with the active-target time projection chamber (AT-TPC). The simulation routines provide a realistic detector response that can be used to understand and benchmark experimental data. Analysis tools and correction routines can be developed and tested from simulations in ATTPCROOT , because they are processed in the same way as the real data. In particular, we study the feasibility of using coincidences with beam-like particles to unambiguously identify the (d, 2 He) reaction channel, and to develop a kinematic fitting routine for future applications. More technically, the impact of space-charge effects in the track reconstruction, and a possible correction method are investigated in detail. Finally, this analysis and simulation package constitutes an essential part of the software development for the fast-beams program with the AT-TPC.

(d,2He)↗