Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spherical geometry”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Anisotropic line emission and the geometry of the broad-line region in active galactic nuclei

Attention is given to an escape probability-based formalism, including the interaction with a background opacity, which is used to predict line-beaming factors for broad-line region (BLR) clouds of AGN. Several computational details necessary to model gas exposed to an intense nonthermal continuum are outlined. It is found that all hydrogen lines are strongly inwardly beamed, and that C IV 1549 is moderately beamed. The implications of this for the geometry of the emission-line regions of AGN are discussed. Beaming is found to cause the source-cloud separation to be overestimated by variability studies (so that the radiation field at the BLR is even more intense than previously thought); it can even cause spherical geometries to have little instantaneous response to continuum changes.

Ferland, G. J.↗
NASA NTRS1978

Design and performance of a modularized NaI/Tl/ detector /The Crystal Ball prototype/

A prototype NaI(Tl) detector (the Cluster of 54) of spherical geometry subtending a solid angle of 7.5% of 4pi at its center, has recently been assembled and tested. This detector consists of 54 close-packed but optically isolated NaI(Tl) modules and the associated electronic circuitry. The Cluster of 54 is the predecessor of an almost complete spherical detector, the Crystal Ball, which will cover 94% of 4pi. The latter detector is now under construction and is especially designed for the study of gamma-rays produced in electron-positron collisions at colliding beam facilities. This article outlines the mechanical, optical, and electronic assembly of the prototype system.

Chan, Y.↗
DOE OSTI2025

Portable Distortion Free Large Solid Angle Coverage Detector

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.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

An implicit particle code with exact energy and charge conservation for studies of dense plasmas in axisymmetric geometries

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

The linear Boltzmann equation in slab geometry - Development and verification of a reliable and efficient solution

The linear Boltzmann equation can be cast in a form mathematically identical to the radiation-transport equation. A multigroup procedure is used to reduce the energy (or velocity) dependence of the transport equation to a series of one-speed problems. Each of these one-speed problems is equivalent to the monochromatic radiative-transfer problem, and existing software is used to solve this problem in slab geometry. The numerical code conserves particles in elastic collisions. Generic examples are provided to illustrate the applicability of this approach. Although this formalism can, in principle, be applied to a variety of test particle or linearized gas dynamics problems, it is particularly well-suited to study the thermalization of suprathermal particles interacting with a background medium when the thermal motion of the background cannot be ignored. Extensions of the formalism to include external forces and spherical geometry are also feasible.

Stamnes, K.↗
NASA NTRS2013

Vertical Distribution of Dust and Water Ice Aerosols from CRISM Limb-geometry Observations

[1] Near-infrared spectra taken in a limb-viewing geometry by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on board the Mars Reconnaissance Orbiter provide a useful tool for probing atmospheric structure. Specifically, the observed radiance as a function of wavelength and height above the limb enables the vertical distribution of both dust and water ice aerosols to be retrieved. More than a dozen sets of CRISM limb observations have been taken so far providing pole-to-pole cross sections, spanning more than a full Martian year. Radiative transfer modeling is used to model the observations taking into account multiple scattering from aerosols and the spherical geometry of the limb observations. Both dust and water ice vertical profiles often show a significant vertical structure for nearly all seasons and latitudes that is not consistent with the well-mixed or Conrath-v assumptions that have often been used in the past for describing aerosol vertical profiles for retrieval and modeling purposes. Significant variations are seen in the retrieved vertical profiles of dust and water ice aerosol as a function of season. Dust typically extends to higher altitudes (approx. 40-50km) during the perihelion season than during the aphelion season (<20km), and the Hellas region consistently shows more dust mixed to higher altitudes than other locations. Detached water ice clouds are common, and water ice aerosols are observed to cap the dust layer in all seasons.

Atmosphere↗
NASA NTRS2002

Defining Top-of-Atmosphere Flux Reference Level for Earth Radiation Budget Studies

To estimate the earth's radiation budget at the top of the atmosphere (TOA) from satellite-measured radiances, it is necessary to account for the finite geometry of the earth and recognize that the earth is a solid body surrounded by a translucent atmosphere of finite thickness that attenuates solar radiation differently at different heights. As a result, in order to account for all of the reflected solar and emitted thermal radiation from the planet by direct integration of satellite-measured radiances, the measurement viewing geometry must be defined at a reference level well above the earth s surface (e.g., 100 km). This ensures that all radiation contributions, including radiation escaping the planet along slant paths above the earth s tangent point, are accounted for. By using a field-of- view (FOV) reference level that is too low (such as the surface reference level), TOA fluxes for most scene types are systematically underestimated by 1-2 W/sq m. In addition, since TOA flux represents a flow of radiant energy per unit area, and varies with distance from the earth according to the inverse-square law, a reference level is also needed to define satellite-based TOA fluxes. From theoretical radiative transfer calculations using a model that accounts for spherical geometry, the optimal reference level for defining TOA fluxes in radiation budget studies for the earth is estimated to be approximately 20 km. At this reference level, there is no need to explicitly account for horizontal transmission of solar radiation through the atmosphere in the earth radiation budget calculation. In this context, therefore, the 20-km reference level corresponds to the effective radiative top of atmosphere for the planet. Although the optimal flux reference level depends slightly on scene type due to differences in effective transmission of solar radiation with cloud height, the difference in flux caused by neglecting the scene-type dependence is less than 0.1%. If an inappropriate TOA flux reference level is used to define satellite TOA fluxes, and horizontal transmission of solar radiation through the planet is not accounted for in the radiation budget equation, systematic errors in net flux of up to 8 W/sq m can result. Since climate models generally use a plane-parallel model approximation to estimate TOA fluxes and the earth radiation budget, they implicitly assume zero horizontal transmission of solar radiation in the radiation budget equation, and do not need to specify a flux reference level. By defining satellite-based TOA flux estimates at a 20-km flux reference level, comparisons with plane-parallel climate model calculations are simplified since there is no need to explicitly correct plane-parallel climate model fluxes for horizontal transmission of solar radiation through a finite earth.

Loeb, N. G.↗
NASA NTRS1997

Spherical Diffusion Flames: Structure and Dynamics

The spherical geometry is the most suitable one to study fundamental issues concerning the structure and the dynamics of diffusion flames. From a theoretical point of view, it is the only geometry that permits the existence of a truly one dimensional stationary diffusion flame. A stationary planar diffusion flame with the fuel supplied upstream, at x approaches -infinity say, and the oxidant downstream, at x approaches +infinity is not possible. For a steady diffusion flame to exist, one must have nonzero fluxes of fuel and oxidant towards the flame. However, in the unlimited region behind the planar flame the only bounded solutions to the reaction-free convective-diffusive operator are constants. Hence the oxidant concentration behind the flame remains constant and there is no mechanism to generate the necessary flux towards the flame. A one-dimensional problem can be formulated if the reactants are supplied at finite locations; but the boundary conditions in this case introduce unnecessary complications and do not appropriately model the physical reality. Indeed, a planar diffusion flame can be established in the stagnation-point flow of two opposed jets but the flame in this case is stretched and the flow is essentially two-dimensional. The only stationary one-dimensional diffusion flame in an unlimited environment is therefore the spherical flame.

Matalon, Moshe↗
NASA NTRS2011

Vertical Distribution of Aersols and Water Vapor Using CRISM Limb Observations

Near-infrared spectra taken in a limb-viewing geometry by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter (MRO) provide a useful tool for probing atmospheric structure. Specifically, the observed radiance as a function of wavelength and height above the limb allows the vertical distribution of both dust and ice aerosols to be retrieved. These data serve as an important supplement to the aerosol profiling provided by the MRO/MCS instrument allowing independent validation and giving additional information on particle physical and scattering properties through multi-wavelength studies. A total of at least ten CRISM limb observations have been taken so far covering a full Martian year. Each set of limb observations nominally contains about four dozen scans across the limb giving pole-to-pole coverage for two orbits at roughly 100 and 290 W longitude over the Tharsis and Syrtis/Hellas regions, respectively. At each longitude, limb scans are spaced roughly 10 degrees apart in latitude, with a vertical spatial resolution on the limb of roughly 800 m. Radiative transfer modeling is used to model the observations. We compute synthetic CRISM limb spectra using a discrete-ordinates radiative transfer code that accounts for multiple scattering from aerosols and accounts for spherical geometry of the limb observations by integrating the source functions along curved paths in that coordinate system. Retrieved are 14-point vertical profiles for dust and water ice aerosols with resolution of 0.4 scale heights between one and six scale heights above the surface. After the aerosol retrieval is completed, the abundances of C02 (or surface pressure) and H20 gas are retrieved by matching the depth of absorption bands at 2000 nm for carbon dioxide and at 2600 run for water vapor. In addition to the column abundance of water vapor, limited information on its vertical structure can also be retrieved depending on the signal available from aerosol scattering.

Smith, Michael D.↗
NASA NTRS1988

CSDUST3 - A radiation transport code for a dusty medium with 1-D planar, spherical or cylindrical geometry

The program solves the radiation transport problem in a dusty medium with one-dimensional planar, spherical or cylindrical geometry. It determines self-consistently the effects of multiple scattering, absorption, and re-emission of photons on the temperature of dust grains and the characteristics of the internal radiation field. The program can treat radiation field anisotropy, linear anisotropic scattering, and multi-grain components. The program output consists of the dust-temperature distribution, flux spectrum, surface brightness at each frequency and the observed intensities (involving a convolution with a telescope beam pattern).

Egan, Michael P.↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

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).

Anisotropy↗
DOE OSTI2026

Effect of a transparent water-cooled cathode on the neutron production rate in the glow discharge-type fusion neutron source

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.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

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.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Suppression and enhancement of electromagnetic pulses from laser–target interactions by strong magnetic fields

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.

electromagnetic pulse↗
DOE OSTI2026

The role of thermal conduction on the size scaling of laser direct-drive implosions

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.

Ablation↗
DOE OSTI2026

Dark Energy Survey Year 3 results: Simulation-based 𝑤CDM inference from weak lensing and galaxy clustering maps with deep learning: Analysis design

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.

Thomsen, A. [Zurich, ETH] (ORCID:0000000203099021)↗
DOE OSTI2025

Frequency shifts in the electron paramagnetic resonance spectrum of 39 K due to spin-exchange collisions with polarized 3 He and precise 3 He polarimetry

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.

Zeeman effect↗