SEARCH · Engineering Papers
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.
Rayleigh is a 3-D convection code designed for the study of dynamo behavior in spherical geometry. It evolves the incompressible and anelastic MHD equations in spherical geometry using a pseudo-spectral approach.
We present a theory of two-component atomic Fermi superfluid with tunable pairing interaction in a thin spherical shell subject to spin-orbit coupling (SOC). By incorporating SOC into the Fermi superfluid in the BCS-Bose Einstein condensation (BEC) crossover, we obtain the energy spectrum and equations of state. While the order parameter and chemical potential are suppressed by SOC on the BCS side, the former remains positive but the latter may be pushed to negative values by SOC.Meanwhile, the compressibility exhibits kinks as the pairing interaction or SOC varies, indicating singular behavior of higher derivatives of the free energy despite the presence of the superfluid order parameter. The minimum of the energy dispersion indicates a decoupling of the energy gap from the order parameter, and the compressibility kinks occur when the energy gap approaches 0. We found the particle-hole mixing in the Fermi superfluid essential for the compressibility kinks since a Fermi gas with the same SOC but no pairing interaction only exhibits complicated dispersions but no singular behavior. Therefore, our results show that a combination of superfluid, SOC, and geometry can lead to interesting phenomena. Here, we also discuss possible experimental realizations and implications.
Flow imaging analysis (FIA), provides composition-agnostic morphological characterization. These measurements of particle size and shape are valuable to mass-based analysis, such as single particle inductively coupled plasma time-of-flight mass spectrometry (sp-ICP-TOFMS), which provides quantitative data on elements within particles. Using these two methods together enables informed use of geometric assumptions required by sp-ICP-TOFMS, as particle mass is typically converted to a particle diameter using assumed-spherical geometry. To validate this concept, parallel measurements to determine particle diameters were performed by FIA and sp-ICP-TOFMS on four particle suspensions: 300 nm polystyrene Eu-doped nanoparticles, 1 μm Fe-rich beads, 3 μm four element calibration polystyrene beads and 5 μm polystyrene beads. The Fe-particles obtained the highest percent difference from the manufacturer’s nominal diameter, as the mean diameter obtained by FIA was overestimated by 21% and sp-ICP-TOFMS underestimated the mean diameter by 20.7%. Two types of particles were selected to test the effect of varying the particle number concentrations (PNC) on sizing accuracy, and both methods accurately sized each particle population at the PNC expected. Single particle analysis of carbon has continued to be a popular research topic, with direct applications to environmental pollutants in terms of nano- and micro- plastics. Real-world plastic particles were studied, and FIA’s measured circularity values demonstrated that the particles deviated from spherical geometries, therefore sp-ICP-TOFMS data should be interpreted as mass-based rather than size-based. Combining these techniques enables improved interpretation of particle populations and evaluation of particle sizes.
We have extended the radiation diffusion model of Hammer and Rosen [Phys. Plasmas 10, 1829 (2003)] to diverging spherical and cylindrical geometries. The effect of curvilinear geometry on the supersonic, expanding wavefront increases as the internal radius of a spherical or cylindrical shell approaches zero. Small spherical geometries are important for modeling systems at the size scale of inertial confinement fusion capsules, at these scales existing quasi-analytic models for planar geometry significantly disagree with the results of simulation. With this method, the benefits of rapid iteration can be applied to common spherical systems at much smaller length scales. We present comparisons between numerical diffusion solutions and the analytic model to give ranges of applicability for the model.
The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.
Neutron single crystal diffractometers require large solid angle coverage for optimum performance. This can be achieved by tiling flat detectors in a cylindrical or spherical geometry, but results in large gaps in detector coverage and parallax distortion. The detector edges exhibit degraded resolution, distortion, and gamma rejection. Since the detector edge regions are a significant fraction of the detector active area, they must be removed from the experimental data set, requiring extra beam time to collect enough analysis data. A spherical detector with a continuous surface would effectively address this issue while eliminating most boundary ‘dead’ areas. Here we report on the development of a novel hemispherical shaped neutron detector using seamlessly tiled readout modules to form the desired shape. The heart of the detector is a specially developed curved neutron scintillator coupled to high resolution silicon photomultiplier (SiPM) Anger cameras via custom made fiber optic tapers (FOTs). The detector has been assembled and initial tests have been conducted at the High Flux Isotope Reactor (HFIR) beamlines at Oak Ridge National Laboratory (ORNL). Here, in this work, we describe details of the scintillator design, fabrication and characterization, evaluation of individual detector modules, the details of the detector design implementation, and evaluation of the assembled detector at ORNL beamlines.
A collisional particle code based on implicit energy- and charge-conserving methods in axisymmetric geometries is presented. A new particle pusher for axisymmetric systems is introduced that is compatible with exact energy and charge conservation and yields improved accuracy compared to other methods. How to appropriately treat all aspects of the algorithm near the r = 0 axis of symmetry is described in detail. Here, the axisymmetric model is verified by simulating the free expansion of a plasma sphere in 2D cylindrical and 1D spherical geometries. The algorithm's ability to study the dynamic compression of a dense plasma is illustrated by simulating the dynamic Z-pinch in 1D cylindrical geometry.
The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).
In glow discharge-type deuterium–deuterium fusion neutron sources, fusion reactions occurring on the cathode surface contribute significantly to the neutron production rate (NPR); therefore, cathode design plays a crucial role in enhancing NPR. Although the NPR generally increases with discharge current, deuterium desorption caused by cathode heating tends to stagnate it in the high-current regime. Water cooling has been shown to mitigate this stagnation for nontransparent cathodes, but its effectiveness for transparent cathodes has not been experimentally clarified. Here, in this study, a transparent cathode with active water cooling and a spherical geometry was designed and fabricated using stainless steel tubing. Neutron production experiments were conducted with and without water cooling: with water cooling, applied voltages of 20–40 kV were investigated, while without water cooling, the applied voltage was limited to 25 kV. The results showed that the NPR increased continuously with discharge current without stagnation when water cooling was applied. A maximum NPR of (1.8 ± 0.02) × 10 6 n/s was achieved at 40 kV and 50 mA. Compared with our previously reported disk cathode and nontransparent water-cooled cathode of the same diameter, the transparent water-cooled cathode exhibited a significantly enhanced NPR. These findings indicate that combining cathode transparency with active cooling is an effective strategy for improving neutron production in glow discharge-type fusion neutron sources.
This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.
Laser–target interactions generate intense electromagnetic pulses (EMPs) that can interfere with measurements and damage equipment. In this paper we show that applying a magnetic field to nanosecond pulse laser–target interactions decreases the magnitude of EMPs. We demonstrate this effect in two experiments with different geometries (spherical versus planar), laser intensities ( ∼10 13 versus ∼10 15 W/cm 2 ) and applied field strength (12 versus 0.1 T) that both observed suppression of EMPs in approximately the 1 GHz band (by factors of 0.65× and 0.32×, respectively). We then observe the opposite effect at high intensities with a picosecond pulse: for planar experiments with laser intensities of approximately 10 19 W/cm 2 and magnetic fields of 6–10 T, the magnitude of EMPs is increased by a factor of 1.75×. These results provide a benchmark for models of EMP generation, but suggest that magnetic fields are not a viable solution for mitigating EMPs in the high-intensity laser facilities where they are most damaging.
Ablation pressure is a key parameter controlling laser drive performance in direct-drive inertial confinement fusion implosions. In this work, we develop a theoretical framework to clarify the role of thermal conduction in determining how ablation pressure scales with laser intensity, wavelength, and target size. We focus on a steady-state, spherically expanding plasma produced when an over-dense pellet is irradiated by a high-intensity laser. Scaling laws for ablation pressure and the mass ablation rate are derived with respect to laser intensity, wavelength, and target size. As a result, the theoretical predictions are found to be in good agreement with one-dimensional hydrodynamics simulations in spherical geometry.
Data-driven approaches using deep learning are emerging as powerful techniques to extract non-Gaussian information from cosmological large-scale structure. Here, this work presents the first simulation-based inference (SBI) pipeline that combines weak lensing and galaxy clustering maps in a realistic Dark Energy Survey Year 3 (DES Y3) configuration and serves as preparation for a forthcoming analysis of the survey data. We develop a scalable forward model based on the CosmoGridV1 suite of N-body simulations to generate over one million self-consistent mock realizations of DES Y3 at the map level. Leveraging this large dataset, we train deep graph convolutional neural networks on the full survey footprint in spherical geometry to learn low-dimensional features that approximately maximize mutual information with target parameters. These learned compressions enable neural density estimation of the implicit likelihood via normalizing flows in a ten-dimensional parameter space spanning cosmological 𝑤CDM, intrinsic alignment, and linear galaxy bias parameters, while marginalizing over baryonic, photometric redshift, and shear bias nuisances. To ensure robustness, we extensively validate our inference pipeline using synthetic observations derived from both systematic contaminations in our forward model and independent Buzzard galaxy catalogs. Our forecasts yield significant improvements in cosmological parameter constraints, achieving 2−3× higher figures of merit in the 𝛺 𝑚 − 𝑆 8 plane relative to our implementation of baseline two-point statistics and effectively breaking parameter degeneracies through probe combination. These results demonstrate the potential of SBI analyses powered by deep learning for upcoming Stage-IV wide-field imaging surveys.
The Zeeman splittings and electron paramagnetic resonance frequencies of alkali-metal atoms are shifted in the presence of a polarized noble gas. For a spherical geometry, the shift is enhanced over what is expected classically by a dimensionless atomic parameter 𝜅 0 that is unique to each alkali-metal-atom–noble-gas pair. Here, we present a precise measurement of 𝜅 0 for the 39 K − 3 He system with a relative accuracy of better than 1%. A critical component of achieving subpercent accuracy involved characterizing the shape of our samples using both MRI and CT medical-imaging techniques. The parameter 𝜅 0 plays an important role in establishing the absolute polarization of 3 He in a variety of contexts, including polarized targets for electron-scattering experiments and MRI of the gas space of the lungs. Our measurement more than doubles the accuracy possible when using 𝜅 0 for polarimetry purposes. Just as important, the work presented here represents a direct measurement of 𝜅 0 for the 39 K − 3 He system; previous values for 𝜅 0 in the 39 K − 3 He system relied on a chain of measurements that were benchmarked by previous measurements of 𝜅 0 in the Rb- 3 He system.
We demonstrate that formulating the composite-fermion theory of the fractional quantum Hall (FQH) effect in terms of quaternions greatly expands its reach and opens the door into many interesting issues that were previously not amenable to quantitative theoretical investigation. As an illustration, we explore the possibility of a nematic or a charge-density wave instability of the composite-fermion Fermi sea at half-filled Landau level and of the nearby FQH states by looking for a gap closing instability of the neutral magneto-roton excitation. As a result, our quaternion formulation of the FQH effect has been inspired by mathematical developments in the theoretical analyses of gravitational wave modes and cosmic microwave background radiation, where an important role is played by spin-weighted spherical harmonics that are nothing but monopole harmonics appearing in the spherical geometry for the FQH effect.
In this paper I attempt to document what I have learned and still remember about temperature effects, particularly Doppler broadening, on nuclear data, as it applies to ENDF formatted data. My focus is on the SIGMA1 method of Doppler broadening that I developed over 50 years ago. However, I do realize that there are currently many different methods used in computer codes, and I have tried to keep my discussion general as it applies to ALL of these methods. In particular, ALL of the MYTHS I describe below apply to all Doppler broadening methods: they all solve the parabolic diffusion equation in spherical geometry in (speed, reaction rate) versus temperature. These methods differ only in how the cross section is represented. So, at 85 years old here is what I still remember.
This project presents the development and evaluation of a general drag model for gas–solid multiphase flows using deep learning techniques. A comprehensive database of more than 4,000 experimental and numerical data points for spherical and non spherical particles was compiled, incorporating geometric features such as sphericity, aspect ratio, and orientation. Several predictive approaches—including traditional em pirical correlations, machine learning, and deep neural networks—were benchmarked, with the proposed Drag Coefficient Correlation-aided Deep Neural Network (DCC DNN) demonstrating superior accuracy. To account for particle–particle interactions, additional drag data were generated using CFD-based simulations of packed and flu idized beds, leading to the development of a retrained model capable of incorporat ing volume fraction effects. Integration of the trained model with the MFiX CFD solver was achieved using FTorch, enabling drag predictions during discrete element method (DEM) simulations. Validation against experimental data for single particles and fluidized beds confirmed the model’s improved predictive ability, particularly for non-spherical geometries. While the model performed strongly under fluidized con ditions, limitations remained in unfluidized regimes, suggesting a need for expanded datasets. Overall, this study demonstrates the feasibility of combining deep learning with physics-informed CFD to improve drag modeling for gas–solid flows, with promis ing implications for scaling multiphase simulations in industrial applications.
is an open-source program to select asymmetric geometries of spherically bent crystal analyzers for X-ray spectroscopy, thus achieving improved energy resolution and operational flexibility.